Dear list,
I tried to apply the logistic regression to different response variables from
a dataframe and would like to store the results using the function
logistic.display from the "epicalc" package in a list, but got an
error message "Error in eval(expr, envir, enclos) : y values must be 0
<= y <= 1". All the response variables have value of 0 or 1. It
worked well for summary and anova but not for logistic.display in the loop.
Thanks for any help,
Tom
logreg1<-vector("list", length(l))
logreg.anov1<-vector("list",length(l))
logreg.summ1<-vector("list",length(l))
logreg.conf1<-vector("list",length(l))
for (i in c(13:16)){
logreg1[[i-x]]<-glm(dat[,i]~group + age, family=binomial, data=Ndat)
logreg.anov1[[i-x]]<-anova(logreg1[[i-x]])
logreg.summ1[[i-x]]<-summary(logreg1[[i-x]])
logreg.conf1[[i-x]]<-logistic.display(logreg1[[i-x]], crude=FALSE)
}
---------------------------------
Går det långsamt? Skaffa dig en snabbare bredbandsuppkoppling.
Sök och jämför priser hos Kelkoo.
[[alternative HTML version deleted]]
