Displaying 1 result from an estimated 1 matches for "salich".
Did you mean:
salic
2010 May 26
1
validation logistic regression
...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<-profilesample$SALIC.H.1
salic.lr<-glm(salich~wetnessindex, profilesample, family=binomial('logit'))
summary(salic.lr)
salichpred<-predict(salic.lr, newdata=profilevalidation)
expsalichpred<-exp(salichpred)
salichprediction<-(expsalichpred/(1+expsalichpred))
So,
table(sali...