У меня тут простой вопрос. Я использую пакет gee
для запуска регрессии 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)
Когда я запускаю тот же анализ в r
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
как с spss?