xingzhaoh 发表于 2013-1-14 08:35 
你可以把程序传上来,然后大家一起修改一下
EQ<-read.table("data.txt",header=TRUE)
EQ.nlmeR002<-nlme(Y)~a+b*X+c*log10(sqrt(R^2+D^2))+e*Fn+f*Fr,data=EQ,fixed=list(a~1,b~1,c~1,D~1,e~1,f~1),random=c~1|X,start=c(-1,0.8,-2,10,-0.2,0.3),method="ML")
summary(EQ.nlmeR002)
程序很短,,X,R, Fn ,Fr为自变量,Y为因变量,
a,b,c,d,e,f是需要拟合的系数。混合效应是把组与组(不同的X之间)的区别用系数的变化来表示了,我这里只把C为可变的,理论上,a,c,D,e,f,都是可变的。