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]] ______________________________________________ R-help@r-project.org mailing list https://stat.ethz.ch/mailman/listinfo/r-help PLEASE do read the posting guide http://www.R-project.org/posting-guide.html and provide commented, minimal, self-contained, reproducible code.