存储嵌套循环的输出

Question

我正在尝试做一个逻辑回归的嵌套循环。

我尝试为离散化值和每个类运行一个循环。这是到目前为止的代码。我无法获得每个不同迭代的输出。

class <- c(1,2,3,4,5)
discretization_value <- seq(0.25, 0.75, by =0.05)
output<-data.frame(matrix(nrow=500, ncol=5))
names(output)=c("discretization_value", "class", "var1_coef", "var2_coef", "var3_coef")

for (i in discretization_value){
  for (j in class) {
  df$discretization_value <- ifelse(df$score >= i,1,0)
  result <- (glm(discretization_value ~ 
                   var1 + var2 + var3, 
                 data = df[df$class == j,], family= "binomial"))
  output[i,1] <- i
  output[i,2] <- j
  output[i,3] <- coef(summary(result))[c("var1"),c("Estimate")]
  output[i,4] <- coef(summary(result))[c("var2"),c("Estimate")]
  output[i,5] <- coef(summary(result))[c("var3"),c("Estimate")]
  }
}

我的df的一小段

class score var1 var2 var3
1 0.3 0.18 0.33 356
1 0.5 0.22 0.55 33
1 0.6 0.77 0.44 35
2 0.9 0.99 0.55 2
3 0 0 0 0
3 0.4 0.5 0.11 5
4 0 0.6 0 7
4 0 0.6 0 9
4 0.6 0.2 0.1 6

Answer 1

这会是问题所在吗？data = dfdf$class == j，，family=“二项式”)

我会试着去掉方括号前的逗号。