Question

我试图在循环中将模型摘要数据导出到excel中。我需要输出2个变量＆＃39;系数（变量gain＆amp; loss）并且我已经成功地编写了截距和变量1的系数，但是R告诉我没有找到第三个变量的对象。

我的代码：按参与者编号PID运行模型，PIDs是PID的列表。

for (i in 1: length(PIDs)) {
   subject<-df[df$PID == PIDs[i],]
   myModel <- glm(gamble~Gain + Loss, data = subject, family=binomial)
   summ <- summary(myModel)

   #save results
   ID[i] <- subject$PID
   intercept_coef[i]<-summ$coefficients[1,1]
   gain_coef[i]<-summ$coefficients[2,1]
   loss_coef[i]<-summ$coefficients[3,1]

}

系数汇总表如下所示，我注意到该表已关闭，因为标题与每列不对应。可能是这个问题吗？

 Estimate   Std. Error     z value  Pr(>|z|)
 (Intercept) 13.4214135 3353.1375049 0.004002643 0.9968064
 Gain         0.2929938    0.1635471 1.791494960 0.0732139
 Loss         8.3144005 1619.8731372 0.005132748 0.9959047

错误：

occurrednumber of items to replace is not a multiple of replacement length 

Error in loss_coef[i] <- summ$coefficients[3, 1] : 
 object 'loss_coef' not found

这是什么问题？我可以Intercept和Gain一切正常。

谢谢！

Answer 1

你有没有机会初始化lostcoef（或者loss_coef，也许可以查看打字）？

Answer 2

在我们将值存储在循环中之前，需要初始化变量。

初始化：

ID <- 0
intercept_coef <- 0
gain_coef <- 0
loss_coef <- 0

#loop

在R中导出模型汇总表时找不到对象

2 个答案: