如何创建一个循环以运行2个变量广义回归模型?

时间:2018-08-30 06:16:46

标签: r

我有19个变量,我想运行19个不同的回归,这些回归由数据集中的2个独立变量组成。

*更新-这是我的数据集的结构:

$ Failure_Response_Var_Yr: num  0 0 0 0 0 0 0 0 0 0 ...
 $ exp_var_nocorr_2       : num  4.61 5.99 6.13 3.17 4.4 ...
 $ exp_var_nocorr_3       : num  4.16 5.46 5.24 2.86 3.72 ...
 $ exp_var_nocorr_4       : num  0.00191 2.23004 0.5613 1.07986 0.99836 ...
 $ exp_var_nocorr_5       : num  0.709 2.79 6.846 15.478 11.418 ...
 $ exp_var_nocorr_6       : num  0.724 0.497 1.782 0.156 2.525 ...
 $ exp_var_nocorr_7       : num  0 168.17 92.041 0.584 265.338 ...
 $ exp_var_nocorr_8       : num  -38.64 4.89 1.5 24.8 16.56 ...
 $ exp_var_nocorr_9       : num  116 88.3 56.4 60.6 57.6 ...
 $ exp_var_nocorr_10      : num  0 10.3 0 93.7 0 ...
 $ exp_var_nocorr_11      : num  1.02 1.23 1.31 2.06 1.33 ...
 $ exp_var_nocorr_12      : num  60 140 124 275 203 ...
 $ exp_var_nocorr_13      : num  10.835 5.175 1.838 0.347 0.783 ...
 $ exp_var_nocorr_14      : num  59 60.2 87.2 42.2 84.2 ...
 $ exp_var_nocorr_15      : num  61.9 68.3 99 50.2 103.9 ...
 $ exp_var_nocorr_16      : num  4.4 11.24 8.23 6.9 8.84 ...
 $ exp_var_nocorr_17      : num  6.43 18.62 10.72 15.62 10.35 ...

我写了这段代码:     col17 <- names(my.sample)[-c(1:9,26:29)]

dput(col17)现在给出的结果是:  c("exp_var_nocorr_2", "exp_var_nocorr_3", "exp_var_nocorr_4", "exp_var_nocorr_5", "exp_var_nocorr_6", "exp_var_nocorr_7", "exp_var_nocorr_8", "exp_var_nocorr_9", "exp_var_nocorr_10", "exp_var_nocorr_11", "exp_var_nocorr_12", "exp_var_nocorr_13", "exp_var_nocorr_14", "exp_var_nocorr_15", "exp_var_nocorr_16", "exp_var_nocorr_17" )

`logit.test2 <- vector("list", length(col17))  
#start of loop #
for(i in seq_along(col17)){
  for(k in seq_along(col17)){
    logit.test2[i] <- glm(reformulate(col17[i]+col17[k], "Failure_Response_Var_Yr"), 
                          family=binomial(link='logit'), data=my.sample)
  }
}`

# end of loop #

但是它打印出了这个问题: "Error in col17[i] + col17[k] : non-numeric argument to binary operator"

有人可以帮我分发可以解决此问题的代码吗?

0 个答案:

没有答案