Hospitals provide essential, vital care to the communities they serve. There are many healthcare facilites spread out all over the country. I am setting out to create a leaflet map of the hospitals in the United States and its territories.
The data being used was compiled by the federal government and put in an interactive map. The data includes all healthcare faclities that receive some type of governemnt funding. The map helps people compare healthcare facilities by location, services offered, vacancy of beds available, wait times, etc…I downloaded it from https://hifld-geoplatform.opendata.arcgis.com/datasets/hospitals/explore?location=12.691987%2C-15.457895%2C2.32&showTable=true after being directed from the health and Human services.
hosp<-read_csv("~/Desktop/stat 4233/HW5/Hospitals.csv")
##
## ── Column specification ────────────────────────────────────────────────────────
## cols(
## .default = col_character(),
## X = col_double(),
## Y = col_double(),
## OBJECTID = col_double(),
## POPULATION = col_double(),
## LATITUDE = col_double(),
## LONGITUDE = col_double(),
## NAICS_CODE = col_double(),
## TTL_STAFF = col_double(),
## BEDS = col_double()
## )
## ℹ Use `spec()` for the full column specifications.
str(hosp)
## spec_tbl_df [7,596 × 34] (S3: spec_tbl_df/tbl_df/tbl/data.frame)
## $ X : num [1:7596] -13318895 -13226507 -13156201 -13171905 -13132075 ...
## $ Y : num [1:7596] 4346975 4049626 4031978 4041752 4037270 ...
## $ OBJECTID : num [1:7596] 1 2 3 4 5 6 7 8 9 10 ...
## $ ID : chr [1:7596] "0005793230" "0053391362" "0011190023" "0017090028" ...
## $ NAME : chr [1:7596] "CENTRAL VALLEY GENERAL HOSPITAL" "LOS ROBLES HOSPITAL & MEDICAL CENTER - EAST CAMPUS" "EAST LOS ANGELES DOCTORS HOSPITAL" "SOUTHERN CALIFORNIA HOSPITAL AT HOLLYWOOD" ...
## $ ADDRESS : chr [1:7596] "1025 NORTH DOUTY STREET" "150 VIA MERIDA" "4060 WHITTIER BOULEVARD" "6245 DE LONGPRE AVENUE" ...
## $ CITY : chr [1:7596] "HANFORD" "WESTLAKE VILAGE" "LOS ANGELES" "HOLLYWOOD" ...
## $ STATE : chr [1:7596] "CA" "CA" "CA" "CA" ...
## $ ZIP : chr [1:7596] "93230" "91362" "90023" "90028" ...
## $ ZIP4 : chr [1:7596] "NOT AVAILABLE" "NOT AVAILABLE" "NOT AVAILABLE" "NOT AVAILABLE" ...
## $ TELEPHONE : chr [1:7596] "NOT AVAILABLE" "NOT AVAILABLE" "NOT AVAILABLE" "(323) 462-2271" ...
## $ TYPE : chr [1:7596] "GENERAL ACUTE CARE" "GENERAL ACUTE CARE" "GENERAL ACUTE CARE" "GENERAL ACUTE CARE" ...
## $ STATUS : chr [1:7596] "CLOSED" "OPEN" "OPEN" "OPEN" ...
## $ POPULATION: num [1:7596] 49 62 127 100 95 172 49 101 16 78 ...
## $ COUNTY : chr [1:7596] "KINGS" "VENTURA" "LOS ANGELES" "LOS ANGELES" ...
## $ COUNTYFIPS: chr [1:7596] "06031" "06111" "06037" "06037" ...
## $ COUNTRY : chr [1:7596] "USA" "USA" "USA" "USA" ...
## $ LATITUDE : num [1:7596] 36.3 34.2 34 34.1 34.1 ...
## $ LONGITUDE : num [1:7596] -120 -119 -118 -118 -118 ...
## $ NAICS_CODE: num [1:7596] 622110 622110 622110 622110 622110 ...
## $ NAICS_DESC: chr [1:7596] "GENERAL MEDICAL AND SURGICAL HOSPITALS" "GENERAL MEDICAL AND SURGICAL HOSPITALS" "GENERAL MEDICAL AND SURGICAL HOSPITALS" "GENERAL MEDICAL AND SURGICAL HOSPITALS" ...
## $ SOURCE : chr [1:7596] "http://www.oshpd.ca.gov/hid/facility-listing.html" "http://www.oshpd.ca.gov/hid/facility-listing.html" "http://www.oshpd.ca.gov/hid/facility-listing.html" "http://www.oshpd.ca.gov/hid/facility-listing.html" ...
## $ SOURCEDATE: chr [1:7596] "2016/02/04 00:00:00" "2019/08/02 00:00:00" "2019/08/02 00:00:00" "2019/08/02 00:00:00" ...
## $ VAL_METHOD: chr [1:7596] "IMAGERY/OTHER" "IMAGERY/OTHER" "IMAGERY/OTHER" "IMAGERY/OTHER" ...
## $ VAL_DATE : chr [1:7596] "2014/02/10 00:00:00" "2014/02/10 00:00:00" "2014/02/10 00:00:00" "2014/02/10 00:00:00" ...
## $ WEBSITE : chr [1:7596] "http://www.hanfordhealth.com" "http://www.losrobleshospital.com" "http://www.elalax.com" "http://sch-hollywood.com/" ...
## $ STATE_ID : chr [1:7596] "NOT AVAILABLE" "NOT AVAILABLE" "NOT AVAILABLE" "NOT AVAILABLE" ...
## $ ALT_NAME : chr [1:7596] "NOT AVAILABLE" "NOT AVAILABLE" "NOT AVAILABLE" "HOLLYWOOD COMMUNITY HOSPITAL OF HOLLYWOOD" ...
## $ ST_FIPS : chr [1:7596] "06" "06" "06" "06" ...
## $ OWNER : chr [1:7596] "PROPRIETARY" "PROPRIETARY" "PROPRIETARY" "PROPRIETARY" ...
## $ TTL_STAFF : num [1:7596] -999 -999 -999 -999 -999 -999 -999 -999 -999 -999 ...
## $ BEDS : num [1:7596] 49 62 127 100 95 172 49 101 16 78 ...
## $ TRAUMA : chr [1:7596] "NOT AVAILABLE" "NOT AVAILABLE" "NOT AVAILABLE" "NOT AVAILABLE" ...
## $ HELIPAD : chr [1:7596] "N" "N" "N" "N" ...
## - attr(*, "spec")=
## .. cols(
## .. X = col_double(),
## .. Y = col_double(),
## .. OBJECTID = col_double(),
## .. ID = col_character(),
## .. NAME = col_character(),
## .. ADDRESS = col_character(),
## .. CITY = col_character(),
## .. STATE = col_character(),
## .. ZIP = col_character(),
## .. ZIP4 = col_character(),
## .. TELEPHONE = col_character(),
## .. TYPE = col_character(),
## .. STATUS = col_character(),
## .. POPULATION = col_double(),
## .. COUNTY = col_character(),
## .. COUNTYFIPS = col_character(),
## .. COUNTRY = col_character(),
## .. LATITUDE = col_double(),
## .. LONGITUDE = col_double(),
## .. NAICS_CODE = col_double(),
## .. NAICS_DESC = col_character(),
## .. SOURCE = col_character(),
## .. SOURCEDATE = col_character(),
## .. VAL_METHOD = col_character(),
## .. VAL_DATE = col_character(),
## .. WEBSITE = col_character(),
## .. STATE_ID = col_character(),
## .. ALT_NAME = col_character(),
## .. ST_FIPS = col_character(),
## .. OWNER = col_character(),
## .. TTL_STAFF = col_double(),
## .. BEDS = col_double(),
## .. TRAUMA = col_character(),
## .. HELIPAD = col_character()
## .. )
Hospitals <-hosp[!is.na(hosp$LATITUDE)&! is.na(hosp$LONGITUDE)&!is.na(hosp$NAME),]
Hospitals %>%
leaflet() %>%
addTiles() %>%
addMarkers(popup = hosp$LONGITUDE, clusterOptions = markerClusterOptions())
## Assuming "LONGITUDE" and "LATITUDE" are longitude and latitude, respectively
As you can see on the map, healthcare facilities outside the United States are marked as they are either located in US territories and/or they receive some type of federal funding.