R

时间:2015-10-09 16:07:58

标签: r regression bandwidth confidence-interval

所以当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。

1 个答案:

答案 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)