背景:
使用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,即认为存在第一种生活方式的减肥效果弱于第二种生活方式,第二种生活方式的减肥效果弱于第三种生活方式。