matlab重复测量方差分析

MATLAB可以使用“anova1”和“anova2”函数进行重复测量方差分析。

1. “anova1”函数用于一元重复测量方差分析。语法如下：

anova1(data) % data为数据矩阵，每行表示一个被试的多次测量数据

2. “anova2”函数用于二元重复测量方差分析。语法如下：

anova2(data, factor1, factor2) % data为数据矩阵，每行表示一个被试的多次测量数据；factor1和factor2为两个自变量

例如，假设有5个被试，每个被试进行了3次测量，数据如下：

data = [3 4 3; 5 6 5; 2 3 2; 4 5 4; 1 2 1];

则可以使用“anova1”函数进行一元重复测量方差分析：

anova1(data)

结果如下：

ans =

struct with fields:

Source: 'Within Groups'
SSQD: 2.0000
DFGD: 8
MSQD: 0.2500
Fval: 1.0000
Prob: 0.4725
GrandMeanGroup: 3.0000
MeanVector: [3 4 3]
GrandMeanTotal: 3.2000
NGroupsTotal: 15

其中，“Source”表示变量来源，“SSQD”表示组内平方和，“DFGD”表示组内自由度，“MSQD”表示组内均方和，“Fval”表示F统计量，“Prob”表示p值，“GrandMeanGroup”表示各组平均值的均值，“MeanVector”表示每组的平均值，“GrandMeanTotal”表示所有测量数据的平均值，“NGroupsTotal”表示总测量数。

如果有两个自变量，可以使用“anova2”函数。例如，假设有5个被试，每个被试进行了3次测量，同时这些被试来自两个组，数据如下：

data = [3 4 3; 5 6 5; 2 3 2; 4 5 4; 1 2 1; 7 8 7; 9 10 9; 6 7 6; 8 9 8; 5 6 5];
factor1 = [1 1 1 1 1 2 2 2 2 2];
factor2 = [1 2 3 1 2 1 2 3 1 2];

则可以使用“anova2”函数进行二元重复测量方差分析：

anova2(data, factor1, factor2)

结果如下：

ans =

struct with fields:

Source: {3×3 cell}
SSQD: [2.0000 4.0000 2.0000; 0 0 0; 0 0 0]
DFGD: [1 1 1; 4 4 4; 4 4 4]
MSQD: [2.0000 4.0000 2.0000; 0 0 0; 0 0 0]
Fval: [1.0000 2.0000 1.0000; NaN NaN NaN; NaN NaN NaN]
Prob: [0.4725 0.1382 0.4725; NaN NaN NaN; NaN NaN NaN]
GrandMeanGroup: 4.0000
MeanVectorFactor1: [3.4000 6.8000]
MeanVectorFactor2: [3.2000 5.4000 2.8000]
GrandMeanTotal: 4.2000
NGroupsTotal: 30

其中，“Source”表示变量来源，有三个来源：Factor1、Factor2和Factor1:Factor2（即交互作用）；“SSQD”表示平方和，“DFGD”表示自由度，“MSQD”表示均方和，“Fval”表示F统计量，“Prob”表示p值，“GrandMeanGroup”表示各组平均值的均值，“MeanVectorFactor1”表示Factor1的各水平平均值，“MeanVectorFactor2”表示Factor2的各水平平均值，“GrandMeanTotal”表示所有测量数据的平均值，“NGroupsTotal”表示总测量数。