kk1.m

来自「多元统计程序包」· M 代码 · 共 32 行

function [bzh,  remarix, dwhpcs, var, per_ex] = kk1(X)
[m,n]=size(X);
% 绘制第一个图形
subplot(1,2,1);
boxplot(X,0,'+',0); 
xlabel('变量值的范围');
ylabel('五个指标（变量）') ;

%计算样本的标准差
stdr=std(X,0,1);
%计算样本各列的均值
meanr=mean(X);
%标准化原始数据 
sr=(X-meanr(ones(m,1),:))./stdr(ones(m,1),:);
%样本相关矩阵
R=sr'*sr/(m-1);
 
%进行主成分分析
[pcs,newdata,variances,t2]=princomp(sr) ;
%方差贡献率
percent_explained=100*variances/sum(variances); 
 
%绘制第二个图形
subplot(1,2,2); 
pareto(percent_explained);
xlabel('主成分');
ylabel('方差占的比重(%)');
bzh = sr;
remarix = R;
dwhpcs = pcs;
var = variances;
per_ex = percent_explained;