data-raw/datacreate_bgtest.R

In ejanalysis/ejanalysis: Tools for Environmental Justice (EJ) Analysis

## code that could recreate a `bgtest` dataset

varnames = c("FIPS", "FIPS.TRACT", "FIPS.COUNTY", "FIPS.ST", "ST", "statename",
             "REGION", "OBJECTID", "pop", "povknownratio", "age25up", "hhlds",
             "builtunits", "unemployedbase", "mins", "pctmin", "lowinc", "pctlowinc",
             "lths", "pctlths", "lingiso", "pctlingiso", "under5", "pctunder5",
             "over64", "pctover64", "pre1960", "pctpre1960", "unemployed",
             "pctunemployed", "VSI.eo", "VNI.eo", "VDI.eo", "dpm", "cancer",
             "resp", "traffic.score", "proximity.npdes", "proximity.npl",
             "proximity.rmp", "proximity.tsdf", "o3", "pm", "ust", "EJ.DISPARITY.pctpre1960.eo",
             "EJ.DISPARITY.dpm.eo", "EJ.DISPARITY.cancer.eo", "EJ.DISPARITY.resp.eo",
             "EJ.DISPARITY.traffic.score.eo", "EJ.DISPARITY.proximity.npdes.eo",
             "EJ.DISPARITY.proximity.npl.eo", "EJ.DISPARITY.proximity.rmp.eo",
             "EJ.DISPARITY.proximity.tsdf.eo", "EJ.DISPARITY.o3.eo", "EJ.DISPARITY.pm.eo",
             "EJ.DISPARITY.ust.eo", "pctile.pctmin", "pctile.pctlowinc", "pctile.pctlths",
             "pctile.pctlingiso", "pctile.pctunder5", "pctile.pctover64",
             "pctile.pctpre1960", "pctile.pctunemployed", "pctile.VSI.eo",
             "pctile.dpm", "pctile.cancer", "pctile.resp", "pctile.traffic.score",
             "pctile.proximity.npdes", "pctile.proximity.npl", "pctile.proximity.rmp",
             "pctile.proximity.tsdf", "pctile.o3", "pctile.pm", "pctile.ust",
             "pctile.EJ.DISPARITY.pctpre1960.eo", "pctile.EJ.DISPARITY.dpm.eo",
             "pctile.EJ.DISPARITY.cancer.eo", "pctile.EJ.DISPARITY.resp.eo",
             "pctile.EJ.DISPARITY.traffic.score.eo", "pctile.EJ.DISPARITY.proximity.npdes.eo",
             "pctile.EJ.DISPARITY.proximity.npl.eo", "pctile.EJ.DISPARITY.proximity.rmp.eo",
             "pctile.EJ.DISPARITY.proximity.tsdf.eo", "pctile.EJ.DISPARITY.o3.eo",
             "pctile.EJ.DISPARITY.pm.eo", "pctile.EJ.DISPARITY.ust.eo", "bin.pctmin",
             "bin.pctlowinc", "bin.pctlths", "bin.pctlingiso", "bin.pctunder5",
             "bin.pctover64", "bin.pctpre1960", "bin.unemployed", "bin.VSI.eo",
             "bin.dpm", "bin.cancer", "bin.resp", "bin.traffic.score", "bin.proximity.npdes",
             "bin.proximity.npl", "bin.proximity.rmp", "bin.proximity.tsdf",
             "bin.o3", "bin.pm", "bin.ust", "bin.EJ.DISPARITY.pctpre1960.eo",
             "bin.EJ.DISPARITY.dpm.eo", "bin.EJ.DISPARITY.cancer.eo", "bin.EJ.DISPARITY.resp.eo",
             "bin.EJ.DISPARITY.traffic.score.eo", "bin.EJ.DISPARITY.proximity.npdes.eo",
             "bin.EJ.DISPARITY.proximity.npl.eo", "bin.EJ.DISPARITY.proximity.rmp.eo",
             "bin.EJ.DISPARITY.proximity.tsdf.eo", "bin.EJ.DISPARITY.o3.eo",
             "bin.EJ.DISPARITY.pm.eo", "bin.EJ.DISPARITY.ust.eo", "arealand",
             "areawater", "count.NPL", "count.TSDF", "Shape_Length", "area",
             "countyname", "flagged", "hisp", "nhwa", "nhba", "nhaiana", "nhaa",
             "nhnhpia", "nhotheralone", "nhmulti", "nonmins", "pcthisp", "pctnhwa",
             "pctnhba", "pctnhaiana", "pctnhaa", "pctnhnhpia", "pctnhotheralone",
             "pctnhmulti")

  bgtest <- ejscreen::bg22[sample(1:NROW(bg22), size = 32698), varnames]

  usethis::use_data(bgtest, overwrite = TRUE)


### obsolete/ unused:

 #    bgtest_small <- structure(list(state = structure(c(1L, 2L, 3L, 4L, 5L, 6L, 7L,
 # 8L, 10L, 11L, 12L, 13L, 14L, 15L, 16L, 17L, 18L, 19L, 20L, 21L,
 # 22L, 23L, 24L, 25L, 26L, 27L, 28L, 29L, 30L, 31L, 32L, 33L, 34L,
 # 35L, 36L, 37L, 38L, 39L, 41L, 42L, 43L, 44L, 45L, 46L, 47L, 48L,
 # 49L, 50L, 51L, 52L), .Label = c("Alabama", "Alaska", "Arizona",
 # "Arkansas", "California", "Colorado", "Connecticut", "Delaware",
 # "District of Columbia", "Florida", "Georgia", "Hawaii", "Idaho",
 # "Illinois", "Indiana", "Iowa", "Kansas", "Kentucky", "Louisiana",
 # "Maine", "Maryland", "Massachusetts", "Michigan", "Minnesota",
 # "Mississippi", "Missouri", "Montana", "Nebraska", "Nevada", "New Hampshire",
 # "New Jersey", "New Mexico", "New York", "North Carolina", "North Dakota",
 # "Ohio", "Oklahoma", "Oregon", "Pennsylvania", "Puerto Rico",
 # "Rhode Island", "South Carolina", "South Dakota", "Tennessee",
 # "Texas", "Utah", "Vermont", "Virginia", "Washington", "West Virginia",
 # "Wisconsin", "Wyoming"), class = "factor"), pcthisp = c(0.0381527239296627,
 # 0.056769492321473, 0.296826116572835, 0.0635169313105461, 0.375728960493789,
 # 0.206327251656949, 0.134422275491411, 0.0813548250199138, 0.2249082771481,
 # 0.0878682317249484, 0.0901734019211436, 0.111947738331921, 0.158094676641822,
 # 0.0599941716405598, 0.0495552961203499, 0.104741084665558, 0.0301921562004411,
 # 0.0425162900517816, 0.0129720920573869, 0.0816325860392955, 0.0960897601513277,
 # 0.0442948677533508, 0.0470583828855611, 0.0264973278249911, 0.0354626134972627,
 # 0.0292535716628734, 0.0914761950105784, 0.265451497002445, 0.0283535007142456,
 # 0.177132117215957, 0.463472498001496, 0.176607017430808, 0.0834317084560556,
 # 0.0209906647364226, 0.0307719359181436, 0.0883052970054721, 0.117261303415395,
 # 0.0568442805511265, 0.124769233546578, 0.0503041778042313, 0.0279186292931113,
 # 0.0454229079840698, 0.376044616311455, 0.129379195461843, 0.0151111594281676,
 # 0.0788112971314249, 0.111945098129999, 0.0119028512046327, 0.0589405823830593,
 # 0.0893971780534219), pop = c(3615, 365, 2212, 2110, 21198, 2541,
 # 3100, 579, 8277, 4931, 868, 813, 11197, 5313, 2861, 2280, 3387,
 # 3806, 1058, 4122, 5814, 9111, 3921, 2341, 4767, 746, 1544, 590,
 # 812, 7333, 1144, 18076, 5441, 637, 10735, 2715, 2284, 11860,
 # 931, 2816, 681, 4173, 12237, 1203, 472, 4981, 3559, 1799, 4589,
 # 376), murder = c(15.1, 11.3, 7.8, 10.1, 10.3, 6.8, 3.1, 6.2,
 # 10.7, 13.9, 6.2, 5.3, 10.3, 7.1, 2.3, 4.5, 10.6, 13.2, 2.7, 8.5,
 # 3.3, 11.1, 2.3, 12.5, 9.3, 5, 2.9, 11.5, 3.3, 5.2, 9.7, 10.9,
 # 11.1, 1.4, 7.4, 6.4, 4.2, 6.1, 2.4, 11.6, 1.7, 11, 12.2, 4.5,
 # 5.5, 9.5, 4.3, 6.7, 3, 6.9), area = c(50708, 566432, 113417,
 # 51945, 156361, 103766, 4862, 1982, 54090, 58073, 6425, 82677,
 # 55748, 36097, 55941, 81787, 39650, 44930, 30920, 9891, 7826,
 # 56817, 79289, 47296, 68995, 145587, 76483, 109889, 9027, 7521,
 # 121412, 47831, 48798, 69273, 40975, 68782, 96184, 44966, 1049,
 # 30225, 75955, 41328, 262134, 82096, 9267, 39780, 66570, 24070,
 # 54464, 97203), temp = c(62.8, 26.6, 60.3, 60.4, 59.4, 45.1, 49,
 # 55.3, 70.7, 63.5, 70, 44.4, 51.8, 51.7, 47.8, 54.3, 55.6, 66.4,
 # 41, 54.2, 47.9, 44.4, 41.2, 63.4, 54.5, 42.7, 48.8, 49.9, 43.8,
 # 52.7, 53.4, 45.4, 59, 40.4, 50.7, 59.6, 48.4, 48.8, 50.1, 62.4,
 # 45.2, 57.6, 64.8, 48.6, 42.9, 55.1, 48.3, 51.8, 43.1, 42)), .Names = c("state",
 # "pcthisp", "pop", "murder", "area", "temp"), class = "data.frame", row.names = c(NA,
 # -50L))