方差檢驗可以評估組間的差異。依據檢驗的結果,雖然你可以拒絕不存在差異的原假設,但方差檢驗並沒有告訴你哪些組顯著地與其他組有不同。Robert 在 《R in Action》一書中推薦了一個包-npmc: 該包提供了一種非參數多組比較程序。在控制犯第一類錯誤的概率(發現一個事實上並不存在的差異的概率)的前提下,執行可以同步進行的多組比較,這樣可以直接完成所有組之間的成對比較.
由於該包所依賴的mvtnorm包的一些更新,npmc包的最后更新版本(1.0.7)無法正確運行在R 3.0以上的版本。我在Rstudio中做了調試和編譯,對npmc.R的原代碼做了一些細微的更改以支持最新的mvtnorm包。npmc包最后編譯運行成功。
最新原程序共享在GitHub:
https://github.com/HK-Zhang/Rice/tree/master/npmc
也可以通過下面鏈接直接下載這個包
http://pan.baidu.com/s/1nuqdXcX
測試代碼和結果:
> class<-state.region > var<-state.x77[,c('Illiteracy')] > mydata<-as.data.frame(cbind(class,var)) > rm(class,var) > library(npmc) > summary(npmc(mydata),type='BF')
結果如下,和Robert在書中給出的結果基本一致(不完全一致,是因為隨機數的問題)