%计算累积频率Рp(end)=p(end)-1/(4*n); %修正最后一个累积频率РProbit=norminv(p,0,1)+5; %计算标准正态分布的p分位数+5РProbit=[ones(n,1), Probit, Probit.^2, Probit.^3];Р[b,bint,r,rint,stats]=regress(WRSR,Probit); Р%三次多项式回归, 一次或二次多项式回归R方和p值较差Рb РstatsРWRSRfit=Probit*b %计算WRSR的估计值Р[s,ind]=sort(WRSRfit,'descend') %对WRSR估计值从大到小排序Р运行结果：Рb = 10.7950 -6.4953 1.3137 -0.0851РРstats = 0.7932 7.6735 0.0178 0.0066РРWRSRfit = 0.4335Р 0.3851Р 0.4133Р 0.4655Р 0.5280Р 0.5953Р 0.6634Р 0.7243Р 0.7463Р 0.5454Рs = 0.7463Р 0.7243Р 0.6634Р 0.5953Р 0.5454Р 0.5280Р 0.4655Р 0.4335Р 0.4133Р 0.3851РРind = 9 8 7 6 10 5 4 1 3 2Р说明：回归方程为РWRSR=10.795-6.4953* Probit +1.3137* Probit2-0.0851* Probit3РРР 评价从高到低，依次是：第9样本（1991年，），第8样本(1990年)，……
全文预览
