Error in contrasts adding bonferroni adjustment

classic Classic list List threaded Threaded
2 messages Options
Reply | Threaded
Open this post in threaded view
|

Error in contrasts adding bonferroni adjustment

Laura Coco
Hello

I am trying to add a Bonferroni adjustment to two comparisons, and I'm
getting an error msg.

The contrasts run just fine; I used a contrast matrix with dummy
variables: Experimental
vs Control at post1 (first contrast) and Experimental vs Control at post2
(second contrast).

contrast.outcome <-contrast(emmeans.outcome, method =
list("experimental_post1 - control_post1" = c(1, 0, -1, 0),
"experimental_post2 - control_post2" = c(0, 1, 0, -1), adjust = "bonf"))

The code runs except when I add the adjust = "bonf" - The error is "Error
in tcmat %*% linfct : requires numeric/complex matrix/vector arguments"

I see many tutorials online about adjusting in pairwise comparisons and
ANOVA but not in this situation. Any guidance you can provide would be much
appreciated.

Thank you!

Laura

        [[alternative HTML version deleted]]

______________________________________________
[hidden email] 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.
Reply | Threaded
Open this post in threaded view
|

Re: Error in contrasts adding bonferroni adjustment

Patrick (Malone Quantitative)
Looks like a misplaced right parenthesis. Try

contrast.outcome <-contrast(emmeans.outcome, method =
list("experimental_post1 - control_post1" = c(1, 0, -1, 0),
"experimental_post2 - control_post2" = c(0, 1, 0, -1)), adjust = "bonf")

Pat


On Sat, Jan 9, 2021 at 3:10 AM Laura Coco <[hidden email]> wrote:

>
> Hello
>
> I am trying to add a Bonferroni adjustment to two comparisons, and I'm
> getting an error msg.
>
> The contrasts run just fine; I used a contrast matrix with dummy
> variables: Experimental
> vs Control at post1 (first contrast) and Experimental vs Control at post2
> (second contrast).
>
> contrast.outcome <-contrast(emmeans.outcome, method =
> list("experimental_post1 - control_post1" = c(1, 0, -1, 0),
> "experimental_post2 - control_post2" = c(0, 1, 0, -1), adjust = "bonf"))
>
> The code runs except when I add the adjust = "bonf" - The error is "Error
> in tcmat %*% linfct : requires numeric/complex matrix/vector arguments"
>
> I see many tutorials online about adjusting in pairwise comparisons and
> ANOVA but not in this situation. Any guidance you can provide would be much
> appreciated.
>
> Thank you!
>
> Laura
>
>         [[alternative HTML version deleted]]
>
> ______________________________________________
> [hidden email] 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.



--
Patrick S. Malone, Ph.D., Malone Quantitative
NEW Service Models: http://malonequantitative.com

He/Him/His

______________________________________________
[hidden email] 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.