ANLY 512 - Problem Set 5

Mapping Trials

Objective

In this report, we Create a series of 5 maps that highlight spatial differences in hospital service coverage for the state of Pennsylvania. We will explore the data set on Pennsylvania hospitals (http://www.arcgis.com/home/item.html?id=eccee5dfe01e4c4283c9be0cfc596882) for our analysis.

Analysis

Map 1: Kidney Trasnplant in Philadelphia

pennsylv_data <- read.dbf("/Users/Himanshu/Downloads/pennsylv/pennsylv.dbf")
pennsylv_mapKT <- qmap('Philadelphia', zoom = 11, source = "google", maptype = "roadmap")

## Map from URL : http://maps.googleapis.com/maps/api/staticmap?center=Philadelphia&zoom=11&size=640x640&scale=2&maptype=roadmap&language=en-EN&sensor=false

## Information from URL : http://maps.googleapis.com/maps/api/geocode/json?address=Philadelphia&sensor=false

## Warning: `panel.margin` is deprecated. Please use `panel.spacing` property
## instead

Kidney_Transplant <- subset(pennsylv_data, kidney_tra =="Y")
Kidney_TransplantPA <- subset(Kidney_Transplant, county == "Philadelphia")
pennsylv_mapKT + geom_point(aes(x, y), size = 8, alpha = 1/2, colour = "blue", data = Kidney_TransplantPA) + labs(title="Kidney Trasnplant in Philadelphia")

Map 2: Heart Trasnplant in Philadelphia

pennsylv_mapHT <- qmap('Philadelphia', zoom = 11, source = "google", maptype = "hybrid")

## Map from URL : http://maps.googleapis.com/maps/api/staticmap?center=Philadelphia&zoom=11&size=640x640&scale=2&maptype=hybrid&language=en-EN&sensor=false

## Information from URL : http://maps.googleapis.com/maps/api/geocode/json?address=Philadelphia&sensor=false

## Warning: `panel.margin` is deprecated. Please use `panel.spacing` property
## instead

Heart_Transplant <- subset(pennsylv_data, heart_tran =="Y")
Heart_TransplantPA <- subset(Heart_Transplant, county == "Philadelphia")
pennsylv_mapHT + geom_point(aes(x, y), size = 8, alpha = 1/2, colour = "red", data = Heart_TransplantPA) + labs(title="Heart Trasnplant in Philadelphia")

Map 3: Chemo Therapy in Philadelphia

pennsylv_mapCT <- qmap('Philadelphia', zoom = 11, source = "stamen", maptype = "watercolor")

## Map from URL : http://maps.googleapis.com/maps/api/staticmap?center=Philadelphia&zoom=11&size=640x640&scale=2&maptype=terrain&sensor=false

## Information from URL : http://maps.googleapis.com/maps/api/geocode/json?address=Philadelphia&sensor=false

## Map from URL : http://tile.stamen.com/watercolor/11/595/774.jpg

## Map from URL : http://tile.stamen.com/watercolor/11/596/774.jpg

## Map from URL : http://tile.stamen.com/watercolor/11/597/774.jpg

## Map from URL : http://tile.stamen.com/watercolor/11/595/775.jpg

## Map from URL : http://tile.stamen.com/watercolor/11/596/775.jpg

## Map from URL : http://tile.stamen.com/watercolor/11/597/775.jpg

## Map from URL : http://tile.stamen.com/watercolor/11/595/776.jpg

## Map from URL : http://tile.stamen.com/watercolor/11/596/776.jpg

## Map from URL : http://tile.stamen.com/watercolor/11/597/776.jpg

## Warning: `panel.margin` is deprecated. Please use `panel.spacing` property
## instead

chemo_therapy <- subset(pennsylv_data, chemo =="Y")
chemo_therapyPA <- subset(chemo_therapy, county == "Philadelphia")
pennsylv_mapCT + geom_point(aes(x, y), size = 8, alpha = 1/2, colour = "black", data = chemo_therapyPA) + labs(title="Chemo Therapy in Philadelphia")

Map 4: Hospitals with ICU Beds in Pennsylvania

qmplot(x,y, data = pennsylv_data, size = icu_beds, maptype = "toner", color = I("dark green"), extent = "panel", main = "Hospitals with ICU Beds in Pennsylvania")

## Using zoom = 8...

## Map from URL : http://tile.stamen.com/toner/8/70/94.png

## Map from URL : http://tile.stamen.com/toner/8/71/94.png

## Map from URL : http://tile.stamen.com/toner/8/72/94.png

## Map from URL : http://tile.stamen.com/toner/8/73/94.png

## Map from URL : http://tile.stamen.com/toner/8/74/94.png

## Map from URL : http://tile.stamen.com/toner/8/70/95.png

## Map from URL : http://tile.stamen.com/toner/8/71/95.png

## Map from URL : http://tile.stamen.com/toner/8/72/95.png

## Map from URL : http://tile.stamen.com/toner/8/73/95.png

## Map from URL : http://tile.stamen.com/toner/8/74/95.png

## Map from URL : http://tile.stamen.com/toner/8/70/96.png

## Map from URL : http://tile.stamen.com/toner/8/71/96.png

## Map from URL : http://tile.stamen.com/toner/8/72/96.png

## Map from URL : http://tile.stamen.com/toner/8/73/96.png

## Map from URL : http://tile.stamen.com/toner/8/74/96.png

## Map from URL : http://tile.stamen.com/toner/8/70/97.png

## Map from URL : http://tile.stamen.com/toner/8/71/97.png

## Map from URL : http://tile.stamen.com/toner/8/72/97.png

## Map from URL : http://tile.stamen.com/toner/8/73/97.png

## Map from URL : http://tile.stamen.com/toner/8/74/97.png

## Warning: Removed 208 rows containing missing values (geom_point).

Map 5: Hospitals with Pharmacy in Pennsylvania

pennsylv_mapPA <- qmap('Pennsylvania', zoom = 7, source = "google", maptype = "terrain")

## Map from URL : http://maps.googleapis.com/maps/api/staticmap?center=Pennsylvania&zoom=7&size=640x640&scale=2&maptype=terrain&language=en-EN&sensor=false

## Information from URL : http://maps.googleapis.com/maps/api/geocode/json?address=Pennsylvania&sensor=false

## Warning: `panel.margin` is deprecated. Please use `panel.spacing` property
## instead

pennsylv_mapPA + stat_bin2d(aes(x,y,color = pharmacy, fill = pharmacy), size = 0.5, bins = 20, alpha = 1/2, data = pennsylv_data) + labs(title="Hospitals with Pharmacy in Pennsylvania")