作者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