Problem Set 5:
Hospital Service Coverage in Pennsylvania

---

1. Spread of hospitals across Pennsylvania

qmplot(x,y,data = dat)

2. Philadelphia Hospitals with Life Support facilities

ggmap(get_map(location="Philadelphia"))+geom_point(aes(x, y, colour = als, fill = bas_ls), data = dat)

3. Pennsylvania Hospitals with Pediatric Services

qmplot(x,y,data = dat, col = pediatric)

4. Philadelphia Hospitals with Physical Therapy and Physical Medicine facilities

ggmap(get_map(geocode("Philadelphia"))) + 
  stat_bin2d(aes(x, y, fill = phys_med, colour = phys_ther), data = dat)

5. Hospitals with Cardiac facilities in Philadelphia

ggmap(get_map(location="Philadelphia", zoom=11)) +
  stat_density2d(aes(x, y, colour = cardiac), data = dat, bins = 3)