Перевод кода Stata в R: веса опроса

Я пытаюсь взвесить данные опроса в R. Я использую код Stata в качестве ссылки.

Код, использованный в этой программе для взвешивания моих данных, был:

svyset [iw=wgtp], vce(sdr), sdrweight(wgtp1-wgtp80) mse

Я перевел этот код в код R следующим образом:

svydesign(ids = ~1, data = house.data, weights = (~wgtp1-wgtp80))

Однако я не уверен, правильный ли это перевод, и правильно ли я применяю его в R. Любые идеи и предложения будут очень благодарны!


r survey stata weighting
person Matt    schedule 13.10.2017    source источник
comment
привет, это выглядит неправильно ... возможно, см. таблицу 2 в journal.r-project.org/archive/2009-2/RJournal_2009-2_Damico.pdf, а также примечание на странице 63 pdf в cran.r-project.org/web/packages/survey/survey.pdf?   -  person Anthony Damico    schedule 15.10.2017

Ответы (1)


arrow_upward
0
arrow_downward

С текущей версией "опроса" вы хотите

svrepdesign(data=house.data, weights=~wgtp, repweights="wgtp[1-9]+",type="successive-difference")

Это svrepdesign для обследования с повторными весами, весами выборки wgtp, повторными весами wgtp1-wgtp80 и весами "последовательной разницы".

person Thomas Lumley    schedule 16.05.2020