similar to: glm.cv delta value?

Dear list members, I am using a cross validation of a generalised linear model (glm). The cv.glm function (from boot package) returns an error as so-called ?delta? value. I would like to get to a (cross-validated) squared q, because I want to directly compare it to the squared correlation coefficient r. I tried to find an an equation for the raw and/or adjusted cross-validation estimate of

Hi, I've done some searching where others have had trouble with this error (or "warning" actually), but I'm unable to solve my problem. I have a data sheet with 13 columns and 36 rows. Each column has exactly the same number of rows. I've created glms and now want to do cross-validation on 2 of them. Please be gentle-- I'm new to R (and statistics, too, for that

I have a couple quick questions about the use of cv.glm for cross-validation. 1. If we have a Poisson GLM with counts Y~Poisson(mu) and ln(mu)=beta0+beta1*x1+..., is the prediction error (delta) that is output from cv.glm provided in terms of the counts (y) or the (mu)? 2. Can cv.glm be used for negative binomial models fit using glm.nb? It appears to work, but since NB models aren't

I have a logistic model fitted with the following R function: glmfit<-glm(formula, data, family=binomial) A reasonable cutoff value in order to get a good data classification (or confusion matrix) with the fitted model is 0.2 instead of the mostly used 0.5. And I want to use the `cv.glm` function with the fitted model: cv.glm(data, glmfit, cost, K) Since the response in the fitted

Dear list, I am modelling fish distributions using the glm-function followed by the step-function, and then want to cross-validate the model via the cv.glm-function from the {boot} package. I am working on fish distributions on coral reefs. The code I have works for one fish species. Since I have 227 fishes, I wrote a loop. Now the cv.glm-function comes up with an error message: "Error in

Hi, everyone, I have a question regarding function cv.glm in library 'boot'. Basically cv.glm can calculate the estimated K-fold cross-validation prediction error for generalized linear models. My question is this: if I am fitting a logit model, what kind of threshold will it use to calculate the prediction error (saved in 'delta')? It will use 0.5 as the threshold or pick a

Hi, I have a question regarding the cv.glm function in the package boot. What is exactly the cost? Is it the threshold value for an estimated value to be classified as either 0 or 1? I have troubles understanding the explanation in R. Lets say I want all estimations >0.65 to be classified as 1s and <0.35 as 0s, how do I write that? And if the cost is something else, how do I set the

Hey all, I need some help with a cross validation. I'm new with R and as well with statistics. I had a group work to create a tool for remote sensing class that extracts the best bands of hyperspectral satellite images that describe vegetation. Its a regression between a linear function of using a normalized differenced index (i-j)/(i+j) while i and j are the bands (in the data these are the

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

Dear list members. I have problems with the usage of cv.glm from the boot package. Here are some parts of the script I wanted to use: data <- read.table("selected_2D.csv", header=TRUE, sep=",") ? glm.fitted <- glm("ydata$ y ~ 1 + density + vsurf_ID6 + vsurf_S ", data=data) error <- cv.glm(data=data, glm.fitted, K=6) ydata$y is a separate data set, where

Hi. Is there a way to get the values predicted from (leave-one-out) cv.glm? It seems like a useful diagnostic to plot observed vs. predicted values. 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

I am trying to write a wrapper for the last example in help(cv.glm) that deals with leave-one-out-cross-validation (LOOCV) for a logistic model. This wrapper will be used as part of a bigger program. Here is my wrapper funtion : logistic.LOOCV.err <- function( formu=NULL, data=NULL ){ cost.fn <- function(cl, pred) mean( abs(cl-pred) > 0.5 ) glmfit <- glm(

Greetings R users, I have been hoping someone would be familiar with this topic. I understand fully everything on this list is from the good graces of those who wish to help. Thanks to those who have helped in multiple circumstances. However, I wanted to post this question once more. I hope someone will read it and perhaps have some additional insight into what I may be able to do. Thank

Greetings R users, I have been working on running plsda and I would like to have the R2 and Q2 values. I know the function R2 from the 'pls' package will generate both R2 and Q2 but they are for each separate class. Is there a way to get the cumulative R2 and Q2 for the whole model? R2(pls.new, estimate="all") Response: B (Intercept) 1 comps 2 comps train

Hello, I would like to use cv.tree to run a 10-fold cross-validation experiment on a tree object to help me choose a tree size. Many users seem to allow their cases to be assigned to CV groups randomly, but I have assigned each case to one of 10 cv groups, such that the data from each of my experimental units is included in only one cv-group. According to the manual for the tree Package

Dear All, I'm having trouble tweaking a forest plot made using the R meta-analysis package metafor. I did the analysis based upon the correlation coeff from studies and plotted the corresponding forest plot easily > q2<-rma(yi,vi,mods=cbind(grupo),data=qim) > q2 > forest (q2,transf=transf.ztor,digits=3, ... ,alim=c(0,1),refline=.5) > text(-1.55,42,"Esp?cie

Georges Monette and I would like to announce a new package, cv, now on CRAN, which implements cross-validation of regression models. Some of the functions supplied by the package: - cv() is a generic function with a default method and computationally efficient "lm" and "glm" methods, along with a method for a list of competing models. There are also experimental

Georges Monette and I would like to announce a new package, cv, now on CRAN, which implements cross-validation of regression models. Some of the functions supplied by the package: - cv() is a generic function with a default method and computationally efficient "lm" and "glm" methods, along with a method for a list of competing models. There are also experimental

Hi I am using clang-6 to compile this C++ code and I see a strange temporary variable allocated at expression address 0x7ff1131536e8. If I change the ternary operator at line 483 to an if-else, the temporary is not allocated. Thanks Variables: ========= FFMPEGVideoCapture ffmpeg_video_capture_; cv::VideoCapture opencv_video_capture_; Function: ======== bool

Hi there, I am trying to use predict() with an object returned by cv.glmnet(), and get the following error: no applicable method for 'predict' applied to an object of class "cv.glmnet" What's wrong? my code: x=matrix(rnorm(100*20),100,20) y=rnorm(100) cv.fit=cv.glmnet(x,y) predict(cv.fit,newx=x[1:5,]) coef(cv.fit) Thanks so much, Asaf -- View this message in context: