Matlab中多组(22,22)位置单因素方差分析(ANOVA1)实现求助

核心问题点 #

• 变量名冲突：循环变量使用了x，但此前已用x存储加载的数据集，这会直接覆盖原变量，导致循环执行到第二次就因变量类型错误中断。
• 未保存P值结果：p_all初始化为空矩阵，且循环中未将anova1返回的p值赋值给p_all的对应位置（注释掉的赋值语句正是关键步骤），最终无法生成目标P值矩阵。
• 内存未预分配（优化点）：动态扩展空矩阵会降低运行效率，建议提前预分配固定维度的矩阵。

修正后的代码 #

% 加载数据，重命名变量避免冲突
data = load("plv_8_12.mat");
data = data.plv;
size(data)

% 提取4个类的数据集，每个类72个样本
S1C1 = data(1:72,:,:);
S1C2 = data(649:720,:,:);
S1C3 = data(1297:1368,:,:);
S1C4 = data(1945:2016,:,:);

% 预分配22×22的P值矩阵内存
p_all = zeros(22,22);

% 遍历所有22×22的位置
for i = 1:22
    for j = 1:22
        % 提取对应位置的4组样本
        tc1 = S1C1(:,i,j);
        tc2 = S1C2(:,i,j);
        tc3 = S1C3(:,i,j);
        tc4 = S1C4(:,i,j);
        
        % 组织成anova1要求的格式：每列对应一个组，共4列72行
        temp = [tc1, tc2, tc3, tc4];
        
        % 执行单因素方差分析，仅保留P值并存入矩阵
        [p,~,~] = anova1(temp);
        p_all(i,j) = p;
    end
end

% 关闭所有anova1生成的统计图表（可选）
close all

修正说明 #

1. 将原变量x重命名为data，循环变量改用i/j，彻底解决变量冲突问题。
2. 预分配p_all矩阵，既提升运行效率，也确保矩阵维度符合预期。
3. 恢复并修正P值赋值语句，将每个位置的分析结果存入对应矩阵位置。
4. 使用~忽略不需要的返回值，简化代码结构。
5. 将close all移至循环结束后，避免重复关闭窗口影响运行效率。

内容的提问来源于stack exchange，提问作者Abdul