在数量生态学-R语言应用 第一二版中，有个做方差分析（ANOVA)的LSD多重比较的函数boxplert和boxplerk，此函数由Numerical Ecology with R 的原作者提供。这两个函数在做多重比较作图过程，字母标号出现错位的情况。我在讲课过程虽感觉用这两个函数的的结果有点奇怪，但是从来没有怀疑是原函数出错，还觉得是多重比较p值校正后引起的箱线图与字母标号不匹配的情况，直至有位学员发现此错位的问题并告知，我认真检查了函数才发现函数的bug后，并及时给原书作者反映。原书作者承认此错误是由于两个函数调用了agricolae包，但是并没有及时跟上agricolae的更改导致的。原书作者及时解决函数bug，而且还加了基于ggplot2的相同功能的函数，现在新的函数做出的多重比较的图是准确而且更美观。

所以，在此对于我的培训中"一元统计”部分的多重比较的代码应该重新更改如附件代码或是百度网盘

多重比较代码更新.zip

链接: https://pan.baidu.com/s/1z4CXjH_KyYoOZoxhrhhgbw 提取码: h5yx 

现在正确的结果如下

boxplert的结果（双因素方差分析结果）

ggplert的结果（双因素方差分析结果）