search for: recordyear

I am try to cross validate some logistic regressions. cv.glm allows me to do this randomly but since I have data over a number of years and over a number of distince areas, I would like to cross-validate temporarly and spatially. I've already attached the temporal and spatial attributes to my data, but I'm unsure as to how to achieve this, as it seems I can't use cv.glm for this

...purposes, but since I'm not deeply familiar with R programming, I don't know why it doesn't work. Now checking the r-help list faq it seems this is an appropriate question. I've included my attempted function below. The error I get is: logcv(basp.data, form, 'basp', 'recordyear') Error in order(na.last, decreasing, ...) : Argument 1 is not a vector My questions are, why doesn't this work, and how do I fix it. I'm using the formula function to create the formula that I'm sending to my function. And the mdata is a data.frame. I'm assumed that...

...onse') < 0.5, 0, 1)) return(ret) } # logerrorest takes four arguments # mdata is a data frame holding the data to be modeled # form is the output of the is.formula function # rvar is the response variable as a string, for example 'birdx' # fvar is the fold variable, for example 'recordyear' logerrorest <- function(mdata, form, rvar, fvar) { require(Hmisc) require(ipred) # determine index of variables rpos <- match(rvar, names(mdata)) fpos <- match(fvar, names(mdata)) # get fold values and count for each group vardesc <- describe(mdata[[fpos]])$values fvarlist <-...

...the model. Thanks in advance. T -- Trevor Wiens twiens at interbaun.com The significant problems that we face cannot be solved at the same level of thinking we were at when we created them. (Albert Einstein) ---------------------------- summary(datamatrix) siteid block recordyear cclo 564-125: 5 Min. :1.000 Min. :2000 Min. :0.0000 564-130: 5 1st Qu.:2.000 1st Qu.:2001 1st Qu.:1.0000 564-135: 5 Median :3.000 Median :2002 Median :1.0000 564-140: 5 Mean :3.042 Mean :2002 Mean :0.7509 564-145: 5 3rd Qu.:...