Pruebas post-hoc sencillas al realizar un meta-análisis con el paquete `metafor` en r

Question

Estoy realizando un meta-análisis utilizando metafor en R. Me gustaría comparar entre 7 niveles en un factor (es decir, diferentes tipos de tratamientos).

fit <- rma (yi, vi, mods = type_of_treatment - 1, data = dat)
fit

He encontrado varios sitios web que explican cómo comparar entre niveles utilizando anova() (por ejemplo, Análisis post hoc tras la metarregresión & http://www.metafor-project.org/doku.php/tips:testing_factors_lincoms?s[]=anova ); sin embargo, como tengo siete niveles (es decir, siete tipos diferentes de tratamientos), hacer comparaciones por mí mismo es bastante problemático.

He intentado glht() en el paquete multcomp, pero el siguiente comando da error.

summary(glht(fit, linfct = mcp(~type_of_treatment = "Tukey"))

Les agradecería mucho que me dijeran si hay alguna forma fácil de realizar una comparación por pares entre niveles, como la glht() o emmeans() para lm.

Answer 1

Puede utilizar el contrMat() función. Algo como esto debería funcionar:

summary(glht(fit, linfct=cbind(contrMat(rep(1,7), type="Tukey"))), test=adjusted("none"))

Sin embargo, podría considerar un ajuste para las pruebas múltiples. Véase help(summary.glht) para algunas opciones.