2015-10-29 182 views
-1

給定一個迴歸模型從線性模型值:估計預測中的R

Y = B0 + B1(x)的

其中x和y是連續的。

擬合模型後,我希望估算y的預測平均值和95%CI當x是在一定的值,比方說,100

在Stata,它可以與邊緣來實現:

reg y x 
margins, at (x = 100) 

在SAS,它可以用做估算:

proc glm; 
model y = x/clparm solution; 
estimate "Test x = 100" intercept 1 x 100; 
run; 

我的問題是:如何實現中的R相同的動作?我嘗試了lsmeans包,但如果我的模型中沒有任何分類變量,它似乎不起作用。

+3

'predict.lm'嗎?只要通過'interval =「confidence」' – jenesaisquoi

+2

@bunk - 如果數據超出了現有數據,那麼'預測'間隔可能更合適。 – thelatemail

+0

@thelatemail確實!謝謝 – jenesaisquoi

回答

1
predict(fit,newdata=data.frame(x=100),interval="confidence") 

(我不@ thelatemail的意見,預測區間爲優選同意;如果你想允許殘差預測區間被指定...)

+0

謝謝!這非常有幫助。 –