作者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/cn.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