Estoy intentando extraer la fórmula de un modelo lineal generalizado mixto (GLMM).
El modelo lo he realizado a partir de este dataset:
https://drive.google.com/file/d/1ZqhBev5dFATvyWmCGdgU4KEi9Fsrj5lz/view?usp=sharing
Este dataset contiene datos de concentración (Conc: variable respuesta), maceración (Grinding : variable predictora 1) y kit (Kit: variable predictora 2).
El GLMM lo realizo a partir de este código:
library(lme4)
library(readxl)
library(equatiomatic)
glmer.conc.Grinding_Kit <- glmer(Conc ~ Grinding * Kit + (1|Especie),
family = Gamma(link = "log"),data = DATA)
He visto que existe un paquete llamado equatomatic que permite extraer la fórmula del modelo. He ejecutado el código de dicho paquete:
equatiomatic::extract_eq(glmer.conc.Grinding_Kit)
Obteniendo el siguiente error:
Error in UseMethod("extract_lhs", model) :
no applicable method for 'extract_lhs' applied to an object of class "c('glmerMod', 'merMod')"
¿Dónde está mi error? Por favor, me gustaría saber como obtener la fórmula del modelo. Si no es con esta función que sea con otra, pero lo necesito.
Muchas gracias
lme4
. Por ahora acepta objetollm
y algunas de las regresiones implementadas enMASS
yordinal
. cran.r-project.org/web/packages/equatiomatic/vignettes/… la referencia