|
在数量生态学-R语言应用 第一二版中,有个做方差分析(ANOVA)的LSD多重比较的函数boxplert和boxplerk,此函数由Numerical Ecology with R 的原作者提供。这两个函数在做多重比较作图过程,字母标号出现错位的情况。我在讲课过程虽感觉用这两个函数的的结果有点奇怪,但是从来没有怀疑是原函数出错,还觉得是多重比较p值校正后引起的箱线图与字母标号不匹配的情况,直至有位学员发现此错位的问题并告知,我认真检查了函数才发现函数的bug后,并及时给原书作者反映。原书作者承认此错误是由于两个函数调用了agricolae包,但是并没有及时跟上agricolae的更改导致的。原书作者及时解决函数bug,而且还加了基于ggplot2的相同功能的函数,现在新的函数做出的多重比较的图是准确而且更美观。
所以,在此对于我的培训中"一元统计”部分的多重比较的代码应该重新更改如附件代码或是百度网盘
链接: https://pan.baidu.com/s/1z4CXjH_KyYoOZoxhrhhgbw 提取码: h5yx
现在正确的结果如下
boxplert的结果(双因素方差分析结果)
ggplert的结果(双因素方差分析结果)
Archiver|手机版|科学网 ( 京ICP备07017567号-12 )
GMT+8, 2024-9-23 11:18
Powered by ScienceNet.cn
Copyright © 2007- 中国科学报社