similar to: Nested for loop

Dear Lister, I understand merge() can be used to join 2 data frames based on 1 variable. But how about merge based on more than 2 variables? Thank you so much! -- WenSui Liu (http://spaces.msn.com/statcompute/blog) Senior Decision Support Analyst Health Policy and Clinical Effectiveness Cincinnati Children Hospital Medical Center [[alternative HTML version deleted]]

Hi all, We've just had a paper accepted for publication in Nature. We used R for 95% of our analyses (one of my co-authors sneaked in some GenStat when I wasn't looking.). The preprint is available from the Nature web site, in the open peer-review trial section. I searched Nature for previous references to "R Development Core Team", and I received no hits. So I tentatively

I had a question about scan in R. For better code readability, I would like to have lines in the block of data to be scanned that are commented - not just lines that have a comment at the end. For example #age, weight, height 33,128,65 34,56,155 instead of having to do something like 33,128,65 #age, weight, height 34,56,155 Is this at all possible?

Hi, I'm not sure if this is an R problem, a Debian problem, a GNU Emacs problem or an ESS problem. I get the following warning the first time I call help on a command in a session, for example: > ?sessionInfo Warning message: use of NULL environment is deprecated This only occurs in GNU Emacs, not when running R from bash, and not with Xemacs. It started happening when I updated R

>From time to time people request symbolic computations beyond what D() and deriv() etc can provide. A brief look at the internet shows that there are many more or less developed computer algebra packages freely available. Therefore, I wondered if it would be an idea to try to 'integrate' one of these packages in R, which I guess can be done in more or less elegant ways... I do not know

Dear there, I am wondering how to test whether a simple linear regression model (e.g. y=1.05x) is significantly different from a 1 to 1 line (i.e. y=x). Thanks. Regards, Jin [[alternative HTML version deleted]]

Hello Group, What is the meaning of the error. is there any place to look for this. I guess 'atomic' seems to be OOP related concept. thank you srini

I need to specify a model similar to this lme.formula(fixed = sqrt(lbPerAc) ~ y + season + y:season, data = cy, random = ~y | observer/set, correlation = corARMA(q = 6)) except that observer and set are actually crossed instead of nested. observer and set are factors y and lbPerAc are numeric If you know how to do it or have suggestions for reading I will be grateful. eal ps I have

Hello. I am looking for some help using anova's in RGui. My experiment ie data, has a fixed burning treatment (Factor A) 2 levels, unburnt/burnt. Nested within each level of Factor A are 2 random sites (Factor B). All sites are crossed with a fixed temperature treatment (Factor C) 2 levels, 0 degreesC/2 degreesC, with many replicates of these temperature treatments randomly located at

Does someone have code, or point to a source to get it, to model repeated measures zero-inflated poisson data. The data come from a replicated field trial comparing two treatments - a control and a test treatment. Thanks in advance ------------------------------------ Subhash Chandra, DSc Senior Biometrician Primary Industries Research Victoria Department of Primary Industries 1 Ferguson Road

Hi netters, This is probably a silly question,but I can't find the answer after searching the R-help archives online. ok, I have a matrix. I know there is a "10" somewhere in it. Now I want to know the index of the element "10" in this matrix. That is, if X[i,j]=10, I want to know i and j. Is there a R function to do this? Just like the "find" function in

Dear List, Is there a way to extract variance components from lmeObjects or summary.lme objects without using intervals()? For my purposes I don't need the confidence intervals which I'm obtaining using parametric bootstrap. Thanks, Mike [[alternative HTML version deleted]]

Hello, I am trying to create a stripchart for my data, where y axis labels are characters (ie,names of cities). I would like to change the orientation of the y - axis labels, ie perpendicular to y axis. Below is the code i am using: par(srt=90) with(ozone.ne.trim, stripchart(Median~City,main = "stripchart(ozone)")) The par option doesn't seem to working. Kindly help. Thanks

Hello! I just noticed new link on R wiki on R galleries and wanted to share this info with YOU! - R graphical manuals (this is awesome page as there are all help pages of all packages on CRAN and probably even more and all graphics examples are displayed! - more than 8000 images!) http://bg9.imslab.co.jp/Rhelp/ This is a very nice addition to already existing R graph and movies galleries - R

Dear All: I would like to center the data in 'x' by 'group'. The following code scale the data and I have not been able to figure out how to change it so I get the centered data. x <- c(1, 2, 3, 4, 5, 6, 7, 8) group <- c(1,1,1,2,2,2,2,2) unsplit(lapply(split(x,group),scale),group) I would appreciate your help. Ashraf

Hi, How do I get R to output the 95% cutoff from a distribution of a continous variable? summary() only displays a few statistics.... Thanks!

I am using R in a Windows environment. I have a basic question regarding lm(). I have a dataframe ?data1? with ncol=w. I know that my dependent variable is in column1. Is there a way to write the regression formula so that I can use columns 2 thru w as my independent variables? e.g. something like: ? lm(data1[,1] ~ data1[,2:w] ) ? Thanks

Hello, I've just located the illuminating explanation by Douglas Bates on degrees of freedom in mixed models. The take-home message appears to be: don't trust the p-values from lme. Questions: Should I give up hypothesis testing for fixed effects terms in mixed models? Has my time spent reading Pinheiro & Bates been in vain? Is there a publication on this issue? Thanks, Dan Bebber

Can someone please recommend a novice/beginner's reading list for non-programmers learning R? --------------------------------- 8:00? 8:25? 8:40? Find a flick in no time [[alternative HTML version deleted]]

Hi, I am trying to find out a collinearity in explanatory variables with alias(). I creat a dataframe: dat <- ds[,sapply(ds,nlevels)>=2] dat$Y <- Response Explanatory variables are factor and response is continuous random variable. When I run a regression, I have the following error: fit <- aov( Y ~ . , data = dat) Error in "contrasts<-"(`*tmp*`, value =