Confidence Band в R

Итак, я пытаюсь найти 95% -ную доверительную полосу для линии регрессии, когда Xh = 28, но я не уверен, как заставить ее работать в R. X представляют ACT-оценку, а Y представляют GPA. вот мой сценарий до сих порConfidence Band в R

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

Я не знаю, как использовать metline в этом случае как для группы прогнозирования и доверительной полосы при х = 28.

источник

2015-10-09 april510

** Какой результат вы получите, что вы уже пробовали? ** Вы не предоставили никаких данных, поэтому ваши вызовы на печать не являются информативными для читателя. –

Эта ссылка может быть использована для вас, http://stackoverflow.com/questions/15180008/how-to-calculate-the-95-confidence-interval-for-the-slope-in-a-linear-regressio? – user1945827

#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 значения, которые в два раза диафрагменное распределения и квадратный корень его

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)

источник

2016-02-15 06:58:59

ps Мы оба в регрессионном анализе, по-видимому, ха-ха, потому что я просто работал над той же проблемой –

ответ

Смежные вопросы