Access the information in help(predict) help(fitted)
With your reproducible code, run the following predict(model1) predict(model1, type="response") predict(model1, type="response", se.fit=TRUE) require(lattice) xyplot(value+fitted(model1)~treatments, data=mydata1, distribute.type=TRUE, type=c("p","a")) Explore the documentation and examples showed there. At your disposal. Walmes. ----- ..oooO .................................................................................................. ..(....)... 0ooo... Walmes Zeviani ...\..(.....(.....)... Master in Statistics and Agricultural Experimentation ....\_)..... )../.... walmeszevi...@hotmail.com, Lavras - MG, Brasil ............ (_/............................................................................................ -- View this message in context: http://n4.nabble.com/confidence-intervals-for-mean-GLM-tp1100091p1100132.html Sent from the R help mailing list archive at Nabble.com. ______________________________________________ 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.