******************************************************************************* *Paper Title: Practice Specialization and Post-acute Outcomes of Patients with *Dementia in Skilled Nursing Facilities ******************************************************************************* ******************************************************************************* *(1) Clean up data ******************************************************************************* clear use "/project/Ryskina_SNFists/Seiyoun_working/ADRD/Stata_code_and_data/group2_after_match_12_20_23.dta" drop if weight==. gen group2=1 drop _treated save "/project/Ryskina_SNFists/Seiyoun_working/ADRD/Stata_code_and_data/group2_12_20_23.dta", replace use "/project/Ryskina_SNFists/Seiyoun_working/ADRD/Stata_code_and_data/group3_after_match_12_20_23.dta" drop if weight==. gen group2=0 save "/project/Ryskina_SNFists/Seiyoun_working/ADRD/Stata_code_and_data/group3_12_20_23.dta", replace append using "/project/Ryskina_SNFists/Seiyoun_working/ADRD/Stata_code_and_data/group2_12_20_23.dta" save "/project/Ryskina_SNFists/Seiyoun_working/ADRD/Stata_code_and_data/regression_final_data_12_20_23.dta", replace *Declare panel egen provider_id = group(SNF_CCN) xtset provider_id gen post=0 replace post=1 if SNF_ADMIT_YEAR==2019 *deflate 2019 cost to 2012 gen adj_total_cost_30_days=total_cost_30_days replace adj_total_cost_30_days=total_cost_30_days*(414.924/498.413) if SNF_ADMIT_YEAR==2019 gen adj_total_cost_60_days=total_cost_60_days replace adj_total_cost_60_days=total_cost_60_days*(414.924/498.413) if SNF_ADMIT_YEAR==2019 *winsorize winsor2 adj_total_cost_30_days, cuts(1 99) by(SNF_ADMIT_YEAR) winsor2 adj_total_cost_60_days, cuts(1 99) by(SNF_ADMIT_YEAR) save "/project/Ryskina_SNFists/Seiyoun_working/ADRD/Stata_code_and_data/regression_final_data_2_12_20_23.dta" ******************************************************************************* *(2) Regression ******************************************************************************* use "/project/Ryskina_SNFists/Seiyoun_working/ADRD/Stata_code_and_data/regression_final_data_2_12_20_23.dta" *Create global variables global risk_read_adj_var b0700 intact_cognitive missing_cognitive c1600 e0800 e0900 g0110c g0110d wander_walk two_person_support g0110h k0100c j1100a j1100b j1400 j1550d m1030 m1040e h0100c k0500a k0500b n0400f o0100a o0100b o0100c o0100f o0100h o0100i O0400d2 i0100 i0200 i1300 i2400 i4200 i4800 i5400 return_to_nursinghome first_assess global risk_ed_adj_var b0700 intact_cognitive missing_cognitive c1600 e0800 e0900 g0110c g0110d wander_walk two_person_support g0110h k0100c j1100a j1100b j1400 j1550d m1030 m1040e m1040f h0100c k0500a k0500b n0400e n0400f o0100c o0100e o0100f o0100h o0100i O0400d2 i0100 i2400 return_to_nursinghome global dc_risk_adj_var b0700 b0800 b1000 cognitive_impairment c1600 cps psychosis_bhvrl e0900 transfer_depend room_walk_depend corridor_walk_depend locomotion_depend dress_depend eat_depend toilet_depend hygine_depend bath_depend moving_balance walking_balance id_dd urinary_incontinence j1100a j1100b j1100c swallow_disorder weight_loss i2500 hemiplegia paraplegia quadriplegia multi_sclerosis huntinton i5400 m1040a m1040b m1040e k0500 max_injection o0100a o0100b o0100c o0100d o0100e o0100f o0100i o0100j i0200 i0600 i0700 i2000 i2400 i2900 i3200 i3900 i4000 i4200 i4800 i5600 i5700 i5900 i5950 i6000 i6200 *Improvement in ADL score sum one_point_adl_improve [aweight=_weight] reg one_point_adl_improve i.group2##i.post [pw=weight] margins i.group2##i.post margins r.group2##r.post reg one_point_adl_improve i.group2##i.post $risk_read_adj_var MA_pene VBP i.HHCI_cat paymedicare paymedicaid rn_hours lpn_hours cna_hours EPISODE_LOS_DAY_CNT [pw=weight] *Unplanned readmission rate sum UNPLANNED_READMIT30 [aweight=_weight] reg UNPLANNED_READMIT30 i.group2##i.post [pw=weight] margins i.group2##i.post margins r.group2##r.post reg UNPLANNED_READMIT30 i.group2##i.post $risk_read_adj_var MA_pene VBP i.HHCI_cat paymedicare paymedicaid rn_hours lpn_hours cna_hours EPISODE_LOS_DAY_CNT [pw=weight] *Discharge to community sum DISCHARGE_TO_COMMUNITY [aweight=_weight] reg DISCHARGE_TO_COMMUNITY i.group2##i.post [pw=weight] margins i.group2##i.post margins r.group2##r.post reg DISCHARGE_TO_COMMUNITY i.group2##i.post $dc_risk_adj_var MA_pene VBP i.HHCI_cat paymedicare paymedicaid rn_hours lpn_hours cna_hours EPISODE_LOS_DAY_CNT [pw=weight] *ED_FLAG sum ED_FLAG [aweight=_weight] reg ED_FLAG i.group2##i.post [pw=weight] margins i.group2##i.post margins r.group2##r.post reg ED_FLAG i.group2##i.post $risk_ed_adj_var MA_pene VBP i.HHCI_cat paymedicare paymedicaid rn_hours lpn_hours cna_hours EPISODE_LOS_DAY_CNT [pw=weight] *Observational care visits replace Obs_visit_flag_revised=0 if Obs_visit_flag_revised==. sum Obs_visit_flag_revised [aweight=_weight] reg Obs_visit_flag_revised i.group2##i.post [pw=weight] margins i.group2##i.post margins r.group2##r.post reg Obs_visit_flag_revised i.group2##i.post $risk_read_adj_var MA_pene VBP i.HHCI_cat paymedicare paymedicaid rn_hours lpn_hours cna_hours EPISODE_LOS_DAY_CNT [pw=weight] *Became long stay after SNF discharge within 6 month sum BECAME_LONG_STAY_WITHIN6MONTHS [aweight=_weight] reg BECAME_LONG_STAY_WITHIN6MONTHS i.group2##i.post [pw=weight] margins i.group2##i.post margins r.group2##r.post reg BECAME_LONG_STAY_WITHIN6MONTHS i.group2##i.post $risk_read_adj_var MA_pene VBP i.HHCI_cat paymedicare paymedicaid rn_hours lpn_hours cna_hours EPISODE_LOS_DAY_CNT [pw=weight] *SNF LOS sum EPISODE_LOS_DAY_CNT [aweight=_weight] reg EPISODE_LOS_DAY_CNT i.group2##i.post [pw=weight] margins i.group2##i.post margins r.group2##r.post reg EPISODE_LOS_DAY_CNT i.group2##i.post $risk_read_adj_var MA_pene VBP i.HHCI_cat paymedicare paymedicaid rn_hours lpn_hours cna_hours [pw=weight] *30 day Cost sum adj_total_cost_30_days_w [aweight=_weight] reg adj_total_cost_30_days_w i.group2##i.post [pw=weight] margins i.group2##i.post margins r.group2##r.post sum adj_total_cost_60_days_w [aweight=_weight] reg adj_total_cost_60_days_w i.group2##i.post [pw=weight] margins i.group2##i.post margins r.group2##r.post gen dd=group2*post glm adj_total_cost_30_days_w i.group2 i.post dd $risk_read_adj_var MA_pene VBP i.HHCI_cat paymedicare paymedicaid rn_hours lpn_hours cna_hours EPISODE_LOS_DAY_CNT [pw=weight], family(gamma) link(log) margins, dydx(dd) glm adj_total_cost_60_days_w i.group2 i.post dd $risk_read_adj_var MA_pene VBP i.HHCI_cat paymedicare paymedicaid rn_hours lpn_hours cna_hours EPISODE_LOS_DAY_CNT [pw=weight], family(gamma) link(log) margins, dydx(dd)