Please look at the documentation for emmeans::contrast
; it provides for such comparisons:
library("emmeans")
emm <- emmeans(m1, "Condicion")
contrast(emm, list(
`DB-DF` = c(-1,1,0,0,0,0),
`NDB-NDF` = c(0,0,-1,0,0,0),
`BB-BF` = c(0,0,0,0,-1,1),
`DB-BB` = c(-1,0,0,0,1,0)
))
Note: This is based on the assumption that the levels of Condicion
are in the order given in the question. If they are not, then they are probably in alphabetical order (the default) and you need to either permute the 1s and -1s accordingly, or create
Datos$Condicion <- factor(as.character(Datos$Condicion),
levels = c("DB", "DF", "NDB", "NDF", "BB", "BF"))
and then re-fit the model before using the above code to obtain the comparisons.