Violinplot 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. Proporciona informacion sobre la distribución de cada grupo, así como la media y los quantiles.
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
Violinplot de los contrastes posthoc en el modelo todas las variables explicativas del modelo.
Examples
library(tlamatini)
data(iris)
iris$Species <- as.factor(iris$Species)
mod <- glm(Petal.Width ~ Species, family = gaussian("log"), data=iris)
violinplot_emmeans(formula = Petal.Width ~ Species, modelo=mod, data= iris)
#> Función en desarrollo. Agregar valores de p manualmente.
#> # 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 ****