Getting p-values from leave-one-out in R
问题 I have a data frame of 96 observations (patients) and 1098 variables (genes). The response is binary (Y and N) and the predictors are numeric. I am trying to perform leave-one-out cross validation, but my interest is not standard error, but the p-values for each variable from each of the 95 logistic regression models created from LOOCV. These are my attempts thus far: #Data frame 96 observations 1098 variables DF2 fit <- list() for (i in 1:96){ df <- DF2[-i,] fit[[i]] <- glm (response ~.,