我正在尝试复制论文。本文计算了一个线性混合模型,我需要使用R中的lme函数进行重现。该模型包含一个四向交互的“ T:P:Race:Sex”(T和P为连续变量和“种族”和“性别”具有两个因子水平,即黑人和白人以及男性和女性)。他们的结果模型只显示了T:P:“男性和白人”,因此对于这种互动,只有一个估算值,而在R中(使用lme)
model <- lme(A ~ Race + Sex + T + P + P:Sex:Race:T,
random = ~ P|S,
data = data)
输出显示此交互的三个交互估计:
T:P:Male:Black
T:P:Male:White
T:P:Female:White
我已经尝试将变量Race和Sex结合起来,因此“ Male and White”和“ Not Male and White”之间只有一个区别,但是不幸的是,这仍然会导致三向互动,因此我们仍然见
T:P:not Male and White
T:P:Male and White
但是我们只需要Male和White以及另一个作为“参考类别”。 因此,我们最终想知道是否有可能将“ Not Male and Black”作为参考类别,因此仅显示带有T和P的“ Male and Black”估算值。