As I was suspecting, the answer is in contrast::contrast(). I'm going to
take this over to r-sig-mixed-models since that mailing list might be more
appropriate than ecology. Thanks Thierry, I appreciate the help.
Cheers
MVS
=
Matthew Van Scoyoc
https://sites.google.com/site/s
tely I would like to create a dataframe so I can plot the contrasts,
something like this...
> x = summary(glht(math.lmm, linfct = cc$X))
> # Contrast data frame
> math.contr = data.frame(Effect.Interaction = ..., Estimate =
x[["test"]]$coefficients, Std.Error = x[["test&quo
MVS
=
Matthew Van Scoyoc
https://sites.google.com/site/scoyoc/
=
Think SNOW!
On Sun, Nov 30, 2014 at 11:44 PM, Chris Howden wrote:
> Hi Mathew,
>
> Are you sure you want to 'weight' by elevation? That would imply higher
> elevations have a greater weight on the
tes = lsms[[1]]$Estimate
At this point I would be plotting the estimates of the significant response
variables and interactions to look at the differences. As I understand, the
ecosystem effect and the interaction between ecosystem and quality are with
an average elevation for all plots, and this could
Scoyoc
Graduate Research Assistant, Ecology
Wildland Resources Department & Ecology Center
Quinney College of Natural Resources
Utah State University
Logan, UT
=
Think SNOW!
--
View this message in context:
http://r-sig-ecology.471788.n2.nabble.com/Identification-of-outliers-using-NMDS-good