Boxplot con pvalues obtenidos con emmeans

Boxplot con los valores de p obtenidos con emmeans de un LM, LMM, GLM, GLMM, etc. Gráficos generados con paqueteria "ggplot", por ahora solo permite una variable explicativa sin interacciones.

Usage

boxplot_emmeans(formula, modelo, data, p.adjust.method = NULL)

Arguments

formula

Formula de los contrastes que nos interesan obtener.

modelo

Modelo LM, GLM, GLMM, etc.

data

Base de datos con que se ajusto el modelo.

p.adjust.method

método para ajustar los valores p para comparaciones múltiples. Se utiliza cuando se realizan comparaciones por pares. Los valores permitidos incluyen "holm", "hochberg", "hommel", "bonferroni", "BH", "BY", "fdr", "none". Si no desea ajustar el valor p (no recomendado), utilice p.adjust.method = "none".

Value

Boxplot de los contrastes posthoc en el modelo todas las variables explicativas del modelo.

Examples

library(tlamatini)
data(iris)
mod <- glm(Petal.Width ~ Species, family = gaussian("log"), data=iris)
boxplot_emmeans(formula = Petal.Width ~ Species, modelo=mod, data= iris)

#> # A tibble: 3 × 7
#>   group1     group2        df statistic        p    p.adj p.adj.signif
#>   <chr>      <chr>      <dbl>     <dbl>    <dbl>    <dbl> <chr>       
#> 1 setosa     versicolor   147     -14.1 4.84e-29 1.45e-28 ****        
#> 2 setosa     virginica    147     -17.8 1.74e-38 5.22e-38 ****        
#> 3 versicolor virginica    147     -16.3 1.22e-34 3.67e-34 ****