R中的geeglm：关于可变长度的错误消息

我正在尝试使用以下代码运行一个带有logit结果的GEE模型

mod.gee <- geeglm(general_elec~activity_outside_home+econ_scale_6_pt, 
                              data=D_work_small, id="pidlink", 
                              family=binomial(link="logit"), corstr="ar1")

我的数据是来自两次调查的纵向数据。我尝试过省略NAs的响应，并更改变量的数据类型，但没有任何效果。有什么建议吗，或者以前有人遇到过这个问题吗

我的数据结构如下：

> head(D_work_small)
     pidlink econ_scale_6_pt activity_outside_home general_elec wave age female java educ_level   pid_unit
1  001220001               3                     1            1    3  48      0    0          1 0012200013
10 001220002               3                     1            1    3  47      1    0          1 0012200023
19 001220003               2                     1            1    3  27      0    0          4 0012200033
77 001250003               2                     1            1    3  27      0    0          1 0012500033
79 001290001               2                     1            1    3  52      0    0          1 0012900013
88 001290002               2                     1            1    3  49      1    0          1 0012900023

> summary(D_work_small)
   pidlink          econ_scale_6_pt activity_outside_home  general_elec         wave      
 Length:44106       Min.   :1.000   Min.   :0.0000        Min.   :0.0000   Min.   :3.000  
 Class :character   1st Qu.:2.000   1st Qu.:0.0000        1st Qu.:1.0000   1st Qu.:3.000  
 Mode  :character   Median :3.000   Median :1.0000        Median :1.0000   Median :4.000  
                    Mean   :2.894   Mean   :0.7048        Mean   :0.8304   Mean   :3.608  
                    3rd Qu.:3.000   3rd Qu.:1.0000        3rd Qu.:1.0000   3rd Qu.:4.000  
                    Max.   :6.000   Max.   :1.0000        Max.   :1.0000   Max.   :4.000  
      age             female            java          educ_level    pid_unit        
 Min.   : 14.00   Min.   :0.0000   Min.   :0.0000   Min.   :1.0   Length:44106      
 1st Qu.: 26.00   1st Qu.:0.0000   1st Qu.:0.0000   1st Qu.:1.0   Class :character  
 Median : 35.00   Median :1.0000   Median :0.0000   Median :2.0   Mode  :character  
 Mean   : 37.35   Mean   :0.5118   Mean   :0.4171   Mean   :2.2                     
 3rd Qu.: 47.00   3rd Qu.:1.0000   3rd Qu.:1.0000   3rd Qu.:3.0                     
 Max.   :999.00   Max.   :1.0000   Max.   :1.0000   Max.   :5.0        

---

您的

id

变量应该不带引号，如

id=pidlink

中所示。此外，由于在计算中将其转换为数值（通过

as.numeric

），您应该将其从字符转换为

factor

ie

id=factor（pidlink）

。[除此之外，

activity\u outside\u home

和

econ\u scale\u 6\u pt

似乎也应该是因子变量，您应该通过

pidlink

对数据进行排序。]。

> head(D_work_small)
     pidlink econ_scale_6_pt activity_outside_home general_elec wave age female java educ_level   pid_unit
1  001220001               3                     1            1    3  48      0    0          1 0012200013
10 001220002               3                     1            1    3  47      1    0          1 0012200023
19 001220003               2                     1            1    3  27      0    0          4 0012200033
77 001250003               2                     1            1    3  27      0    0          1 0012500033
79 001290001               2                     1            1    3  52      0    0          1 0012900013
88 001290002               2                     1            1    3  49      1    0          1 0012900023

> summary(D_work_small)
   pidlink          econ_scale_6_pt activity_outside_home  general_elec         wave      
 Length:44106       Min.   :1.000   Min.   :0.0000        Min.   :0.0000   Min.   :3.000  
 Class :character   1st Qu.:2.000   1st Qu.:0.0000        1st Qu.:1.0000   1st Qu.:3.000  
 Mode  :character   Median :3.000   Median :1.0000        Median :1.0000   Median :4.000  
                    Mean   :2.894   Mean   :0.7048        Mean   :0.8304   Mean   :3.608  
                    3rd Qu.:3.000   3rd Qu.:1.0000        3rd Qu.:1.0000   3rd Qu.:4.000  
                    Max.   :6.000   Max.   :1.0000        Max.   :1.0000   Max.   :4.000  
      age             female            java          educ_level    pid_unit        
 Min.   : 14.00   Min.   :0.0000   Min.   :0.0000   Min.   :1.0   Length:44106      
 1st Qu.: 26.00   1st Qu.:0.0000   1st Qu.:0.0000   1st Qu.:1.0   Class :character  
 Median : 35.00   Median :1.0000   Median :0.0000   Median :2.0   Mode  :character  
 Mean   : 37.35   Mean   :0.5118   Mean   :0.4171   Mean   :2.2                     
 3rd Qu.: 47.00   3rd Qu.:1.0000   3rd Qu.:1.0000   3rd Qu.:3.0                     
 Max.   :999.00   Max.   :1.0000   Max.   :1.0000   Max.   :5.0