私はRで助けを求めて書いています。
次のスクリプトを使用して簡単なRCBD分析を行い、特性「X」の遺伝子型(名前)を比較しています。
library(stats)
data_1=read.table(file="test.txt", head=TRUE)
result_X =aov(X~Block+Name, data=data_1)
sink("result_X.txt")
summary(result_X)
sink()
データに欠測データ(「NA」)があります。そこで、LSDを計算した後、遺伝子型を降順で比較したいと思います。一部の「名前」の「ブロック」が欠落しているため、平均は適切ではないと思います。
したがって、問題は、「最小二乗平均」を出力するためのスクリプトは何であるかということです。これは、単純な平均よりもLSDと比較するのに最適だと思います。ご協力ありがとうございました、
オズワルド