Re: [R] Mac woes

2009-01-08 Thread Prof Brian Ripley
On Thu, 8 Jan 2009, Stuart Vyse wrote: I am new to R but have spent a little time messing around with it. I am running version 2.8.1 under Leopard OS 10.5 on a power PC and also on my laptop, which is Intel. Unless, I am wrong, this version is just a bit buggy. 1. In the R Commander, there is n

Re: [R] longtable example

2009-01-08 Thread Dieter Menne
Felipe Carrillo yahoo.com> writes: > > Could someone try to run this example? I can't make it work. It would be better if you would show a complete sample; the error could be in the header part, and chances are better to get a reply when other do not have to guess what to add. > \begin{smal

[R] shake rattle() and roll

2009-01-08 Thread Dr Eberhard W Lisse
Graham, nice package this :-)-O. It was easy to set up MySQL but I struggled immensely to get the stock ODBC working with PostgreSQL (8.3), which has been discussed on the mailing list before (Prof Brian): Warning messages: 1: In odbcDriverConnect(st, ...) : [RODBC] ERROR: Could not SQLDriverC

Re: [R] AT&T Researchers and the New York Times

2009-01-08 Thread Dr Eberhard W Lisse
Robert, go ahead, fix whatever bothers you, this is Open Sauce, not Jet-Engine Science :-)-O el On 09 Jan 2009, at 07:41 , Johannes Huesing wrote: stephen sefick [Fri, Jan 09, 2009 at 03:16:11AM CET]: It has worked wonders for me over the last years. On Thu, Jan 8, 2009 at 8:58 PM, Ro

Re: [R] Mac woes

2009-01-08 Thread John Fox
Dear Stuart, > -Original Message- > From: r-help-boun...@r-project.org [mailto:r-help-boun...@r-project.org] On > Behalf Of Stuart Vyse > Sent: January-08-09 10:41 PM > To: r-help@r-project.org > Subject: [R] Mac woes > > I am new to R but have spent a little time messing around with it.

Re: [R] AT&T Researchers and the New York Times

2009-01-08 Thread Johannes Huesing
stephen sefick [Fri, Jan 09, 2009 at 03:16:11AM CET]: > It has worked wonders for me over the last years. > On Thu, Jan 8, 2009 at 8:58 PM, Robert Wilkins wrote: > > And by the way, ARE YOU GUYS EVER GOING TO FIX your mailing list > > platform? It is extremely user-unfriendly and a technologic

[R] Programming Question (setting ylim generally)

2009-01-08 Thread stephen sefick
library(StreamMetabolism) day <- (structure(c(10.08, 10.08, 10.05, 10.03, 10, 9.98, 9.96, 9.95, 9.95, 9.96, 9.96, 9.98, 10.01, 10.05, 10.06, 10.09, 10.11, 10.11, 10.13, 10.13, 10.15, 10.15, 10.13, 10.14, 10.11, 10.13, 10.14, 10.13, 10.12, 10.13, 10.14, 10.16, 10.18, 10.19, 10.23, 10.27, 10.33, 10.

Re: [R] VCOV Source Code

2009-01-08 Thread Henrik Bengtsson
> stats::vcov.lm Error: 'vcov.lm' is not an exported object from 'namespace:stats' > stats:::vcov.lm function (object, ...) { so <- summary.lm(object, corr = FALSE) so$sigma^2 * so$cov.unscaled } /Henrik On Thu, Jan 8, 2009 at 1:44 AM, Carlos J. Gil Bellosta wrote: > Hello, > > You can

[R] VCOV Source Code

2009-01-08 Thread Murray Cooper
Mostly I am interested in using R for statistics. I am also interested in being able to look at source code. I hope to be able to write extensions. I tried the suggestion below but was unable to access vcov.lm. methods(vcov) [1] vcov.Arima* vcov.glm* vcov.lm*vcov.mlm* vcov.nls* Non-v

[R] Mac woes

2009-01-08 Thread Stuart Vyse
I am new to R but have spent a little time messing around with it. I am running version 2.8.1 under Leopard OS 10.5 on a power PC and also on my laptop, which is Intel. Unless, I am wrong, this version is just a bit buggy. 1. In the R Commander, there is no direct entry of data. The Edit men

Re: [R] survey statistics, rate/proportions with standard errors

2009-01-08 Thread Stas Kolenikov
On 1/8/09, Robert Wilkins wrote: > what does R have to compare with , say , proc surveymeans, estimate survey > means/proportions with standard errors, using Taylor methods? survey package by Thomas Lumley has pretty much everything you would need. Way more than SAS, at any rate. There are means

Re: [R] Problem with aovlmer.fnc in languageR

2009-01-08 Thread Chuck Clifton
I have been corresponding with Harald Baayen about this problem, and it appears that something is broken with the mcmcsamp function that pvals.fnc calls. Running pvals.fnc on an lmer object actually changes the contents of the object (get a summary before and after running pvals.fnc and you'll see

Re: [R] Extracting File Basename without Extension

2009-01-08 Thread jim holtman
You can use 'sub' to get rid of the extensions: > sub("^([^.]*).*", "\\1", 'filename.extension') [1] "filename" > sub("^([^.]*).*", "\\1", 'filename.extension.and.more') [1] "filename" > sub("^([^.]*).*", "\\1", 'filename without extension') [1] "filename without extension" On Thu, Jan 8, 2009 a

Re: [R] AT&T Researchers and the New York Times

2009-01-08 Thread stephen sefick
It has worked wonders for me over the last years. It clunkily and reliably delivers messages to my inbox that are either someone elses question about R or an answer to one of my questions. clunkily yours Stephen Sefick On Thu, Jan 8, 2009 at 8:58 PM, Robert Wilkins wrote: > Is anyone in the lea

[R] Extracting File Basename without Extension

2009-01-08 Thread Gundala Viswanath
Dear all, The basename() function returns the extension also: > myfile <- "path1/path2/myoutput.txt" > basename(myfile) [1] "myoutput.txt" Is there any other function where it just returns plain base: "myoutput" i.e. without 'txt' - Gundala Viswanath Jakarta - Indonesia

[R] What does R have for age-adjusted survey analysis?

2009-01-08 Thread Robert Wilkins
A procedure that , after adjusting for sampling weights, also explicitly does an age adjustment to conform with an age distribution of an older census? __ R-help@r-project.org mailing list https://stat.ethz.ch/mailman/listinfo/r-help PLEASE do read the p

[R] survey statistics, rate/proportions with standard errors

2009-01-08 Thread Robert Wilkins
what does R have to compare with , say , proc surveymeans, estimate survey means/proportions with standard errors, using Taylor methods? [[alternative HTML version deleted]] __ R-help@r-project.org mailing list https://stat.ethz.ch/mailman/listi

Re: [R] Creating file names.

2009-01-08 Thread Tony Breyal
Hi, if i understood correctly, you probably want the something like: ## R format(Sys.time(), "LMSMOD_%b_%d_%y_%H%M") (works on windows, don't know about Mac) Hope that helps a little, Tony Breyal On 8 Jan, 23:47, greggal...@gmail.com wrote: > I have a function that is called from a "sourced" s

[R] Using apply for two datasets

2009-01-08 Thread Christos Argyropoulos
It depends on how the data are arranged ## x<-matrix(c(1,2,3,2,8,2,4,5,6),nrow=3) y<-matrix(c(10,2,13,0,8,4,4.2,5.2,6.2),nrow=3) q<-mapply(t.test,as.data.frame(x),as.data.frame(y)) q ## The ith column of q contain the results of applying t.te

[R] AT&T Researchers and the New York Times

2009-01-08 Thread Robert Wilkins
Is anyone in the leadership of the R-project going to contact the New York Times and clarify that the article gave remarkably short shrift to the people who designed the user interface for R, to a large extent AT&T researchers from an earlier generation? It would be the appropriate thing to do. Th

[R] The AT&T researchers and the New York Times

2009-01-08 Thread Robert Wilkins
Is anyone in the leadership of the R-project going to contact the New York Times and clarify that the article gave remarkably short shrift to the people who designed the user interface for R, to a large extent AT&T researchers from an earlier generation? It would be the appropriate thing to do. Th

Re: [R] Creating file names.

2009-01-08 Thread Don MacQueen
Your easiest starting point is probably file = format(Sys.time()) and then look at the online help for the strptime function to learn about formatting options. For example format( Sys.time() , '%m%d%y' ) Use paste() and some formatting of Sys.time() to construct the form you are looking

[R] Graph Script/X-axis label Problem

2009-01-08 Thread jimdare
Hi, I created the script below to plot multiple graphs. It seems to work pretty well (see BYX7) except when the max value for catch or TACC is very small. In this case the x-axis labels are removed from the x-axis (see BYX10). I can't figure out why given that y co-ordinate for the text functi

Re: [R] Creating file names.

2009-01-08 Thread Peter Alspach
Gregg You could either manipulate the character returned by today(), or try using Sys.time() instead. For example, on Windows, paste(strsplit(date(), ' ')[[1]][c(2,3,5,4)], collapse='_') gives Jan_09_2009_13:40:57 which approaches what you want. Alternatively format(Sys.time(), "%b_%d_%y_%H

Re: [R] help

2009-01-08 Thread Duncan Mackay
Unfortunately a little long Try this (I have used Miktex 2.7 on win XP and R 2.8 for this) I have not gone into the detail of the first chunk I just created the header strings in R and copied to make the Sweave file lines With a little thought it can be done I suppose ## Sweave script \documen

[R] Creating file names.

2009-01-08 Thread greggallen
I have a function that is called from a "sourced" script, whenever certain conditions are met. I have this inside it: if(P > 0.9)save( file = today() , myLMSmodel ); This gives me files that look like this: Thu Jan 8 14:04:43 2009 That works fine on my Mac, but it would be much more convenien

Re: [R] Bar Plot with Connected Points on 1 Y-Axis

2009-01-08 Thread jimdare
Cheers Marc, it works perfectly now. Thanks for your help! Marc Schwartz wrote: > > You can generalize the approach by using something like: > > ... > ylim = c(0, max(DF$TACC, DF$Catch) * 1.1) > ... > > > That would allow you to use the max value of the two columns, multiplied > by a fud

Re: [R] Alignment of image plot overlay

2009-01-08 Thread Greg Snow
The updateusr function in the TeachingDemos package can be used to change the user coordinates of the first plot to match what you think they should before adding to the plot. Hope this helps, -- Gregory (Greg) L. Snow Ph.D. Statistical Data Center Intermountain Healthcare greg.s...@imail.org

[R] Alignment of image plot overlay

2009-01-08 Thread tsippel
I'm having trouble with alignment of a trend line overlayed onto an image plot. The two should be plotted on the same x-axis (time-series). However, the trend line begins about an inch into the image plot x-axis and ends about an inch off of end of the image plot. Once I have the alignment sort

[R] Want to post

2009-01-08 Thread Alon Ben-Ari
Dont be so American Cut the paper work lets get going [[alternative HTML version deleted]] __ R-help@r-project.org mailing list https://stat.ethz.ch/mailman/listinfo/r-help PLEASE do read the posting guide http://www.R-project.org/posting-guide.

Re: [R] 'equivalent' sign in plot title

2009-01-08 Thread Peter Alspach
Try plot(1:10, main=expression(x%==%y)) For more details see ?plotmath or demo(plotmath) HTH ... Peter Alspach > -Original Message- > From: r-help-boun...@r-project.org > [mailto:r-help-boun...@r-project.org] On Behalf Of jjh21 > Sent: Friday, 9 January 2009 10:07 a.m. > To: r-hel

[R] 'equivalent' sign in plot title

2009-01-08 Thread jjh21
Hello, I would like to add the 'equivalent' sign (i.e., three horizontal lines, almost like an equal sign) to a plot. Does R have anything like the LaTeX command '\equiv' that I could put into a plot title? Thanks -- View this message in context: http://www.nabble.com/%27equivalent%27-sign-in-

Re: [R] Fw: Re: R Stacked Histogram

2009-01-08 Thread Jason Rupert
Thank you again for your response.   This worked great.   Quick question about the legend for qplot.  Instead of being outside the plot, is it possible to move the location of the legend to the upper left or right corner of the plot?  Could you possibly provide an example.   Thank you again for

[R] help

2009-01-08 Thread Felipe Carrillo
Hi: I am going through some of the xtable examples and I can't make the one below work. I need to create a longtable on the fly keeping the column headers for all the pages and I thought this example could give some ideas on how to do it. I am using Sweave and xtable to create my tables and gra

Re: [R] Another newbie question

2009-01-08 Thread AllenL
Thank you all! In future I will include examples of my code to make things simpler for you. This is what I settled on: Sp.presence<-Data[,14:31] ##The subset of my data set I'm interested in (the presence/absence data) Sp.presence$Species<-apply(Sp.presence,1,function(x) {c(paste(names(Sp.prese

Re: [R] R in the NY Times

2009-01-08 Thread Carlos J. Gil Bellosta
On Thu, 2009-01-08 at 13:52 -0600, Marc Schwartz wrote: > Reading the posts on SAS-L since yesterday via Google RSS, where the > NYT > article was also posted, some have noted that SAS itself offers online > support forums (http://support.sas.com/forums/index.jspa). From a > quick > review, it look

Re: [R] R in the NY Times

2009-01-08 Thread Marc Schwartz
on 01/08/2009 01:12 PM Andrew Choens wrote: > On Thu, 2009-01-08 at 10:42 -0600, Stas Kolenikov wrote: >> A really good measure for R will be the total # of the downloads of >> r-base for all platforms from all CRAN mirrors (and I would expect >> that # can be found from the servers' logs). Given t

Re: [R] Convert to as.Date

2009-01-08 Thread Rowe, Brian Lee Yung (Portfolio Analytics)
You can use zoo to do this: library(zoo) DA <- c("1991q1", "1993q2") as.Date(as.yearqtr(DA)) -Original Message- From: r-help-boun...@r-project.org [mailto:r-help-boun...@r-project.org] On Behalf Of Shruthi Jayaram Sent: Thursday, January 08, 2009 1:46 PM To: r-help@r-project.org Subject:

Re: [R] R in the NY Times

2009-01-08 Thread Rahul-A.Agarwal
I believe R as a package has everything people with little knowledge of programming can handle quite easily. Moreover even if someone has no programming knowledge can learn R without much effort. I also believe if people in corporate world start using R instead of other complex software which are

Re: [R] R in the NY Times

2009-01-08 Thread ohri2007
Yes I think R as a package can really learn from SAS and SPSS in making GUI more user friendly , even at the risk of dumbing down some complexity.. also as a consultant I know that selling software requires a lot of marketing follow ups..which is why R has lagged behind in actual implementation an

Re: [R] R in the NY Times

2009-01-08 Thread Carlos J. Gil Bellosta
On Thu, 2009-01-08 at 10:42 -0600, Stas Kolenikov wrote: > A really good measure for R will be the total # of the downloads of > r-base for all platforms from all CRAN mirrors (and I would expect > that # can be found from the servers' logs). Hello, You obviate here that many of us are downloadi

Re: [R] R in the NY Times

2009-01-08 Thread Louis Bajuk-Yorgan
As the product manager for S+, I'd like to comment as well. I think the burgeoning interest in R demonstrates that there's demand for analytics to solve real, business-critical problems in a broad spectrum of companies and roles, and that some of the incumbent analytics offerings, in particular SA

[R] F-test

2009-01-08 Thread Swanton0822
hi, i was wondering if i have been giving ANOVA table in R: Response:MPG DF Sum SqMean Sq F vaule Pr(>F) Model 1 216.750216.7506.1272 0.04811* Model.Mixture4 167.44

Re: [R] R in the NY Times

2009-01-08 Thread Andrew Choens
On Thu, 2009-01-08 at 10:42 -0600, Stas Kolenikov wrote: > A really good measure for R will be the total # of the downloads of > r-base for all platforms from all CRAN mirrors (and I would expect > that # can be found from the servers' logs). Given that it is so easy > to download everything nice a

Re: [R] Excluding data with apply

2009-01-08 Thread David Winsemius
Well, testing would have been wise.. The last variable in the sample dataframe had 5 elements and the use "negation" of a character vector is not proper. > df[ , -c("a","d")] Error in -c("a", "d") : invalid argument to unary operator So if you limit yourself to negative indexing of numeric r

Re: [R] cosinor analysis

2009-01-08 Thread Charles C. Berry
On Thu, 8 Jan 2009, Anne Berger wrote: Hallo, I didn?t found any facilities for Halbergs cosinor analysis in R. This analysis is well known in the Chronobiology as the least square approximation of time series using cosine function of known period (in my case of 24hours-period). I trie

Re: [R] how to plot histogram plot and fitted distributions on the same graph

2009-01-08 Thread Greg Snow
Here is a direct quote from the help page for hist: "if 'FALSE', probability densities, component 'density', are plotted (so that the histogram has a total area of one)." That looks like it may apply to your problem, just figure out what the above quote is referring to. The last example on that

Re: [R] R and Excel

2009-01-08 Thread Erich Neuwirth
Can't you just set a link to our server? I do not know how to embed it. I used Camtasia to produce this video, and it does all the embedding. Ajay ohri wrote: > Hi Erich, > > I would like to share and embed the RExcel Training video (just like > youtube allows me to) . How can I do that ? > > R

Re: [R] R and Excel

2009-01-08 Thread Ajay ohri
Hi Erich, I would like to share and embed the RExcel Training video (just like youtube allows me to) . How can I do that ? Regards, Ajay www.decisionstats.com On Fri, Jan 9, 2009 at 12:20 AM, Ajay ohri wrote: > Hi Erich, > I saw that it uses a remote server ( which can be the same machine )

Re: [R] R and Excel

2009-01-08 Thread Ajay ohri
Hi Erich, I saw that it uses a remote server ( which can be the same machine ) to compute. Here is the question- What is the remote server is Amazon EC2 which has upscalaing and downscaling facillity for RAM and CPU... Will it work ? is there a SaaS version of this? Regards, Ajay On Thu, Jan

[R] Convert to as.Date

2009-01-08 Thread Shruthi Jayaram
Hi, I have an vector object that looks like DA <- c("1991q1", "1993q2") (first quarter of 1991 etc) and I want to convert it into a date object using as.Date(). I did this for montly data but am stumped when it comes to dealing with quarterly data and as.Date. Would anyone be able to help?

Re: [R] Excluding data with apply

2009-01-08 Thread Stavros Macrakis
On Thu, Jan 8, 2009 at 5:52 AM, Christian Kamenik < christian.kame...@giub.unibe.ch> wrote: > > 'Apply' is a great thing for running functions on rows or columns of a > matrix: > > X <- rnorm(20, mean = 0, sd = 1) > dim(X) <- c(5,4) > apply(X,2,sum) > > Is there a way to use apply for excluding ro

Re: [R] R Stacked Histogram

2009-01-08 Thread hadley wickham
> # Unfortunately due to a small bug in the current version of ggplot, > you need to use > # jpeg with an e to get jpegs > ggsave(file = "ggsave_qplot_output.jpeg") > > # But I'd recommend using png instead - much better quality > ggsave(file = "ggsave_qplot_output.png") Here's a cute illustration

Re: [R] interpolation to abscissa

2009-01-08 Thread Greg Snow
If you want to just linearly interpolate, then use the functions approx or approxfun from the stats package (one of those that is loaded by default). See the function TkApprox in the TeachingDemos package for an interactive way to plot the approximations with the interpolations plotted. Hope

Re: [R] R and Excel

2009-01-08 Thread Erich Neuwirth
There is RExcel (available by downloading the CRAN package RExcelInstaller. It allows to transfer data between R and Excel, and run R code from within Excel. So you can start with your data in Excel, let R do an analysis, and transfer the results back to Excel. You can write VBA macros which do thi

Re: [R] ftp connections for uploading files

2009-01-08 Thread Duncan Temple Lang
Prof Brian Ripley wrote: Try system() with curl or a decent ftp client (I don't see that package RCurl covers this, but it might despite its description only mentioning HTTP). It does support FTP, and all of the protocols that are supported in the installed libcurl, so it depends the configu

Re: [R] VaR-Monte carlo Simulation, Historic simulation, Variance-Covariance Simulation

2009-01-08 Thread Liviu Andronic
On 1/8/09, Maithili Shiva wrote: > Is there any fucntion in R that will help me calculate Value at Risk (VaR) > using Monte carlo Simulation , Historic simulation and Variance - Covariance > Simulation. > There are some "un-published" Crystal Ball functions for R [1], which may help in running

[R] R and Excel

2009-01-08 Thread ohri2007
Even using the VBA back of Excel to create interfaces with R would make a lot of sense. Suppose I could have access to VBA macros that import and export data into R , it would be great. The R GUI series like Rattle come even closer to Excel...so a VBA _R_ExCel package might be useful to ordinary

[R] lattice question: independent per-row or per-column scaling?

2009-01-08 Thread René J.V. Bertin
Hello - and happy newyear to all of you! I've got some data that I'm plotting with bwplot, a 3x2x3 design where the observable decreases with the principle independent factor, but at different rates. I'd like to get lattice to impose not a single set of axes ranges identical for all panels, but r

Re: [R] Regarding Books on R

2009-01-08 Thread David M Smith
In addition to the other excellent examples already given, I maintain a list of resources for newcomers to R here: http://www.revolution-computing.com/community/resources.php#TOC-Resources-for-Newcomers-to-R I can also suggest this blog category (which will grow further with time): http://blog.r

Re: [R] how to plot histogram plot and fitted distributions on the same graph

2009-01-08 Thread Albyn Jones
You are plotting the histogram in the frequency scale. A quick look at the doc page for hist() would reveal the freq option: hist(x,freq=FALSE) then you can add the densities with lines() albyn Quoting Xin Shi : Dear: I am trying to plot the histogram graph for my observed data.

Re: [R] NAs are not allowed in subscripted assignments

2009-01-08 Thread Prof Brian Ripley
On Thu, 8 Jan 2009, Rune Schjellerup Philosof wrote: Browse[1]> j <- c(1,2,NA) Browse[1]> j[j==1][-1] [1] NA Browse[1]> j[j==1][-2] [1] 1 Browse[1]> j[j==1][-2] <- 2 Error during wrapup: NAs are not allowed in subscripted assignments As far as I can see, I have no NA in the lhs (not after the s

[R] longtable example

2009-01-08 Thread Felipe Carrillo
Hi: Could someone try to run this example? I can't make it work. ## Not run: \begin{small} \setlongtables \begin{longtable}{ <>= cat(paste(c('c', rep('cc', 34/2-1), 'c'), collapse='@{\hspace{2pt}}')) @ } \hline \endhead \hline \endfoot <>= library(xtable) set.seed(2345) res <- matr

Re: [R] Letter-based representation of pairwise comparisons

2009-01-08 Thread hadley wickham
On Thu, Jan 8, 2009 at 5:40 AM, retama wrote: > > Hi! > > I have been working several years with R but it's my first public question. > I hope I'll be clear :) . > > This question is related to obtaining letter-based representation of > non-parametric pairwise comparisons. > > I have a dataframe w

Re: [R] ftp connections for uploading files

2009-01-08 Thread Prof Brian Ripley
On Thu, 8 Jan 2009, Gábor Csárdi wrote: On Thu, Jan 8, 2009 at 5:42 PM, Thomas Loridan wrote: Many thanks for that. I am running R under windows and I am not sure what you mean by this curl? http://curl.haxx.se/ And it comes as standard on Mac OS X and many Linux distros. (Generally people

Re: [R] Excluding data with apply

2009-01-08 Thread David Winsemius
Using indexing and putting a minus sign in front of a vector of column names that you want to exclude would be a typical approach: df <- data.frame(a=LETTERS[1:4], b= rnorm(4), c=rnorm(4), d= letters[5:9]) apply(df[ , -c("a","d")], 2, sum) (Pretty sure this will run properly but don't have

Re: [R] ftp connections for uploading files

2009-01-08 Thread Gábor Csárdi
On Thu, Jan 8, 2009 at 5:42 PM, Thomas Loridan wrote: > Many thanks for that. > I am running R under windows and I am not sure what you mean by this curl? http://curl.haxx.se/ > Is'nt there a way to ftp via an intrinsic R function? That I don't know, perhaps some else does. You might want to do

[R] Why Excel is our friend

2009-01-08 Thread Stavros Macrakis
> > "Some people familiar with R describe it as a supercharged version of > Microsoft's Excel spreadsheet software..." > It is easy to ridicule this line from the NYT article. But this is not only a very sensible comment by a smart reporter, but also one that is good for R: It is good for R beca

Re: [R] Generating GUI for r-scripts

2009-01-08 Thread Bert Gunter
If you are on Windows (you didn't say!) see also ?winMenuAdd, ?winDialog ?choose.files etc. -- Bert Gunter Genentech -Original Message- From: r-help-boun...@r-project.org [mailto:r-help-boun...@r-project.org] On Behalf Of Mike Prager Sent: Thursday, January 08, 2009 8:29 AM To: r-h...@st

[R] ftp connections for uploading files

2009-01-08 Thread Thomas Loridan
Hi all, I would like to upload some plots I create wth R via ftp or something similar but I don t really understand which command/syntax I should use: should I go for make.socket + write.socket or try and create environment variables like frp_proxy_user and then ftp my files? how? many thanks fo

Re: [R] Numerical Integration Problems

2009-01-08 Thread Christos Argyropoulos
Hi, You may want to try the double exponential transformation on the numerator and the denominator on this one. The method is described in detail here: http://projecteuclid.org/DPubS?service=UI&version=1.0&verb=Display&handle=euclid.prims/1145474600 If you want to give it a shot outside

[R] Excluding data with apply

2009-01-08 Thread Christian Kamenik
Dear all, 'Apply' is a great thing for running functions on rows or columns of a matrix: X <- rnorm(20, mean = 0, sd = 1) dim(X) <- c(5,4) apply(X,2,sum) Is there a way to use apply for excluding rows or columns from a matrix to run functions on the remaining rows or columns? I know, I could

[R] cosinor analysis

2009-01-08 Thread Anne Berger
Hallo, I didn´t found any facilities for Halbergs cosinor analysis in R. This analysis is well known in the Chronobiology as the least square approximation of time series using cosine function of known period (in my case of 24hours-period). I tried to write a script but crashed... Can you give m

Re: [R] Converting data frame to symmetric matrix

2009-01-08 Thread poastpd
Thank you. It worked! Quoting Henrique Dallazuanna : Try this: reshape(d, timevar = "x", idvar = "y", direction = "wide") On Tue, Jan 6, 2009 at 6:29 AM, wrote: Dear Sir or Madam, I have the following data frame (which is just a toy example of my larger dataset) L3 <- LETTERS[1:3] x=c(1

[R] Letter-based representation of pairwise comparisons

2009-01-08 Thread retama
Hi! I have been working several years with R but it's my first public question. I hope I'll be clear :) . This question is related to obtaining letter-based representation of non-parametric pairwise comparisons. I have a dataframe with this structure (but with quite more rows and cols): A B C

[R] how to plot histogram plot and fitted distributions on the same graph

2009-01-08 Thread Xin Shi
Dear: I am trying to plot the histogram graph for my observed data. Then plot fitted distribution on the same graph of histogram plot in R. 1. histogram plot y. 2. based on "1", plotting y1 v. x; 3. based on "1", plotting y2 v. x; 4. based on "1", plotting y3 v. x; A

Re: [R] ftp connections for uploading files

2009-01-08 Thread Thomas Loridan
Many thanks for that. I am running R under windows and I am not sure what you mean by this curl? Is'nt there a way to ftp via an intrinsic R function? Thanks again Thomas 2009/1/8 Prof Brian Ripley : > Try system() with curl or a decent ftp client (I don't see that package > RCurl covers this, b

Re: [R] R in the NY Times

2009-01-08 Thread Stas Kolenikov
On 1/7/09, Gabor Grothendieck wrote: > Here is the same number of messages/posts data > for each of S, SAS, R: > - reworked into a 3 column ts class time series > - with Jan 2009 removed since its not complete > - leading and trailing NA rows removed My software of choice is Stata, so here ar

Re: [R] Generating GUI for r-scripts

2009-01-08 Thread Mike Prager
Daren Tan wrote: >>> I have developed some scripts that basically ask for input tab-limited format files, do some processing, and output several pictures or csv. Now I need to have some gui to wrap on top of the scripts, [...] Please advice me if there is any tools or project suitable for my task

[R] interpolation to abscissa

2009-01-08 Thread e-letter
Readers, I have looked at various documents hosted on the web site; I couldn't find anything on interpolation. So I started r and accessed the help (help.start()). (by the way is it possible to configure r to open help in opera instead of firefox?) Initially I read the help for the akima package b

Re: [R] ftp connections for uploading files

2009-01-08 Thread Prof Brian Ripley
Try system() with curl or a decent ftp client (I don't see that package RCurl covers this, but it might despite its description only mentioning HTTP). From 'man curl' curl offers a busload of useful tricks like proxy support, user authen- tication, ftp upload, HTTP post, SSL connections

Re: [R] R in the NY Times

2009-01-08 Thread Doran, Harold
The open-source mentality is invaluable, as most on this list know. That is what keeps the R evolution progressing at a pace that SAS cannot keep up with. On a side note (a very side note), I am a zealot for an exercise program called Crossfit. Crossfit has adopted the same open-source mentality

Re: [R] xtable-longtable question

2009-01-08 Thread Felipe Carrillo
Thanks Duncan, I'll try that with my data. --- On Wed, 1/7/09, Duncan Mackay wrote: > From: Duncan Mackay > Subject: Re: [R] xtable-longtable question > To: mazatlanmex...@yahoo.com > Cc: r-help@r-project.org > Date: Wednesday, January 7, 2009, 1:15 PM > I use Sweave and longtable but have not

[R] ftp connections for uploading files

2009-01-08 Thread Thomas Loridan
Hi all, I would like to upload some plots I create wth R via ftp or something similar but I don t really understand which command/syntax I should use: should I go for make.socket + write.socket or try and create environment variables like frp_proxy_user and then ftp my files? how? many thanks fo

Re: [R] R in the NY Times

2009-01-08 Thread Max Kuhn
More commentary on Slashdot: http://developers.slashdot.org/article.pl?sid=09/01/07/2316227 __ R-help@r-project.org mailing list https://stat.ethz.ch/mailman/listinfo/r-help PLEASE do read the posting guide http://www.R-project.org/posting-guide.html a

[R] NAs are not allowed in subscripted assignments

2009-01-08 Thread Rune Schjellerup Philosof
Browse[1]> j <- c(1,2,NA) Browse[1]> j[j==1][-1] [1] NA Browse[1]> j[j==1][-2] [1] 1 Browse[1]> j[j==1][-2] <- 2 Error during wrapup: NAs are not allowed in subscripted assignments As far as I can see, I have no NA in the lhs (not after the second subscript anyway). Besides, I have a single value

[R] Reading UTF-8 from MySQL in Windows

2009-01-08 Thread Tomas Krilavičius
Dear R-help, I am reading strings (varchar) from MySQL database and instead of some characters get garbled chars. Strings are encoded in UTF-8. I have tried to change encoding in MySQL to latin7, but it does not help. I was not able to find solution, how to get strings in proper encoding/form.

Re: [R] Faster Printing Alternatives to 'cat'

2009-01-08 Thread jim holtman
Here is one way of doing it. To write out 1 million rows on my system took 21 seconds. > # create some data > dataSize <- 1e6 > foo <- runif(dataSize) > bar <- runif(dataSize) > n <- 1000 # number of items to write out each time > output <- file('/output.txt', 'w') > # now split the indices into

Re: [R] Dataframe with unequal rows

2009-01-08 Thread David Winsemius
I doubt that you have a dataframe with those features, since R would not allow such an event; more likely you have data in a file. If your goal is to determine the number of items, then you should definitely look at: ?count.fields count.fields(filename, sep=",") # would give you 1 + the n

Re: [R] R in the NY Times

2009-01-08 Thread Marc Schwartz
on 01/07/2009 09:47 PM Gabor Grothendieck wrote: > On Wed, Jan 7, 2009 at 10:26 PM, Dirk Eddelbuettel wrote: >> On 7 January 2009 at 18:24, Gabor Grothendieck wrote: >> | By running the code below we see that the: >> | - sum of the three seems to be rising at a constant rate >> | - S is declining

Re: [R] Faster Printing Alternatives to 'cat'

2009-01-08 Thread Petr PIKAL
Hi r-help-boun...@r-project.org napsal dne 08.01.2009 14:26:32: > Dear Jim and Henrik, > > > What exactly is the problem you are trying to solve. > > Is it going to be read by some other program? > > I simply want to print the data out. Surely, this data > will be manipulated (with Excel or ot

Re: [R] Faster Printing Alternatives to 'cat'

2009-01-08 Thread ONKELINX, Thierry
Yoo should have a look at ?write.table dataset <- cbind(foo, bar, qux) #write to a file write.table(x = dataset, file = "your_file_name.txt", sep = "\t", row.names = FALSE) #wrfite to the console write.table(x = dataset, file = "", sep = "\t", row.names = FALSE) HTH, Thierry -

Re: [R] Faster Printing Alternatives to 'cat'

2009-01-08 Thread Gundala Viswanath
Dear Jim and Henrik, > What exactly is the problem you are trying to solve. > Is it going to be read by some other program? I simply want to print the data out. Surely, this data will be manipulated (with Excel or other programming languages) by other people suit to their purpose. Typically the

Re: [R] Faster Printing Alternatives to 'cat'

2009-01-08 Thread jim holtman
What exactly is the problem you are trying to solve. What is going to be done with the data? Is it going to be read by some other program? How much physical memory do you have on your machine? Is there paging occuring due to the size of the objects? Have you consider creating a structure with 1

Re: [R] troubles performing Moran.I test

2009-01-08 Thread Roger Bivand
agrar.uni-giessen.de> writes: > > dear R users, > > I have troubles performing Moran.I test as suggested on > http://www.ats.ucla.edu/stat/r/faq/morans_i.htm Using the ape package which is written for a different application area is not necessarily a good idea - the web page you refer to a

Re: [R] Faster Printing Alternatives to 'cat'

2009-01-08 Thread Henrik Bengtsson
On Thu, Jan 8, 2009 at 3:12 AM, Gundala Viswanath wrote: > Dear all, > > I found that printing with 'cat' is very slow. > > For example in my machine this snippet > > __BEGIN__ > > # I need to resolve to use this type of loop. > # because using write(), I need to create a matrix which > # consume

[R] Faster Printing Alternatives to 'cat'

2009-01-08 Thread Gundala Viswanath
Dear all, I found that printing with 'cat' is very slow. For example in my machine this snippet __BEGIN__ # I need to resolve to use this type of loop. # because using write(), I need to create a matrix which # consumes so much memory. Note that "foo, bar, qux" object # is already very large (

Re: [R] resolution (dpi) problem

2009-01-08 Thread Gabor Grothendieck
On Tue, Apr 29, 2008 at 7:56 AM, Ted Harding wrote: > Certainly for Windows users, the suggested change could be > wise, since Windows (for reasons best known to someone who > should have known better) conceals the extension from the I think the reason for hiding extensions is that its normally a

Re: [R] VaR-Monte carlo Simulation, Historic simulation, Variance-Covariance Simulation

2009-01-08 Thread Carlos J. Gil Bellosta
Yes, there are: replicate and quantile are your friends. You will find better support in the R-Finance list, though. Best regards, Carlos J. Gil Bellosta http://www.datanalytics.com On Thu, 2009-01-08 at 01:36 -0800, Maithili Shiva wrote: > Dear R helpers > > Suppose I have a portfolio of sec

Re: [R] library error on start up, error on update packages

2009-01-08 Thread Prof Brian Ripley
Please do read the rw-FAQ (as requested in the posting guide): it will enlighten you. On Wed, 7 Jan 2009, Greg Blevins wrote: Hello, I am unable to update packages. I am running R version 2.8.1 (2008-12-22) under VISTA. When I start R I get the following error: I presume that is *Windows*

  1   2   >