Hello - new to StataList service. I am trying to get trends in surgical mortality estimates at a hospital referral region to use as an an outcome predictor for another model. This is a patient level model, I am also risk adjusting with demographic and comorbidity data.
local elix "htn_c chf valve pulmcirc perivasc para neuro chrnlung dm dmcx hypothy renlfail liver ulcer aids lymph mets tumor arth coag obese wghtloss lytes bldloss anemdef psych depress"
regress death30 year i.hrrnum i.hrrnum#c.year i.race1 i.sex1 i.elective i.`elix', vce(cluster hrrnum)
lincom year + 345.hrrnum#c.year
I want to sum the year + i.hhrnum#c.year coeffecients to get the risk-adjusted sloped for each hrr leve (306 levels).
It seems that I would have to write a loop that would calculate lincom for each hrr level and store the estimates? Is there an easier way to do this?
local elix "htn_c chf valve pulmcirc perivasc para neuro chrnlung dm dmcx hypothy renlfail liver ulcer aids lymph mets tumor arth coag obese wghtloss lytes bldloss anemdef psych depress"
regress death30 year i.hrrnum i.hrrnum#c.year i.race1 i.sex1 i.elective i.`elix', vce(cluster hrrnum)
lincom year + 345.hrrnum#c.year
I want to sum the year + i.hhrnum#c.year coeffecients to get the risk-adjusted sloped for each hrr leve (306 levels).
It seems that I would have to write a loop that would calculate lincom for each hrr level and store the estimates? Is there an easier way to do this?
Comment