First, package the new values of your predictor variables into a data frame:
 
     data.frame(VAR1=8, VAR2=64)
 
Then, use the predict() method to apply the model to the predictor values:
 
     predict(model, newdata=data.frame(VAR1=8, VAR2=64))
 
Finally, tell the predict() method that these are new observations, not
original observations, and you want the confidence interval of the
prediction:
 
     predict(model, newdata=data.frame(VAR1=8, VAR2=64), interval="pred",
level=0.95)
 
That should print the fitted value for OUTCOME, the lower bound, and the
upper bound.
 
To see the standard error, include "predict(..., st.fit=TRUE)" in your call
to predict().
 
Hope that helps.
 
Paul

 
Paul Teetor
Elgin, IL   USA
 
 

	[[alternative HTML version deleted]]

