matlab anova2

在MATLAB中进行双因素方差分析（Two-Way ANOVA），可以使用函数`anova2`。该函数的语法如下：

[p, tbl, stats] = anova2(y, group, replicates, displayopt)

其中，`y`是一个包含所有观测值的向量或矩阵，`group`是一个包含两个因素水平的矩阵，`replicates`是一个包含每个组合下的重复次数的矩阵，`displayopt`是一个控制输出结果格式的字符串。

例如，假设我们有一个数据集`data`，包含两个因素`A`和`B`，每个因素各有三个水平，每个组合下有5个重复观测值，可以使用以下代码进行方差分析：

A = repmat([1 2 3], 1, 9)';
B = repmat([1 2 3]', 3, 1);
data = [ 55    76    78    90    95    94    90    88    85
         66    79    80    91    98    97    92    89    86
         64    77    81    92    99    98    93    90    87
         68    81    82    93   100    99    94    91    88
         70    83    84    95   102   101    96    93    90];
group = [A B];
replicates = 5;
[p,tbl,stats] = anova2(data, group, replicates)

`p`为一个包含三个元素的向量，分别表示因素A、因素B和交互作用的显著性水平。`tbl`是一个包含方差分析表的结构体，其中包括每个因素和交互作用的均方、自由度和F值等信息。`stats`是一个包含组间均值、标准误差和置信区间等信息的结构体。