How to do specific, custom contrasts in EMMEANs with multiple nested factor levels but without subsetting data

问题

Here is my data frame (my real DF has way more data points):

    rearing.temp<-c("15", "15", "15", "15", "19", "19", "19", "19")
    source<-c("field", "field", "woods", "woods", "field", "field", "woods", "woods")
    runway.temp<-c("40","20","40","20","40","20","40","20")
    velocity<-c("2.3", "2.1", "1.9", "1.9", "2.3", "2.2", "2.3", "2.0")
    snail<-data.frame(rearing.temp, source, runway.temp, velocity)

Here is my model:

mod <- lmer(velocity ~ runway.temp*source*rearing.temp + (1|family) + (1|collection site) + (1|individual.plus.family.id.combined), data=snail)

When I do an emmeans contrast:

emmeans(mod, pairwise~runway.temp*source*rearing.temp)

I get 28 different comparisons, but I am only interested in looking at the difference between the velocity of field snails reared at 15° tested at the 40° runway temperature compared to woods snails reared at 15° tested at the 40° runway temperature. I just want to do one comparison between the snails that are both reared at the same temp, are tested at the same temp, but are sourced from different habitats. How can I do this?

Thank you,

Ricardo

回答1:

The contrast function provides for specifying a named list of desired contrast coefficients. In your case, apparently you have 8 means (2x2x2 array), and they are arranged in order with the first index alternating fastest and the last factor the slowest - (1,1,1), (2,1,1), (1,2,1), ..., (2,2,2). So do something like this:

emm <- emmeans(mod, ~runway.temp*source*rearing.temp)
custom <- list(`111vs211` = c(1,0,0,0,-1,0,0,0),
               `211vs121` = c(0,1,-1,0,0,0,0,0))
contrast(emm, custom)

There is abundant documentation and examples. See for example ? contrast and vignette("contrasts", "emmeans")

来源：https://stackoverflow.com/questions/65329821/how-to-do-specific-custom-contrasts-in-emmeans-with-multiple-nested-factor-leve