0

我已经设法从一组数据中制作了 2 的不同组合,以便我得到数据集中的名称对。我想在使用时有一个函数,mapply以便我可以使用每对中的每个名称来引用它们对应的数据集。现在我有:

myPairs <- combn(names(iris[1:4]), 2)

f <- function(x,y)
{
#Want to make a lm(x ~ y) and other potential calculations 
}

a <- myPairs[1,]
b <- myPairs[2,]

mapply(f, a, b)

换句话说,我想这样计算:

data <- iris
attach(data)
lm(Sepal.Length ~ Sepal.Width)
lm(Sepal.Length ~ Petal.Length)
lm(Sepal.Length ~ Petal.Width)
lm(Sepal.Width ~ Petal.Length)
lm(Petal.Length ~ Petal.Width)

但使用组合中的名称作为对数据集的引用

4

1 回答 1

2

首先,不要附加您的数据。使用data参数 forlm代替。

让我们构建公式:

myFormula = apply(myPairs, MARGIN = 2, FUN = paste, collapse = " ~ ")

然后制作模型

myModels = lapply(myFormula, function(x) lm(formula = x, data = iris))

然后我们可以检查一些结果:

myModels[[1]]
# Call:
# lm(formula = x, data = iris)
#
# Coefficients:
# (Intercept)  Sepal.Width  
#      6.5262      -0.2234  
于 2016-06-15T00:31:11.303 回答