Skip to content

lsmeans for multiple response variables

2 messages · Paul Tanger, Ben Bolker

#
On 14-03-30 12:47 AM, Paul Tanger wrote:
The short answers are that (1) this is indeed not a mixed-model question;
try r-help at r-project.org or Stack Overflow (2) you may not be aware that
you're fitting two separate linear models here --  from ?lm

   If ?response? is a matrix a linear model is fitted separately by
     least-squares to each column of the matrix.

This is a relatively less-used feature, so presumably lsmeans()
doesn't expect it.  You could fit the two models separately, or
use lapply() ... or request of the lsmeans maintainer that they
change lsmeans to allow such support.

  Ben Bolker