Acredito que efeitos como esses são freqüentemente causados por colinearidade (veja esta pergunta ). Eu acho que o livro sobre modelagem multinível de Gelman e Hill fala sobre isso. O problema é que IV1
está correlacionado com um ou mais dos outros preditores e, quando todos são incluídos no modelo, sua estimativa se torna irregular.
Se o inversão do coeficiente é devido à colinearidade, não é realmente interessante relatar, porque não é devido à relação entre seus preditores e o resultado, mas realmente devido à relação entre preditores.
O que vi sugerido para resolver esse problema é a residualização. Primeiro, você ajusta um modelo e IV2 ~ IV1
, em seguida, toma os resíduos desse modelo como rIV2
. Se todas as suas variáveis estiverem correlacionadas, você deve realmente residualizar todas elas. Você pode optar por fazer isso assim
rIV2 <- resid(IV2 ~ IV1)
rIV3 <- resid(IV3 ~ IV1 + rIV2)
rIV4 <- resid(IV4 ~ IV1 + rIV2 + rIV3)
Agora, ajuste o modelo final com
DV ~ IV1 + rIV2 + rIV3 + rIV4
Agora, o coeficiente para rIV2
representa o efeito independente de IV2
dada sua correlação com IV1
. Ouvi dizer que você não obterá o mesmo resultado se você se restabelecer em uma ordem diferente, e que escolher a ordem de residualização é realmente um julgamento na sua pesquisa.