This won't make you an R aficionado, but depending on your needs
library(broom)
tidy(SumTukey)
might be useful. This converts the output to a familiar data.frame, making
it much easier to work with.
Best,
Ista
On Jun 24, 2016 9:48 AM, "John Sorkin" <jsor...@grecc.umaryland.edu> wrote:
>
> I am trying to make the leap from an R users to an R aficionado . . .
>
> I am trying to understand how add a column to the output of summary (and
> to understand how summary() works).
>
> I have run a glmer
> fit0 <- glmer(Fall ~
> Group+(1|PID),family=poisson(link="log"),data=data[data[,"Group"]!=0,])
>
> and I want to perform adjusted multiple comparisons:
>
> SumTukey <- summary(glht(fit0, linfct= mcp(Group="Tukey")))
>
> which gives beautiful output:
>
> > SumTukey
>
> Simultaneous Tests for General Linear Hypotheses
>
> Multiple Comparisons of Means: Tukey Contrasts
>
>
> Fit: glmer(formula = Fall ~ Group + (1 | PID), data = data[data[,
> "Group"] != 0, ], family = poisson(link = "log"))
>
> Linear Hypotheses:
> Estimate Std. Error z value Pr(>|z|)
> 2 - 1 == 0 0.5320 0.5075 1.048 0.717
> 3 - 1 == 0 0.6554 0.5000 1.311 0.551
> 4 - 1 == 0 0.9357 0.4655 2.010 0.181
> 3 - 2 == 0 0.1234 0.4174 0.296 0.991
> 4 - 2 == 0 0.4037 0.3754 1.075 0.700
> 4 - 3 == 0 0.2803 0.3651 0.768 0.867
> (Adjusted p values reported -- single-step method)
>
> I want to add a column to the output (unadjusted p-values), but I don't
> see how this might be done. The output
> is not a dataframe, nor is it a matix.
> > class(SumTukey)
> [1] "summary.glht" "glht"
>
> It is some class of objects that I don't understand and know nothing
> about. How can I add a column to the output of SumTukey [a.k.a.
> summary(glht(fit0, linfct= mcp(Group="Tukey")))] ?
>
> Thank you
> John
>
>
> John David Sorkin M.D., Ph.D.
> Professor of Medicine
> Chief, Biostatistics and Informatics
> University of Maryland School of Medicine Division of Gerontology and
> Geriatric Medicine
> Baltimore VA Medical Center
> 10 North Greene Street
> GRECC (BT/18/GR)
> Baltimore, MD 21201-1524
> (Phone) 410-605-7119
> (Fax) 410-605-7913 (Please call phone number above prior to faxing)
>
> Confidentiality Statement:
> This email message, including any attachments, is for ...{{dropped:16}}
______________________________________________
R-help@r-project.org mailing list -- To UNSUBSCRIBE and more, see
https://stat.ethz.ch/mailman/listinfo/r-help
PLEASE do read the posting guide http://www.R-project.org/posting-guide.html
and provide commented, minimal, self-contained, reproducible code.
