Hello,
I am examining health care satisfaction with particular interest on those with neurological conditions. I am having trouble identifying the appropriate syntax for gologit2 with mi and svy/gsvy commands.
The commands for gologit2 work well when I do complete cases analysis using STATA 13. The proportional odds assumption is violated so I use gologit2 to create partial proportional odds model following steps highlighted in a paper by Richard Williams:
Mymodel runs: gsvy: gologit2 RecdCareSatisfaction AGECAT Gender MARSTAT EDUC NeuCon SelfPerceivedHealth SPerMHealth GLifeSat UnmetCare qualhealthcarerecd, autofit lrf
OR
gsvy: gologit2 RecdCareSatisfaction AGECAT Gender MARSTAT EDUC NeuCon SelfPerceivedHealth SPerMHealth GLifeSat UnmetCare qualhealthcarerecd, pl(NeuCon Gender MARSTAT GLifeSat EDUC) lrf
I do receive a warning with the output that I don't yet know if it is troubling because I'm concentrating on checking that the commands work. WARNING! 92 in-sample cases have an outcome with a predicted probability that is less than 0. See the gologit2 help section on Warning Messages for more information.
However, I can't seem to find the appropriate syntax when I apply multiple imputation to my survey data.
The code for my model works: mi estimate, dots or: svy: ologit RecdCareSatisfaction ib0.AGECAT 1.Gender ib0.MARSTAT ib0.EDUC ib0.INCOME 1.NeuCon ib0.SelfPerceivedHealth ib0.SPerMHealth ib0.GLifeSat 1.UnmetCare ib0.qualhealthcarerecd
My proportional odds assumption is violated and I am trying to use gologit2.
I have tried the following and received the output:
1. mi estimate, dots or: gsvy: gologit2 RecdCareSatisfaction AGECAT MARSTAT EDUC INCOME NeuCon SelfPerceivedHealth SPerMHealth GLifeSat UnmetCare qualhealthcarerecd
command prefix gsvy: not allowed
r(198);
2. mi estimate, dots or: cmdok: gsvy: gologit2 RecdCareSatisfaction AGECAT MARSTAT EDUC INCOME NeuCon SelfPerceivedHealth SPerMHealth GLifeSa
> t UnmetCare qualhealthcarerecd
command prefix cmdok: not allowed
r(198);
I have spent the past week reading every article by Richard Williams and the svy and mi commands seem to be supported for STATA. Is there something I am missing with the code for gologit2? Is there something lacking in my method? Any suggestions would be greatly appreciated.
I am examining health care satisfaction with particular interest on those with neurological conditions. I am having trouble identifying the appropriate syntax for gologit2 with mi and svy/gsvy commands.
The commands for gologit2 work well when I do complete cases analysis using STATA 13. The proportional odds assumption is violated so I use gologit2 to create partial proportional odds model following steps highlighted in a paper by Richard Williams:
Mymodel runs: gsvy: gologit2 RecdCareSatisfaction AGECAT Gender MARSTAT EDUC NeuCon SelfPerceivedHealth SPerMHealth GLifeSat UnmetCare qualhealthcarerecd, autofit lrf
OR
gsvy: gologit2 RecdCareSatisfaction AGECAT Gender MARSTAT EDUC NeuCon SelfPerceivedHealth SPerMHealth GLifeSat UnmetCare qualhealthcarerecd, pl(NeuCon Gender MARSTAT GLifeSat EDUC) lrf
I do receive a warning with the output that I don't yet know if it is troubling because I'm concentrating on checking that the commands work. WARNING! 92 in-sample cases have an outcome with a predicted probability that is less than 0. See the gologit2 help section on Warning Messages for more information.
However, I can't seem to find the appropriate syntax when I apply multiple imputation to my survey data.
The code for my model works: mi estimate, dots or: svy: ologit RecdCareSatisfaction ib0.AGECAT 1.Gender ib0.MARSTAT ib0.EDUC ib0.INCOME 1.NeuCon ib0.SelfPerceivedHealth ib0.SPerMHealth ib0.GLifeSat 1.UnmetCare ib0.qualhealthcarerecd
My proportional odds assumption is violated and I am trying to use gologit2.
I have tried the following and received the output:
1. mi estimate, dots or: gsvy: gologit2 RecdCareSatisfaction AGECAT MARSTAT EDUC INCOME NeuCon SelfPerceivedHealth SPerMHealth GLifeSat UnmetCare qualhealthcarerecd
command prefix gsvy: not allowed
r(198);
2. mi estimate, dots or: cmdok: gsvy: gologit2 RecdCareSatisfaction AGECAT MARSTAT EDUC INCOME NeuCon SelfPerceivedHealth SPerMHealth GLifeSa
> t UnmetCare qualhealthcarerecd
command prefix cmdok: not allowed
r(198);
I have spent the past week reading every article by Richard Williams and the svy and mi commands seem to be supported for STATA. Is there something I am missing with the code for gologit2? Is there something lacking in my method? Any suggestions would be greatly appreciated.
Comment