Hi,
Not sure this is the best place to ask since I'm fitting a fixed effect
model.. please advise if there is a better forum..
I'm trying to get lsmeans for a set of response variables, but I'm getting
this error and I can't figure out how to get it working.. it just spits out
the lsmeans for the first response.
Here is an example:

data(iris)
mydata = iris
# add another fixed effect
mydata$location = unlist(rep(list("US", "USSR"),
length(mydata$Sepal.Length)/2))
mydata$location = as.factor(mydata$location)
# run model with two response variables
modellist = lm(cbind(Sepal.Length,Sepal.Width) ~ Species*location,
data=mydata)
summary(modellist)
# try to lsmeans for each response
all.lsmeans = lsmeans(modellist, ~ Species*location)
summary(all.lsmeans)

Error in make.link(misc$tran) : 'cbind' link not recognised
 Species    location lsmean         SE  df lower.CL upper.CL
 setosa     US        4.252 0.07546007 144 4.102848 4.401152
 versicolor US        4.384 0.07546007 144 4.234848 4.533152
 virginica  US        4.720 0.07546007 144 4.570848 4.869152
 setosa     USSR      4.182 0.07546007 144 4.032848 4.331152
 versicolor USSR      4.322 0.07546007 144 4.172848 4.471152
 virginica  USSR      4.842 0.07546007 144 4.692848 4.991152

Thanks for any ideas!

	[[alternative HTML version deleted]]

