如果有人能向我解释为什么我在从一个从线性模型中提取系数的函数中导出一些数据时遇到问题,我会很高兴也很感激。我有数百个这样做,我希望建立一个循环来处理它,但已经陷入了早期的障碍。 我正在使用从这个东西更聪明的人那里借来的方法: https://stat.ethz.ch/pipermail/r-sig-ecology/2008-May/000062.html 相关位(数据创建,功能以及最后我尝试导出数据的内容如下)但首先我要提到的是,“export”数据作为单个COLUMN导出到experiment.csv文件中。我被告知write.table函数的Append属性仅适用于行。因此,它会覆盖先前运行的相同命令集,而不是成功附加它们 错误消息的形式如下:(它们都是相同的,每条信息都有一个)。 警告讯息: 1:在write.csv中(export,file =“experiment.csv”,append = TRUE,quote = TRUE ,: 尝试将'append'设置为忽略
#DATA CREATION
# create an empty list
mod <- list()
# start a loop for create 5 objects of class 'lm'
for (i in 1:5) {
x <- rnorm(i*10)
y <- rnorm(i*10)
mod[[paste("run",i,sep="")]] <- lm(y ~ x)
}
# FUNCTION TO EXTRACT DATA
myFun <-
function(lm)
{
out <- c(lm$coefficients[1],
lm$coefficients[2],
length(lm$model$y),
summary(lm)$coefficients[2,2],
pf(summary(lm)$fstatistic[1], summary(lm)$fstatistic[2],
summary(lm)$fstatistic[3], lower.tail = FALSE),
summary(lm)$r.squared)
names(out) <- c("intercept","slope","n","slope.SE","p.value","r.squared")
return(out)}
# FAILED ATTEMPT TO EXPORT
export <-myFun(mod$run1)
write.csv(export, file = "experiment.csv", append = TRUE, quote = TRUE, sep = " ",
eol = "\n", na = "NA", dec = ".", row.names = FALSE,
col.names = FALSE, qmethod = c("escape", "double"),
fileEncoding = "")