Я использую пакет svydesign в R для выполнения взвешенной логит-регрессии, как показано ниже:
sdobj <- svydesign(id = ~0, weights = ~chweight, strata = ~strata, data = svdat)
model1 <- svyglm(formula=formula1,design=sdobj,family = quasibinomial)
Однако в документации содержится предостережение относительно регрессий без указания поправок на конечную популяцию (FPC):
If fpc is not specified then sampling is assumed to be
with replacement at the top level and only the first stage of
cluster is used in computing variances.
К сожалению, у меня нет достаточной информации, чтобы указать мои популяции на каждом уровне (из которых я отбираю очень мало). Любая информация о том, как указать веса опроса без информации FPC, была бы очень полезной.