***********************************************************************************
*Paper Title: Nursing Home–Based Care for Children Younger Than 18 Years: 2012–2019

*Authors: Seiyoun Kim, PhD, Nurah Koney-Laryea, BA,Ziwei Pan, MS, Hye-Young Jung, PhD, 
*Jaya Aysola, MD, MPH, Sophia Jan, MD, MPH, Kira L. Ryskina, MD, MSHP
***********************************************************************************

*Import the values from the figures in Excel into Stata.

*Short Stay 
import excel "C:\Users\kimsei\Documents\Work\Kira\HCBS\HCBS\Pediatrics\short_stay_under_18_correct.xlsx", sheet("Sheet2") firstrow clear
gen Total_pop=NonHispanicWhite+Black+Hispanic+Asian+OtherRace
gen pct_nhw=(NonHispanicWhite/Total_pop)*100
gen pct_black=(Black/Total_pop)*100
gen pct_hispanic=(Hispanic/Total_pop)*100
gen pct_asian=(Asian/Total_pop)*100
gen pct_other=(OtherRace/Total_pop)*100

rename * *_short
rename Year_short Year

save "C:\Users\kimsei\Documents\Work\Kira\HCBS\HCBS\Pediatrics\short_stay_under_18_correct", replace

*Long Stay 
import excel "C:\Users\kimsei\Documents\Work\Kira\HCBS\HCBS\Pediatrics\long_stay_under_18_correct.xlsx", sheet("Sheet2") firstrow clear
gen Total_pop=NonHispanicWhite+Black+Hispanic+Asian+OtherRace
gen pct_nhw=(NonHispanicWhite/Total_pop)*100
gen pct_black=(Black/Total_pop)*100
gen pct_hispanic=(Hispanic/Total_pop)*100
gen pct_asian=(Asian/Total_pop)*100
gen pct_other=(OtherRace/Total_pop)*100

rename * *_long
rename Year_long Year

drop if Year==.

save "C:\Users\kimsei\Documents\Work\Kira\HCBS\HCBS\Pediatrics\long_stay_under_18_correct", replace

merge 1:1 Year using "C:\Users\kimsei\Documents\Work\Kira\HCBS\HCBS\Pediatrics\short_stay_under_18_correct"

gen log_pop_short=log(Total_at_risk_short)
gen log_pop_long=log(Total_at_risk_long)

*Negative binomial regression
nbreg Total_pop_short Year, offset(log_pop_short) vce(robust) 
est store NB

nbreg Total_pop_long Year, offset(log_pop_long) vce(robust) 

*Linear regression 
reg pct_nhw_short Year, vce(robust)
reg pct_black_short Year, vce(robust)
reg pct_hispanic_short Year, vce(robust)
reg pct_asian_short Year, vce(robust)
reg pct_other_short Year, vce(robust)

reg pct_nhw_long Year, vce(robust)
reg pct_black_long Year, vce(robust)
reg pct_hispanic_long Year, vce(robust)
reg pct_asian_long Year, vce(robust)
reg pct_other_long Year, vce(robust)