我有一个属于不同物种的植物性状数据集,我正在将其放入 GLMM 中。对于一个变量(0 < Z < 1),我通过首先对变量进行 logit 转换,然后使用 Species 作为分组因子对其进行缩放,从而显着提高了拟合度
data$Z_logit <- logit(data$Z, percents=max(data$Z, na.rm = TRUE) > 1)
data = data %>% group_by(Species) %>% mutate(Z_ls = scale(Z_logit, center = T))
m_Z <– glmmTMB(Z_ls~ Predictor 1* Predictor 2 +
(1|Species), family=gaussian, data=data)
但是,我现在无法对变量进行反向转换,这对于使用 emmeans 计算估计均值并绘制它们很有用。
我想我首先必须缩小比例,然后反转 logit 转换;但是,由于它是一个dplyr管道命令,我还没有找到任何解决方案。