作者meg97 (meg)
看板Statistics
標題[程式] 用R求信賴區間
時間Fri Oct 6 15:32:07 2017
[程式問題]:
有一data S,
ID Age FEV Height Sex
1 9 1.95 58 M
2 14 2.23 66 F
3 13 3.20 61 F
.
.
.
幫他建立了一個簡單線性回歸模型model.S
code:
model.S<-lm(S$FEV~S$Age)
題目:用這個模型來預測新個案K的FEV的95%信賴區間
(mean response及predicted response的95%信賴區間皆要)
K:
Age Height Sex
15 57 M
目前想法是
先建立一個矩陣K:ID Age Height Sex
66 15 57 M
再A<-rbind(S,K),
再model.A<-lm(A$FEV~A$Age)
再confidence<-predict(model.A,interval='confidence')
prediction<-predict(model.A,interval='prediction')
然後各看第66項就可以得到答案,
不過因為K缺FEV項沒有辦法合併起來
或者直接使用
confint(model.S,level=0.95)
得到的結果:
2.5% 97.5%
(Intercept) 1.124340408 3.2695921
S$Age 0.001680815 0.1580307
又代表什麼意思呢?
不知道版上的前輩們有沒有比較好的方法?
這是我的生統課作業,寫了前三題,最後一題想了一天無法解決,希望大家可以幫忙
R初學者,請多包涵~謝謝大家
[軟體熟悉度]:
新手
--
※ 發信站: 批踢踢實業坊(ptt.cc), 來自: 140.112.206.184
※ 文章網址: https://webptt.com/m.aspx?n=bbs/Statistics/M.1507275130.A.30F.html
1F：推 asdiy: 你要把 回歸係數成 k的係數相加才對 10/09 00:45
2F：推 asdiy: Prediction(Model.S,K,interval="confident") 10/09 00:51
3F：→ asdiy: 選字錯誤 自行修改 10/09 00:51
4F：→ meg97: 謝謝a大!! 10/11 11:15