你试一下library(multcomp)
用simint()
这有一个巨详细的例子：
http://math.furman.edu/~dcs/courses/math47/R/library/multcomp/doc/Rmc.pdf

金币虽不多，还请大家帮忙啊

asreml是计算混线性模型的，一般数据不平衡是用的较多。
多重比较，比如lsd，用的是sed*t，前提是sed是一样的，但是不平衡数据中，sed是不同的。
因此，可以用mean sed *t(0.05,df)，来计算lsd，然后用lsd进行多重比较。