我一直在尝试使用二进制结果(outcome.csv文件)运行原始文件(预测变量)中所有列的逻辑回归。
raw <- data.frame(matrix(nrow=500, ncol=10))
out <- read.csv(file="outcome.csv", header=T)
models <- list()
res <- list()
a <- colnames(raw)
for(i in 1:length(raw){
models[[i]] <- summary(glm(out$blue ~ raw[,i] + out$sex, data=raw , family= binomial ) )
res[[i]] <- paste("logistic", a[[i]], ".txt", sep="")
write.table(models, res, row.names=FALSE, quote=FALSE, sep="\t")
}
它一直在说这个
Error in model.frame.default(formula = out$blue ~ raw[, i] + out$MALE, :
variable lengths differ (found for 'raw[, i]')
任何与循环或申请相关的建议都将非常感谢。 感谢。