library(foreign)
library(ggmap)
## Loading required package: ggplot2
library(ggplot2)
dat <- read.dbf("pennsylv/pennsylv.dbf") # change this to connect to your version
All the hospitals in the dataset
qmplot(x, y, data = dat, size = I(3),alpha = I(0.8),extent = "device", geom = "point",main = "Locations of all the hospitals in the dataset", xlab = "Longitude", ylab = "Latitude",source = "google",maptype = "roadmap",zoom = 7)
## Map from URL : http://maps.googleapis.com/maps/api/staticmap?center=40.942826,-77.681596&zoom=7&size=640x640&scale=2&maptype=roadmap&language=en-EN&sensor=false
## Warning: `panel.margin` is deprecated. Please use `panel.spacing` property
## instead
Hospitals and the number of beds
qmplot(x, y, data = dat, size = beds_sus,alpha = I(0.8),extent = "device", geom = "point",main = "Hospitals and the number of beds", xlab = "Longitude", ylab = "Latitude",source = "google",maptype = "roadmap",zoom = 7)
## Map from URL : http://maps.googleapis.com/maps/api/staticmap?center=40.942826,-77.681596&zoom=7&size=640x640&scale=2&maptype=roadmap&language=en-EN&sensor=false
## Warning: `panel.margin` is deprecated. Please use `panel.spacing` property
## instead
## Warning: Removed 31 rows containing missing values (geom_point).
Hospitals that offer air ambulance service
amb<-subset(dat, air_amb == "Y")
qmplot(x, y, data = amb, size = I(3),alpha = I(0.8),,extent = "device", geom = "point",main = "Hospitals that offers air ambulance service", xlab = "Longitude", ylab = "Latitude",source = "google",maptype = "roadmap",zoom = 7)
## Map from URL : http://maps.googleapis.com/maps/api/staticmap?center=41.040851,-77.606184&zoom=7&size=640x640&scale=2&maptype=roadmap&language=en-EN&sensor=false
## Warning: `panel.margin` is deprecated. Please use `panel.spacing` property
## instead
## Warning: Ignoring unknown parameters: NA
Hospitals with have emergency department and air ambulance service
ED<-subset(dat, emer_dept == "Y")
qmplot(x, y, data = ED, size = I(3),alpha = I(0.8),colour=air_amb,extent = "device", geom = "point",main = "Hospitals with have emergency department and air ambulance service", xlab = "Longitude", ylab = "Latitude",source = "google",maptype = "roadmap",zoom = 7)
## Map from URL : http://maps.googleapis.com/maps/api/staticmap?center=40.942648,-77.681596&zoom=7&size=640x640&scale=2&maptype=roadmap&language=en-EN&sensor=false
## Warning: `panel.margin` is deprecated. Please use `panel.spacing` property
## instead
Hospitals that take care of pets
pet_hos<-subset(dat, pet == "Y")
qmplot(x, y, data = pet_hos, size = I(3),alpha = I(0.8),extent = "device", geom = "point",main = "Hospitals that take care of pets", xlab = "Longitude", ylab = "Latitude",source = "google",maptype = "roadmap",zoom = 7)
## Map from URL : http://maps.googleapis.com/maps/api/staticmap?center=40.960082,-77.71003&zoom=7&size=640x640&scale=2&maptype=roadmap&language=en-EN&sensor=false
## Warning: `panel.margin` is deprecated. Please use `panel.spacing` property
## instead
