MANOVA 中的交互项与 R

如何解决MANOVA 中的交互项与 R

我想根据 R 中交互项的多变量分析计算 Pillai's Trace。

我的数据的一个子集是：

structure(list(id = structure(c(1L,3L,4L,5L,6L,7L,8L,10L,11L,12L,13L,14L,15L,16L,17L,18L,19L,20L,21L,22L,23L,24L,25L,26L,27L,28L,29L,30L,31L,32L,33L,34L,36L,37L,38L,39L),.Label = c("5356","5357","5358","5359","5360","5363","5373","5381","5383","5386","5395","5397","5400","5401","5444","5445","5446","5448","5450","5451","5454","5472","5473","5475","5476","5477","5478","5480","5481","5483","5487","5494","5495","5504","5505","5506","5507","5508","5509","5513","5514","5515","5516","5517","5518","5519","5521","5523","5524","5526","5527","5528","5544","5545","5546","5547","5551","5552","5553","5554","5555","5558","5559","5560","5562","5564","5566","5573","5574","5575","5576","5577","5578","5579","5584","5585","5587","5588","5589","5591","5594","5595","5604","5611","5612","5613","5615","5616","5619","5620","5621","5622","5626","5627","5628","5631","5632","5634","5635","5643","5652","5653","5654","5655","5656","5657","5659","5660","5661","5664","5665","5666","5669","5671","5672","5673","5678","5680","5688","5689","5690","5691","5692","5698","5699","5700","5702","5703","5704","5706","5708","5709","5710","5730","5731","5732","5733","5734","5735","5739","5740","5741","5742","5743","5744","5745","5746","5747","5748","5749","5750","5753","5754","5755","5766","5767","5776","5777","5778","5779","5780","5781","5787","5788","5789","5790","5791","5792","5793","5797","5798","5799","5800","5801","5810","5811","5812","5813","5814","5819","5820","5821","5822","5823","5824","5825","5827","5828","5829","5830","5857","5859","5874","5875","5876","5877","5878","5879","5883","5884","5886","5887","5888","5889","5890","5892","5893","5896","5899","5900","5909","5910","5918","5919","5920","5921","5922","5923","5927","5929","5931","5932","5933","5934","5936","5937","5941","5943","5944","5949","5950","5951","5952","5956","5957","5958","5959","5971","5972","5973","5976","5979","5980","5981","6001","6002","6003","6004","6005","6009","6027","6028","6033","6042","6054","6063","6067","6073","6076","6077","6078","6079","6080","6081","6082","6083","6098","6102","6103","6104","6105","6106","6107","6111","6119","6133","6146","6147","6157","6158","6160","6161","6162","6163","6164","6165","6166","6167","6168","6169","6170","6171","6172","6173","6174","6175","6190","6193","6195","6196","6197","6208","6228","6229","6232","6255","6268","6269","6270","6275"),class = "factor"),Age = structure(c(2L,2L,1L,2L
),.Label = c("Young","Old"),Gender = structure(c(1L,1L),.Label = c("Male","Female"),ROI_Coded = structure(c(2L,2L),.Label = c("LC","Noradrenergic","FPT"),leftright = structure(c(1L,.Label = c("left","right"),FA = c(0.335363597,0.390655279,0.402167588,0.372361958,0.411948681,0.454060167,0.416042715,0.356336266,0.380758196,0.363425195,0.435529858,0.369446278,0.357858658,0.434102625,0.410618424,0.409390539,0.435052991,0.404053301,0.445090771,0.430938751,0.427978307,0.346751571,0.394940138,0.392188877,0.382936746,0.397471696,0.412624329,0.363124639,0.372451693,0.377948016,0.374344319,0.378421336,0.367845714,0.358090758,0.409444004,0.435775071)),row.names = c(NA,36L),class = "data.frame")

子集的样子：

     id   Age Gender     ROI_Coded leftright        FA
1  5356   Old   Male Noradrenergic      left 0.3353636
2  5358   Old Female Noradrenergic      left 0.3906553
3  5359   Old   Male Noradrenergic      left 0.4021676
4  5360   Old   Male Noradrenergic      left 0.3723620
5  5363   Old   Male Noradrenergic      left 0.4119487

我可以像这样运行 MANOVA：

manova(cbind(Age,ROI_Coded,leftright,Gender) ~ FA,data = df3) %>% 
  summary()

但我想要 Age*ROI_Coded*leftright 的交互，但不知道如何将其适当地输入到模型中。它似乎不像线性模型那样采用交互术语。