求助关于nlme程序中如何在固定作用图形上加随机作用曲线

是 否 +2 论坛币

k人 参与回答

经管之家送您一份

应届毕业生专属福利!

求职就业群

赵安豆老师微信：zhaoandou666

经管之家联合CDA

送您一个全额奖学金名额~ !

立即领取

感谢您参与论坛问题回答

经管之家送您两个论坛币！

+2 论坛币

求助高人：

非线性混合模型拟合完后如何在原图形中加入随机作用曲线。

扫码加我 拉你入群

请注明：姓名-公司-职位

以便审核进群资格，未注明则拒绝

分享0 收藏0 回帖

关键词：nlme nlm LME 混合模型 模型拟合 程序 如何 模型

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
这是我拟合的混合模型程序和结果，就是想画既有原来普通拟合的曲线还能有含有混合作用拟合出来的曲线。求高人指点。谢谢