I would appreciate some help writing R code to plot predicted species richness vs. observed species richness (nat_est) with 95% CI lines. I'm using glm to get model coefficients and remove-one cross validation to get predicted (cv.glm). Here is what I have for the code so far: SRCOUNT <- read.table(file.choose(),header=T) library(boot) library(MASS) quk.native <- glm.nb(nat_est ~ UK + I(UK^2), SRCOUNT, link=log) cv.quk <- cv.glm(SRCOUNT, quk) Many thanks, Jeff **************************************** Jeffrey A. Stratford, Ph.D. Postdoctoral Associate 331 Funchess Hall Department of Biological Sciences Auburn University Auburn, AL 36849 334-329-9198 FAX 334-844-9234 http://www.auburn.edu/~stratja