2

我对两个重复测量的固定效应和随机效应的建模有疑问。我的数据集如下所示:

subject id      condition      trial      DV
    1                1           1        2
    1                1           2        5
    1                1           3        6    
    1                2           1        3    
    1                2           2        8   
    1                2           3        9    
    2                1           1        5    
    2                1           2        4    
    2                1           3        3    
    2                2           1        7    
    2                2           2        8    
    2                2           3        6    
    3                1           1        4    
    3                1           2        5   
    3                1           3        7    
    3                2           1        9    
    3                2           2        10    
    3                2           3        2    

因此,有一个具有两个条件的任务,每个参与者都执行两个条件(重复测量),并且每个条件都存在 3 次试验(重复测量)。然后是连续的DV。

(1)首先,我想知道变量subject id,condition和trial是嵌套还是交叉。所有科目都在做两个级别的条件,所以我认为科目和条件是交叉的。那是对的吗?
尽管条件 1 和条件 2 都具有三个级别的试验,但它们在时间上当然不是相同的试验(因为参与者不是同时玩任务的条件,而是彼此相继)。这是否意味着条件和试验是嵌套的?还是因为他们都是三重试炼,他们也越过了?
此外,我认为主题和试验是交叉的,因为所有参与者都经历了 6 个级别的试验。这个对吗?

(2) 其次,我喜欢运行多级模型,但我不确定如何指定固定效应和随机效应。最低级别将是试验/时间,因为我想对两种条件下 DV 随时间(试验)的曲线/增长轨迹进行建模。实际上,每个条件有 20 次试验。我现在有以下代码来为每个参与者建模随机截距:

m1 = gls(DV~1, data=dataset,method = "ML", na.action = na.exclude) 
m2 = lme(DV~ 1 , data = dataset, random = ~1|subject id, method = "ML", na.action = na.exclude)

M2 比 M1 更适合数据。但这是否意味着在每个参与者的所有 6 次试验中都建立了随机截距模型?我实际上想知道每个条件。那么我应该只为每个条件制作一个多级模型(因此将两种条件的所有分析分开)?或者我可以在代码中包含条件吗?代码应该是:

m3 = lme(DV~ 1 , data = dataset, random = ~1|subject id/condition, method = "ML", na.action = na.exclude)

如果 M3 比 M1 更适合数据,这是否意味着这两种情况都存在随机截距?

(3)如果我想知道每个条件是否有试验效果,我应该如何建模?像这样?

m4 = lme(DV~ trial , data = dataset, random = ~1|subject id/condition, method = "ML", na.action = na.exclude)

试验的显着效果是否意味着在两种情况下试验都存在线性效应?

(4)最后,我想知道每个条件(所以不是整个 6 次试验)DV 和试验之间的关系是否存在随机斜率。我不想知道整个 6 次试验但每个条件(实际上每个条件 20 次试验)。模型是:

m5 = lme(DV~ trial, random = ~condition/trial|subject id/condition, method = "ML", data = dataset) 

总的来说,我很困惑如何处理模型中的两个重复变量(条件和试验)。我希望有人能帮帮忙!

4

0 回答 0