我正在运行 lme4 并且我正在处理一个数据集,其中变量 a、d、e 是对数转换的。我想知道如何处理有关“对”功能的反向转换。我如何获得对比度?它被设置为一个函数,因为实际上有很多列我将循环遍历。
此外,如果我在 lmer 模型中运行数据集之前对数据集进行 logtransform,我能说数据是 logtransformed,即使它不在模型中作为 log(df[,i]) 吗? 非常感谢!
a <- runif(21)
d <- runif(21)
e <- runif(21)
trt <- rep(LETTERS[1:3],7)
id <- rep(1:7, each=3)
df <- data.frame(id,trt,a,d,e)
func <- function(i){
dflmer <- lmer(formula= log(df[,i]) ~trt+(1|id), data=df)
emm <- emmeans(dflmer,"trt")
emm_df <- data.frame(emm)# I could get the data back-transformed like this
exp(emm[,c(2,5,6)])
emm_pairs<- pairs(emm, adjust="fdr") # but how about the pairs function?
contrast(regrid(e, type="log")) #?
}
func(3) # 示例