MANOVA.RM R 中的语法问题
MANOVA.RM syntax issues in R
我有一个由三个连续响应变量、两个分类变量和一个主题 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.RM 的帮助文件使用以下语法作为示例
EEG_MANOVA <- MANOVA(resp ~ sex * diagnosis,
data = EEG, subject = "id", resampling = "paramBS",
iter = 1000, alpha = 0.05)
然而,公式的左侧(即 'resp')仅指单个响应变量,而我想指的是变量矩阵(在本例中为 y1, y2, and y3
)
如何更改语法以包含三个响应变量?
如果您像这样绑定响应列(这在 manova 包中有效):
test <- manova(cbind(y1,y2,y3) ~ x1, data = data, subject = "id")
我有一个由三个连续响应变量、两个分类变量和一个主题 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.RM 的帮助文件使用以下语法作为示例
EEG_MANOVA <- MANOVA(resp ~ sex * diagnosis,
data = EEG, subject = "id", resampling = "paramBS",
iter = 1000, alpha = 0.05)
然而,公式的左侧(即 'resp')仅指单个响应变量,而我想指的是变量矩阵(在本例中为 y1, y2, and y3
)
如何更改语法以包含三个响应变量?
如果您像这样绑定响应列(这在 manova 包中有效):
test <- manova(cbind(y1,y2,y3) ~ x1, data = data, subject = "id")