假设模型是幂函数形式的模型Cden~a*DBH^b*NDR^c 以TH和PART进行分组，代码要怎么写？
nlme1<-nlme(Cden~a*DBH^b*NDR^c,data=data, method='ML', fixed=a+b+c~1, random= a ~ 1 | PART/TH ,start=c( a=1,b=1,c=1 ) )
这个感觉出来的结果不对。求大神指导一下。