[R] Nolinear mixed-effects models (nlme)

2005-07-01 Thread Alex Bach
Hello, I am trying to fit a nonlinear model of the form of: A*x^b*exp(-c*x) This represents a lactation curve. I have a bunch of cows, so I want COW to be a random effect. I have been trying the following code with very littel success: > fm1 <- nlme(yield ~ A*(DIM^B)*(exp(-C*DIM)), +

[R] Installing sna package on MacOSX

2005-04-28 Thread Alex Bach
Hi folks, I am trying to install the sna package on a MacOS X 10.3.9 and 10.3.4. In neither case I have been able to get the R CMD Install command to succeed from the Terminal window. With 10.3.9 I get an error with the Make command, and with 10.3.4 I get an error on the Description(???) Has

[R] Obtaining Least Square Means with lme (mixed models)

2004-08-31 Thread Alex Bach
Hello, I was wondering how one could obtain LSM for fixed factors in a mixed-effect model using lme. In other words, if I have model such as: model<-lme(yield~period+treatment+period*treatment, data=data,random=~1|cow) how can I get LSM for period, treatment and their interaction? Thanks a lot!

[R] Repeated measures

2004-06-06 Thread Alex Bach
Dear R-gurus, I am pretty much new on R. I am trying to to do a repeated analysis of a linear mixed model with R, and I consistently fail... The problem is: Cow is the random factor, treatment is the fixed factor. The dependent variable is milk yield, which is measured several times (repeatedly

[R] Repeated measures

2004-06-06 Thread Alex Bach
Dear R-gurus, I am pretty much new on R. I am trying to to do a repeated analysis of a linear mixed model with R, and I consistently fail... The problem is: Cow is the random factor, treatment is the fixed factor. The dependent variable is milk yield, which is measured several times (repeatedly