sta12_2.m

几个用于统计方法的matlab程序

X=[13.0 9.7 1.5 6.4
   1.0 7.5 1.5 6.5
   20.6 12.5 2.3 7.0
   33.3 19.0 2.8 5.8
   20.5 14.2 1.9 6.9
   10.0 6.7 2.2 7.0
   12.7 5.7 2.9 6.7
   36.5 15.7 2.3 7.2
   37.1 14.3 2.1 7.2
   25.5 12.9 1.9 7.3
   26.5 14.9 2.4 6.7
   22.3 8.4 4.0 7.0
   30.8 7.4 2.7 6.4
   25.3 7.0 4.8 7.3
   31.2 11.6 2.4 6.3
   22.7 10.1 3.3 6.2
   31.2 9.6 2.4 6.0
   13.2 6.6 2.0 5.8
   11.1 6.7 2.2 7.2
   20.7 9.6 3.1 5.9];
[pcs,newdata,variances,t2]=princomp(X)
percent_explained=100*variances/sum(variances);
pareto(percent_explained)
xlabel('Principal Component')
ylabel('Variance Explained(%)')