Re: [R] validation logistic regression

2010-05-26 Thread Frank E Harrell Jr
00), how I tell R that omit one by one data? Is validationsize=100? Thanks alot Azam --- On Wed, 5/26/10, Joris Meys wrote: From: Joris Meys Subject: Re: [R] validation logistic regression To: "azam jaafari" Cc: r-help@r-project.org Date: Wednesday, May 26, 2010, 5:00 AM Hi

Re: [R] validation logistic regression

2010-05-26 Thread azam jaafari
ct: Re: [R] validation logistic regression To: "azam jaafari" Cc: r-help@r-project.org Date: Wednesday, May 26, 2010, 5:00 AM Hi, first of all, you shouldn't backtransform your prediction, use the option type=response instead : salichpred<-predict(salic.lr, newdata=profilevalida

Re: [R] validation logistic regression

2010-05-26 Thread Frank E Harrell Jr
On 05/26/2010 07:00 AM, Joris Meys wrote: Hi, first of all, you shouldn't backtransform your prediction, use the option type=response instead : salichpred<-predict(salic.lr, newdata=profilevalidation,type="response") limit<- 0.5 salichpredcat<- ifelse(salichpred If you want to use split-sampl

Re: [R] validation logistic regression

2010-05-26 Thread Joris Meys
Hi, first of all, you shouldn't backtransform your prediction, use the option type=response instead : salichpred<-predict(salic.lr, newdata=profilevalidation,type="response") limit <- 0.5 salichpredcat <- ifelse(salichpredwrote: > Hi > > I did validation for prediction by logistic regression ac

[R] validation logistic regression

2010-05-26 Thread azam jaafari
Hi   I did validation for prediction by logistic regression according to following:   validationsize <- 23 set.seed(1) random<-runif(123) order(random) nrprofilesinsample<-sort(order(random)[1:100]) profilesample <- data[nrprofilesinsample,] profilevalidation <- data[-nrprofilesinsample,] salich<-p