The goal of addPmData is to add PM estimates to geocoded data.
More information on the development of the spatial model can be found at the st_pm_hex repository and in the scientific manuscript describing model development: Brokamp, C. A High Resolution Spatiotemporal Fine Particulate Matter Exposure Assessment Model for the Contiguous United States. Environmental Advances. In Press. 2021. (accepted preprint available online: https://doi.org/10.1016/j.envadv.2021.100155).
You can install the development version from GitHub with:
# install.packages("remotes")
remotes::install_github("geomarker-io/addPmData")
library(addPmData)
d <- tibble::tribble(
~id, ~lat, ~lon, ~start_date, ~end_date,
'55000100280', 39.2, -84.6, '2008-09-09', '2008-09-11',
'55000100281', 39.2, -84.6, '2007-08-05', '2007-08-08',
'55000100282', 39.2, -84.6, '2015-08-31', '2015-09-02')
add_pm(d)
#> matching lat/lon to h3 cells...
#> downloading PM chunk files...
#> ℹ all files already exist
#> Reading in and joining PM data...
#> # A tibble: 10 × 11
#> id lat lon start_date end_date date year h3 h3_3 pm_pred
#> <chr> <dbl> <dbl> <date> <date> <date> <dbl> <chr> <chr> <dbl>
#> 1 55000… 39.2 -84.6 2008-09-09 2008-09-11 2008-09-09 2008 882a… 832a… 8.23
#> 2 55000… 39.2 -84.6 2008-09-09 2008-09-11 2008-09-10 2008 882a… 832a… 9.58
#> 3 55000… 39.2 -84.6 2008-09-09 2008-09-11 2008-09-11 2008 882a… 832a… 13.0
#> 4 55000… 39.2 -84.6 2007-08-05 2007-08-08 2007-08-05 2007 882a… 832a… 29.3
#> 5 55000… 39.2 -84.6 2007-08-05 2007-08-08 2007-08-06 2007 882a… 832a… 28.5
#> 6 55000… 39.2 -84.6 2007-08-05 2007-08-08 2007-08-07 2007 882a… 832a… 22.3
#> 7 55000… 39.2 -84.6 2007-08-05 2007-08-08 2007-08-08 2007 882a… 832a… 20.7
#> 8 55000… 39.2 -84.6 2015-08-31 2015-09-02 2015-08-31 2015 882a… 832a… 12.7
#> 9 55000… 39.2 -84.6 2015-08-31 2015-09-02 2015-09-01 2015 882a… 832a… 17.2
#> 10 55000… 39.2 -84.6 2015-08-31 2015-09-02 2015-09-02 2015 882a… 832a… 19.4
#> # … with 1 more variable: pm_se <dbl>