0

我正在对按国家和年份嵌套的调查数据进行分析。被调查的受访者从来都不是相同的，但被调查的国家却是重复的。

数据看起来像这样，其中yDVx是 IV，g是我有兴趣与 DV 交互的组变量x。数据按国家co和年份嵌套t。

dat <- data.frame(y = c(1, 2, 3, 4, 5, 4, 3, 2, 1, 2, 5, 2, 4, 3, 1),
                  x = c(1, 6, 3, 9, 3, 6, 4, 4, 9, 2, 8, 2, 5, 3, 7),
                  g = c(1, 1, 2, 2, 1, 1, 2, 2, 1, 1, 2, 2, 1, 1, 2),
                  t = c(0, 1, 0, 1, 0, 1, 0, 1, 0, 1, 0, 1, 0, 1, 0),
                  co = c("A", "A", "A", "A", "A", "A", "A", "B", "B", "B", "B", "B", "B", "B", "B"))

基本上，我想在x*g预测的地方进行纵向分析y。鉴于这是纵向分析，我认为我需要将效果与 year 进行交互，对t吗？另外，我认为我需要控制 country 的随机效应和斜率co。所以这就是我所做的：

model1 <- glmer(y ~ x*t*g + (1+x|co) + (1|co), data = dat)

stargazer(model1, type = "text")

===============================================
                        Dependent variable:    
                    ---------------------------
                                 y             
-----------------------------------------------
x                             -0.390           
                              (0.818)          

t                             -3.317           
                              (6.415)          

g                             -1.971           
                              (3.012)          

x:t                            0.471           
                              (1.218)          

x:g                            0.379           
                              (0.514)          

t:g                            1.566           
                              (4.312)          

x:t:g                         -0.359           
                              (0.798)          

Constant                       5.517           
                              (4.367)          

-----------------------------------------------
Observations                    15             
Log Likelihood                -22.627          
Akaike Inf. Crit.             71.253           
Bayesian Inf. Crit.           80.458           
===============================================
Note:               *p<0.1; **p<0.05; ***p<0.01

我不确定这是否是进行这样的纵向分析的正确方法，所以想问问是否有人可以确认或纠正这一点。谢谢。

4

0 回答 0