《非参数统计》与MATLAB编程 第六章 多样本问题

第六章 多样本问题§6.1Kruskal-Wallis检验p=kruskalwallis(X)p=kruskalwallis(X,group)p=kruskalwallis(X,group,'on')%显示方差分析表和箱线图p=kruskalwallis(X,group,'off')%不显示方差分析表和箱线图[p,table]=kruskalwallis(...)[p,table,stats]=kruskalwallis(...)详细解说请看MATLAB的帮助文件score=[8472759572907575658055956958788062657242];college—{a,a,a,a,a,a,a,b,b,b,b,b‘ >i> '一' '一' '一' '一' '一' '一' '一' 〕_,b,c,c,c,c,c,c,c,}；[p,anvotab,stats]—kruskalwallis(score,college)p—0.1294anvotab—,Source,,SS,,df,,MS’,Chi-sq,，Prob>Chi-sq，,Groups,[141.9821][2][70.9911][4.0905][ 0.1294]'Error'[517.5179][17][30.4422][][],Total'[659.5000][19][][][]stats—gnames:{3x1cell}n:[767]source:kruskalwallismeanranks:[13.785710.25007.4286]sumt:66Clickonthegroupyouwanttotest2 4 6 8 10 12 14 16 18Nogroupshavemeanrankssignificantlydifferentfromc并得到方差分析表和下面的箱线图:SourceSSdfMSChi-sqProb>Chi-sqGroups141.98270.9914.090.1294Error517.521730.442Total659.519ABC60ABC60[p,table,stats]=kruskalwallis(score,college,'off')p=0.1294table='Source''SS''df''MS''Chi-sq''Prob>Chi-sq''Groups'[141.9821][2][70.9911][4.0905][ 0.1294]'Error'[517.5179][17][30.4422][][]'Total'[659.5000][19][][][]stats=gnames:{3x1cell}n:[767]source:'kruskalwallis'meanranks:[13.785710.25007.4286]sumt:66此题有结，这里是按结的处理方法对待。对于普通的单因素方差分析：[p,table,stats]=anova1(score,college,'off')0.0928table='Source''SS''df''MS''F''Prob>F',Groups,[802.9738][2][401.4869][2.7433][0.0928]'Error'[2.4880e+003][17][146.3515][][]'Total'[3.2909e+003][19][][][]stats=gnames:{3x1cell}n:[767]source:'anova1'means:[80.428673.166765.2857]df:17s:12.0976§6.2han=[260200240170270205190200250200310310190225170210280210280240320260360310270380240295260250];zu={，，a,'a','a'，，,a,'一,'一',一,a,a,a,'a','a','a','b','b','b','b','b','b','b','b','b','b','c','c','c','c','c','c','c','c','c','c'};[p,table,stats]=kruskalwallis(han,zu)p=0.0062table=,Source' 'SS' 'df' 'MS' 'Chi-sq' 'Prob>Chi-sq',Groups, [ 783.8000] [2] [391.9000] [10.1633] [ 0.0062]'Error' [1.4527e+003] [27] [53.8037] :] :]'Total' [2.2365e+003] [29] :] :] :]stats=gnames:{3x1cell}n:[101010]source:'kruskalwallis'meanranks:[9.800014.500022.2000]sumt:132T=stats.meanranks*10*[10;30;50]T=16430N=30ET=N"2*(N+1)/2ET=13950g=11其中长度为2的结有6个，长度为3的结有4个，DT=(N*(N”2-1)-6*6-4*24)*([103050].”2*[10;1