Can someone help me calculate predicted probability after running a sub-population analysis using logistic regression. For example,
I just ran the following:
svy, subpop (hispanicorigin): logistic mammogramadherent i.legalstatus
svy, subpop (hispanicorigin): logistic mammogramadherent i.legalstatus i.age i.education i.employment i.married i.income
svy, subpop (hispanicorigin): logistic mammogramadherent i.legalstatus i.age i.education i.married i.employment i.income i.healthcare i.insurance
margins legalstatus, at (age =(0 1 2 3 4 5 )) vsquish post vce(uncond)
marginsplot, ///
legend(rows(1) symxsize (5) region(fcolor(none)lcolor(none)) position(6)) ///
ytitle (Predicted probabilities with 95% CI, size(medlarge)justification(center) color (black)) ////
title (Citizenship Status , size(large)justification(center) color (black))
However, I got the following error message:
. margins legalstatus, at (age =(0 1 2 3 4 5 )) vsquish post vce(uncond)
at values for factor age do not sum to 1
r(198);
Can someone please tell me the proper syntax?
I just ran the following:
svy, subpop (hispanicorigin): logistic mammogramadherent i.legalstatus
svy, subpop (hispanicorigin): logistic mammogramadherent i.legalstatus i.age i.education i.employment i.married i.income
svy, subpop (hispanicorigin): logistic mammogramadherent i.legalstatus i.age i.education i.married i.employment i.income i.healthcare i.insurance
margins legalstatus, at (age =(0 1 2 3 4 5 )) vsquish post vce(uncond)
marginsplot, ///
legend(rows(1) symxsize (5) region(fcolor(none)lcolor(none)) position(6)) ///
ytitle (Predicted probabilities with 95% CI, size(medlarge)justification(center) color (black)) ////
title (Citizenship Status , size(large)justification(center) color (black))
However, I got the following error message:
. margins legalstatus, at (age =(0 1 2 3 4 5 )) vsquish post vce(uncond)
at values for factor age do not sum to 1
r(198);
Can someone please tell me the proper syntax?
Comment