[程式] 用R求信賴區間
[程式問題]:
有一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://www.ptt.cc/bbs/Statistics/M.1507275130.A.30F.html
推
10/09 00:45,
8年前
, 1F
10/09 00:45, 1F
推
10/09 00:51,
8年前
, 2F
10/09 00:51, 2F
→
10/09 00:51,
8年前
, 3F
10/09 00:51, 3F
→
10/11 11:15,
8年前
, 4F
10/11 11:15, 4F