我正在使用R中的for语句对不同的个体进行条件逻辑回归分析。这个代码非常简单:
for(ID in unique(Hour168Fin$BAND)){
modelone = clogit(Hour168Fin$OBSERVED ~ Hour168Fin$LNSTEPLENG + Hour168Fin$PowCross + Shrub +
strata(Hour168Fin$STEPID), data=Hour168Fin, subset = which(ID==Hour168Fin$BAND))
我对输出的非常特定的部分感兴趣,所以我构造了输出,以便使用它来准确地给出我需要的系数:
x1beta = as.numeric(summary(modelone)$coef[1,1])
x2beta = as.numeric(summary(modelone)$coef[2,1])
x3beta = as.numeric(summary(modelone)$coef[3,1])
x1SE = as.numeric(summary(modelone)$coef[1,3])
x2SE = as.numeric(summary(modelone)$coef[2,3])
x3SE = as.numeric(summary(modelone)$coef[3,3])
x1pvalue = as.numeric(summary(modelone)$coef[1,5])
x2pvalue = as.numeric(summary(modelone)$coef[2,5])
x3pvalue = as.numeric(summary(modelone)$coef[3,5])
modelAIC = AIC(modelone)
results = table(x1beta, x1SE, x1pvalue, x2beta, x2SE, x2pvalue, x2beta, x2SE, x2pvalue, modelAIC, rownames = ID)}
在R中,我可以看到我正在寻找的格式的所有结果,但是当我使用它来将这些结果变成csv时:
write.csv = (results, file = "TrialOut.csv")
我只获得1个唯一ID的结果。我已经尝试在write语句中嵌入write.csv语句,并在其外部使用它具有相同的结果。有什么建议?我真的很困惑,因为我可以在R中看到结果,但似乎无法将其转换为csv。
谢谢你的时间!
答案 0 :(得分:0)
尝试在循环中包含write.csv
调用,并使用append = TRUE
:
for (...) {
# ...
# ...
write.csv(results, file = "someFile.csv", append = TRUE)
}