多指标单因素方差分析程序.doc

多指标单因素方差分析程序————————————————————————————————作者:————————————————————————————————日期: DimnindexAsLongDimindexNAME(20)AsStringDimindexMAX(20)AsLongDimindexORDER(20,20)AsStringDimnsampleAsLongDimindextable(20,500)AsIntegerDimcurrentROWAsLong'''迭代的最大次数ConstEULER==1E-30'''为防止除0使用的常数''',不能应用复杂统计程序求解,以减少计算时间及程序复杂性Subbasestat(index()AsInteger,X()AsDouble,ByValNAsInteger,ByValNKAsInteger,_AVEI()AsDouble,NI()AsInteger,SEI2()AsDouble,mini()AsDouble,maxi()AsDouble,_aveAsDouble,se2AsDouble,minAsDouble,maxAsDouble)DimIAsInteger,JAsIntegerReDimSUMi2(NK)AsDouble,SUMI(NK)AsDoubleDimSUMAsDouble,SUM2AsDoublemax=-1E-99'极值初值min=1E+100ForI=1ToNKmaxi(I)=-1E-99mini(I)=1E+100NI(I)=0NextForI=1ToN''基本数据计算SUM=SUM+X(I)SUM2=SUM2+X(I)*X(I)'MsgBox"index"&index(I)SUMI(index(I))=SUMI(index(I))+X(I)SUMi2(index(I))=SUMi2(index(I))+X(I)*X(I)NI(index(I))=NI(index(I))+1Ifmini(index(I))>X(I)Thenmini(index(I))=X(I)Ifmaxi(index(I))>X(I)Thenmaxi(index(I))=X(I)Ifmin>X(I)Thenmin=X(I)Ifmax>X(I)Thenmax=X(I)NextForI=1ToNK''二级数据计算'MsgBoxI&""&NI(I)AVEI(I)=SUMI(I)/NI(I)SEI2(I)=(SUMi2(I)-SUMI(I)*SUMI(I)/NI(I))/NI(I)Nextave=SUM/Nse2=(SUM2-SUM*SUM/N)/NEndSub'''()'''测试用20110307晚DimIAsInteger,JAsInteger,NAsIntegerDimNKAsInteger,nxyAsInteger:DimtoolAsNewstat_tool_OK_20120414DimanovaAsNewssanddf_A_R_okNK=8nxy=4N=244ReDima(N,3)AsString,X(N,11)AsDouble,indexa(N,N)AsInt