Como desescalar os coeficientes de um modelo lmer () equipado com uma resposta em escala
Eu coloquei um modelo emR
com olmer()
-função dolme4
pacote. Eu dimensionei a variável dependente:
mod <- lmer(scale(Y)
~ X
+ (X | Z),
data = df,
REML = FALSE)
Eu olho para os coeficientes de efeito fixo comfixef(mod)
:
> fixef(mod)
(Intercept) X1 X2 X3 X4
0.08577525 -0.16450047 -0.15040043 -0.25380073 0.02350007
É muito fácil calcular os meios manualmente a partir dos coeficientes de efeitos fixos. No entanto, quero que eles não sejam dimensionados e não tenho certeza de como fazer isso exatamente. Estou ciente de que escalar significa subtrair a média de todos osY
e desviar pelo desvio padrão. Mas ambos, média e desvio padrão, foram calculados a partir dos dados originais. Posso simplesmente reverter esse processo depois de instalar umlmer()
-model usando a média e o desvio padrão dos dados originais?
Obrigado por qualquer ajuda!
Atualização: A maneira como apresentei o modelo acima parece implicar que a variável dependente é escalada, assumindo a média de todas as respostas e dividindo pelo desvio padrão de todas as respostas. Geralmente, isso é feito de maneira diferente. Em vez de calcular a média geral e o desvio padrão, as respostas são padronizadas por sujeito, usando a média e o desvio padrão das respostas desse sujeito. (Isso é estranho em umlmer()
Eu acho que a interceptação aleatória deve cuidar disso ... Sem mencionar o fato de que estamos falando de calcular meios em uma escala ordinal ...) O problema, no entanto, permanece o mesmo: uma vez que eu instalei esse modelo, existe uma maneira limpa de redimensionar os coeficientes do modelo ajustado?