R中的MANOVA.RM语法问题

时间:2017-03-13 05:21:19

标签: r multivariate-testing manova

我有一个由三个连续响应变量,两个分类变量和一个主题ID组成的数据框

library(MANOVA.RM)
y1<-sample(1:150, 100, replace=T)
y2<-sample(1:150, 100, replace=T)
y3<-sample(1:150, 100, replace=T)
x1<-sample(x=c("S1", "S2"), size=100, replace=T, prob=rep(1/2,2))
x2<-sample(x=c("T1", "T2"), size=100, replace=T, prob=rep(1/2,2))
id <- seq (from=1, to=100, by=1)

data <- data.frame(y1, y2, y3, x1, x2, id)

我正在尝试运行MANOVA但是在分配响应变量方面存在语法问题。

MANOVA.RM的帮助文件使用以下语法作为示例

EEG_MANOVA <- MANOVA(resp ~ sex * diagnosis, data = EEG, subject = "id", resampling = "paramBS", iter = 1000, alpha = 0.05)

然而,公式的左侧(即&#39; resp&#39;)仅指单个响应变量,而我想引用变量矩阵(在本例中为y1, y2, and y3

如何更改语法以包含三个响应变量?

1 个答案:

答案 0 :(得分:0)

如果您绑定响应列(例如,在man man包中),它应该可以工作:

leg.get_frame().set_linewidth(7.0)