[R] cv.glm help

Jeffrey Stratford stratja at auburn.edu
Thu Nov 17 18:44:06 CET 2005


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




More information about the R-help mailing list