Question

我正在运行一个线性混合模型，以查看任务的反应时间是否在受试者，实验条件或目标之间有所不同。但是，当我运行lme时，它会警告我有关奇异配合的信息。

我知道单数拟合可能表示模型过度拟合，但我不明白为什么我的模型因我拥有的数据量而过度拟合。

有关更多信息，该实验涉及对象为一系列照片命名，我们记录反应时间（RT）。每个参与者看到所有图片（目标）并具有所有4个条件。有440个目标，每种条件下有110个目标。

我的第一个模型没有奇异拟合问题：

model1 = lmer(log(RTs300ms)~Condition+(1|Targets)+(1|Subject),data=beh_acc2)


REML criterion at convergence: -6030.481
Random effects:
 Groups   Name        Std.Dev.
 Targets  (Intercept) 0.07918 
 Subject  (Intercept) 0.13678 
 Residual             0.17972 
Number of obs: 10985, groups:  Targets, 110; Subject, 27
Fixed Effects:
  (Intercept)    ConditionTh  ConditionUnTa  ConditionUnTh  
      6.67960       -0.03549       -0.01475       -0.01700

但是从第二个模型开始，我开始遇到问题：

model2 = lmer(log(RTs300ms)~Condition+(1|Targets)+(1+Condition|Subject),data=beh_acc2)

REML criterion at convergence: -6037.6

Scaled residuals: 
    Min      1Q  Median      3Q     Max 
-7.3985 -0.6456 -0.1593  0.4551  5.3105 

Random effects:
Groups   Name          Variance  Std.Dev. Corr             
Targets  (Intercept)   6.270e-03 0.079184                  
 Subject  (Intercept)   1.947e-02 0.139546                  
          ConditionTh   2.482e-05 0.004982  0.58            
          ConditionUnTa 2.273e-04 0.015078 -0.50  0.41      
          ConditionUnTh 1.200e-04 0.010956 -0.64  0.26  0.99
 Residual               3.226e-02 0.179597                  
Number of obs: 10985, groups:  Targets, 110; Subject, 27

Fixed effects:
               Estimate Std. Error t value
(Intercept)    6.679624   0.028108 237.641
ConditionTh   -0.035441   0.004949  -7.162
ConditionUnTa -0.014807   0.005648  -2.621
ConditionUnTh -0.017071   0.005293  -3.225

Correlation of Fixed Effects:
            (Intr) CndtnT CondtnUnT
ConditionTh  0.022                 
ConditinUnT -0.321  0.463          
ConditnUnTh -0.321  0.471  0.597   
convergence code: 0
boundary (singular) fit: see ?isSingular

理想情况下，我希望最终模型能够正常工作，但这是三个模型中最复杂的一个：

model3 = lmer(log(RTs300ms)~Condition+(1+Condition|Targets)+(1+Condition|Subject),data=beh_acc2)

REML criterion at convergence: -6068.3

Scaled residuals: 
    Min      1Q  Median      3Q     Max 
-7.3684 -0.6393 -0.1575  0.4541  5.2834 

Random effects:
 Groups   Name          Variance  Std.Dev. Corr             
 Targets  (Intercept)   6.931e-03 0.083252                  
          ConditionTh   1.008e-03 0.031748 -0.19            
          ConditionUnTa 1.419e-03 0.037676 -0.35  0.58      
          ConditionUnTh 1.727e-03 0.041553 -0.29  0.72  0.98
 Subject  (Intercept)   1.951e-02 0.139662                  
          ConditionTh   2.575e-05 0.005074  0.57            
          ConditionUnTa 2.371e-04 0.015399 -0.50  0.43      
          ConditionUnTh 1.237e-04 0.011124 -0.63  0.28  0.99
 Residual               3.187e-02 0.178527                  
Number of obs: 10985, groups:  Targets, 110; Subject, 27

Fixed effects:
               Estimate Std. Error t value
(Intercept)    6.679466   0.028234 236.574
ConditionTh   -0.035074   0.005785  -6.063
ConditionUnTa -0.014748   0.006706  -2.199
ConditionUnTh -0.016823   0.006607  -2.546

Correlation of Fixed Effects:
            (Intr) CndtnT CondtnUnT
ConditionTh -0.008                 
ConditinUnT -0.325  0.497          
ConditnUnTh -0.306  0.548  0.722   
convergence code: 0
boundary (singular) fit: see ?isSingular

我不确定如何解决在model2和model3上遇到的奇异拟合问题。我已经阅读了一些尝试不熟悉的贝叶斯模型的建议。

在此问题上的任何建议或进一步的建议，将不胜感激！

Answer 1

在第一个模型中，您不包含任何随机效果。但是，在返回错误的模型中，您现在包括了条件的随机效应。我的猜测是，您对病情的影响实际上没有随机性。

您可以使用以下代码使用rePCA中的lme4函数进行检查：结果将提供一个表格，该表格提供随机效应结构中解释的方差比例。如果您有任何一栏说明接近0的方差比例，则可能是问题所在，并导致了奇异拟合误差。确认这一点的另一种方法是由参与者绘制每种条件的平均RT，以创建所谓的“意大利面条图”。添加行并将其按参加者分组，以便您可以查看参加者针对每种情况的RT的变化。如果这些线全部平行（或几乎平行），则表明您的条件影响可能没有随机性。

如何在线性混合模型中固定奇异拟合？

1 个答案: