library(tidyverse)
library(tidycensus)
library(zipcodeR)
KEY <- '30e13ab22563318ff59286e433099f4174d4edd4'
TRACTS <- get_tracts(search_city('Diamondville','WY')$zipcode) %>% full_join(get_tracts(search_city('Kemmerer','WY')$zipcode))
TRACTS 
census_api_key(KEY, install = TRUE)
vars <- c("P2_002N")
get_decennial(geography="tract",variables=vars,state='WY',county='lincoln',GEOID='56023978400')
vars <- c("B01001_004","B01001_005")

get_acs(geography="tract",variables=vars,state='WY',county='lincoln',GEOID='56023978400')
######MALE

#Variable names from: https://api.census.gov/data/2019/acs/acs1/variables.html
VARS <-c(
#Under_Five
'B01001_003E',
#5-9
'B01001_004E',
#10-14
'B01001_005E',
#15-17
'B01001_006E',
#18-19
'B01001_007E',
#20
'B01001_008E',
#21
'B01001_009E',
#22-24
'B01001_010E',
#25-29
'B01001_011E',
#85+
'B01001_025E')
#Testing age Comparison between the two
TEMP <- get_acs(geography="tract",variables=VARS,state='WY',county='lincoln')
KEM <- TEMP %>% filter(GEOID=='56023978400')
KEM
19/sum(KEM$estimate)

LN <- TEMP %>% filter(GEOID!='56023978400')
LN
19/sum(KEM$estimate)