Northen Virginia 2015 Family Income By County

library(leaflet)
set.seed(2018-01-01)

fairfax_cities <- data.frame(name = c("Fairfax", "Arlington", "Fauquier", "Loudoun", "Prince William", "Spotsylvania", "Stafford", "Alexandria"),
income = c(129538, 143843, 107479, 141487, 108243, 84401, 103705, 112337), 
 lat = c(38.846224, 38.8816211, 38.720879, 39.076767, 38.783894, 38.204165, 38.422069, 38.804836),
 lng = c(-77.306373, -77.090981,-77.790978, -77.653613, -77.607786, -77.607786, -77.408309, -77.046921))
 
fairfax_cities %>%
leaflet() %>%
addTiles() %>%
addCircles(weight = 3, radius = sqrt(fairfax_cities$income) * 30)

## Assuming 'lng' and 'lat' are longitude and latitude, respectively