Segue um script que pode ser útil.

http://www.leg.ufpr.br/~walmes/cursoR/cnpaf2/bibi.R

Onde estiver popMeans() e popMatrix(), troque por LSmeans() e LSmatrix().

À disposição.

Walmes.