Displaying 1 result from an estimated 1 matches for "gpen".
Did you mean:
gen
2010 Dec 13
0
SAS "/slice" and "/diff" equivalent? - Automated multiple comparisons in nlme package?
...l y = g h yr yr*h;
random b;
lsmeans g / diff;
The codes above give equivalent outputs (analysis of a 5(g) by 2(h)
factorial, RCBD(b)), except for the "lsmeans g / diff" line, which I
cannot reproduce in R!
The below gives one of the comparisons I am looking for.
anova(model, L=c("gPEN"=1, "gC5"=-1))
Is there a command to automate such thing? I understand from
https://stat.ethz.ch/pipermail/r-help/2007-March/128151.html
and
https://stat.ethz.ch/pipermail/r-help/2004-February/046684.html
the issues with "lsmeans" (this is *not* what I am trying to get h...