所以当Xh = 28时,我试图找到回归线的95%置信区间,但我不确定如何在R中使用它.X表示ACT得分而Y表示GPA。到目前为止这是我的剧本
attach(GPA1)
plot(GPA ~ ACT,
xlab="ACT scores", ylab="GPA",
main="GPA vs. ACT scores")
abline(score$coef,col="red")
score<-lm(GPA~ACT)
GPA.clim <- predict(score, new<-data.frame(ACT=28),se.fit=TRUE,interval="confidence")
GPA.clim
GPA.clim<-predict(score,interval="confidence")
GPA.clim
GPA.plim<-predict(score,interval="predict")
GPA plim
我不确定如何在x = 28的情况下对预测带和置信带使用metline。
答案 0 :(得分:0)
#Here is some code for it to work
#first create a data frame with the only value being a act score of 28
#and name it act28, and in this case gpafit is lm(GPA ~ ACT).
你必须确保计算W值是f分布的两倍并且它是平方根
ci.wh <- function(gpafit, act28, alpha = 0.05)
{
df <- 118
W <- sqrt( 2 * qf(1 - alpha, 2, df) ) #2.43
ci <- predict(gpafit, act28, se.fit = TRUE)
x <- cbind(
'x' = act28,
's' = ci$se.fit,
' fit' = ci$fit,
'lwr' = ci$fit - W * ci$se.fit,
'upr' = ci$fit + W * ci$se.fit)
return(x)
}
ci.wh(gpafit, act28, alpha=0.05)