time_subset
time_subset.Rmd
library(cityClimateHealth)
library(data.table)
#>
#> Attaching package: 'data.table'
#> The following object is masked from 'package:base':
#>
#> %notin%To create an exposure matrix for analysis, the time period must be specified in years, months, or days. For example, we can make an exposure matrix for Massachusetts maximum daily heat exposure between 2012 and 2015:
exposure_columns <- list(
"date" = "date",
"exposure" = "tmax_C",
"geo_unit" = "TOWN20",
"geo_unit_grp" = "COUNTY20"
)
ma_exposure_matrix <- make_exposure_matrix(
ma_exposure,
exposure_columns,
time_subset = list(year = 2012:2015)
)
#> -- NA values automatically removed
#> > grp_level == FALSE, so using geo_unit as strata
#> strata dt_by = 'day', setting strata as geo_unit:yr:mn:dow
ma_exposure_matrix
#> date TOWN20 COUNTY20 tmax_C strata
#> <IDat> <char> <char> <num> <char>
#> 1: 2012-01-01 ABINGTON PLYMOUTH 4.0592 ABINGTON:yr2012:mn01:dow01
#> 2: 2012-01-02 ABINGTON PLYMOUTH 9.6909 ABINGTON:yr2012:mn01:dow02
#> 3: 2012-01-03 ABINGTON PLYMOUTH 4.5612 ABINGTON:yr2012:mn01:dow03
#> 4: 2012-01-04 ABINGTON PLYMOUTH -4.9586 ABINGTON:yr2012:mn01:dow04
#> 5: 2012-01-05 ABINGTON PLYMOUTH -4.5164 ABINGTON:yr2012:mn01:dow05
#> ---
#> 511346: 2015-12-27 YARMOUTH BARNSTABLE 9.7547 YARMOUTH:yr2015:mn12:dow01
#> 511347: 2015-12-28 YARMOUTH BARNSTABLE 10.6210 YARMOUTH:yr2015:mn12:dow02
#> 511348: 2015-12-29 YARMOUTH BARNSTABLE -0.4467 YARMOUTH:yr2015:mn12:dow03
#> 511349: 2015-12-30 YARMOUTH BARNSTABLE -0.5062 YARMOUTH:yr2015:mn12:dow04
#> 511350: 2015-12-31 YARMOUTH BARNSTABLE 1.1314 YARMOUTH:yr2015:mn12:dow05
#> match_strata explag1 explag2 explag3 explag4 explag5
#> <char> <num> <num> <num> <num> <num>
#> 1: ABINGTON:2012-01-01 3.0964 0.1248 5.5685 12.8328 4.4747
#> 2: ABINGTON:2012-01-02 4.0592 3.0964 0.1248 5.5685 12.8328
#> 3: ABINGTON:2012-01-03 9.6909 4.0592 3.0964 0.1248 5.5685
#> 4: ABINGTON:2012-01-04 4.5612 9.6909 4.0592 3.0964 0.1248
#> 5: ABINGTON:2012-01-05 -4.9586 4.5612 9.6909 4.0592 3.0964
#> ---
#> 511346: YARMOUTH:2015-12-27 16.5056 18.9899 11.3674 13.1147 9.1770
#> 511347: YARMOUTH:2015-12-28 9.7547 16.5056 18.9899 11.3674 13.1147
#> 511348: YARMOUTH:2015-12-29 10.6210 9.7547 16.5056 18.9899 11.3674
#> 511349: YARMOUTH:2015-12-30 -0.4467 10.6210 9.7547 16.5056 18.9899
#> 511350: YARMOUTH:2015-12-31 -0.5062 -0.4467 10.6210 9.7547 16.5056We can also specify the exposure matrix for months. For example, we can create a matrix that uses time exclusively from June through August, the highest-heat months in Massachusetts. Note that months should be specified as a plain integer, not a named list.
exposure_columns <- list(
"date" = "date",
"exposure" = "tmax_C",
"geo_unit" = "TOWN20",
"geo_unit_grp" = "COUNTY20"
)
ma_exposure_matrix <- make_exposure_matrix(
ma_exposure,
exposure_columns,
time_subset = list(month = 6:8) # 6 = June, 7 = July, 8 = August
)
#> -- NA values automatically removed
#> > grp_level == FALSE, so using geo_unit as strata
#> strata dt_by = 'day', setting strata as geo_unit:yr:mn:dow
ma_exposure_matrix
#> date TOWN20 COUNTY20 tmax_C strata
#> <IDat> <char> <char> <num> <char>
#> 1: 2010-06-01 ABINGTON PLYMOUTH 28.5766 ABINGTON:yr2010:mn06:dow03
#> 2: 2010-06-02 ABINGTON PLYMOUTH 26.1360 ABINGTON:yr2010:mn06:dow04
#> 3: 2010-06-03 ABINGTON PLYMOUTH 29.0638 ABINGTON:yr2010:mn06:dow05
#> 4: 2010-06-04 ABINGTON PLYMOUTH 27.8157 ABINGTON:yr2010:mn06:dow06
#> 5: 2010-06-05 ABINGTON PLYMOUTH 29.4976 ABINGTON:yr2010:mn06:dow07
#> ---
#> 354196: 2020-08-27 YARMOUTH BARNSTABLE 23.1338 YARMOUTH:yr2020:mn08:dow05
#> 354197: 2020-08-28 YARMOUTH BARNSTABLE 18.9940 YARMOUTH:yr2020:mn08:dow06
#> 354198: 2020-08-29 YARMOUTH BARNSTABLE 26.0081 YARMOUTH:yr2020:mn08:dow07
#> 354199: 2020-08-30 YARMOUTH BARNSTABLE 24.1151 YARMOUTH:yr2020:mn08:dow01
#> 354200: 2020-08-31 YARMOUTH BARNSTABLE 22.4953 YARMOUTH:yr2020:mn08:dow02
#> match_strata explag1 explag2 explag3 explag4 explag5
#> <char> <num> <num> <num> <num> <num>
#> 1: ABINGTON:2010-06-01 26.9500 25.3234 26.3356 26.0529 35.4982
#> 2: ABINGTON:2010-06-02 28.5766 26.9500 25.3234 26.3356 26.0529
#> 3: ABINGTON:2010-06-03 26.1360 28.5766 26.9500 25.3234 26.3356
#> 4: ABINGTON:2010-06-04 29.0638 26.1360 28.5766 26.9500 25.3234
#> 5: ABINGTON:2010-06-05 27.8157 29.0638 26.1360 28.5766 26.9500
#> ---
#> 354196: YARMOUTH:2020-08-27 29.2008 31.6397 32.1062 30.4607 30.6083
#> 354197: YARMOUTH:2020-08-28 23.1338 29.2008 31.6397 32.1062 30.4607
#> 354198: YARMOUTH:2020-08-29 18.9940 23.1338 29.2008 31.6397 32.1062
#> 354199: YARMOUTH:2020-08-30 26.0081 18.9940 23.1338 29.2008 31.6397
#> 354200: YARMOUTH:2020-08-31 24.1151 26.0081 18.9940 23.1338 29.2008