如何使用存储在向量中的列名调用 cbind 的列（尝试在 R 中使用 manova 对列进行成对比较）

Question

我正在尝试执行成对 manova 分析，在其中循环遍历所有可能的列对。我认为最好通过示例进行交流：

varList <- colnames(iris)

m1 <- manova(cbind(varList[1], varList[2]) ~ Species, data = iris)
# Error in model.frame.default(formula = cbind(varList[1], varList[2]) ~  :
# variable lengths differ (found for 'Species')

m2 <- manova(cbind(noquote(varList[1]), noquote(varList[2])) ~ Species,
             data = iris)

# Error in model.frame.default(formula = cbind(noquote(varList[1]), noquote(varList[2])) ~  : 
# variable lengths differ (found for 'Species')

m3 <- manova(cbind(Sepal.Length, Petal.Length) ~ Species, data = iris)
m4 <- manova(cbind(iris[ ,1], iris[ ,3]) ~ Species, data = iris)

summary(m3)

#            Df Pillai approx F num Df den Df    Pr(>F)    
# Species     2 0.9885   71.829      4    294 < 2.2e-16 ***
# Residuals 147                                            
# ---
# Signif. codes:  0 ‘***’ 0.001 ‘**’ 0.01 ‘*’ 0.05 ‘.’ 0.1 ‘ ’ 1

R.version.string
# [1] "R version 3.4.2 (2017-09-28)"

RStudio.Version()$version
# [1] ‘1.1.383’

我认为这与在我的 cbind() 函数中从向量引用 colnames 更相关。我从这个问题中看到了有关使用括号的一些信息，但无法让它对我的情况起作用。我可以按编号调用列（请参阅 m4），但如果可能，我更愿意使用列名。

Answer 1

您需要包装您使用 eval(as.symbol()) 调用的向量中的每个条目。

所以： m1 <- manova(cbind(eval(as.symbol(varList[1])), eval(as.symbol(varList[2]))) ~ Species, data = iris) 应该可以。

如何使用存储在向量中的列名调用 cbind 的列（尝试在 R 中使用 manova 对列进行成对比较）

How to call columns for cbind using column names stored in a vector (attempting to perform pairwise comparison of columns with manova in R)

r

dataframe

cbind

manova