成对函数的反向转换(lmer 分析)

时间:2021-01-19 10:54:02

标签: r

我正在运行 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) # 示例

0 个答案:

没有答案