ANLY 512 Problem Set-5

#install.packages("ggplot2")
library(ggplot2)
#install.packages("ggmap")
library(ggmap)

## Google's Terms of Service: https://cloud.google.com/maps-platform/terms/.

## Please cite ggmap if you use it! See citation("ggmap") for details.

#install.packages("foreign")
library(foreign)
dat <- read.dbf("C:/Users/gmutya048/Downloads/pennsylv.dbf")
names(dat)

##   [1] "acc_trauma" "air_amb"    "als"        "arc_street" "arc_zone"  
##   [6] "bas_ls"     "bassinets"  "bb_id"      "bc_beds"    "bc_sus_bed"
##  [11] "beds_sus"   "birthing_r" "bone_marro" "burn_car"   "burn_care" 
##  [16] "card_beds"  "card_surge" "card_sus_b" "cardiac"    "cardiac_ca"
##  [21] "cardio_reh" "chemo"      "city"       "clin_lab"   "clin_psych"
##  [26] "county"     "countyname" "ct_scan"    "cty_key"    "cystoscopi"
##  [31] "deliv_rms"  "dental"     "detox_alc_" "diag_radio" "diag_xray" 
##  [36] "doh_hosp"   "doh_phone"  "emer_dept"  "endoscopie" "fac_id"    
##  [41] "facility"   "flu_old"    "fred_con_1" "fred_conta" "fred_email"
##  [46] "fred_fax"   "fred_hosp"  "fred_pager" "fred_phone" "gamma_knif"
##  [51] "gen_outpat" "gene_counc" "heart_tran" "helipad"    "hemodial_c"
##  [56] "hemodial_m" "hosp_id"    "hospice"    "hyper_cham" "icu"       
##  [61] "icu_beds"   "icu_sus_be" "inpat_flu_" "inpat_pneu" "kidney_tra"
##  [66] "labor_rms"  "lic_beds"   "lic_dent"   "lic_dos"    "lic_mds"   
##  [71] "lic_pod"    "linear_acc" "lithotrips" "liver_tran" "loc_method"
##  [76] "ltc"        "mcd"        "mcd_key"    "mcd_name"   "medical"   
##  [81] "mob_ccu"    "mob_icu"    "mri"        "ms1"        "neo2_beds" 
##  [86] "neo2_sus_b" "neo3_beds"  "neo3_sus_b" "neuro_surg" "neurology" 
##  [91] "obs_gyn"    "occ_ther"   "optometry"  "organ_bank" "ped_trauma"
##  [96] "pediatric"  "pet"        "pharmacy"   "phys_med"   "phys_ther" 
## [101] "podiatry"   "providerid" "psych"      "psych_inpa" "reg_trauma"
## [106] "resp_ther"  "so_flu_65u" "social_wor" "speech_pat" "street"    
## [111] "surgical"   "surgical_s" "thera_radi" "typ_org"    "typ_serv"  
## [116] "ultrasound" "x"          "y"          "zip"

Number of Hospitals in My County (Delaware County)

county = subset(dat, county=="Delaware")
qmplot(x,y, data = county, colour = I('darkblue'), size = I(4), zoom = 11, extent = "panel", source = "stamen", maptype = "terrain-background", main = "Hospitals in Delaware county ", xlab = "Longitude", ylab = "Latitude")

## Source : http://tile.stamen.com/terrain/11/594/775.png

## Source : http://tile.stamen.com/terrain/11/595/775.png

## Source : http://tile.stamen.com/terrain/11/594/776.png

## Source : http://tile.stamen.com/terrain/11/595/776.png

Hospitals in PA with Organ Banks

organbank = subset(dat, organ_bank=="Y")
qmplot(x,y, data = organbank, colour = I('red'), size = I(4), zoom = 8, extent = "panel", source = "osm", maptype = "toner", main = "Hospitals in Pennsylvania with Organ Bank", xlab = "Longitude", ylab = "Latitude")

## Source : http://tile.stamen.com/terrain/8/70/95.png

## Source : http://tile.stamen.com/terrain/8/71/95.png

## Source : http://tile.stamen.com/terrain/8/72/95.png

## Source : http://tile.stamen.com/terrain/8/73/95.png

## Source : http://tile.stamen.com/terrain/8/74/95.png

## Source : http://tile.stamen.com/terrain/8/70/96.png

## Source : http://tile.stamen.com/terrain/8/71/96.png

## Source : http://tile.stamen.com/terrain/8/72/96.png

## Source : http://tile.stamen.com/terrain/8/73/96.png

## Source : http://tile.stamen.com/terrain/8/74/96.png

## Source : http://tile.stamen.com/terrain/8/70/97.png

## Source : http://tile.stamen.com/terrain/8/71/97.png

## Source : http://tile.stamen.com/terrain/8/72/97.png

## Source : http://tile.stamen.com/terrain/8/73/97.png

## Source : http://tile.stamen.com/terrain/8/74/97.png

Hospitals in Philadelphia with Ultrasound

ULT = subset(dat, ultrasound == "Y")
ULTP = subset(ULT, city == "Philadelphia")
qmplot(x,y, data = ULTP , colour = I('seagreen1'), size = I(4), zoom = 11, extent = "panel", source = "stamen", maptype = "terrain-lines", main = "Hospitals in Philadelphia with Ultrasound", xlab = "Longitude", ylab = "Latitude")

## Source : http://tile.stamen.com/terrain/11/595/774.png

## Source : http://tile.stamen.com/terrain/11/596/774.png

## Source : http://tile.stamen.com/terrain/11/597/774.png

## Source : http://tile.stamen.com/terrain/11/596/775.png

## Source : http://tile.stamen.com/terrain/11/597/775.png

Non Profitable Pediatric Hospitals in PA

pediatric = subset(dat, pediatric == "Y")
pediatricNP = subset(pediatric, typ_org=="Non-profit Corp.")
qmplot(x,y, data = pediatricNP, colour = I('tomato'), size = I(4), zoom = 7, extent = "panel", source = "osm", maptype = "hybrid", main = "Non Profit Pediatric Hospitals in PA", xlab = "Longitude", ylab = "Latitude")

## Source : http://tile.stamen.com/terrain/7/35/47.png

## Source : http://tile.stamen.com/terrain/7/36/47.png

## Source : http://tile.stamen.com/terrain/7/37/47.png

## Source : http://tile.stamen.com/terrain/7/35/48.png

## Source : http://tile.stamen.com/terrain/7/36/48.png

## Source : http://tile.stamen.com/terrain/7/37/48.png

Hospitals in PA with Heart Transplantation

heart_transplant = subset(dat, heart_tran == "Y")
qmplot(x, y, data = heart_transplant, legend = "none", colour= I('blue'), mapcolor = "color", main = "Hospitals in PA with Heart Transplant Surgery Capability", xlab = "Longitude", ylab = "Latitude", source = 'stamen', maptype = 'toner-labels', zoom= 8,size = I(1.5))