Re: [R] logistic regression model + Cross-Validation

From: Frank E Harrell Jr <>
Date: Sun 21 Jan 2007 - 23:17:32 GMT

nitin jindal wrote:

> If validate.lrm does not has this option, do any other function has it.
> I will certainly look into your advice on cross validation. Thnx.
> nitin

Not that I know of, but easy to program. Frank

> On 1/21/07, Frank E Harrell Jr <> wrote:

>> nitin jindal wrote:
>>> Hi,
>>> I am trying to cross-validate a logistic regression model.
>>> I am using logistic regression model (lrm) of package Design.
>>> f <- lrm( cy ~ x1 + x2, x=TRUE, y=TRUE)
>>> val <- validate.lrm(f, method="cross", B=5)
>> val <- validate(f, ...) # .lrm not needed
>>> My class cy has values 0 and 1.
>>> "val" variable will give me indicators like slope and AUC. But, I also
>> need
>>> the vector of predicted values of class variable "cy" for each record
>> while
>>> cross-validation, so that I can manually look at the results. So, is
>> there
>>> any way to get those probabilities assigned to each class.
>>> regards,
>>> Nitin
>> No, validate.lrm does not have that option. Manually looking at the
>> results will not be easy when you do enough cross-validations. A single
>> 5-fold cross-validation does not provide accurate estimates. Either use
>> the bootstrap or repeat k-fold cross-validation between 20 and 50 times.
>> k is often 10 but the optimum value may not be 10. Code for averaging
>> repeated cross-validations is in
>> along with simulations of bootstrap vs. a few cross-validation methods
>> for binary logistic models.
>> Frank
>> --
>> Frank E Harrell Jr Professor and Chair School of Medicine
>> Department of Biostatistics Vanderbilt University mailing list PLEASE do read the posting guide and provide commented, minimal, self-contained, reproducible code. Received on Mon Jan 22 10:20:58 2007

Archive maintained by Robert King, hosted by the discipline of statistics at the University of Newcastle, Australia.
Archive generated by hypermail 2.1.8, at Sun 21 Jan 2007 - 23:30:34 GMT.

Mailing list information is available at Please read the posting guide before posting to the list.