Перевод кода Stata в R: веса обследования

Question

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

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

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

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

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

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

Answer 1

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

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

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