library(ggplot2)

library(sjPlot)

library(sjmisc)

modelo2=glm(crecimiento~fitofago+trat,family=gaussian(link='log'))

plot_model(modelo2, type = "pred", terms = c("fitofago","trat"))

theme_set(theme_sjplot())

library(emmeans) #####LO VEMOS EL MARTES  CON EL EJERCICIO######

comparaciones_post1<- emmeans(modelo2, "fitofago")

pairs(comparaciones_post1)

summary(comparaciones_post1,type="response")