R语言 为进行JT检验,安装SAGx包时报错当前版本没有适合的包的解决方法

背景:

使用R语言进行JT检验时,教材上提示需要安装相应的包SAGx.

安装时报错:

Warning in install.packages :
  package ‘SAGx’ is not available for this version of R

解决:

后续在网上找到了两种解决方法:

一是这一篇文章方差分析Dunn函数&JT检验_dunn检验-优快云博客,是博主“leechee___”自己写的JT检验函数

二是在网上找到了一个Github项目GitHub - cran/SAGx: Statistical Analysis of the GeneChip,里面包含有JT检验的r语言程序(对应文件名为JT.test.r)。

运用:

示例:运行JT.test.r内包含的代码后,进行教材144页4.5题的JT检验:

x = c(3.7, 3.7, 3.0, 3.9, 2.7)
y = c(7.3, 5.2, 5.3, 5.7, 6.5)
z = c(9.0, 4.9, 7.1, 8.7)
g = c(rep(c(1,2,3), c(5,5,4)))
tapply(c(x, y, z), g, median)
JT.test(data = t(c(x, y, z)), class = g)

运行结果:

     trend p-value   1   2   3 rank correlation
[1,]   0.001913353 3.7 5.7 7.9        0.8185622

p值小于0.05,拒绝H0,即认为存在第一种生活方式的减肥效果弱于第二种生活方式,第二种生活方式的减肥效果弱于第三种生活方式。

评论
添加红包

请填写红包祝福语或标题

红包个数最小为10个

红包金额最低5元

当前余额3.43前往充值 >
需支付:10.00
成就一亿技术人!
领取后你会自动成为博主和红包主的粉丝 规则
hope_wisdom
发出的红包
实付
使用余额支付
点击重新获取
扫码支付
钱包余额 0

抵扣说明:

1.余额是钱包充值的虚拟货币,按照1:1的比例进行支付金额的抵扣。
2.余额无法直接购买下载,可以购买VIP、付费专栏及课程。

余额充值