[R] Error Message Help

Hi... Can someone please explain to me what this error message means? It only appears if I include the correlation= part, but I have no idea why. Thanks, Joe > ARModel <-lme(Life_Satisfaction ~Time, data = restructuredData, random = > ~Time|Person, method = "ML", na.action = na.exclude, contro

[R] Executing a R-string

Hello... In other languages (e.g. php, perl), you have the ability to create a valid string and execute the string to get the result. For example (in pseudo-R): S<-"which(m==4)" R<-exec(S) I know this does not work, but was wondering if there was an equivalent mechanism that I cannot fin

[R] GAM plots

Hello all... I'm attempting to write my own GAM plot function, so I can overlay it on top of an already existing plot. Problem is that after I do the gam, e.g. m<-gam(...), I cannot match the graph that gam.plot outputs when I attempt to plot the values from m$residuals, m$linear.predict

[R] 2-Y-axes on same plot

Hi... Is this possible in R? I have 2-sets of data, that were collected simultaneously using 2- different data acquisition schemes. The x-values are the same for both. The y-values have different ranges (16.4-37.5 using one method, 557-634 using another). In theory, if you plot both plots

[R] Removing rows from a data frame

Hi... I have a rather large dataframe that I'm trying to remove rows from. I'm issuing the command: dtx[-which(dtx$rdate > "2008-06-16"),] and it tries to print out over 170,000 lines of output. So...I did: options(max.print=1e6) and ran it again. It worked, but when I did a: which(dtx$rdate

[R] R question for Stata users

Hi... In Stata, there is the ability to display scatter plots with data points at the same (x,y) location, using the 'jitter' command of the twoway scatter stata command. Anyone know of a way that I can do the equivalent thing in R? For non-Stata readers, if jitter is enabled in stata, and n

[R] Creating a data.frame

OK...newbie question here. Either I'm reading the docs wrong, or I'm totally confused. Given the following: x<-c("aaa","bbb","ccc") y<-rep(0,3) z<-rep(0,3) is.character(x) [1] TRUE is.numeric(y) [1] TRUE Now...I want to create a data frame, but keep the data types. In reading the docs, I assum

[R] Newbie on plotting

Hi... I am issuing a: plot(x,y), where x is a factor. i expect a box-and-whisker plot and I do get it. The data is 'bucketed' into 10 buckets on the x-axis. When I process some data, I get data in buckets 1,7,8,9 and 10. Problem is when I plot it, it shows only those buckets on the x-axis. if I su

[R] Newbie question on subsets

Hi... I'm working through the book, A Handbook of Statistical Analyses using R by Everitt, and I'm trying to do the following (p. 19 of his book): boxplot(log(marketvalue)~country, data = subset(Forbes2000, country %in% c("United Kingdom","Germany","India","T

[R] Newbie syntax question

Hi... I'm trying to understand the following syntax: cor.test(~mortality + hardness,data=water,method="pearson") which is the same as: cor.test(water$mortality,water$hardness,data=water,method="pearson") Can anyone point me to the correct doc or explain to me how to interpret "~mortality + h