source('../../Data/get_hospital_costs_data.R') DF = get_hospital_costs_data() DF$mortality_rate = DF$n_in_hospital_deaths / DF$cases DF$avg_paid = DF$avg_paid_medicare + DF$avg_paid_insurance m1 = lm(avg_paid ~ postsurgical_length, data=DF) print(summary(m1)) m2 = lm(postsurgical_length ~ avg_paid, data=DF) print(summary(m2)) DF$medicare_to_private = DF$avg_paid_medicare / DF$avg_paid_insurance m3 = lm(medicare_to_private ~ postsurgical_length + mortality_rate, data=DF) print(summary(m3)) m4 = lm(mortality_rate ~ medicare_to_private + postsurgical_length, data=DF) print(summary(m4)) plot(m4, which=1) print(DF[c(4, 11, 8), ])