MATLAB实战——方差分析

例1

例2

例3

代码

% 读取数据
[x,y] = xlsread('examp8_1_1.xls');

score = x(:,2);   % 提取得分
college = y(2:end,4);   % 所在学院

% 正态性检验
for i = 1:4
    scorei = score(x(:,1)==i);
    [h(i),p(i)] = lillietest(scorei);
end

% 方差齐性检验
[p,stats] = vartestn(score,college)

% 单因素一元方差分析
[p,table,stats] = anova1(score,college);
c = multcompare(stats);

% X = [125 60 338 210;
%     63 51 260 203;
%     130 65 403 205;
%     94 33 260 280;
%     55 40 390 295;
%     120 56 416 280;
%     114 63 395 380;
%     110 69 377 260;
%     60 57 273 260];
% 
% data = xlsread('例2.xlsx');
% X = data(:,2:end);
% group = num2str(data(:,1));
% [d,p,stats] = manova1(X,group)    % d=1，不能拒绝共线性，d=2则存在一个共面关系


% [x,y] = xlsread('例3.xlsx');
% num = length(y(2:end,1));
% x2 = x(:);
% group = num2str([ones(num,1);2*ones(num,1);3*ones(num,1);4*ones(num,1)]);
% [p,table,stats] = anova1(x2,group);
% c = multcompare(stats);