Таблица коэффициентов не отображает результаты для обеих независимых переменных после запуска GEE

Question

У меня есть простой вопрос здесь. Я использую gee Пакет для запуска Джи регрессии на данных ниже. Я использовал тот же набор данных в spss и взял эту таблицу в качестве результата.

Набор данных

round<-c( 0.125150,  0.045800, -0.955299, -0.232007,  0.120880, -0.041525,  0.290473, -0.648752,  0.113264, -0.403685)
    square<-c(-0.634753,  0.000492, -0.178591, -0.202462, -0.592054, -0.583173, -0.632375, -0.176673, -0.680557, -0.062127)
    ideo<-c(0,1,0,1,0,1,0,0,1,1)
    ex<-data.frame(round,square,ideo)

Когда я запускаю тот же анализ в г

library(gee)  
    exmen<-summary(gee(round ~ square,
                        data = ex, id = ideo,
                        corstr = "independence"))
    exmen

Я получил:

Coefficients:
            Estimate Naive S.E. Naive z Robust S.E. Robust z
(Intercept)   -0.510      0.181   -2.82       0.229    -2.23
square        -0.939      0.399   -2.36       0.404    -2.32

То же самое происходит с geepack пакет

library(geepack)
    exmen<-summary(geeglm(round ~ square,
                        data = ex, id = ideo,
                        corstr = "independence"))
    exmen

 Coefficients:
            Estimate Std.err Wald Pr(>|W|)  
(Intercept)   -0.510   0.229 4.95    0.026 *
square        -0.939   0.404 5.40    0.020 *
---

Так что мне интересно 2 вещи. 1. Почему я получаю результаты только для квадрата? 2.(необязательно) Можно ли воссоздать точно такую ​​же таблицу с gee или же geepack как со спсс?

Answer 1

Потому что вы указали квадрат как единственную объясняющую переменную. Зависимая переменная LeftvsRight не существует в R. При ее добавлении укажите модель следующим образом:

gee(LeftvsRight ~ square + round,data = ex, id = ideo,
                    corstr = "independence")