0

也许是微不足道的，但我正在尝试解决这个问题：

我必须数据框，一个有 25 列，另一个有 9 列。现在，我需要做的是拟合多项式方程，其中我的因变量在 25 列的数据框中，而我的自变量在 9 列的数据框中。目前，我将这些列组合在一起并创建了一个名为“my.data”的数据框，因此我当时使用一个自变量遍历因变量。但是，我想自动执行循环中的函数 25 * 9 次。有没有办法做到这一点？

setwd("C:\\......")

 my.data <- read.table("MyData.txt", header = TRUE, sep = "\t")


for(i in seq_along(my.data))
 {

    fit1b <- lm(my.data[ ,i] ~ my.data$V1)
    fit2b <- lm(my.data[ ,i] ~ poly(my.data$V1, 2, raw=TRUE))
    fit3b <- lm(my.data[ ,i] ~ poly(my.data$V1, 3, raw=TRUE))
    poly1 <-capture.output(summary(fit1b))
    poly2 <-capture.output(summary(fit2b))
    poly3 <-capture.output(summary(fit3b))


con = file(description = "MyResults.txt", open="a")
write.table(poly1, file= con, append = TRUE, quote=F, col.names=FALSE, row.names= F)
write.table(poly2, file= con, append = TRUE, quote=F, col.names=FALSE, row.names= F)
write.table(poly3, file= con, append = TRUE, quote=F, col.names=FALSE, row.names= F)
close(con)
 }
4

2 回答 2

1

mapply这是一个使用和使用的绝佳机会expand.grid

例如。

# some dummy data
xx <- data.frame(replicate(5, runif(50)))
yy <- setNames(data.frame(replicate(3, runif(50))), paste0('Y',1:3))
# all combinations
cs <- expand.grid(list(pred = names(xx), resp = names(yy)), stringsAsFactors= FALSE)

# a function to do the fitting
fitting <- function(pred, resp, dd){
  # fit linear model
  ff <- reformulate(pred, resp)
  lmf <- lm(ff, data =dd)
  # create a formula for poly(,2)
  ff.poly2 <- update(ff, .~poly(.,2, raw=TRUE))
  # and poly(,3)
  ff.poly3 <- update(ff, .~poly(.,3, raw=TRUE))
  # fit these models
  lmp2 <- lm(ff.poly2, data = dd)
  lmp3 <- lm(ff.poly3, data = dd)
  # return a list with these three models
  list(linear = lmf, poly2 = lmp2, poly3 = lmp3)
}

biglist <- mapply('fitting', pred = as.list(cs[['pred']]), 
        resp = as.list(cs[['resp']]),
       MoreArgs = list(dd = cbind(xx,yy)), SIMPLIFY = FALSE)

# give this list meaningful names

names(biglist) <- do.call(paste, c(cs, sep = ':'))

lapply 然后，您可以使用一些嵌套语句提取事物/总结事物

例如所有线性模型的总结

lapply(lapply(biglist, `[[`,'linear'), summary)

二次模型的

lapply(lapply(biglist, `[[`,'poly2'), summary)

如果要从print(summary(lm))单个文件中提取信息，例如

capture.output(lapply(biglist, function(x) lapply(x, summary)), file = 'results.txt')

将创建一个名为的文件results.txt，其中打印所有结果。

于 2013-03-28T00:23:16.213 回答
0
于 2013-03-28T01:42:12.430 回答