fmdbh.nlm<-nlme(zdbh~a*(1-exp(-b*age))^c,data=rd,fixed = list(a+b+c ~ 1), random = list(zl = pdDiag(list(a+b~1)),jxm=pdDiag(list(a+b~1))),start=c(a=25,b=0.1,c=1.2))
> summary(fmdbh.nlm)
Nonlinear mixed-effects model fit by maximum likelihood
  Model: zdbh ~ a * (1 - exp(-b * age))^c
Data: rd
       AIC      BIC    logLik
  1632.281 1668.656 -808.1407

Random effects:
Formula: list(a ~ 1, b ~ 1)
Level: zl
Structure: Diagonal
               a         b
StdDev: 0.730029 0.0177251

Formula: list(a ~ 1, b ~ 1)
Level: jxm %in% zl
Structure: Diagonal
               a         b  Residual
StdDev: 1.603547 0.0136651 0.6796938

Fixed effects: list(a + b + c ~ 1)
      Value Std.Error  DF  t-value p-value
a 19.849257 0.5123396 672 38.74238       0
b  0.081901 0.0075680 672 10.82196       0
c  2.044081 0.0591179 672 34.57636       0
Correlation:
  a      b     
b -0.160      
c -0.397  0.325

Standardized Within-Group Residuals:
        Min          Q1         Med          Q3         Max
-3.89657295 -0.64339021 -0.06845265  0.65748673  3.16987259

Number of Observations: 697
Number of Groups:
         zl jxm %in% zl
8        23
这是我拟合的混合模型程序和结果，就是想画既有原来普通拟合的曲线还能有含有混合作用拟合出来的曲线。求高人指点。谢谢