作者avocadosweet (avocadosweet)
看板Statistics
标题[问题] R NLME的model
时间Sat Sep 27 17:24:52 2008
用NLME做的logistic model(参数为asymp、b、rc)
最後可以用summary(model's name)跑出random effect及fixed effect的结果
例如
Random effects:
Formula: asymp ~ 1 | plot
asymp.(Intercept) Residual
StdDev: 0.3480761 0.2534289
Fixed effects: list(asymp + b ~ spacing, rc ~ 1)
Value Std.Error DF t-value p-value
asymp.(Intercept) 7.286464 0.326402 76 22.323603 0.0000
asymp.spacing3m -0.937608 0.406871 76 -2.304433 0.0239
b.(Intercept) 9.311714 2.061776 76 4.516355 0.0000
b.spacing3m 3.752211 2.593882 76 1.446562 0.1521
rc 0.133469 0.008858 76 15.067969 0.0000
请问有办法算出logistic curve不同x值(假设x为时间,y为asymp)下的b值吗?
该用什麽指令?
因为用augPred(model's name)可以看到每年model的预测值
可是我不知道要用什麽指令才可以让R跑出每年的b值
--
※ 发信站: 批踢踢实业坊(ptt.cc)
◆ From: 140.112.82.134
1F:→ clickhere:b是解释变数, 为已知了. 就是你给的资料. 09/29 05:08
2F:→ avocadosweet:谢谢你~不过那是最後一年的b,我指的是每一年的b 10/02 09:11
3F:→ avocadosweet:抱歉我写的可能不太清楚 10/02 09:11