方差检验可以评估组间的差异。依据检验的结果,虽然你可以拒绝不存在差异的原假设,但方差检验并没有告诉你哪些组显著地与其他组有不同。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在书中给出的结果基本一致(不完全一致,是因为随机数的问题)

转载于:https://www.cnblogs.com/1zhk/p/4676624.html

R 中同步进行的多组比较的包:npmc相关推荐

  1. R语言转换并保存json文件--使用jsonlite包

    R语言转换并保存json文件--使用jsonlite包 钱亦欣 发表于 2017-07-03 00:10   2477 阅读 http://www.ituring.com.cn/article/468 ...

  2. r中如何求变量的对数转换_对数转换以求阳性。

    r中如何求变量的对数转换 In Simple terms, log transform squashes or compresses range of large numbers and expand ...

  3. 简单介绍一下R中的几种统计分布及常用模型

    统计学上分布有很多,在R中基本都有描述.因能力有限,我们就挑选几个常用的.比较重要的简单介绍一下每种分布的定义,公式,以及在R中的展示. 统计分布每一种分布有四个函数:d――density(密度函数) ...

  4. 《机器学习与数据科学(基于R的统计学习方法)》——2.11 R中的SQL等价表述...

    本节书摘来异步社区<机器学习与数据科学(基于R的统计学习方法)>一书中的第2章,第2.11节,作者:[美]Daniel D. Gutierrez(古铁雷斯),更多章节内容可以访问云栖社区& ...

  5. R语言笔记6:在R中写一些简单的函数、functions基础和作用域

    R语言基础系列: 1数据类型(向量.数组.矩阵. 列表和数据框) 2读写数据所需的主要函数.与外部环境交互 3数据筛选--提取对象的子集 4向量.矩阵的数学运算 5控制结构 Your first R ...

  6. 你知道R中的赋值符号箭头(-)和等号(=)的区别吗?

    作为一门高级语言,R语言拥有独特的语法,比如今天说道的赋值符号.在其他语言里,赋值符合通常用一个等号(=)表示,而在R语言里,承担这个任务的可以是箭头(<-)符号,也可以是等号(=).这就导致许 ...

  7. 你知道R中的赋值符号箭头和等号的区别吗?

    作者:陈亮 单位:中国科学院微生物研究所 你知道R中的赋值符号箭头(<-)和等号(=)的区别吗? 作为一门高级语言,R语言拥有独特的语法,比如今天说道的赋值符号.在其他语言里,赋值符合通常用一个 ...

  8. Co-occurrence网络图在R中的实现

    作者:陈亮 单位:中科院微生物所 编者按:上个月菌群月坛,在军科院听取王军组陈亮博士分享网络分析的经验,不仅使我对网络的背景知识有了更全面的认识,更使我手上一个关于菌根的课题有极大的启示.这么好的知识 ...

  9. 独家 | 在R中使用LIME解释机器学习模型

    作者:PURVAHUILGOL 翻译:陈丹 校对:欧阳锦 本文约3200字,建议阅读15分钟 本文为大家介绍如何在R中使用LIME来解释机器学习模型,并提供了相关代码. 关键词:机器学习模型解释.R语 ...

最新文章

  1. 你不知道的18个Python高效编程技巧
  2. runtime无法执行grep_让你的 Shell 命令执行可视化和告警
  3. 生成Base58格式的UUID(Hibernate Base64格式的UUID续)
  4. java当前时间转化毫秒_Java学习笔记-时间相关类
  5. error LNK2038: 检测到“RuntimeLibrary”的不匹配项: 值“MDd_DynamicDebug”不匹配值“MTd_StaticDebug”...
  6. 基于spring-redis发布订阅模式的实现
  7. 通过伙伴系统申请内核内存的函数有哪些?
  8. 如何使用JSZip(How to use JSZip)
  9. Hi3519A调试记录
  10. 理论知识:电感感应电压公式,纯电感正弦电路感应电压公式,纯电容正弦电路电流公式的推导
  11. 数据库---[复习2]---数据查询---设有一个SPJ数据库,包括S、P、J及SPJ4个关系模式··· ···
  12. Python中的立方根
  13. 前端项目发版后用户需要手动刷新页面获取的解决方案
  14. 致远项目管理SPM系统案例:中关村创客小镇数据采集云管理
  15. 电子元器件B2B商城系统开发:赋能企业构建进销存标准化流程实例
  16. 小明左、右手中分别拿两张纸牌(比如:黑桃10和红桃8,数字10和8可通过键盘录入),要求编写代码交换小明手中的牌
  17. KubeEdge环境搭建(支持网络插件flannel)
  18. 【学习笔记+习题集】字符相关(输入输出流,字典树,AC自动机,后缀自动机)(4598字)(更新至2022.12.28)
  19. Java_Map_Map详解
  20. windows server 2016磁盘安全与管理_磁盘管理工具哪一款好用?

热门文章

  1. android 沙箱 逆向,【原创】沙箱Sandboxie v3.40 逆向完整源码
  2. 搭建oracle运行环境,搭建 Oracle Warehouse Builder 11
  3. pbl和sbl_PBL与SBL联合教学模式在妇产科教学中的应用研究
  4. 塔菲尔曲线斜率的大小_中国第一塔,与埃菲尔铁塔齐名,到底有多强悍?
  5. 对比原生Node封装的Express路由 和 express框架路由
  6. php7 电子书 下载,PHP 7: Real World Application Development
  7. linux压缩文件命令_Linux基础篇(二)--Linux常用命令
  8. 调整样式_“寒来暑往”美国海军陆战队应季节调整迷彩样式的规定变化
  9. java如何对foo bar调用方法_关于java:如何测试工厂方法传递给构造函数的参数?...
  10. 中兴手机数据通道打不开_换了个新手机