我有一个练习说
找到当 x2=2300、x7=56 和 x8=2100 时球队赢得的平均比赛次数的 95% 的置信区间。
R 中有一个函数可以直接给出这样的置信区间吗?
我已经考虑过使用该功能confint(f)
,但是当它是关于一个或多个参数时,这个函数给出结果,据我所知,我没有一个参数,而是一个像 beta0+beta1xi 这样的函数,其中已经估计了参数 beta 并且点 xi 将是x2、x7 和 x8。
另一种方法是“手动”进行,但这很复杂,因为我必须计算标准误差、方差、t 值等。
你能帮忙吗?
先感谢您
最佳答案
你需要看的不是 confint
而是 predict.lm
:
您需要为 newdata
参数设置一个与模型拟合中使用的列名称相同的数据框,其中包含您想要预测的设置值。
这是一个显示如何使用 newdata
的示例:
x1<-c(1,2,5,6); x2<-c(3,2,4,1); x3<-c(5,4,3,4); y<-c(21,21,27,23)
res<-lm(y~x1+x2+x3)
predict.lm(res,newdata=data.frame(x1=4,x2=4,x3=2),
interval="confidence")
(即,您将需要 data.frame(x2= ..., x7=...
等形式的内容,但您可以在其中填写所需的值)但是,您还需要告诉它您需要的间隔类型。
(
predict
是通用的;如果您在 predict
对象上调用 lm
,它将调用 predict.lm
,但要获得正确的帮助,您需要直接查看特定函数)关于r - 计算均值的 95% 置信区间,我们在Stack Overflow上找到一个类似的问题:https://stackoverflow.com/questions/53689298/