library(data.table) #fread
library(sp) #used by rgdal
library(rgdal) #readOGR
library(ggplot2) #fortify
library(plyr) #join
library(dplyr) #left_join
library(scales) #percent
library(viridisLite)
library(viridis)
library(openintro) #state2abbr
library(rgeos) #gCentroid
library(cdlTools) #fips
library(ggmap) #theme_nothing

ARIS 2008

ARIS (American Religious Identification Survey)

The Association of Religion Data Archives (ARDA)

American Religious Identification Survey, 2008

The American Religious Identification Survey (ARIS) 2008 is the third in a landmark time series of large, nationally representative surveys that track changes in the religious loyalties of the U.S. adult population within the 28 contiguous states from 1990 to 2008. The 2001 and 2008 surveys are replicas of the 1990 survey, and are led by the same academic research team using an identical methodology of random-digit-dialed telephone interviews (RDD) and the same unprompted, open-ended key question “What is your religion, if any?” Interviewers did not prompt or offer a suggested list of potential answers. Moreover, the self-description of respondents was not based on whether established religious bodies or institutions considered them to be members. To the contrary, the surveys sought to determine whether the respondents regarded themselves as adherents of a religious community. The surveys tap subjective, rather than objective, standards of religious identification. The value of this unique series of national surveys, which allows scientific monitoring of change over time, has been recognized by the U.S. Bureau of the Census The Bureau itself is constitutionally precluded from such an inquiry into religion, and so has incorporated NSRI/ARIS findings into its official publication the Statistical Abstract of the United States since 2003.

Cases: 54,461
Variables: 33
Weight Variable: POPWGHT
Date Collected: February - November 2008

cf. Religious Population by State

ARIS2008 <- fread("ARIS2008.csv")
str(ARIS2008)
Classes ‘data.table’ and 'data.frame':  54461 obs. of  38 variables:
 $ CASEID    : int  102 166 184 199 228 255 283 285 297 316 ...
 $ WEEK      : int  203 203 203 203 203 203 203 203 203 203 ...
 $ POPWGHT   : num  0.412 1.884 0.915 0.927 0.062 ...
 $ STATE     : chr  "GA" "TX" "AZ" "CA" ...
 $ DMA       : int  524 709 753 807 524 501 537 521 810 807 ...
 $ NEILSEN   : chr  "D" "D" "B" "A" ...
 $ COUNTYC   : int  137 401 15 81 59 25 19 5 77 85 ...
 $ ADI       : int  197 323 275 65 197 9 357 47 339 65 ...
 $ METRO     : chr  "Non-Metro" "Non-Metro" "Suburban (Metro)" "Suburban (Metro)" ...
 $ REGION    : chr  "South" "South" "West" "West" ...
 $ DIVISION  : chr  "South Atlantic" "West South Central" "Mountain" "Pacific" ...
 $ OWNHOME   : chr  "Rented" "Rented" "Rented" "Owned" ...
 $ MARITAL   : chr  "Single, living with a partner" "Single, living with a partner" "Married" "Married" ...
 $ HEADHH    : chr  "Don't know" "No" "No" "Yes" ...
 $ EMPLOY    : chr  "Temporarily unemployed" "Other not employed" "Housewife" "Full-time" ...
 $ HHCOMP    : chr  "Six" "Refused" "Five" "Four" ...
 $ ADULTS    : chr  "Two" "" "Two" "Two" ...
 $ ADULTM    : chr  "One" "" "One" "One" ...
 $ ADULTF    : chr  "One" "" "One" "One" ...
 $ KIDS1217  : chr  "Three" "" "None" "Two" ...
 $ KIDS611   : chr  "None" "" "Three" "None" ...
 $ KIDS06    : chr  "One" "" "None" "None" ...
 $ AGE       : chr  "35" "19" "27" "40" ...
 $ EDUC      : chr  "Less than high school graduate" "Refused" "High school graduate" "High school graduate" ...
 $ INCOME    : chr  "Don't know" "Refused" "Less than $10,000" "$75,000 but less than $100,000" ...
 $ RACE      : chr  "Black Hispanic" "White Hispanic" "White Hispanic" "Unspecified" ...
 $ POLAFFIL  : chr  "Refused" "Refused" "An Independent" "A Democrat" ...
 $ REGVOTE   : chr  "No" "Refused" "No" "Yes" ...
 $ OTHTEL    : chr  "None" "Refused" "None" "None" ...
 $ SEX       : chr  "Male" "Female" "Female" "Female" ...
 $ RELIG     : chr  "Don't know" "Catholic, Roman Catholic" "Apostolic/New Apostolic" "Catholic, Roman Catholic" ...
 $ RELIG2    : chr  "" "" "" "" ...
 $ RELIG3    : chr  "" "" "" "" ...
 $ I_SEX     : chr  "Male" "Female" "Female" "Female" ...
 $ I_MARITAL : chr  "Single" "Single" "Married" "Married" ...
 $ I_EDUC    : chr  "Less than high school graduate" "" "High school graduate" "High school graduate" ...
 $ I_PARTYID : chr  "" "" "An Independent or other affiliation" "A Democrat" ...
 $ I_RELIGION: chr  "" "Catholic" "Other religion" "Catholic" ...
 - attr(*, ".internal.selfref")=<externalptr>

Religions by Frequencies

major <- names(head(sort(prop.table(table(ARIS2008$RELIG)), decreasing=TRUE), 17))

major
 [1] "Catholic, Roman Catholic"                          
 [2] "Christian"                                         
 [3] "No religion/none"                                  
 [4] "Baptist"                                           
 [5] "Protestant"                                        
 [6] "Refused"                                           
 [7] "Lutheran"                                          
 [8] "Methodist/Wesleyan"                                
 [9] "Presbyterian"                                      
[10] "Jewish, Judaism"                                   
[11] "Mormon/Church of Jesus Christ of Latter-day Saints"
[12] "Pentecostal"                                       
[13] "Non-denominational Christian"                      
[14] "United Methodist"                                  
[15] "Episcopalian/Anglican"                             
[16] "Don't know"                                        
[17] "Southern Baptist"                                  
ARIS2008_major <- ARIS2008[RELIG %in% major,]

nrow(ARIS2008_major)
[1] 50287
ggplot(ARIS2008_major, aes(RELIG)) +
  geom_bar(mapping = aes(x = RELIG, y = ..prop.., group = 1), stat = "count") +
  scale_y_continuous(labels = scales::percent) +
  coord_flip()+theme(axis.text = element_text(size = 7))

Contingency Tables

prop.table(table(ARIS2008_major$RELIG, ARIS2008_major$POLAFFIL), margin=1)
                                                    
                                                      A Democrat A Republican An Independent
  Baptist                                            0.409696970  0.284155844    0.230129870
  Catholic, Roman Catholic                           0.349144171  0.262624542    0.302733543
  Christian                                          0.254428537  0.369570493    0.302353798
  Don't know                                         0.300601202  0.140280561    0.364729459
  Episcopalian/Anglican                              0.336092715  0.324503311    0.286423841
  Jewish, Judaism                                    0.527412281  0.143640351    0.279605263
  Lutheran                                           0.300906056  0.346208870    0.271816881
  Methodist/Wesleyan                                 0.345850999  0.336765597    0.252574197
  Mormon/Church of Jesus Christ of Latter-day Saints 0.117496807  0.549169860    0.263090677
  No religion/none                                   0.345891966  0.142532910    0.410198215
  Non-denominational Christian                       0.273087071  0.306068602    0.329815303
  Pentecostal                                        0.373056995  0.278497409    0.237046632
  Presbyterian                                       0.287575150  0.405811623    0.256513026
  Protestant                                         0.266201009  0.383391541    0.289483896
  Refused                                            0.214257109  0.136964357    0.299959952
  Southern Baptist                                   0.295711061  0.415349887    0.250564334
  United Methodist                                   0.333333333  0.344000000    0.264000000
                                                    
                                                      Don't know       Other     Refused
  Baptist                                            0.036190476 0.009696970 0.030129870
  Catholic, Roman Catholic                           0.034403474 0.013625138 0.037469131
  Christian                                          0.029119146 0.016622179 0.027905848
  Don't know                                         0.112224449 0.016032064 0.066132265
  Episcopalian/Anglican                              0.013245033 0.014900662 0.024834437
  Jewish, Judaism                                    0.010964912 0.014254386 0.024122807
  Lutheran                                           0.033381021 0.012398665 0.035288507
  Methodist/Wesleyan                                 0.031496063 0.004239855 0.029073289
  Mormon/Church of Jesus Christ of Latter-day Saints 0.029374202 0.014048531 0.026819923
  No religion/none                                   0.044484793 0.018610985 0.038281132
  Non-denominational Christian                       0.032981530 0.022427441 0.035620053
  Pentecostal                                        0.050518135 0.014248705 0.046632124
  Presbyterian                                       0.016032064 0.013026052 0.021042084
  Protestant                                         0.022700815 0.009507179 0.028715561
  Refused                                            0.022426912 0.019223068 0.307168602
  Southern Baptist                                   0.015801354 0.011286682 0.011286682
  United Methodist                                   0.024000000 0.008000000 0.026666667
prop.table(table(ARIS2008_major$RELIG, ARIS2008_major$REGVOTE), margin=1)
                                                    
                                                      Don't know          No     Refused         Yes
  Baptist                                            0.004675325 0.128658009 0.007099567 0.859567100
  Catholic, Roman Catholic                           0.003235970 0.133866985 0.014221238 0.848675807
  Christian                                          0.005095851 0.140499879 0.007037127 0.847367144
  Don't know                                         0.018036072 0.214428858 0.018036072 0.749498998
  Episcopalian/Anglican                              0.003311258 0.086092715 0.003311258 0.907284768
  Jewish, Judaism                                    0.008771930 0.081140351 0.006578947 0.903508772
  Lutheran                                           0.004768717 0.090605627 0.007153076 0.897472580
  Methodist/Wesleyan                                 0.003634161 0.083585706 0.008479709 0.904300424
  Mormon/Church of Jesus Christ of Latter-day Saints 0.002554278 0.108556833 0.003831418 0.885057471
  No religion/none                                   0.005144500 0.188076865 0.008473294 0.798305341
  Non-denominational Christian                       0.003957784 0.160949868 0.003957784 0.831134565
  Pentecostal                                        0.010362694 0.176165803 0.007772021 0.805699482
  Presbyterian                                       0.000000000 0.067134269 0.007014028 0.925851703
  Protestant                                         0.002716337 0.077609624 0.005820722 0.913853318
  Refused                                            0.004405286 0.082098518 0.148177813 0.765318382
  Southern Baptist                                   0.004514673 0.110609481 0.004514673 0.880361174
  United Methodist                                   0.002666667 0.078666667 0.002666667 0.916000000

State FIPS Codes

sort(unique(ARIS2008$STATE))
 [1] "AL" "AR" "AZ" "CA" "CO" "CT" "DC" "DE" "FL" "GA" "IA" "ID" "IL" "IN" "KS" "KY" "LA" "MA" "MD"
[20] "ME" "MI" "MN" "MO" "MS" "MT" "NC" "ND" "NE" "NH" "NJ" "NM" "NV" "NY" "OH" "OK" "OR" "PA" "RI"
[39] "SC" "SD" "TN" "TX" "UT" "VA" "VT" "WA" "WI" "WV" "WY"

State FIPS Codes

lower_48 <- c("01", "04", "05", "06", "08", "09", "10", "11", "12", "13", "16", "17", "18", "19", "20", "21", "22", "23", "24", "25", "26", "27", "28", "29", "30", "31", "32", "33", "34", "35", "36", "37", "38", "39", "40", "41", "42", "44", "45", "46", "47", "48", "49", "50", "51", "53", "54", "55", "56")

States Shapefile

United States Census Bureau, Cartographic Boundary Files - Shapefile

cb_2018_us_state_500k.zip

states.sh <- readOGR("C:/Users/14165/Desktop/Shapefiles/cb_2018_us_state_500k", "cb_2018_us_state_500k", stringsAsFactors=FALSE)
OGR data source with driver: ESRI Shapefile 
Source: "C:\Users\14165\Desktop\Shapefiles\cb_2018_us_state_500k", layer: "cb_2018_us_state_500k"
with 56 features
It has 9 fields
Integer64 fields read as strings:  ALAND AWATER 
str(states.sh@data)
'data.frame':   56 obs. of  9 variables:
 $ STATEFP : chr  "28" "37" "40" "51" ...
 $ STATENS : chr  "01779790" "01027616" "01102857" "01779803" ...
 $ AFFGEOID: chr  "0400000US28" "0400000US37" "0400000US40" "0400000US51" ...
 $ GEOID   : chr  "28" "37" "40" "51" ...
 $ STUSPS  : chr  "MS" "NC" "OK" "VA" ...
 $ NAME    : chr  "Mississippi" "North Carolina" "Oklahoma" "Virginia" ...
 $ LSAD    : chr  "00" "00" "00" "00" ...
 $ ALAND   : chr  "121533519481" "125923656064" "177662925723" "102257717110" ...
 $ AWATER  : chr  "3926919758" "13466071395" "3374587997" "8528531774" ...
states_lower48.sh <- states.sh[states.sh@data$GEOID %in% lower_48,]
str(states_lower48.sh@data)
'data.frame':   49 obs. of  9 variables:
 $ STATEFP : chr  "28" "37" "40" "51" ...
 $ STATENS : chr  "01779790" "01027616" "01102857" "01779803" ...
 $ AFFGEOID: chr  "0400000US28" "0400000US37" "0400000US40" "0400000US51" ...
 $ GEOID   : chr  "28" "37" "40" "51" ...
 $ STUSPS  : chr  "MS" "NC" "OK" "VA" ...
 $ NAME    : chr  "Mississippi" "North Carolina" "Oklahoma" "Virginia" ...
 $ LSAD    : chr  "00" "00" "00" "00" ...
 $ ALAND   : chr  "121533519481" "125923656064" "177662925723" "102257717110" ...
 $ AWATER  : chr  "3926919758" "13466071395" "3374587997" "8528531774" ...

Add “id” column

states_lower48.sh@data$id <- as.integer(states_lower48.sh@data$GEOID)
sort(states_lower48.sh@data$id)
 [1]  1  4  5  6  8  9 10 11 12 13 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37
[33] 38 39 40 41 42 44 45 46 47 48 49 50 51 53 54 55 56
states_lower48.sh@data$STATE <- as.character(states_lower48.sh@data$STUSPS)

Shapefile processing

states_lower48.points = fortify(states_lower48.sh, region = "id")

states_lower48.df = join(states_lower48.points, states_lower48.sh@data, by = "id")

State Religious Ratios and Choropleths

Totals2008 <- ARIS2008[, .(.N), by = .(STATE)]

Totals2008 <- Totals2008[order(STATE)]

Catholics

Catholics2008 <- ARIS2008[RELIG == "Catholic, Roman Catholic", .(.N), by = .(STATE)]

Catholics2008 <- Catholics2008[order(STATE)]

Catholics2008$ratio <- Catholics2008$N/Totals2008$N

Catholics2008_lower48 <- merge(states_lower48.df, Catholics2008, by = "STATE")

Catholics2008_lower48 <- Catholics2008_lower48[order(Catholics2008_lower48$STATE), ]

Make and plot choropleth

ggplot() +
  geom_polygon(data = Catholics2008_lower48, 
               aes(x = long, y = lat, group = group, fill = ratio), 
               color = "white", size = 0.1) + 
  coord_map() + 
  scale_fill_viridis(option="magma", labels=percent, name="Catholics") +
   labs(title="Ratios of Catholics by State, ARIS2008")

Baptists

Baptists2008 <- ARIS2008[RELIG == "Baptist", .(.N), by = .(STATE)]

Baptists2008 <- Baptists2008[order(STATE)]

Baptists2008$ratio <- Baptists2008$N/Totals2008$N

Baptists2008_lower48 <- merge(states_lower48.df, Baptists2008, by = "STATE")

Baptists2008_lower48 <- Baptists2008_lower48[order(Baptists2008_lower48$STATE), ]

Make and plot choropleth

ggplot() +
  geom_polygon(data = Baptists2008_lower48, 
               aes(x = long, y = lat, group = group, fill = ratio), 
               color = "white", size = 0.1) + 
  coord_map() + 
  scale_fill_viridis(option="magma", labels=percent, name="Baptists") +
   labs(title="Ratios of Baptists by State, ARIS 2008")

No Religion

Noreligion2008 <- ARIS2008[RELIG == "No religion/none", .(.N), by = .(STATE)]

Noreligion2008 <- Noreligion2008[order(STATE)]

Noreligion2008$ratio <- Noreligion2008$N/Totals2008$N

Noreligion2008_lower48 <- merge(states_lower48.df, Noreligion2008, by = "STATE")

Noreligion2008_lower48 <- Noreligion2008_lower48[order(Noreligion2008_lower48$STATE), ]

Make and plot choropleth

ggplot() +
  geom_polygon(data = Noreligion2008_lower48, 
               aes(x = long, y = lat, group = group, fill = ratio), 
               color = "white", size = 0.1) + 
  coord_map() + 
  scale_fill_viridis(option="magma", labels=percent, name="No Religion") +
   labs(title="Ratios of No Religion by State, ARIS 2008")

ARIS 2001

Download ARIS 2001 Data

ARIS2001SPSS.zip

ARIS2001 <- fread("ARIS2001.csv")

ARIS 2001 Codebook

ARIS2001 <- ARIS2001[, .(POPWGHT = wgt_pop, STATE = state, COUNTYFIPS = fips, ADI = adi, DMA = dma, MSA = msa, METRO = metro, NIELSEN = neilsen, DIVISION = cen_div, REGION = region, SEX = gender, OWNHOME = dem2, MARITAL = dem3, EMPLOY = dem4, INCOME = newdem9, AGE = dem7, EDUC = dem8, HISPANIC = dem10, RACE = dem11, RELIG = jq1tot, RELIG_SPOUSE = jq2tot, AGECAT = agecat, POLAFFIL = z11a, REGVOTE = z11b)]
ARIS2001$STATE <- state2abbr(ARIS2001$STATE)
major_2001 <- names(head(sort(prop.table(table(ARIS2001$RELIG)), decreasing=TRUE), 17))

major_2001
 [1] "Catholic"                                "Baptist"                                
 [3] "None"                                    "Methodist/Wesleyan"                     
 [5] "Christian"                               "Refused"                                
 [7] "Lutheran"                                "Presbyterian"                           
 [9] "Protestant"                              "Pentecostal/Charismatics"               
[11] "Southern Baptist"                        "Episcopalian"                           
[13] "Mormon/Latter Day Saints"                "Churches of Christ & Christian Churches"
[15] "Non-denominational"                      "Jewish/Judaism"                         
[17] "United Church of Christ/Congregational" 
ARIS2001_major <- ARIS2001[RELIG %in% major_2001,]

nrow(ARIS2001_major)
[1] 46021
ggplot(ARIS2001_major, aes(RELIG)) +
  geom_bar(mapping = aes(x = RELIG, y = ..prop.., group = 1), stat = "count") +
  scale_y_continuous(labels = scales::percent) +
  coord_flip()+theme(axis.text = element_text(size = 7))

Totals2001 <- ARIS2001[, .(.N), by = .(STATE)]

Totals2001 <- Totals2001[order(STATE)]

Catholics

Catholics2001 <- ARIS2001[RELIG == "Catholic", .(.N), by = .(STATE)]

Catholics2001 <- Catholics2001[order(STATE)]

Catholics2001$ratio <- Catholics2001$N/Totals2001$N

Catholics2001_lower48 <- merge(states_lower48.df, Catholics2001, by = "STATE")

Catholics2001_lower48 <- Catholics2001_lower48[order(Catholics2001_lower48$STATE), ]

Make and plot choropleth

ggplot() +
  geom_polygon(data = Catholics2001_lower48, 
               aes(x = long, y = lat, group = group, fill = ratio), 
               color = "white", size = 0.1) + 
  coord_map() + 
  scale_fill_viridis(option="magma", labels=percent, name="Catholics") +
   labs(title="Ratios of Catholics by State, ARIS 2001")

Baptists

Baptists2001 <- ARIS2001[RELIG == "Baptist", .(.N), by = .(STATE)]

Baptists2001 <- Baptists2001[order(STATE)]

Baptists2001$ratio <- Baptists2001$N/Totals2001$N

Baptists2001_lower48 <- merge(states_lower48.df, Baptists2001, by = "STATE")

Baptists2001_lower48 <- Baptists2001_lower48[order(Baptists2001_lower48$STATE), ]

ggplot() +
  geom_polygon(data = Baptists2001_lower48, 
               aes(x = long, y = lat, group = group, fill = ratio), 
               color = "white", size = 0.1) + 
  coord_map() + 
  scale_fill_viridis(option="magma", labels=percent, name="Baptists") +
   labs(title="Ratios of Baptists by State, ARIS 2001")

No Religion

Noreligion2001 <- ARIS2001[RELIG == "None", .(.N), by = .(STATE)]

Noreligion2001 <- Noreligion2001[order(STATE)]

Noreligion2001$ratio <- Noreligion2001$N/Totals2001$N

Noreligion2001_lower48 <- merge(states_lower48.df, Noreligion2001, by = "STATE")

Noreligion2001_lower48 <- Noreligion2001_lower48[order(Noreligion2001_lower48$STATE), ]

ggplot() +
  geom_polygon(data = Noreligion2001_lower48, 
               aes(x = long, y = lat, group = group, fill = ratio), 
               color = "white", size = 0.1) + 
  coord_map() + 
  scale_fill_viridis(option="magma", labels=percent, name="No Religion") +
   labs(title="Ratios of No Religion by State, ARIS 2001")

Change from 2001 to 2008

Catholics

Catholics2008$change <- (Catholics2008$ratio-Catholics2001$ratio)/Catholics2001$ratio
Catholics2008change_lower48 <- merge(states_lower48.df, Catholics2008, by = "STATE")

Catholics2008change_lower48 <- Catholics2008change_lower48[order(Catholics2008change_lower48$STATE), ]

Make and plot choropleth

ggplot() +
  geom_polygon(data = Catholics2008change_lower48, 
               aes(x = long, y = lat, group = group, fill = change), 
               color = "white", size = 0.1) + 
  coord_map() + 
  scale_fill_viridis(option="magma", labels=percent, name="Change") +
   labs(title="(ARIS2008-ARIS2001)/ARIS2001, Catholics")

Baptists

Baptists2008$change <- (Baptists2008$ratio-Baptists2001$ratio)/Baptists2001$ratio

Baptists2008change_lower48 <- merge(states_lower48.df, Baptists2008, by = "STATE")

Baptists2008change_lower48 <- Baptists2008change_lower48[order(Baptists2008change_lower48$STATE), ]

ggplot() +
  geom_polygon(data = Baptists2008change_lower48, 
               aes(x = long, y = lat, group = group, fill = change), 
               color = "white", size = 0.1) + 
  coord_map() + 
  scale_fill_viridis(option="magma", labels=percent, name="Change") +
   labs(title="(ARIS2008-ARIS2001)/ARIS2001, Baptists")

No Religion

Noreligion2008$change <- (Noreligion2008$ratio-Noreligion2001$ratio)/Noreligion2001$ratio

Noreligion2008change_lower48 <- merge(states_lower48.df, Noreligion2008, by = "STATE")

Noreligion2008change_lower48 <- Noreligion2008change_lower48[order(Noreligion2008change_lower48$STATE), ]

ggplot() +
  geom_polygon(data = Noreligion2008change_lower48, 
               aes(x = long, y = lat, group = group, fill = change), 
               color = "white", size = 0.1) + 
  coord_map() + 
  scale_fill_viridis(option="magma", labels=percent, name="Change") +
   labs(title="(ARIS2008-ARIS2001)/ARIS2001, No Religion")

ARIS 2008 Density Maps

State and County FIPS Codes

State FIPS Codes

County FIPS Codes

ARIS2008_counties <- ARIS2008[,.(COUNTYFIPS = COUNTYC, STATE, RELIG)]

ARIS2008_counties$STATEFIPS <- fips(ARIS2008_counties$STATE, to = "FIPS")

Counties Shapefile

United States Census Bureau, Cartographic Boundary Files - Shapefile

cb_2018_us_county_500k.zip

counties.sh <- readOGR("C:/Users/14165/Desktop/Shapefiles/cb_2018_us_county_500k", "cb_2018_us_county_500k", stringsAsFactors=FALSE)
OGR data source with driver: ESRI Shapefile 
Source: "C:\Users\14165\Desktop\Shapefiles\cb_2018_us_county_500k", layer: "cb_2018_us_county_500k"
with 3233 features
It has 9 fields
Integer64 fields read as strings:  ALAND AWATER 
str(counties.sh@data)
'data.frame':   3233 obs. of  9 variables:
 $ STATEFP : chr  "21" "21" "21" "21" ...
 $ COUNTYFP: chr  "007" "017" "031" "065" ...
 $ COUNTYNS: chr  "00516850" "00516855" "00516862" "00516879" ...
 $ AFFGEOID: chr  "0500000US21007" "0500000US21017" "0500000US21031" "0500000US21065" ...
 $ GEOID   : chr  "21007" "21017" "21031" "21065" ...
 $ NAME    : chr  "Ballard" "Bourbon" "Butler" "Estill" ...
 $ LSAD    : chr  "06" "06" "06" "06" ...
 $ ALAND   : chr  "639387454" "750439351" "1103571974" "655509930" ...
 $ AWATER  : chr  "69473325" "4829777" "13943044" "6516335" ...
counties_lower48.sh <- counties.sh[counties.sh@data$STATEFP %in% lower_48,]
counties_lower48.sh@data$STATEFIPS <- as.integer(counties_lower48.sh@data$STATEFP)

counties_lower48.sh@data$COUNTYFIPS <- as.integer(counties_lower48.sh@data$COUNTYFP)
counties_lower48.centroids  <- as.data.table(gCentroid(counties_lower48.sh, byid = TRUE))

counties_lower48.centroids$STATEFIPS <- counties_lower48.sh@data$STATEFIPS

counties_lower48.centroids$COUNTYFIPS <- counties_lower48.sh@data$COUNTYFIPS
str(counties_lower48.centroids)
Classes ‘data.table’ and 'data.frame':  3108 obs. of  4 variables:
 $ x         : num  -89 -84.2 -86.7 -84 -83.7 ...
 $ y         : num  37.1 38.2 37.2 37.7 38.4 ...
 $ STATEFIPS : int  21 21 21 21 21 21 21 21 21 21 ...
 $ COUNTYFIPS: int  7 17 31 65 69 93 99 131 151 155 ...
 - attr(*, ".internal.selfref")=<externalptr> 
str(ARIS2008_counties)
Classes ‘data.table’ and 'data.frame':  54461 obs. of  4 variables:
 $ COUNTYFIPS: int  137 401 15 81 59 25 19 5 77 85 ...
 $ STATE     : chr  "GA" "TX" "AZ" "CA" ...
 $ RELIG     : chr  "Don't know" "Catholic, Roman Catholic" "Apostolic/New Apostolic" "Catholic, Roman Catholic" ...
 $ STATEFIPS : int  13 48 4 6 13 34 23 25 53 6 ...
 - attr(*, ".internal.selfref")=<externalptr> 
ARIS2008_XY <- left_join(ARIS2008_counties, counties_lower48.centroids, by = c("STATEFIPS", "COUNTYFIPS"))
ARIS2008_XY$COUNTYFIPS <- NULL

ARIS2008_XY$STATE <- NULL

ARIS2008_XY$STATEFIPS <- NULL

str(ARIS2008_XY)
'data.frame':   54461 obs. of  3 variables:
 $ RELIG: chr  "Don't know" "Catholic, Roman Catholic" "Apostolic/New Apostolic" "Catholic, Roman Catholic" ...
 $ x    : num  -83.5 -94.8 -113.8 -122.3 -83.4 ...
 $ y    : num  34.6 32.1 35.7 37.4 34 ...
Catholics2008_XY <- ARIS2008_XY[ARIS2008_XY$RELIG == "Catholic, Roman Catholic",]

Catholics2008_XY$RELIG <- NULL
Baptists2008_XY <- ARIS2008_XY[ARIS2008_XY$RELIG == "Baptist",]

Baptists2008_XY$RELIG <- NULL
Noreligion2008_XY <- ARIS2008_XY[ARIS2008_XY$RELIG == "No religion/none",]

Noreligion2008_XY$RELIG <- NULL

Catholic Density Map

ggplot() + 
  geom_polygon(data = states_lower48.df, aes(x = long, y = lat, group = group),
          color = 'black', size = 1, fill=NA) +
  stat_density2d(data=Catholics2008_XY, aes(x=x, y=y, fill=..level..), alpha=0.2, geom = 'polygon', colour = 'black', contour=TRUE) +
  scale_fill_continuous(low="yellow", high="red")+
  theme_nothing(legend = TRUE) +   
  labs(title="Catholic Denity Map, ARIS 2008")

Baptist Density Map

ggplot() + 
  geom_polygon(data = states_lower48.df, aes(x = long, y = lat, group = group),
          color = 'black', size = 1, fill=NA) +
  stat_density2d(data=Baptists2008_XY, aes(x=x, y=y, fill=..level..), alpha=0.2, geom = 'polygon', colour = 'black', contour=TRUE) +
  scale_fill_continuous(low="yellow",high="red")+
  theme_nothing(legend = TRUE) +  
  labs(title="Baptist Density Map, ARIS 2008")

No Religion Density Map

ggplot() + 
  geom_polygon(data = states_lower48.df, aes(x = long, y = lat, group = group),
          color = 'black', size = 1, fill=NA) +
  stat_density2d(data=Noreligion2008_XY, aes(x=x, y=y, fill=..level..), alpha=0.2, geom = 'polygon', colour = 'black', contour=TRUE) +
  scale_fill_continuous(low="yellow",high="red")+
  theme_nothing(legend = TRUE) +  
  labs(title="No Religion Density Map, ARIS 2008")

LS0tDQp0aXRsZTogIkFtZXJpY2FuIFJlbGlnaW91cyBJZGVudGlmaWNhdGlvbiBTdXJ2ZXkgKEFSSVMpIERhdGFzZXRzLCAyMDA4IGFuZCAyMDAxIg0Kb3V0cHV0OiANCiAgaHRtbF9ub3RlYm9vazogDQogICAgdG9jOiB5ZXMNCi0tLQ0KDQoqIEpvcmRhbiBCZWxsDQoqIEp1bHkgMTAsIDIwMTkNCiogPGh0dHBzOi8vam9yZGFuYmVsbDIzNTcuZ2l0aHViLmlvL0FSSVMubmIuaHRtbD4NCg0KDQpgYGB7cn0NCmxpYnJhcnkoZGF0YS50YWJsZSkgI2ZyZWFkDQpsaWJyYXJ5KHNwKSAjdXNlZCBieSByZ2RhbA0KbGlicmFyeShyZ2RhbCkgI3JlYWRPR1INCmxpYnJhcnkoZ2dwbG90MikgI2ZvcnRpZnkNCmxpYnJhcnkocGx5cikgI2pvaW4NCmxpYnJhcnkoZHBseXIpICNsZWZ0X2pvaW4NCmxpYnJhcnkoc2NhbGVzKSAjcGVyY2VudA0KbGlicmFyeSh2aXJpZGlzTGl0ZSkNCmxpYnJhcnkodmlyaWRpcykNCmxpYnJhcnkob3BlbmludHJvKSAjc3RhdGUyYWJicg0KbGlicmFyeShyZ2VvcykgI2dDZW50cm9pZA0KbGlicmFyeShjZGxUb29scykgI2ZpcHMNCmxpYnJhcnkoZ2dtYXApICN0aGVtZV9ub3RoaW5nDQpgYGANCg0KIyBBUklTIDIwMDgNCg0KW0FSSVMgKEFtZXJpY2FuIFJlbGlnaW91cyBJZGVudGlmaWNhdGlvbiBTdXJ2ZXkpXShodHRwczovL2NvbW1vbnMudHJpbmNvbGwuZWR1L2FyaXMvKQ0KDQpbVGhlIEFzc29jaWF0aW9uIG9mIFJlbGlnaW9uIERhdGEgQXJjaGl2ZXMgKEFSREEpXShodHRwOi8vd3d3LnRoZWFyZGEuY29tL0FyY2hpdmUvYnJvd3NlLmFzcCkNCg0KW0FtZXJpY2FuIFJlbGlnaW91cyBJZGVudGlmaWNhdGlvbiBTdXJ2ZXksIDIwMDhdKGh0dHA6Ly93d3cudGhlYXJkYS5jb20vQXJjaGl2ZS9GaWxlcy9EZXNjcmlwdGlvbnMvQVJJUzIwMDguYXNwKQ0KDQo+IFRoZSBBbWVyaWNhbiBSZWxpZ2lvdXMgSWRlbnRpZmljYXRpb24gU3VydmV5IChBUklTKSAyMDA4IGlzIHRoZSB0aGlyZCBpbiBhIGxhbmRtYXJrIHRpbWUgc2VyaWVzIG9mIGxhcmdlLCBuYXRpb25hbGx5IHJlcHJlc2VudGF0aXZlIHN1cnZleXMgdGhhdCB0cmFjayBjaGFuZ2VzIGluIHRoZSByZWxpZ2lvdXMgbG95YWx0aWVzIG9mIHRoZSBVLlMuIGFkdWx0IHBvcHVsYXRpb24gd2l0aGluIHRoZSAyOCBjb250aWd1b3VzIHN0YXRlcyBmcm9tIDE5OTAgdG8gMjAwOC4gVGhlIDIwMDEgYW5kIDIwMDggc3VydmV5cyBhcmUgcmVwbGljYXMgb2YgdGhlIDE5OTAgc3VydmV5LCBhbmQgYXJlIGxlZCBieSB0aGUgc2FtZSBhY2FkZW1pYyByZXNlYXJjaCB0ZWFtIHVzaW5nIGFuIGlkZW50aWNhbCBtZXRob2RvbG9neSBvZiByYW5kb20tZGlnaXQtZGlhbGVkIHRlbGVwaG9uZSBpbnRlcnZpZXdzIChSREQpIGFuZCB0aGUgc2FtZSB1bnByb21wdGVkLCBvcGVuLWVuZGVkIGtleSBxdWVzdGlvbiAiV2hhdCBpcyB5b3VyIHJlbGlnaW9uLCBpZiBhbnk/IiBJbnRlcnZpZXdlcnMgZGlkIG5vdCBwcm9tcHQgb3Igb2ZmZXIgYSBzdWdnZXN0ZWQgbGlzdCBvZiBwb3RlbnRpYWwgYW5zd2Vycy4gTW9yZW92ZXIsIHRoZSBzZWxmLWRlc2NyaXB0aW9uIG9mIHJlc3BvbmRlbnRzIHdhcyBub3QgYmFzZWQgb24gd2hldGhlciBlc3RhYmxpc2hlZCByZWxpZ2lvdXMgYm9kaWVzIG9yIGluc3RpdHV0aW9ucyBjb25zaWRlcmVkIHRoZW0gdG8gYmUgbWVtYmVycy4gVG8gdGhlIGNvbnRyYXJ5LCB0aGUgc3VydmV5cyBzb3VnaHQgdG8gZGV0ZXJtaW5lIHdoZXRoZXIgdGhlIHJlc3BvbmRlbnRzIHJlZ2FyZGVkIHRoZW1zZWx2ZXMgYXMgYWRoZXJlbnRzIG9mIGEgcmVsaWdpb3VzIGNvbW11bml0eS4gVGhlIHN1cnZleXMgdGFwIHN1YmplY3RpdmUsIHJhdGhlciB0aGFuIG9iamVjdGl2ZSwgc3RhbmRhcmRzIG9mIHJlbGlnaW91cyBpZGVudGlmaWNhdGlvbi4gVGhlIHZhbHVlIG9mIHRoaXMgdW5pcXVlIHNlcmllcyBvZiBuYXRpb25hbCBzdXJ2ZXlzLCB3aGljaCBhbGxvd3Mgc2NpZW50aWZpYyBtb25pdG9yaW5nIG9mIGNoYW5nZSBvdmVyIHRpbWUsIGhhcyBiZWVuIHJlY29nbml6ZWQgYnkgdGhlIFUuUy4gQnVyZWF1IG9mIHRoZSBDZW5zdXMgVGhlIEJ1cmVhdSBpdHNlbGYgaXMgY29uc3RpdHV0aW9uYWxseSBwcmVjbHVkZWQgZnJvbSBzdWNoIGFuIGlucXVpcnkgaW50byByZWxpZ2lvbiwgYW5kIHNvIGhhcyBpbmNvcnBvcmF0ZWQgTlNSSS9BUklTIGZpbmRpbmdzIGludG8gaXRzIG9mZmljaWFsIHB1YmxpY2F0aW9uIHRoZSBTdGF0aXN0aWNhbCBBYnN0cmFjdCBvZiB0aGUgVW5pdGVkIFN0YXRlcyBzaW5jZSAyMDAzLg0KDQo+IENhc2VzOiA1NCw0NjEgIA0KPiBWYXJpYWJsZXM6IDMzICANCj4gV2VpZ2h0IFZhcmlhYmxlOiBQT1BXR0hUICANCj4gRGF0ZSBDb2xsZWN0ZWQ6IEZlYnJ1YXJ5IC0gTm92ZW1iZXIgMjAwOA0KDQpjZi4gW1JlbGlnaW91cyBQb3B1bGF0aW9uIGJ5IFN0YXRlXShodHRwczovL2RhdGEud29ybGQvZGF0YWJlYXRzL3JlbGlnaW91cy1wb3B1bGF0aW9uLWJ5LXN0YXRlKQ0KDQpgYGB7cn0NCkFSSVMyMDA4IDwtIGZyZWFkKCJBUklTMjAwOC5jc3YiKQ0KYGBgDQoNCmBgYHtyfQ0Kc3RyKEFSSVMyMDA4KQ0KYGBgDQoNCiMjIFJlbGlnaW9ucyBieSBGcmVxdWVuY2llcw0KDQpgYGB7cn0NCm1ham9yIDwtIG5hbWVzKGhlYWQoc29ydChwcm9wLnRhYmxlKHRhYmxlKEFSSVMyMDA4JFJFTElHKSksIGRlY3JlYXNpbmc9VFJVRSksIDE3KSkNCg0KbWFqb3INCmBgYA0KDQpgYGB7cn0NCkFSSVMyMDA4X21ham9yIDwtIEFSSVMyMDA4W1JFTElHICVpbiUgbWFqb3IsXQ0KDQpucm93KEFSSVMyMDA4X21ham9yKQ0KYGBgDQoNCmBgYHtyfQ0KZ2dwbG90KEFSSVMyMDA4X21ham9yLCBhZXMoUkVMSUcpKSArDQogIGdlb21fYmFyKG1hcHBpbmcgPSBhZXMoeCA9IFJFTElHLCB5ID0gLi5wcm9wLi4sIGdyb3VwID0gMSksIHN0YXQgPSAiY291bnQiKSArDQogIHNjYWxlX3lfY29udGludW91cyhsYWJlbHMgPSBzY2FsZXM6OnBlcmNlbnQpICsNCiAgY29vcmRfZmxpcCgpK3RoZW1lKGF4aXMudGV4dCA9IGVsZW1lbnRfdGV4dChzaXplID0gNykpDQpgYGANCg0KIyMgQ29udGluZ2VuY3kgVGFibGVzDQoNCmBgYHtyfQ0KcHJvcC50YWJsZSh0YWJsZShBUklTMjAwOF9tYWpvciRSRUxJRywgQVJJUzIwMDhfbWFqb3IkUE9MQUZGSUwpLCBtYXJnaW49MSkNCmBgYA0KDQpgYGB7cn0NCnByb3AudGFibGUodGFibGUoQVJJUzIwMDhfbWFqb3IkUkVMSUcsIEFSSVMyMDA4X21ham9yJFJFR1ZPVEUpLCBtYXJnaW49MSkNCmBgYA0KDQojIyBTdGF0ZSBGSVBTIENvZGVzDQoNCmBgYHtyfQ0Kc29ydCh1bmlxdWUoQVJJUzIwMDgkU1RBVEUpKQ0KYGBgDQoNCltTdGF0ZSBGSVBTIENvZGVzXShodHRwczovL3d3dy5ucmNzLnVzZGEuZ292L3dwcy9wb3J0YWwvbnJjcy9kZXRhaWwvP2NpZD1ucmNzMTQzXzAxMzY5NikNCg0KYGBge3J9DQpsb3dlcl80OCA8LSBjKCIwMSIsICIwNCIsICIwNSIsICIwNiIsICIwOCIsICIwOSIsICIxMCIsICIxMSIsICIxMiIsICIxMyIsICIxNiIsICIxNyIsICIxOCIsICIxOSIsICIyMCIsICIyMSIsICIyMiIsICIyMyIsICIyNCIsICIyNSIsICIyNiIsICIyNyIsICIyOCIsICIyOSIsICIzMCIsICIzMSIsICIzMiIsICIzMyIsICIzNCIsICIzNSIsICIzNiIsICIzNyIsICIzOCIsICIzOSIsICI0MCIsICI0MSIsICI0MiIsICI0NCIsICI0NSIsICI0NiIsICI0NyIsICI0OCIsICI0OSIsICI1MCIsICI1MSIsICI1MyIsICI1NCIsICI1NSIsICI1NiIpDQpgYGANCg0KIyMgU3RhdGVzIFNoYXBlZmlsZQ0KDQpVbml0ZWQgU3RhdGVzIENlbnN1cyBCdXJlYXUsIFtDYXJ0b2dyYXBoaWMgQm91bmRhcnkgRmlsZXMgLSBTaGFwZWZpbGVdKGh0dHBzOi8vd3d3LmNlbnN1cy5nb3YvZ2VvZ3JhcGhpZXMvbWFwcGluZy1maWxlcy90aW1lLXNlcmllcy9nZW8vY2FydG8tYm91bmRhcnktZmlsZS5odG1sKQ0KDQpjYl8yMDE4X3VzX3N0YXRlXzUwMGsuemlwDQoNCmBgYHtyfQ0Kc3RhdGVzLnNoIDwtIHJlYWRPR1IoIkM6L1VzZXJzLzE0MTY1L0Rlc2t0b3AvU2hhcGVmaWxlcy9jYl8yMDE4X3VzX3N0YXRlXzUwMGsiLCAiY2JfMjAxOF91c19zdGF0ZV81MDBrIiwgc3RyaW5nc0FzRmFjdG9ycz1GQUxTRSkNCmBgYA0KDQpgYGB7cn0NCnN0cihzdGF0ZXMuc2hAZGF0YSkNCmBgYA0KDQpgYGB7cn0NCnN0YXRlc19sb3dlcjQ4LnNoIDwtIHN0YXRlcy5zaFtzdGF0ZXMuc2hAZGF0YSRHRU9JRCAlaW4lIGxvd2VyXzQ4LF0NCmBgYA0KDQpgYGB7cn0NCnN0cihzdGF0ZXNfbG93ZXI0OC5zaEBkYXRhKQ0KYGBgDQoNCkFkZCAiaWQiIGNvbHVtbg0KYGBge3J9DQpzdGF0ZXNfbG93ZXI0OC5zaEBkYXRhJGlkIDwtIGFzLmludGVnZXIoc3RhdGVzX2xvd2VyNDguc2hAZGF0YSRHRU9JRCkNCmBgYA0KDQpgYGB7cn0NCnNvcnQoc3RhdGVzX2xvd2VyNDguc2hAZGF0YSRpZCkNCmBgYA0KDQpgYGB7cn0NCnN0YXRlc19sb3dlcjQ4LnNoQGRhdGEkU1RBVEUgPC0gYXMuY2hhcmFjdGVyKHN0YXRlc19sb3dlcjQ4LnNoQGRhdGEkU1RVU1BTKQ0KYGBgDQoNClNoYXBlZmlsZSBwcm9jZXNzaW5nDQpgYGB7cn0NCnN0YXRlc19sb3dlcjQ4LnBvaW50cyA9IGZvcnRpZnkoc3RhdGVzX2xvd2VyNDguc2gsIHJlZ2lvbiA9ICJpZCIpDQoNCnN0YXRlc19sb3dlcjQ4LmRmID0gam9pbihzdGF0ZXNfbG93ZXI0OC5wb2ludHMsIHN0YXRlc19sb3dlcjQ4LnNoQGRhdGEsIGJ5ID0gImlkIikNCmBgYA0KDQojIyBTdGF0ZSBSZWxpZ2lvdXMgUmF0aW9zIGFuZCBDaG9yb3BsZXRocw0KDQpgYGB7cn0NClRvdGFsczIwMDggPC0gQVJJUzIwMDhbLCAuKC5OKSwgYnkgPSAuKFNUQVRFKV0NCg0KVG90YWxzMjAwOCA8LSBUb3RhbHMyMDA4W29yZGVyKFNUQVRFKV0NCmBgYA0KDQojIyMgQ2F0aG9saWNzDQoNCmBgYHtyfQ0KQ2F0aG9saWNzMjAwOCA8LSBBUklTMjAwOFtSRUxJRyA9PSAiQ2F0aG9saWMsIFJvbWFuIENhdGhvbGljIiwgLiguTiksIGJ5ID0gLihTVEFURSldDQoNCkNhdGhvbGljczIwMDggPC0gQ2F0aG9saWNzMjAwOFtvcmRlcihTVEFURSldDQoNCkNhdGhvbGljczIwMDgkcmF0aW8gPC0gQ2F0aG9saWNzMjAwOCROL1RvdGFsczIwMDgkTg0KDQpDYXRob2xpY3MyMDA4X2xvd2VyNDggPC0gbWVyZ2Uoc3RhdGVzX2xvd2VyNDguZGYsIENhdGhvbGljczIwMDgsIGJ5ID0gIlNUQVRFIikNCg0KQ2F0aG9saWNzMjAwOF9sb3dlcjQ4IDwtIENhdGhvbGljczIwMDhfbG93ZXI0OFtvcmRlcihDYXRob2xpY3MyMDA4X2xvd2VyNDgkU1RBVEUpLCBdDQpgYGANCg0KTWFrZSBhbmQgcGxvdCBjaG9yb3BsZXRoDQoNCmBgYHtyfQ0KZ2dwbG90KCkgKw0KICBnZW9tX3BvbHlnb24oZGF0YSA9IENhdGhvbGljczIwMDhfbG93ZXI0OCwgDQogICAgICAgICAgICAgICBhZXMoeCA9IGxvbmcsIHkgPSBsYXQsIGdyb3VwID0gZ3JvdXAsIGZpbGwgPSByYXRpbyksIA0KICAgICAgICAgICAgICAgY29sb3IgPSAid2hpdGUiLCBzaXplID0gMC4xKSArIA0KICBjb29yZF9tYXAoKSArIA0KICBzY2FsZV9maWxsX3ZpcmlkaXMob3B0aW9uPSJtYWdtYSIsIGxhYmVscz1wZXJjZW50LCBuYW1lPSJDYXRob2xpY3MiKSArDQogICBsYWJzKHRpdGxlPSJSYXRpb3Mgb2YgQ2F0aG9saWNzIGJ5IFN0YXRlLCBBUklTMjAwOCIpDQpgYGANCg0KIyMjIEJhcHRpc3RzDQoNCmBgYHtyfQ0KQmFwdGlzdHMyMDA4IDwtIEFSSVMyMDA4W1JFTElHID09ICJCYXB0aXN0IiwgLiguTiksIGJ5ID0gLihTVEFURSldDQoNCkJhcHRpc3RzMjAwOCA8LSBCYXB0aXN0czIwMDhbb3JkZXIoU1RBVEUpXQ0KDQpCYXB0aXN0czIwMDgkcmF0aW8gPC0gQmFwdGlzdHMyMDA4JE4vVG90YWxzMjAwOCRODQoNCkJhcHRpc3RzMjAwOF9sb3dlcjQ4IDwtIG1lcmdlKHN0YXRlc19sb3dlcjQ4LmRmLCBCYXB0aXN0czIwMDgsIGJ5ID0gIlNUQVRFIikNCg0KQmFwdGlzdHMyMDA4X2xvd2VyNDggPC0gQmFwdGlzdHMyMDA4X2xvd2VyNDhbb3JkZXIoQmFwdGlzdHMyMDA4X2xvd2VyNDgkU1RBVEUpLCBdDQpgYGANCg0KTWFrZSBhbmQgcGxvdCBjaG9yb3BsZXRoDQoNCmBgYHtyfQ0KZ2dwbG90KCkgKw0KICBnZW9tX3BvbHlnb24oZGF0YSA9IEJhcHRpc3RzMjAwOF9sb3dlcjQ4LCANCiAgICAgICAgICAgICAgIGFlcyh4ID0gbG9uZywgeSA9IGxhdCwgZ3JvdXAgPSBncm91cCwgZmlsbCA9IHJhdGlvKSwgDQogICAgICAgICAgICAgICBjb2xvciA9ICJ3aGl0ZSIsIHNpemUgPSAwLjEpICsgDQogIGNvb3JkX21hcCgpICsgDQogIHNjYWxlX2ZpbGxfdmlyaWRpcyhvcHRpb249Im1hZ21hIiwgbGFiZWxzPXBlcmNlbnQsIG5hbWU9IkJhcHRpc3RzIikgKw0KICAgbGFicyh0aXRsZT0iUmF0aW9zIG9mIEJhcHRpc3RzIGJ5IFN0YXRlLCBBUklTIDIwMDgiKQ0KYGBgDQoNCiMjIyBObyBSZWxpZ2lvbg0KDQpgYGB7cn0NCk5vcmVsaWdpb24yMDA4IDwtIEFSSVMyMDA4W1JFTElHID09ICJObyByZWxpZ2lvbi9ub25lIiwgLiguTiksIGJ5ID0gLihTVEFURSldDQoNCk5vcmVsaWdpb24yMDA4IDwtIE5vcmVsaWdpb24yMDA4W29yZGVyKFNUQVRFKV0NCg0KTm9yZWxpZ2lvbjIwMDgkcmF0aW8gPC0gTm9yZWxpZ2lvbjIwMDgkTi9Ub3RhbHMyMDA4JE4NCg0KTm9yZWxpZ2lvbjIwMDhfbG93ZXI0OCA8LSBtZXJnZShzdGF0ZXNfbG93ZXI0OC5kZiwgTm9yZWxpZ2lvbjIwMDgsIGJ5ID0gIlNUQVRFIikNCg0KTm9yZWxpZ2lvbjIwMDhfbG93ZXI0OCA8LSBOb3JlbGlnaW9uMjAwOF9sb3dlcjQ4W29yZGVyKE5vcmVsaWdpb24yMDA4X2xvd2VyNDgkU1RBVEUpLCBdDQpgYGANCg0KTWFrZSBhbmQgcGxvdCBjaG9yb3BsZXRoDQoNCmBgYHtyfQ0KZ2dwbG90KCkgKw0KICBnZW9tX3BvbHlnb24oZGF0YSA9IE5vcmVsaWdpb24yMDA4X2xvd2VyNDgsIA0KICAgICAgICAgICAgICAgYWVzKHggPSBsb25nLCB5ID0gbGF0LCBncm91cCA9IGdyb3VwLCBmaWxsID0gcmF0aW8pLCANCiAgICAgICAgICAgICAgIGNvbG9yID0gIndoaXRlIiwgc2l6ZSA9IDAuMSkgKyANCiAgY29vcmRfbWFwKCkgKyANCiAgc2NhbGVfZmlsbF92aXJpZGlzKG9wdGlvbj0ibWFnbWEiLCBsYWJlbHM9cGVyY2VudCwgbmFtZT0iTm8gUmVsaWdpb24iKSArDQogICBsYWJzKHRpdGxlPSJSYXRpb3Mgb2YgTm8gUmVsaWdpb24gYnkgU3RhdGUsIEFSSVMgMjAwOCIpDQpgYGANCg0KDQojIEFSSVMgMjAwMQ0KDQpbRG93bmxvYWQgQVJJUyAyMDAxIERhdGFdKGh0dHBzOi8vY29tbW9ucy50cmluY29sbC5lZHUvYXJpcy9zdXJ2ZXlzL2FyaXMtMjAwMS9kb3dubG9hZC1hcmlzLTIwMDEvKQ0KDQpBUklTMjAwMVNQU1MuemlwDQoNCmBgYHtyfQ0KQVJJUzIwMDEgPC0gZnJlYWQoIkFSSVMyMDAxLmNzdiIpDQpgYGANCg0KW0FSSVMgMjAwMSBDb2RlYm9va10oaHR0cHM6Ly9jb21tb25zLnRyaW5jb2xsLmVkdS9hcmlzL3N1cnZleXMvYXJpcy0yMDAxL2FyaXMtMjAwMS1jb2RlYm9vay8pDQoNCmBgYHtyfQ0KQVJJUzIwMDEgPC0gQVJJUzIwMDFbLCAuKFBPUFdHSFQgPSB3Z3RfcG9wLCBTVEFURSA9IHN0YXRlLCBDT1VOVFlGSVBTID0gZmlwcywgQURJID0gYWRpLCBETUEgPSBkbWEsIE1TQSA9IG1zYSwgTUVUUk8gPSBtZXRybywgTklFTFNFTiA9IG5laWxzZW4sIERJVklTSU9OID0gY2VuX2RpdiwgUkVHSU9OID0gcmVnaW9uLCBTRVggPSBnZW5kZXIsIE9XTkhPTUUgPSBkZW0yLCBNQVJJVEFMID0gZGVtMywgRU1QTE9ZID0gZGVtNCwgSU5DT01FID0gbmV3ZGVtOSwgQUdFID0gZGVtNywgRURVQyA9IGRlbTgsIEhJU1BBTklDID0gZGVtMTAsIFJBQ0UgPSBkZW0xMSwgUkVMSUcgPSBqcTF0b3QsIFJFTElHX1NQT1VTRSA9IGpxMnRvdCwgQUdFQ0FUID0gYWdlY2F0LCBQT0xBRkZJTCA9IHoxMWEsIFJFR1ZPVEUgPSB6MTFiKV0NCmBgYA0KDQpgYGB7cn0NCkFSSVMyMDAxJFNUQVRFIDwtIHN0YXRlMmFiYnIoQVJJUzIwMDEkU1RBVEUpDQpgYGANCg0KYGBge3J9DQptYWpvcl8yMDAxIDwtIG5hbWVzKGhlYWQoc29ydChwcm9wLnRhYmxlKHRhYmxlKEFSSVMyMDAxJFJFTElHKSksIGRlY3JlYXNpbmc9VFJVRSksIDE3KSkNCg0KbWFqb3JfMjAwMQ0KYGBgDQoNCmBgYHtyfQ0KQVJJUzIwMDFfbWFqb3IgPC0gQVJJUzIwMDFbUkVMSUcgJWluJSBtYWpvcl8yMDAxLF0NCg0KbnJvdyhBUklTMjAwMV9tYWpvcikNCmBgYA0KDQpgYGB7cn0NCmdncGxvdChBUklTMjAwMV9tYWpvciwgYWVzKFJFTElHKSkgKw0KICBnZW9tX2JhcihtYXBwaW5nID0gYWVzKHggPSBSRUxJRywgeSA9IC4ucHJvcC4uLCBncm91cCA9IDEpLCBzdGF0ID0gImNvdW50IikgKw0KICBzY2FsZV95X2NvbnRpbnVvdXMobGFiZWxzID0gc2NhbGVzOjpwZXJjZW50KSArDQogIGNvb3JkX2ZsaXAoKSt0aGVtZShheGlzLnRleHQgPSBlbGVtZW50X3RleHQoc2l6ZSA9IDcpKQ0KYGBgDQoNCmBgYHtyfQ0KVG90YWxzMjAwMSA8LSBBUklTMjAwMVssIC4oLk4pLCBieSA9IC4oU1RBVEUpXQ0KDQpUb3RhbHMyMDAxIDwtIFRvdGFsczIwMDFbb3JkZXIoU1RBVEUpXQ0KYGBgDQoNCg0KIyMgQ2F0aG9saWNzDQoNCmBgYHtyfQ0KQ2F0aG9saWNzMjAwMSA8LSBBUklTMjAwMVtSRUxJRyA9PSAiQ2F0aG9saWMiLCAuKC5OKSwgYnkgPSAuKFNUQVRFKV0NCg0KQ2F0aG9saWNzMjAwMSA8LSBDYXRob2xpY3MyMDAxW29yZGVyKFNUQVRFKV0NCg0KQ2F0aG9saWNzMjAwMSRyYXRpbyA8LSBDYXRob2xpY3MyMDAxJE4vVG90YWxzMjAwMSRODQoNCkNhdGhvbGljczIwMDFfbG93ZXI0OCA8LSBtZXJnZShzdGF0ZXNfbG93ZXI0OC5kZiwgQ2F0aG9saWNzMjAwMSwgYnkgPSAiU1RBVEUiKQ0KDQpDYXRob2xpY3MyMDAxX2xvd2VyNDggPC0gQ2F0aG9saWNzMjAwMV9sb3dlcjQ4W29yZGVyKENhdGhvbGljczIwMDFfbG93ZXI0OCRTVEFURSksIF0NCmBgYA0KDQpNYWtlIGFuZCBwbG90IGNob3JvcGxldGgNCg0KYGBge3J9DQpnZ3Bsb3QoKSArDQogIGdlb21fcG9seWdvbihkYXRhID0gQ2F0aG9saWNzMjAwMV9sb3dlcjQ4LCANCiAgICAgICAgICAgICAgIGFlcyh4ID0gbG9uZywgeSA9IGxhdCwgZ3JvdXAgPSBncm91cCwgZmlsbCA9IHJhdGlvKSwgDQogICAgICAgICAgICAgICBjb2xvciA9ICJ3aGl0ZSIsIHNpemUgPSAwLjEpICsgDQogIGNvb3JkX21hcCgpICsgDQogIHNjYWxlX2ZpbGxfdmlyaWRpcyhvcHRpb249Im1hZ21hIiwgbGFiZWxzPXBlcmNlbnQsIG5hbWU9IkNhdGhvbGljcyIpICsNCiAgIGxhYnModGl0bGU9IlJhdGlvcyBvZiBDYXRob2xpY3MgYnkgU3RhdGUsIEFSSVMgMjAwMSIpDQpgYGANCg0KIyMgQmFwdGlzdHMNCg0KYGBge3J9DQpCYXB0aXN0czIwMDEgPC0gQVJJUzIwMDFbUkVMSUcgPT0gIkJhcHRpc3QiLCAuKC5OKSwgYnkgPSAuKFNUQVRFKV0NCg0KQmFwdGlzdHMyMDAxIDwtIEJhcHRpc3RzMjAwMVtvcmRlcihTVEFURSldDQoNCkJhcHRpc3RzMjAwMSRyYXRpbyA8LSBCYXB0aXN0czIwMDEkTi9Ub3RhbHMyMDAxJE4NCg0KQmFwdGlzdHMyMDAxX2xvd2VyNDggPC0gbWVyZ2Uoc3RhdGVzX2xvd2VyNDguZGYsIEJhcHRpc3RzMjAwMSwgYnkgPSAiU1RBVEUiKQ0KDQpCYXB0aXN0czIwMDFfbG93ZXI0OCA8LSBCYXB0aXN0czIwMDFfbG93ZXI0OFtvcmRlcihCYXB0aXN0czIwMDFfbG93ZXI0OCRTVEFURSksIF0NCg0KZ2dwbG90KCkgKw0KICBnZW9tX3BvbHlnb24oZGF0YSA9IEJhcHRpc3RzMjAwMV9sb3dlcjQ4LCANCiAgICAgICAgICAgICAgIGFlcyh4ID0gbG9uZywgeSA9IGxhdCwgZ3JvdXAgPSBncm91cCwgZmlsbCA9IHJhdGlvKSwgDQogICAgICAgICAgICAgICBjb2xvciA9ICJ3aGl0ZSIsIHNpemUgPSAwLjEpICsgDQogIGNvb3JkX21hcCgpICsgDQogIHNjYWxlX2ZpbGxfdmlyaWRpcyhvcHRpb249Im1hZ21hIiwgbGFiZWxzPXBlcmNlbnQsIG5hbWU9IkJhcHRpc3RzIikgKw0KICAgbGFicyh0aXRsZT0iUmF0aW9zIG9mIEJhcHRpc3RzIGJ5IFN0YXRlLCBBUklTIDIwMDEiKQ0KYGBgDQoNCiMjIE5vIFJlbGlnaW9uDQoNCmBgYHtyfQ0KTm9yZWxpZ2lvbjIwMDEgPC0gQVJJUzIwMDFbUkVMSUcgPT0gIk5vbmUiLCAuKC5OKSwgYnkgPSAuKFNUQVRFKV0NCg0KTm9yZWxpZ2lvbjIwMDEgPC0gTm9yZWxpZ2lvbjIwMDFbb3JkZXIoU1RBVEUpXQ0KDQpOb3JlbGlnaW9uMjAwMSRyYXRpbyA8LSBOb3JlbGlnaW9uMjAwMSROL1RvdGFsczIwMDEkTg0KDQpOb3JlbGlnaW9uMjAwMV9sb3dlcjQ4IDwtIG1lcmdlKHN0YXRlc19sb3dlcjQ4LmRmLCBOb3JlbGlnaW9uMjAwMSwgYnkgPSAiU1RBVEUiKQ0KDQpOb3JlbGlnaW9uMjAwMV9sb3dlcjQ4IDwtIE5vcmVsaWdpb24yMDAxX2xvd2VyNDhbb3JkZXIoTm9yZWxpZ2lvbjIwMDFfbG93ZXI0OCRTVEFURSksIF0NCg0KZ2dwbG90KCkgKw0KICBnZW9tX3BvbHlnb24oZGF0YSA9IE5vcmVsaWdpb24yMDAxX2xvd2VyNDgsIA0KICAgICAgICAgICAgICAgYWVzKHggPSBsb25nLCB5ID0gbGF0LCBncm91cCA9IGdyb3VwLCBmaWxsID0gcmF0aW8pLCANCiAgICAgICAgICAgICAgIGNvbG9yID0gIndoaXRlIiwgc2l6ZSA9IDAuMSkgKyANCiAgY29vcmRfbWFwKCkgKyANCiAgc2NhbGVfZmlsbF92aXJpZGlzKG9wdGlvbj0ibWFnbWEiLCBsYWJlbHM9cGVyY2VudCwgbmFtZT0iTm8gUmVsaWdpb24iKSArDQogICBsYWJzKHRpdGxlPSJSYXRpb3Mgb2YgTm8gUmVsaWdpb24gYnkgU3RhdGUsIEFSSVMgMjAwMSIpDQpgYGANCg0KIyBDaGFuZ2UgZnJvbSAyMDAxIHRvIDIwMDgNCg0KIyMgQ2F0aG9saWNzDQoNCmBgYHtyfQ0KQ2F0aG9saWNzMjAwOCRjaGFuZ2UgPC0gKENhdGhvbGljczIwMDgkcmF0aW8tQ2F0aG9saWNzMjAwMSRyYXRpbykvQ2F0aG9saWNzMjAwMSRyYXRpbw0KYGBgDQoNCmBgYHtyfQ0KQ2F0aG9saWNzMjAwOGNoYW5nZV9sb3dlcjQ4IDwtIG1lcmdlKHN0YXRlc19sb3dlcjQ4LmRmLCBDYXRob2xpY3MyMDA4LCBieSA9ICJTVEFURSIpDQoNCkNhdGhvbGljczIwMDhjaGFuZ2VfbG93ZXI0OCA8LSBDYXRob2xpY3MyMDA4Y2hhbmdlX2xvd2VyNDhbb3JkZXIoQ2F0aG9saWNzMjAwOGNoYW5nZV9sb3dlcjQ4JFNUQVRFKSwgXQ0KYGBgDQoNCk1ha2UgYW5kIHBsb3QgY2hvcm9wbGV0aA0KDQpgYGB7cn0NCmdncGxvdCgpICsNCiAgZ2VvbV9wb2x5Z29uKGRhdGEgPSBDYXRob2xpY3MyMDA4Y2hhbmdlX2xvd2VyNDgsIA0KICAgICAgICAgICAgICAgYWVzKHggPSBsb25nLCB5ID0gbGF0LCBncm91cCA9IGdyb3VwLCBmaWxsID0gY2hhbmdlKSwgDQogICAgICAgICAgICAgICBjb2xvciA9ICJ3aGl0ZSIsIHNpemUgPSAwLjEpICsgDQogIGNvb3JkX21hcCgpICsgDQogIHNjYWxlX2ZpbGxfdmlyaWRpcyhvcHRpb249Im1hZ21hIiwgbGFiZWxzPXBlcmNlbnQsIG5hbWU9IkNoYW5nZSIpICsNCiAgIGxhYnModGl0bGU9IihBUklTMjAwOC1BUklTMjAwMSkvQVJJUzIwMDEsIENhdGhvbGljcyIpDQpgYGANCg0KDQojIyBCYXB0aXN0cw0KDQpgYGB7cn0NCkJhcHRpc3RzMjAwOCRjaGFuZ2UgPC0gKEJhcHRpc3RzMjAwOCRyYXRpby1CYXB0aXN0czIwMDEkcmF0aW8pL0JhcHRpc3RzMjAwMSRyYXRpbw0KDQpCYXB0aXN0czIwMDhjaGFuZ2VfbG93ZXI0OCA8LSBtZXJnZShzdGF0ZXNfbG93ZXI0OC5kZiwgQmFwdGlzdHMyMDA4LCBieSA9ICJTVEFURSIpDQoNCkJhcHRpc3RzMjAwOGNoYW5nZV9sb3dlcjQ4IDwtIEJhcHRpc3RzMjAwOGNoYW5nZV9sb3dlcjQ4W29yZGVyKEJhcHRpc3RzMjAwOGNoYW5nZV9sb3dlcjQ4JFNUQVRFKSwgXQ0KDQpnZ3Bsb3QoKSArDQogIGdlb21fcG9seWdvbihkYXRhID0gQmFwdGlzdHMyMDA4Y2hhbmdlX2xvd2VyNDgsIA0KICAgICAgICAgICAgICAgYWVzKHggPSBsb25nLCB5ID0gbGF0LCBncm91cCA9IGdyb3VwLCBmaWxsID0gY2hhbmdlKSwgDQogICAgICAgICAgICAgICBjb2xvciA9ICJ3aGl0ZSIsIHNpemUgPSAwLjEpICsgDQogIGNvb3JkX21hcCgpICsgDQogIHNjYWxlX2ZpbGxfdmlyaWRpcyhvcHRpb249Im1hZ21hIiwgbGFiZWxzPXBlcmNlbnQsIG5hbWU9IkNoYW5nZSIpICsNCiAgIGxhYnModGl0bGU9IihBUklTMjAwOC1BUklTMjAwMSkvQVJJUzIwMDEsIEJhcHRpc3RzIikNCmBgYA0KDQojIyBObyBSZWxpZ2lvbg0KDQpgYGB7cn0NCk5vcmVsaWdpb24yMDA4JGNoYW5nZSA8LSAoTm9yZWxpZ2lvbjIwMDgkcmF0aW8tTm9yZWxpZ2lvbjIwMDEkcmF0aW8pL05vcmVsaWdpb24yMDAxJHJhdGlvDQoNCk5vcmVsaWdpb24yMDA4Y2hhbmdlX2xvd2VyNDggPC0gbWVyZ2Uoc3RhdGVzX2xvd2VyNDguZGYsIE5vcmVsaWdpb24yMDA4LCBieSA9ICJTVEFURSIpDQoNCk5vcmVsaWdpb24yMDA4Y2hhbmdlX2xvd2VyNDggPC0gTm9yZWxpZ2lvbjIwMDhjaGFuZ2VfbG93ZXI0OFtvcmRlcihOb3JlbGlnaW9uMjAwOGNoYW5nZV9sb3dlcjQ4JFNUQVRFKSwgXQ0KDQpnZ3Bsb3QoKSArDQogIGdlb21fcG9seWdvbihkYXRhID0gTm9yZWxpZ2lvbjIwMDhjaGFuZ2VfbG93ZXI0OCwgDQogICAgICAgICAgICAgICBhZXMoeCA9IGxvbmcsIHkgPSBsYXQsIGdyb3VwID0gZ3JvdXAsIGZpbGwgPSBjaGFuZ2UpLCANCiAgICAgICAgICAgICAgIGNvbG9yID0gIndoaXRlIiwgc2l6ZSA9IDAuMSkgKyANCiAgY29vcmRfbWFwKCkgKyANCiAgc2NhbGVfZmlsbF92aXJpZGlzKG9wdGlvbj0ibWFnbWEiLCBsYWJlbHM9cGVyY2VudCwgbmFtZT0iQ2hhbmdlIikgKw0KICAgbGFicyh0aXRsZT0iKEFSSVMyMDA4LUFSSVMyMDAxKS9BUklTMjAwMSwgTm8gUmVsaWdpb24iKQ0KYGBgDQoNCiMgQVJJUyAyMDA4IERlbnNpdHkgTWFwcw0KDQojIyBTdGF0ZSBhbmQgQ291bnR5IEZJUFMgQ29kZXMNCg0KW1N0YXRlIEZJUFMgQ29kZXNdKGh0dHBzOi8vd3d3Lm5yY3MudXNkYS5nb3Yvd3BzL3BvcnRhbC9ucmNzL2RldGFpbC8/Y2lkPW5yY3MxNDNfMDEzNjk2KQ0KDQpbQ291bnR5IEZJUFMgQ29kZXNdKGh0dHBzOi8vd3d3Lm5yY3MudXNkYS5nb3Yvd3BzL3BvcnRhbC9ucmNzL2RldGFpbC9uYXRpb25hbC9ob21lLz9jaWQ9bnJjczE0M18wMTM2OTcpDQoNCmBgYHtyfQ0KQVJJUzIwMDhfY291bnRpZXMgPC0gQVJJUzIwMDhbLC4oQ09VTlRZRklQUyA9IENPVU5UWUMsIFNUQVRFLCBSRUxJRyldDQoNCkFSSVMyMDA4X2NvdW50aWVzJFNUQVRFRklQUyA8LSBmaXBzKEFSSVMyMDA4X2NvdW50aWVzJFNUQVRFLCB0byA9ICJGSVBTIikNCmBgYA0KDQojIyBDb3VudGllcyBTaGFwZWZpbGUNCg0KVW5pdGVkIFN0YXRlcyBDZW5zdXMgQnVyZWF1LCBbQ2FydG9ncmFwaGljIEJvdW5kYXJ5IEZpbGVzIC0gU2hhcGVmaWxlXShodHRwczovL3d3dy5jZW5zdXMuZ292L2dlb2dyYXBoaWVzL21hcHBpbmctZmlsZXMvdGltZS1zZXJpZXMvZ2VvL2NhcnRvLWJvdW5kYXJ5LWZpbGUuaHRtbCkNCg0KY2JfMjAxOF91c19jb3VudHlfNTAway56aXANCg0KYGBge3J9DQpjb3VudGllcy5zaCA8LSByZWFkT0dSKCJDOi9Vc2Vycy8xNDE2NS9EZXNrdG9wL1NoYXBlZmlsZXMvY2JfMjAxOF91c19jb3VudHlfNTAwayIsICJjYl8yMDE4X3VzX2NvdW50eV81MDBrIiwgc3RyaW5nc0FzRmFjdG9ycz1GQUxTRSkNCmBgYA0KDQpgYGB7cn0NCnN0cihjb3VudGllcy5zaEBkYXRhKQ0KYGBgDQoNCmBgYHtyfQ0KY291bnRpZXNfbG93ZXI0OC5zaCA8LSBjb3VudGllcy5zaFtjb3VudGllcy5zaEBkYXRhJFNUQVRFRlAgJWluJSBsb3dlcl80OCxdDQpgYGANCg0KYGBge3J9DQpjb3VudGllc19sb3dlcjQ4LnNoQGRhdGEkU1RBVEVGSVBTIDwtIGFzLmludGVnZXIoY291bnRpZXNfbG93ZXI0OC5zaEBkYXRhJFNUQVRFRlApDQoNCmNvdW50aWVzX2xvd2VyNDguc2hAZGF0YSRDT1VOVFlGSVBTIDwtIGFzLmludGVnZXIoY291bnRpZXNfbG93ZXI0OC5zaEBkYXRhJENPVU5UWUZQKQ0KYGBgDQoNCmBgYHtyfQ0KY291bnRpZXNfbG93ZXI0OC5jZW50cm9pZHMgIDwtIGFzLmRhdGEudGFibGUoZ0NlbnRyb2lkKGNvdW50aWVzX2xvd2VyNDguc2gsIGJ5aWQgPSBUUlVFKSkNCg0KY291bnRpZXNfbG93ZXI0OC5jZW50cm9pZHMkU1RBVEVGSVBTIDwtIGNvdW50aWVzX2xvd2VyNDguc2hAZGF0YSRTVEFURUZJUFMNCg0KY291bnRpZXNfbG93ZXI0OC5jZW50cm9pZHMkQ09VTlRZRklQUyA8LSBjb3VudGllc19sb3dlcjQ4LnNoQGRhdGEkQ09VTlRZRklQUw0KYGBgDQoNCmBgYHtyfQ0Kc3RyKGNvdW50aWVzX2xvd2VyNDguY2VudHJvaWRzKQ0KYGBgDQoNCmBgYHtyfQ0Kc3RyKEFSSVMyMDA4X2NvdW50aWVzKQ0KYGBgDQoNCmBgYHtyfQ0KQVJJUzIwMDhfWFkgPC0gbGVmdF9qb2luKEFSSVMyMDA4X2NvdW50aWVzLCBjb3VudGllc19sb3dlcjQ4LmNlbnRyb2lkcywgYnkgPSBjKCJTVEFURUZJUFMiLCAiQ09VTlRZRklQUyIpKQ0KYGBgDQoNCmBgYHtyfQ0KQVJJUzIwMDhfWFkkQ09VTlRZRklQUyA8LSBOVUxMDQoNCkFSSVMyMDA4X1hZJFNUQVRFIDwtIE5VTEwNCg0KQVJJUzIwMDhfWFkkU1RBVEVGSVBTIDwtIE5VTEwNCg0Kc3RyKEFSSVMyMDA4X1hZKQ0KYGBgDQoNCmBgYHtyfQ0KQ2F0aG9saWNzMjAwOF9YWSA8LSBBUklTMjAwOF9YWVtBUklTMjAwOF9YWSRSRUxJRyA9PSAiQ2F0aG9saWMsIFJvbWFuIENhdGhvbGljIixdDQoNCkNhdGhvbGljczIwMDhfWFkkUkVMSUcgPC0gTlVMTA0KYGBgDQoNCg0KYGBge3J9DQpCYXB0aXN0czIwMDhfWFkgPC0gQVJJUzIwMDhfWFlbQVJJUzIwMDhfWFkkUkVMSUcgPT0gIkJhcHRpc3QiLF0NCg0KQmFwdGlzdHMyMDA4X1hZJFJFTElHIDwtIE5VTEwNCmBgYA0KDQpgYGB7cn0NCk5vcmVsaWdpb24yMDA4X1hZIDwtIEFSSVMyMDA4X1hZW0FSSVMyMDA4X1hZJFJFTElHID09ICJObyByZWxpZ2lvbi9ub25lIixdDQoNCk5vcmVsaWdpb24yMDA4X1hZJFJFTElHIDwtIE5VTEwNCmBgYA0KDQojIyBDYXRob2xpYyBEZW5zaXR5IE1hcA0KDQpgYGB7cn0NCmdncGxvdCgpICsgDQogIGdlb21fcG9seWdvbihkYXRhID0gc3RhdGVzX2xvd2VyNDguZGYsIGFlcyh4ID0gbG9uZywgeSA9IGxhdCwgZ3JvdXAgPSBncm91cCksDQogICAgICAgICAgY29sb3IgPSAnYmxhY2snLCBzaXplID0gMSwgZmlsbD1OQSkgKw0KICBzdGF0X2RlbnNpdHkyZChkYXRhPUNhdGhvbGljczIwMDhfWFksIGFlcyh4PXgsIHk9eSwgZmlsbD0uLmxldmVsLi4pLCBhbHBoYT0wLjIsIGdlb20gPSAncG9seWdvbicsIGNvbG91ciA9ICdibGFjaycsIGNvbnRvdXI9VFJVRSkgKw0KICBzY2FsZV9maWxsX2NvbnRpbnVvdXMobG93PSJ5ZWxsb3ciLCBoaWdoPSJyZWQiKSsNCiAgdGhlbWVfbm90aGluZyhsZWdlbmQgPSBUUlVFKSArICAgDQogIGxhYnModGl0bGU9IkNhdGhvbGljIERlbml0eSBNYXAsIEFSSVMgMjAwOCIpDQpgYGANCg0KIyMgQmFwdGlzdCBEZW5zaXR5IE1hcA0KDQpgYGB7cn0NCmdncGxvdCgpICsgDQogIGdlb21fcG9seWdvbihkYXRhID0gc3RhdGVzX2xvd2VyNDguZGYsIGFlcyh4ID0gbG9uZywgeSA9IGxhdCwgZ3JvdXAgPSBncm91cCksDQogICAgICAgICAgY29sb3IgPSAnYmxhY2snLCBzaXplID0gMSwgZmlsbD1OQSkgKw0KICBzdGF0X2RlbnNpdHkyZChkYXRhPUJhcHRpc3RzMjAwOF9YWSwgYWVzKHg9eCwgeT15LCBmaWxsPS4ubGV2ZWwuLiksIGFscGhhPTAuMiwgZ2VvbSA9ICdwb2x5Z29uJywgY29sb3VyID0gJ2JsYWNrJywgY29udG91cj1UUlVFKSArDQogIHNjYWxlX2ZpbGxfY29udGludW91cyhsb3c9InllbGxvdyIsaGlnaD0icmVkIikrDQogIHRoZW1lX25vdGhpbmcobGVnZW5kID0gVFJVRSkgKyAgDQogIGxhYnModGl0bGU9IkJhcHRpc3QgRGVuc2l0eSBNYXAsIEFSSVMgMjAwOCIpDQpgYGANCg0KIyMgTm8gUmVsaWdpb24gRGVuc2l0eSBNYXANCg0KYGBge3J9DQpnZ3Bsb3QoKSArIA0KICBnZW9tX3BvbHlnb24oZGF0YSA9IHN0YXRlc19sb3dlcjQ4LmRmLCBhZXMoeCA9IGxvbmcsIHkgPSBsYXQsIGdyb3VwID0gZ3JvdXApLA0KICAgICAgICAgIGNvbG9yID0gJ2JsYWNrJywgc2l6ZSA9IDEsIGZpbGw9TkEpICsNCiAgc3RhdF9kZW5zaXR5MmQoZGF0YT1Ob3JlbGlnaW9uMjAwOF9YWSwgYWVzKHg9eCwgeT15LCBmaWxsPS4ubGV2ZWwuLiksIGFscGhhPTAuMiwgZ2VvbSA9ICdwb2x5Z29uJywgY29sb3VyID0gJ2JsYWNrJywgY29udG91cj1UUlVFKSArDQogIHNjYWxlX2ZpbGxfY29udGludW91cyhsb3c9InllbGxvdyIsaGlnaD0icmVkIikrDQogIHRoZW1lX25vdGhpbmcobGVnZW5kID0gVFJVRSkgKyAgDQogIGxhYnModGl0bGU9Ik5vIFJlbGlnaW9uIERlbnNpdHkgTWFwLCBBUklTIDIwMDgiKQ0KYGBgDQoNCg==