> We accounted for patient characteristics, physician characteristics, and hospital fixed effects. Patient characteristics included patient age in 5-year increments (the oldest group was categorized as ≥95 years), sex, race/ethnicity (non-Hispanic white, non-Hispanic black, Hispanic, and other), primary diagnosis (Medicare Severity Diagnosis Related Group), 27 coexisting conditions (determined using the Elixhauser comorbidity index28), median annual household income estimated from residential zip codes (in deciles), an indicator variable for Medicaid coverage, and indicator variables for year. Physician characteristics included physician age in 5-year increments (the oldest group was categorized as ≥70 years), indicator variables for the medical schools from which the physicians graduated, and type of medical training (ie, allopathic vs osteopathic29 training).
They don't say _how_ they controlled for those characteristics though. Presumably, they divided each patient's calculated 30-day risk of death [1, table 10] by the relative ratios in risk of death of every category, calculated from the same data set. That should probably cut it for controlling for this particular difference, although I am not a statistician.
It's a regression model, so "accounted for" usually means that the factors were included in the models.
Ideally, we'd ask a group of male doctors and a group of female doctors to independently diagnose and treat the same set of patients. We can't actually do that--in addition to the cost, you obviously can't treat the same patient twice. However, we can try to estimate this effect statistically.
First, they build a model describing the probability of a patient dying within 30 days. They used a linear probability model, which essentially means that the probability of someone dying is the sum of the "weights" related to the patient, doctor, and hospital. These weights are estimated from the data using ordinary least squares (the same way you may have learned to fit a line to points at school).
Having built the model, they then ask what the marginal effect of the doctor's sex is. In other words, if you hold everything but that constant, how does the probability of dying change? They cite a nice Stata guide (#32 in the paper here: http://www.stata-journal.com/sjpdf.html?articlenum=st0260 ) which gives some background and examples.
The rest of the paper looks at different variations (only hospitalists, different diseases, etc), using a pretty similar approach.
> We accounted for patient characteristics, physician characteristics, and hospital fixed effects. Patient characteristics included patient age in 5-year increments (the oldest group was categorized as ≥95 years), sex, race/ethnicity (non-Hispanic white, non-Hispanic black, Hispanic, and other), primary diagnosis (Medicare Severity Diagnosis Related Group), 27 coexisting conditions (determined using the Elixhauser comorbidity index28), median annual household income estimated from residential zip codes (in deciles), an indicator variable for Medicaid coverage, and indicator variables for year. Physician characteristics included physician age in 5-year increments (the oldest group was categorized as ≥70 years), indicator variables for the medical schools from which the physicians graduated, and type of medical training (ie, allopathic vs osteopathic29 training).
They don't say _how_ they controlled for those characteristics though. Presumably, they divided each patient's calculated 30-day risk of death [1, table 10] by the relative ratios in risk of death of every category, calculated from the same data set. That should probably cut it for controlling for this particular difference, although I am not a statistician.
Paper link:
https://jamanetwork.com/journals/jamainternalmedicine/fullar...
[1] Supplemental material https://jamanetwork.com/data/Journals/INTEMED/0/IOI160102sup...