R 中同步進行的多組比較的包:npmc


方差檢驗可以評估組間的差異。依據檢驗的結果,雖然你可以拒絕不存在差異的原假設,但方差檢驗並沒有告訴你哪些組顯著地與其他組有不同。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在書中給出的結果基本一致(不完全一致,是因為隨機數的問題)


免責聲明!

本站轉載的文章為個人學習借鑒使用,本站對版權不負任何法律責任。如果侵犯了您的隱私權益,請聯系本站郵箱yoyou2525@163.com刪除。



 
粵ICP備18138465號   © 2018-2025 CODEPRJ.COM