У меня есть большой опрос, который я хотел бы проанализировать в Stata 16. Я заинтересован в определении прогнозируемого процента смертности для детей с определенным заболеванием на основе их страховки.Я хотел бы получить эти значения для нескольких отдельных лет.Я установил параметры своего опроса и провел регрессию Пуассона для своей подгруппы интересов.
svyset [w=weight], psu(primarysamplingunit) strata(stratum) vce(linearized) singleunit(certainty)
svy, subpop(if child==1 & missing==0): poisson died i.disease age female i.race i.lowincome i.insurance i.year i.disease#i.year i.insurance#i.disease, irr
Используя эту модель, относительный риск смерти для моих переменных, представляющих интерес, показан ниже:
refinsurance#disease 1
medicare#disease 1.24
medicaid#disease 1.66
private#disease 1.4
Теперь я хотел бы получить наценки для каждой страховки по годам.Итак, я хочу знать процент смертности для детей с заболеванием, которые имеют Medicare и которые посещают в 2010 году. Я также хочу знать процент смертности для детей с заболеванием, которые имеют Medicare и наблюдаются в 2011 году, и так далее для разныхвиды страхования в разные годы.
Как я могу точно получить эту информацию?
Я пытался использовать приведенный ниже код, но результаты не имеют смысла.
margins i.year, subpop(if child==1 & missing==0 & disease==1 & insurance==1)
margins i.year, subpop(if child==1 & missing==0 & disease==1 & insurance==2)
Поля показывают, что процентная смертность является самой низкой среди программ Medicare и Medicaid, что противоречит ожидаемым результатам и результатам регрессии Пуассона.
Тогда я попробовал следующее, но я не уверен, что это действительно так.
margins i.insurance, subpop(if child==1 & missing==0 & disease==1) over(year)
Это правильный код для ответа на мой вопрос?