In using the NLME package (R 2.6.1 for Windows), I am having a problem in running an R script that used to run with no problems using a Linux OS in 2004. So I am wondering if during these last ~3 yrs we had major changes in the syntax of the NLME package that I am not aware. 
 
This is the R script:
 
library(nlme)
 
treat=as.factor(c(1,2,1,2,1,2,1,2))
mouse=as.factor(c(1,1,2,2,3,3,4,4))
time=as.factor(c(1,1,1,1,2,2,2,2))

mouse.anova=function(vec){lme(vec~treat+time+treat*time,random=~1|mouse/time)} 
 
lme.out=apply(exprmat,1,mouse.anova)
 
This is the ERROR message:
 
Error in lme.formula(vec ~ treat + time + treat * time, random = ~1 |  : 
  nlminb problem, convergence error code = 1
  message = singular convergence (7)
 
I will appreciate your help in pointing mistakes in this R script. 
Thank you very much.
Roger

Roger L. Vallejo, Ph.D.
Computational Biologist & Geneticist
U.S. Department of Agriculture, ARS
National Center for Cool & Cold Water Aquaculture
11861 Leetown Road
Kearneysville, WV 25430
Voice:  (304) 724-8340 Ext. 2141
Email:  roger.vallejo@ars.usda.gov
http://www.ars.usda.gov/pandp/people/people.htm?personid=37662


	[[alternative HTML version deleted]]

