如何解决如何遍历列表并根据 LSD/HSD 值分配字母/组
我有一组 14 个番茄品种的产量数据,每个品种 2 次。我使用空间分析 (SpAT) 工具创建了一个混合模型(产量 ~ 品种 * 行 * 列)来估计每个基因型的拟合产量值,现在我有了这个拟合产量值 + 标准误差。
我的问题是如何计算每个基因型的平均值之间的显着差异,并相应地对组进行分组。
我已经估计了 HSD 的值(基于此方法 https://rpubs.com/aaronsc32/post-hoc-analysis-tukey)。但是现在,我正在寻找代码来循环遍历 14 个基因型的列表,比较它们的均值,并相应地对组进行分配。请帮忙!
我在 RStudio 工作。
这是一些模拟数据:
genotypes = c("gen1","gen2","gen3","gen4","gen5","gen6","gen7","gen8","gen9","gen10","gen11","gen12","gen13","gen14")
yields = c(348.07811,67.05801,136.71615,212.15367,115.54952,332.44170,113.59584,261.3755,122.36931,199.45903,140.82236,392.07310,178.56068,263.57455)
standarderrors= c(22.41150,23.09072,26.22475,22.69071,25.49619,25.45282,23.64242,25.65258,26.70335,26.35686,32.14079,51.14266,30.07858,26.18179)
mydata= cbind(genotypes,yields,standarderrors)
HSD.value = 50
版权声明:本文内容由互联网用户自发贡献,该文观点与技术仅代表作者本人。本站仅提供信息存储空间服务,不拥有所有权,不承担相关法律责任。如发现本站有涉嫌侵权/违法违规的内容, 请发送邮件至 dio@foxmail.com 举报,一经查实,本站将立刻删除。