1
library(lme4)
library(emmeans)
m <- lmer(angle ~ recipe*temp + (1|replicate), data=cake)

各レシピの単純な効果(ペアワイズ比較ではない)に興味があるので、この投稿に基づいて

out.emtrends <- emtrends(m, pairwise~recipe, var="temp")
emtr <- as.data.frame(out.emtrends$emtrends)
tvalues <- emtr$temp.trend/emtr$SE
dfs <- emtr$df
pvalues <- 2*pt(-abs(tvalues), dfs)

ここで、p 値を修正するために Tukey の検定を使用したいと思います。これは可能ですか、それともTukeyのテストは治療間の違いにのみ適していますか?

4

1 に答える 1