1

我正在尝试在 R 中运行重复测量 MANOVA，它还包含许多因变量（行为任务的关键结果变量）。重复测量是由于交叉设计，其中个人服用药物和安慰剂（按随机顺序）。

我正在运行的代码如下所示：

imatrix <- matrix(c(
  1, 0, 0, 0, 0, 0, 1, 
  1, 0, 0, 0, 0, 0, -1, 
  0, 1, 0, 0, 0, 0, 1, 
  0, 1, 0, 0, 0, 0, -1,
  0, 0, 1, 0, 0, 0, 1, 
  0, 0, 1, 0, 0, 0, -1,
  0, 0, 0, 1, 0, 0, 1, 
  0, 0, 0, 1, 0, 0, -1,
  0, 0, 0, 0, 1, 0, 1, 
  0, 0, 0, 0, 1, 0, -1,
  0, 0, 0, 0, 0, 1, 1, 
  0, 0, 0, 0, 0, 1, -1
), 12, 7, byrow=TRUE)
colnames(imatrix) <- c("BCST", "CGT", "AST", "AGN", "DDT", "FERT", "NAC")
(imatrix <- list(measure=imatrix[,1:6], condition=imatrix[,7]))
contrasts(condition_factor) <- matrix(c(-1,1,1, -1), ncol=2)

doubly.mod<-lm(cbind(bcst_nac$totPersErr,bcst_placebo$totPersErr,cantab_nac$CGT.Delay.aversion,cantab_placebo$CGT.Delay.aversion,cantab_nac$AST.Switching.cost..Mean..correct.,cantab_placebo$AST.Switching.cost..Mean..correct.,cantab_nac$AGN.Affective.response.bias..Mean.,cantab_placebo$AGN.Affective.response.bias..Mean.,aucs_NAC,aucs_placebo,fert_nac$FERTACCHA,fert_placebo$FERTACCHA)~1))
Manova(doubly.mod, imatrix=imatrix, type =3)

结果是这个错误： Anova.III.mlm 中的错误（mod，SSPE，error.df，idata，idesign，icontrasts，：（列表）对象不能被强制输入'double'但是，当我将 imatrix 从列表到一个矩阵，我得到这个错误响应：

do.call(cbind, imatrix) 中的错误：第二个参数必须是一个列表

我基于 car::Anova 包中关于双多变量分析的示例。请让我知道您是否可以提供帮助，或者我是否可以添加任何内容以使这个问题更清楚。

4

0 回答 0