#r #survey #weighted
#r #опрос #взвешенные
Вопрос:
Я создал взвешенные данные в своем опросе, используя anesrake weights пакет and . Однако я не уверен, как использовать веса впоследствии, помимо wpct функции в пакете. Как я могу вычислить, скажем, описательную статистику и интегрировать взвешенные данные с другими функциями / пакетами?
Воспроизводимые данные из anesrake пакета:
data("anes04")
anes04$caseid <- 1:length(anes04$age)
anes04$agecats <- cut(anes04$age, c(0, 25,35,45,55,65,99))
levels(anes04$agecats) <- c("age1824", "age2534", "age3544",
"age4554", "age5564", "age6599")
marriedtarget <- c(.4, .6)
agetarg <- c(.10, .15, .17, .23, .22, .13)
names(agetarg) <- c("age1824", "age2534", "age3544",
"age4554", "age5564", "age6599")
targets <- list(marriedtarget, agetarg)
names(targets) <- c("married", "agecats")
outsave <- anesrake(targets, anes04, caseid=anes04$caseid,
verbose=TRUE)
caseweights <- data.frame(cases=outsave$caseid, weights=outsave$weightvec)
Это даст мне новый вектор с весами для фрейма данных. Итак, мой вопрос в том, как я могу проанализировать данные, которые знаю? Как я могу включить эти веса в сводную статистику?
Комментарии:
1. Пожалуйста, покажите небольшой воспроизводимый пример
2. Извините, у меня нет REPEX, но я подумал, что у меня нет, потому что это вопрос о конкретном пакете / методе.
Ответ №1:
Вы могли бы указать веса в качестве weights= аргумента survey::svydesign . В идеале вы должны выполнить сгребание в пакете опроса, чтобы вы могли учитывать уменьшение дисперсии при сгребании, но довольно стандартно (по крайней мере, в общедоступных данных) анализировать сгребаемые веса, как если бы они были просто весами выборки.
Или, если полученная вами спецификация сгребания была достаточно простой для воспроизведения в survey::rake или survey::calibrate , вы можете повторить сгребание в пакете опроса.
Причиной использования пакета опроса является очень широкий спектр других анализов, которые он позволяет (и даже больше svyVGAM ).