R預測區間

Question

47個obs和5個變量的數據集（男性編碼爲0，女性爲1）試圖預測男性的平均狀態，收入和口頭會花在95％CI上。

我用我的lm<-spending ~ status + income + verbal + sex, teenspend來獲得平均值。 我發現我的係數爲：

mdl$coefficient 
(Intercept) sexfemale  status  income 
22.55565063 -22.11833009 0.05223384 4.96197922 
     verbal 
-2.95949350 

predict(mdl, sex=0, interval='confidence', level=0.90) 

一些問題：我用上面的預測，但我得到的所有的意見，我怎麼找到我的預測？

 fit   lwr  upr 
1 -10.6507430 -21.4372267 0.1357407 
2 -9.3711318 -21.9428731 3.2006095 
3 -5.4630298 -15.0782882 4.1522286 
4 24.7957487 12.5630143 37.0284831 

請澄清一下？

Answer 1

查看predict.lm的文檔，你會發現在這裏不能使用參數sex=0。預測方法會忽略該參數，因此您會得到數據中所有觀測值的擬合值加置信區間。您可以指定下列方式預測：
predict(mdl, newdata=teenspend[teenspend$sex==0,], interval="confidence")
如果你確實需要一個預測區間 - 你在你的問題的標題中使用它 - 你應該選擇interval="prediction"。