R/HouseVal.R
In ColoradoDemography/codemogLib:
Defines functions
HouseVal
#' HouseVal Table Showing the Median House value for Owner-Occupied housing,
#' the Median Gross Rent, and Median Costs as a Percentage of Income for
#' Owners and Renters for a place and the State of Colorado
#'
#' Restructured 5/4/2018 to remove codemog_api calls
#'
#' @param listID the list containing place id and Place names
#' @param ACS Specifies the ACS data set to be used, reads curACS from Shiny program
#' @param oType output type html table or latex table
#' @return kable formatted table and data file
#' @export
#'
HouseVal <- function(listID, ACS, oType, state="08"){
# Collecting place ids from idList, setting default values
ctyfips <- listID$ctyNum
ctyname <- listID$ctyName
placefips <- listID$plNum
placename <- listID$plName
# if(listID$PlFilter == "T") {
# placefips <- ""
# placename <- ""
# }
#Prepping SQL calls
#County SQL Calls
CTYSQL01 <- paste0("SELECT * FROM data.houseval_cty WHERE (geonum = ",paste0("1", state, ctyfips)," AND acs = '",ACS,"');")
if(nchar(placefips) == 0) {
#State SQL Calls
STSQL01 <- paste0("SELECT * FROM data.houseval_st WHERE acs = '",ACS,"';")
} else {
# Place SQL calls
PLSQL01 <- paste0("SELECT * FROM data.houseval_muni WHERE (geonum = ",paste0("1", state, placefips)," AND acs = '",ACS,"');")
}
# Call to Postrgres
pw <- {
"demography"
}
# loads the PostgreSQL driver
drv <- dbDriver("PostgreSQL")
# creates a connection to the postgres database
# note that "con" will be used later in each connection to the database
con <- dbConnect(drv, dbname = 'dola',
host = "104.197.26.248", port = 5433,
user = "codemog", password = pw)
rm(pw) # removes the password
# Read data files
# Raw County data
f.ACSCTY <- dbGetQuery(con, CTYSQL01)
if(nchar(placefips) == 0) {
# Raw State data
f.ACSST <- dbGetQuery(con, STSQL01)
} else {
# Raw Place data
f.ACSPL <- dbGetQuery(con, PLSQL01)
}
#closing the connections
dbDisconnect(con)
dbUnloadDriver(drv)
rm(con)
rm(drv)
#Accounting for Missing value
f.ACSCTY[f.ACSCTY == "-9999"] <- NA
if(nchar(placefips) == 0) {
f.ACSST[f.ACSST == "-9999"] <- NA
} else {
f.ACSPL[f.ACSPL == "-9999"] <- NA
}
# Assemble County Data
tot_OO <- as.numeric(f.ACSCTY$b25095001)
tot_RT <- as.numeric(f.ACSCTY$b25074001)
f.ACSCTY <- f.ACSCTY %>% mutate(
#Values
Med_val = b25077001,
OO_3049 = b25095006 + b25095007 + b25095008 +
b25095015 + b25095016 + b25095017 +
b25095024 + b25095025 + b25095026 +
b25095033 + b25095034 + b25095035 +
b25095042 + b25095043 + b25095044 +
b25095051 + b25095052 + b25095053 +
b25095060 + b25095061 + b25095062 +
b25095069 + b25095070 + b25095071,
OO_50 = b25095009 + b25095018 + b25095027 + b25095036 + b25095045 + b25095054 + b25095063 + b25095072,
PROP_OO_3049 = OO_3049/tot_OO,
PCT_OO_3049 = percent(PROP_OO_3049*100),
PROP_OO_50 = OO_50/tot_OO,
PCT_OO_50 = percent(PROP_OO_50*100),
Med_Rent = b25064001,
RT_3049 = b25074006 + b25074007 + b25074008 +
b25074015 + b25074016 + b25074017 +
b25074024 + b25074025 + b25074026 +
b25074033 + b25074034 + b25074035 +
b25074042 + b25074043 + b25074044 +
b25074051 + b25074052 + b25074053 +
b25074060 + b25074061 + b25074062,
RT_50 = b25074009 + b25074018 + b25074027 + b25074036 + b25074045 + b25074054 + b25074063,
PROP_RT_3049 = RT_3049/tot_RT,
PCT_RT_3049 = percent(PROP_RT_3049*100),
PROP_RT_50 = RT_50/tot_RT,
PCT_RT_50 = percent(PROP_RT_50*100),
#MOE
Med_val_moe = b25077_moe001,
OO_3049_moe = sqrt(b25095_moe006^2 + b25095_moe007^2 + b25095_moe008^2 +
b25095_moe015^2 + b25095_moe016^2 + b25095_moe017^2 +
b25095_moe024^2 + b25095_moe025^2 + b25095_moe026^2 +
b25095_moe033^2 + b25095_moe034^2 + b25095_moe035^2 +
b25095_moe042^2 + b25095_moe043^2 + b25095_moe044^2 +
b25095_moe051^2 + b25095_moe052^2 + b25095_moe053^2 +
b25095_moe060^2 + b25095_moe061^2 + b25095_moe062^2 +
b25095_moe069^2 + b25095_moe070^2 + b25095_moe071^2),
OO_50_moe = sqrt(b25095_moe009^2 + b25095_moe018^2 + b25095_moe027^2 + b25095_moe036^2 +
b25095_moe045^2 + b25095_moe054^2 + b25095_moe063^2 + b25095_moe072^2),
PROP_OO_3049_moe = OO_3049_moe/tot_OO,
PCT_OO_3049_moe = percent(PROP_OO_3049_moe*100),
PROP_OO_50_moe = OO_50_moe/tot_OO,
PCT_OO_50_moe = percent(PROP_OO_50_moe*100),
Med_Rent_moe = b25064_moe001,
RT_3049_moe = sqrt(b25074_moe006^2 + b25074_moe007^2 + b25074_moe008^2 +
b25074_moe015^2 + b25074_moe016^2 + b25074_moe017^2 +
b25074_moe024^2 + b25074_moe025^2 + b25074_moe026^2 +
b25074_moe033^2 + b25074_moe034^2 + b25074_moe035^2 +
b25074_moe042^2 + b25074_moe043^2 + b25074_moe044^2 +
b25074_moe051^2 + b25074_moe052^2 + b25074_moe053^2 +
b25074_moe060^2 + b25074_moe061^2 + b25074_moe062^2),
RT_50_moe = sqrt(b25074_moe009^2 + b25074_moe018^2 + b25074_moe027^2 +
b25074_moe036^2 + b25074_moe045^2 + b25074_moe054^2 + b25074_moe063^2),
PROP_RT_3049_moe = RT_3049_moe/tot_RT,
PCT_RT_3049_moe = percent(PROP_RT_3049_moe*100),
PROP_RT_50_moe = RT_50_moe/tot_RT,
PCT_RT_50_moe = percent(PROP_RT_50_moe*100)
)
f.ACSCTY_val <- f.ACSCTY[,c(2,129,132:135,136,139:142)]
f.ACSCTY_moe <- f.ACSCTY[,c(2,143,146:149,150,153:156)]
f.ACSCTY_valL <- as.data.frame(t(f.ACSCTY_val))
names(f.ACSCTY_valL)[1] <- "CTY_VAL"
f.ACSCTY_valL <- rownames_to_column( f.ACSCTY_valL,"value")
f.ACSCTY_moeL <- as.data.frame(t(f.ACSCTY_moe))
names(f.ACSCTY_moeL)[1] <- "CTY_MOE"
f.ACSCTY_moeL <- rownames_to_column( f.ACSCTY_moeL,"value")
f.ACSCTY_moeL$value <- gsub("_moe","",f.ACSCTY_moeL$value)
f.ACSCTY_Fin <- merge(f.ACSCTY_valL,f.ACSCTY_moeL, by="value")
if(nchar(placefips) == 0) {
# Assemble State Data
tot_OO <- as.numeric(f.ACSST$b25095001)
tot_RT <- as.numeric(f.ACSST$b25074001)
f.ACSST <- f.ACSST %>% mutate(
#Values
Med_val = b25077001,
OO_3049 = b25095006 + b25095007 + b25095008 +
b25095015 + b25095016 + b25095017 +
b25095024 + b25095025 + b25095026 +
b25095033 + b25095034 + b25095035 +
b25095042 + b25095043 + b25095044 +
b25095051 + b25095052 + b25095053 +
b25095060 + b25095061 + b25095062 +
b25095069 + b25095070 + b25095071,
OO_50 = b25095009 + b25095018 + b25095027 + b25095036 + b25095045 + b25095054 + b25095063 + b25095072,
PROP_OO_3049 = OO_3049/tot_OO,
PCT_OO_3049 = percent(PROP_OO_3049*100),
PROP_OO_50 = OO_50/tot_OO,
PCT_OO_50 = percent(PROP_OO_50*100),
Med_Rent = b25064001,
RT_3049 = b25074006 + b25074007 + b25074008 +
b25074015 + b25074016 + b25074017 +
b25074024 + b25074025 + b25074026 +
b25074033 + b25074034 + b25074035 +
b25074042 + b25074043 + b25074044 +
b25074051 + b25074052 + b25074053 +
b25074060 + b25074061 + b25074062,
RT_50 = b25074009 + b25074018 + b25074027 + b25074036 + b25074045 + b25074054 + b25074063,
PROP_RT_3049 = RT_3049/tot_RT,
PCT_RT_3049 = percent(PROP_RT_3049*100),
PROP_RT_50 = RT_50/tot_RT,
PCT_RT_50 = percent(PROP_RT_50*100),
#MOE
Med_val_moe = b25077_moe001,
OO_3049_moe = sqrt(b25095_moe006^2 + b25095_moe007^2 + b25095_moe008^2 +
b25095_moe015^2 + b25095_moe016^2 + b25095_moe017^2 +
b25095_moe024^2 + b25095_moe025^2 + b25095_moe026^2 +
b25095_moe033^2 + b25095_moe034^2 + b25095_moe035^2 +
b25095_moe042^2 + b25095_moe043^2 + b25095_moe044^2 +
b25095_moe051^2 + b25095_moe052^2 + b25095_moe053^2 +
b25095_moe060^2 + b25095_moe061^2 + b25095_moe062^2 +
b25095_moe069^2 + b25095_moe070^2 + b25095_moe071^2),
OO_50_moe = sqrt(b25095_moe009^2 + b25095_moe018^2 + b25095_moe027^2 + b25095_moe036^2 +
b25095_moe045^2 + b25095_moe054^2 + b25095_moe063^2 + b25095_moe072^2),
PROP_OO_3049_moe = OO_3049_moe/tot_OO,
PCT_OO_3049_moe = percent(PROP_OO_3049_moe*100),
PROP_OO_50_moe = OO_50_moe/tot_OO,
PCT_OO_50_moe = percent(PROP_OO_50_moe*100),
Med_Rent_moe = b25064_moe001,
RT_3049_moe = sqrt(b25074_moe006^2 + b25074_moe007^2 + b25074_moe008^2 +
b25074_moe015^2 + b25074_moe016^2 + b25074_moe017^2 +
b25074_moe024^2 + b25074_moe025^2 + b25074_moe026^2 +
b25074_moe033^2 + b25074_moe034^2 + b25074_moe035^2 +
b25074_moe042^2 + b25074_moe043^2 + b25074_moe044^2 +
b25074_moe051^2 + b25074_moe052^2 + b25074_moe053^2 +
b25074_moe060^2 + b25074_moe061^2 + b25074_moe062^2),
RT_50_moe = sqrt(b25074_moe009^2 + b25074_moe018^2 + b25074_moe027^2 +
b25074_moe036^2 + b25074_moe045^2 + b25074_moe054^2 + b25074_moe063^2),
PROP_RT_3049_moe = RT_3049_moe/tot_RT,
PCT_RT_3049_moe = percent(PROP_RT_3049_moe*100),
PROP_RT_50_moe = RT_50_moe/tot_RT,
PCT_RT_50_moe = percent(PROP_RT_50_moe*100)
)
f.ACSST_val <- f.ACSST[,c(2,129,132:135,136,139:142)]
f.ACSST_moe <- f.ACSST[,c(2,143,146:149,150,153:156)]
f.ACSST_valL <- as.data.frame(t(f.ACSST_val))
names(f.ACSST_valL)[1] <- "ST_VAL"
f.ACSST_valL <- rownames_to_column( f.ACSST_valL,"value")
f.ACSST_moeL <- as.data.frame(t(f.ACSST_moe))
names(f.ACSST_moeL)[1] <- "ST_MOE"
f.ACSST_moeL <- rownames_to_column( f.ACSST_moeL,"value")
f.ACSST_moeL$value <- gsub("_moe","",f.ACSST_moeL$value)
f.ACSST_Fin <- merge(f.ACSST_valL,f.ACSST_moeL, by="value")
} else {
# Assemble Place Data
#Raw values
tot_OO <- as.numeric(f.ACSPL$b25095001)
tot_RT <- as.numeric(f.ACSPL$b25074001)
f.ACSPL <- f.ACSPL %>% mutate(
#Values
Med_val = b25077001,
OO_3049 = b25095006 + b25095007 + b25095008 +
b25095015 + b25095016 + b25095017 +
b25095024 + b25095025 + b25095026 +
b25095033 + b25095034 + b25095035 +
b25095042 + b25095043 + b25095044 +
b25095051 + b25095052 + b25095053 +
b25095060 + b25095061 + b25095062 +
b25095069 + b25095070 + b25095071,
OO_50 = b25095009 + b25095018 + b25095027 + b25095036 + b25095045 + b25095054 + b25095063 + b25095072,
PROP_OO_3049 = OO_3049/tot_OO,
PCT_OO_3049 = percent(PROP_OO_3049*100),
PROP_OO_50 = OO_50/tot_OO,
PCT_OO_50 = percent(PROP_OO_50*100),
Med_Rent = b25064001,
RT_3049 = b25074006 + b25074007 + b25074008 +
b25074015 + b25074016 + b25074017 +
b25074024 + b25074025 + b25074026 +
b25074033 + b25074034 + b25074035 +
b25074042 + b25074043 + b25074044 +
b25074051 + b25074052 + b25074053 +
b25074060 + b25074061 + b25074062,
RT_50 = b25074009 + b25074018 + b25074027 + b25074036 + b25074045 + b25074054 + b25074063,
PROP_RT_3049 = RT_3049/tot_RT,
PCT_RT_3049 = percent(PROP_RT_3049*100),
PROP_RT_50 = RT_50/tot_RT,
PCT_RT_50 = percent(PROP_RT_50*100),
#MOE
Med_val_moe = b25077_moe001,
OO_3049_moe = sqrt(b25095_moe006^2 + b25095_moe007^2 + b25095_moe008^2 +
b25095_moe015^2 + b25095_moe016^2 + b25095_moe017^2 +
b25095_moe024^2 + b25095_moe025^2 + b25095_moe026^2 +
b25095_moe033^2 + b25095_moe034^2 + b25095_moe035^2 +
b25095_moe042^2 + b25095_moe043^2 + b25095_moe044^2 +
b25095_moe051^2 + b25095_moe052^2 + b25095_moe053^2 +
b25095_moe060^2 + b25095_moe061^2 + b25095_moe062^2 +
b25095_moe069^2 + b25095_moe070^2 + b25095_moe071^2),
OO_50_moe = sqrt(b25095_moe009^2 + b25095_moe018^2 + b25095_moe027^2 + b25095_moe036^2 +
b25095_moe045^2 + b25095_moe054^2 + b25095_moe063^2 + b25095_moe072^2),
PROP_OO_3049_moe = OO_3049_moe/tot_OO,
PCT_OO_3049_moe = percent(PROP_OO_3049_moe*100),
PROP_OO_50_moe = OO_50_moe/tot_OO,
PCT_OO_50_moe = percent(PROP_OO_50_moe*100),
Med_Rent_moe = b25064_moe001,
RT_3049_moe = sqrt(b25074_moe006^2 + b25074_moe007^2 + b25074_moe008^2 +
b25074_moe015^2 + b25074_moe016^2 + b25074_moe017^2 +
b25074_moe024^2 + b25074_moe025^2 + b25074_moe026^2 +
b25074_moe033^2 + b25074_moe034^2 + b25074_moe035^2 +
b25074_moe042^2 + b25074_moe043^2 + b25074_moe044^2 +
b25074_moe051^2 + b25074_moe052^2 + b25074_moe053^2 +
b25074_moe060^2 + b25074_moe061^2 + b25074_moe062^2),
RT_50_moe = sqrt(b25074_moe009^2 + b25074_moe018^2 + b25074_moe027^2 +
b25074_moe036^2 + b25074_moe045^2 + b25074_moe054^2 + b25074_moe063^2),
PROP_RT_3049_moe = RT_3049_moe/tot_RT,
PCT_RT_3049_moe = percent(PROP_RT_3049_moe*100),
PROP_RT_50_moe = RT_50_moe/tot_RT,
PCT_RT_50_moe = percent(PROP_RT_50_moe*100)
)
f.ACSPL_val <- f.ACSPL[,c(2,129,132:135,136,139:142)]
f.ACSPL_moe <- f.ACSPL[,c(2,143,146:149,150,153:156)]
f.ACSPL_valL <- as.data.frame(t(f.ACSPL_val))
names(f.ACSPL_valL)[1] <- "PL_VAL"
f.ACSPL_valL <- rownames_to_column( f.ACSPL_valL,"value")
f.ACSPL_moeL <- as.data.frame(t(f.ACSPL_moe))
names(f.ACSPL_moeL)[1] <- "PL_MOE"
f.ACSPL_moeL <- rownames_to_column( f.ACSPL_moeL,"value")
f.ACSPL_moeL$value <- gsub("_moe","",f.ACSPL_moeL$value)
f.ACSPL_Fin <- merge(f.ACSPL_valL,f.ACSPL_moeL, by="value")
}
# Joining Fles and calculating tests
if(nchar(placefips) == 0) {
f.HouseVal <- merge(f.ACSCTY_Fin,f.ACSST_Fin, by="value")
} else {
f.HouseVal <- merge(f.ACSPL_Fin,f.ACSCTY_Fin, by="value")
}
m.HouseVal <- as.matrix(f.HouseVal)
m.test <- matrix(nrow=11, ncol=2)
#calculating statistical test
for(i in c(2,3,8:11)) {
m.test[i,1] <-abs(as.numeric(m.HouseVal[i,2]) - as.numeric(m.HouseVal[i,4]))/sqrt((as.numeric(m.HouseVal[i,3])^2) + (as.numeric(m.HouseVal[i,5])^2))
m.test[i,2] <- ifelse(m.test[i,1] < 1,"No","Yes")
}
for(i in 4:7) {
m.test[i,1] = m.test[i+4,1]
m.test[i,2] = m.test[i+4,2]
}
m.HouseVal <- cbind(m.HouseVal,m.test)
# Bulding Renter and owner-occupied mataices for table
m.rental <- m.HouseVal[c(2,6,7),c(1:5,7)]
m.oocc <- m.HouseVal[c(3:5),c(1:5,7)]
m.rental[1,1] <- "Median Gross Rent of Rental Households (Current Dollars)"
m.rental[2,1] <- "Percentage of Rental Households paying 30-49% of income on housing"
m.rental[3,1] <- "Percentage of Rental Households paying 50% or more of income on housing"
m.oocc[1.1] <- "Median Value of Owner-Occupied Households (Current Dollars)"
m.oocc[2,1] <- "Percentage of Owner-Occupied Households paying 30-49% of income on housing"
m.oocc[3,1] <- "Percentage of Owner-Occupied Households paying 50% or more of income on housing"
#formatting values
m.rental[1,2:5] <- paste0("$",formatC(as.numeric(m.rental[1,2:5]),format="f",digits=0,big.mark=","))
m.oocc[1,2:5] <- paste0("$",formatC(as.numeric(m.oocc[1,2:5]),format="f",digits=0,big.mark=","))
#Creating output data set
m.FinTab <- rbind(m.oocc,m.rental)
f.HouseVal_Fin <- as.data.frame(m.FinTab)
if(nchar(placefips) == 0) {
names(f.HouseVal_Fin) <- c("Variable",paste0("Value: ",ctyname), paste0("MOE: ",ctyname),
paste0("Value: Colorado"), paste0("MOE: Colorado"), "Siginficant Difference?")
} else {
names(f.HouseVal_Fin) <- c("Variable",paste0("Value: ",placename), paste0("MOE: ",placename),
paste0("Value: ",ctyname), paste0("MOE: ",ctyname), "Siginficant Difference?")
}
# Setting up table
#Column Names
names_spaced <- c("Variable","Value","MOE","Value","MOE","Sig. Diff.?")
#Span Header
# create vector with colspan
if(nchar(placefips) == 0) {
tblHead1 <- c(" " = 1, ctyname = 2, "Colorado" = 2, " " = 1)
# set vector names
names(tblHead1) <- c(" ", ctyname, "Colorado", " ")
} else {
tblHead1 <- c(" " = 1, placename = 2, ctyname = 2, " " = 1)
# set vector names
names(tblHead1) <- c(" ", placename, ctyname, " ")
}
if(oType == "html") {
Housing_tab1 <- m.oocc %>%
kable(format='html', table.attr='class="cleanTable"',
row.names=FALSE,
align='lrrrrr',
caption="Comparative Owner-Occupied Housing Values",
col.names = names_spaced,
escape = FALSE) %>%
kable_styling(bootstrap_options = "condensed",full_width = F,font_size=11) %>%
row_spec(0, align = "c") %>%
column_spec(1, width = "3in") %>%
column_spec(column=2:6, width = "0.33in") %>%
add_header_above(header=tblHead1) %>%
footnote(captionSrc("ACS",ACS))
Housing_tab2 <- m.rental %>%
kable(format='html', table.attr='class="cleanTable"',
row.names=FALSE,
align='lrrrrr',
caption="Comparative Rental Housing Values",
col.names = names_spaced,
escape = FALSE) %>%
kable_styling(bootstrap_options = "condensed",full_width = F,font_size=11) %>%
row_spec(0, align = "c") %>%
column_spec(1, width = "3in") %>%
column_spec(2:6, width = "0.33in") %>%
add_header_above(header=tblHead1) %>%
footnote(captionSrc("ACS",ACS))
#Building FlexTables
# Owner-Occupied
OO <- data.frame(m.oocc)
names(OO) <- c("V1","V2","V3","V4","V5","V6")
FTOO <- regulartable(OO)
FTOO <- set_header_labels(FTOO, V1 = "Variable",
V2="Value", V3="Margin of Error",
V4="Value",V5="Margin of Error",
V6= "Significant Difference?"
)
if(nchar(placefips) == 0) {
FTOO <- add_header(FTOO,V2=ctyname,V4="Colorado",top=TRUE)
} else {
FTOO <- add_header(FTOO,V2=placename,V4=ctyname,top=TRUE)
}
FTOO <- add_header(FTOO,V1="Comparative Owner-Occupied Housing Values",top=TRUE)
FTOO <- add_footer(FTOO,V1=captionSrc("ACS",ACS))
FTOO <- merge_at(FTOO,i=1, j = 1:6, part = "header")
FTOO <- merge_at(FTOO,i=2,j=2:3,part="header")
FTOO <- merge_at(FTOO,i=2,j=4:5,part="header")
FTOO <- merge_at(FTOO,i=1, j = 1:6, part = "footer")
FTOO <- align(FTOO,i=1, j=1, align="left",part="header")
FTOO <- align(FTOO,i=2:3, j=2:6, align="center",part="header")
FTOO <- align(FTOO,i=1, align="left",part="footer")
FTOO <- align(FTOO, j=1, align="left", part="body")
FTOO <- autofit(FTOO)
FTOO <- width(FTOO,j=1, width=3)
FTOO <- width(FTOO,j=2:6, width=1)
# Rental
RT <- data.frame(m.rental)
names(RT) <- c("V1","V2","V3","V4","V5","V6")
FTRT <- regulartable(RT)
FTRT <- set_header_labels(FTRT, V1 = "Variable",
V2="Value", V3="Margin of Error",
V4="Value",V5="Margin of Error",
V6= "Significant Difference?"
)
if(nchar(placefips) == 0) {
FTRT <- add_header(FTRT,V2=ctyname,V4="Colorado",top=TRUE)
} else {
FTRT <- add_header(FTRT,V2=placename,V4=ctyname,top=TRUE)
}
FTRT <- add_header(FTRT,V1="Comparative Rental Housing Values",top=TRUE)
FTRT <- add_footer(FTRT,V1=captionSrc("ACS",ACS))
FTRT <- merge_at(FTRT,i=1, j = 1:6, part = "header")
FTRT <- merge_at(FTRT,i=2,j=2:3,part="header")
FTRT <- merge_at(FTRT,i=2,j=4:5,part="header")
FTRT <- merge_at(FTRT,i=1, j = 1:6, part = "footer")
FTRT <- align(FTRT,i=1, j=1, align="left",part="header")
FTRT <- align(FTRT,i=2:3, j=2:6, align="center",part="header")
FTRT <- align(FTRT,i=1, align="left",part="footer")
FTRT <- align(FTRT, j=1, align="left", part="body")
FTRT <- autofit(FTRT)
FTRT <- width(FTRT,j=1, width=3)
FTRT <- width(FTRT,j=2:6, width=1)
outList <- list("OOTab" = Housing_tab1, "RTTab" = Housing_tab2, "data" = f.HouseVal_Fin,
"FlexTableOO" = FTOO,"FlexTableRT" = FTRT)
return(outList)
}
if(oType == "latex") {
Housing_tab1 <- kable(m.oocc,
col.names = names_spaced,
align=c("lrrrrr"),
caption="Comparison of Owner-Occupied Housing Values", row.names=FALSE,
format="latex", booktabs=TRUE) %>%
kable_styling(latex_options=c("scale_down","HOLD_position"),font_size=10) %>%
row_spec(0, align = "c") %>%
column_spec(1,width="3in") %>%
add_header_above(header=tblHead1) %>%
footnote(captionSrc("ACS",ACS))
Housing_tab2 <- kable(m.rental,
col.names = names_spaced,
align=c("lrrrrr"),
caption="Comparison of Rental Housing Values", row.names=FALSE,
format="latex", booktabs=TRUE) %>%
kable_styling(latex_options=c("scale_down","HOLD_position"),font_size=10) %>%
row_spec(0, align = "c") %>%
column_spec(1,width="3in") %>%
add_header_above(header=tblHead1) %>%
footnote(captionSrc("ACS",ACS))
outList <- list("OOTab" = Housing_tab1, "RTTab" = Housing_tab2)
return(outList)
}
}
ColoradoDemography/codemogLib documentation built on Dec. 15, 2020, 2:42 a.m.
R Package Documentation
Browse R Packages
We want your feedback!
Note that we can't provide technical support on individual packages. You should contact the package authors for that.
Defines functions HouseVal
#' HouseVal Table Showing the Median House value for Owner-Occupied housing,
#' the Median Gross Rent, and Median Costs as a Percentage of Income for
#' Owners and Renters for a place and the State of Colorado
#'
#' Restructured 5/4/2018 to remove codemog_api calls
#'
#' @param listID the list containing place id and Place names
#' @param ACS Specifies the ACS data set to be used, reads curACS from Shiny program
#' @param oType output type html table or latex table
#' @return kable formatted table and data file
#' @export
#'
HouseVal <- function(listID, ACS, oType, state="08"){
# Collecting place ids from idList, setting default values
ctyfips <- listID$ctyNum
ctyname <- listID$ctyName
placefips <- listID$plNum
placename <- listID$plName
# if(listID$PlFilter == "T") {
# placefips <- ""
# placename <- ""
# }
#Prepping SQL calls
#County SQL Calls
CTYSQL01 <- paste0("SELECT * FROM data.houseval_cty WHERE (geonum = ",paste0("1", state, ctyfips)," AND acs = '",ACS,"');")
if(nchar(placefips) == 0) {
#State SQL Calls
STSQL01 <- paste0("SELECT * FROM data.houseval_st WHERE acs = '",ACS,"';")
} else {
# Place SQL calls
PLSQL01 <- paste0("SELECT * FROM data.houseval_muni WHERE (geonum = ",paste0("1", state, placefips)," AND acs = '",ACS,"');")
}
# Call to Postrgres
pw <- {
"demography"
}
# loads the PostgreSQL driver
drv <- dbDriver("PostgreSQL")
# creates a connection to the postgres database
# note that "con" will be used later in each connection to the database
con <- dbConnect(drv, dbname = 'dola',
host = "104.197.26.248", port = 5433,
user = "codemog", password = pw)
rm(pw) # removes the password
# Read data files
# Raw County data
f.ACSCTY <- dbGetQuery(con, CTYSQL01)
if(nchar(placefips) == 0) {
# Raw State data
f.ACSST <- dbGetQuery(con, STSQL01)
} else {
# Raw Place data
f.ACSPL <- dbGetQuery(con, PLSQL01)
}
#closing the connections
dbDisconnect(con)
dbUnloadDriver(drv)
rm(con)
rm(drv)
#Accounting for Missing value
f.ACSCTY[f.ACSCTY == "-9999"] <- NA
if(nchar(placefips) == 0) {
f.ACSST[f.ACSST == "-9999"] <- NA
} else {
f.ACSPL[f.ACSPL == "-9999"] <- NA
}
# Assemble County Data
tot_OO <- as.numeric(f.ACSCTY$b25095001)
tot_RT <- as.numeric(f.ACSCTY$b25074001)
f.ACSCTY <- f.ACSCTY %>% mutate(
#Values
Med_val = b25077001,
OO_3049 = b25095006 + b25095007 + b25095008 +
b25095015 + b25095016 + b25095017 +
b25095024 + b25095025 + b25095026 +
b25095033 + b25095034 + b25095035 +
b25095042 + b25095043 + b25095044 +
b25095051 + b25095052 + b25095053 +
b25095060 + b25095061 + b25095062 +
b25095069 + b25095070 + b25095071,
OO_50 = b25095009 + b25095018 + b25095027 + b25095036 + b25095045 + b25095054 + b25095063 + b25095072,
PROP_OO_3049 = OO_3049/tot_OO,
PCT_OO_3049 = percent(PROP_OO_3049*100),
PROP_OO_50 = OO_50/tot_OO,
PCT_OO_50 = percent(PROP_OO_50*100),
Med_Rent = b25064001,
RT_3049 = b25074006 + b25074007 + b25074008 +
b25074015 + b25074016 + b25074017 +
b25074024 + b25074025 + b25074026 +
b25074033 + b25074034 + b25074035 +
b25074042 + b25074043 + b25074044 +
b25074051 + b25074052 + b25074053 +
b25074060 + b25074061 + b25074062,
RT_50 = b25074009 + b25074018 + b25074027 + b25074036 + b25074045 + b25074054 + b25074063,
PROP_RT_3049 = RT_3049/tot_RT,
PCT_RT_3049 = percent(PROP_RT_3049*100),
PROP_RT_50 = RT_50/tot_RT,
PCT_RT_50 = percent(PROP_RT_50*100),
#MOE
Med_val_moe = b25077_moe001,
OO_3049_moe = sqrt(b25095_moe006^2 + b25095_moe007^2 + b25095_moe008^2 +
b25095_moe015^2 + b25095_moe016^2 + b25095_moe017^2 +
b25095_moe024^2 + b25095_moe025^2 + b25095_moe026^2 +
b25095_moe033^2 + b25095_moe034^2 + b25095_moe035^2 +
b25095_moe042^2 + b25095_moe043^2 + b25095_moe044^2 +
b25095_moe051^2 + b25095_moe052^2 + b25095_moe053^2 +
b25095_moe060^2 + b25095_moe061^2 + b25095_moe062^2 +
b25095_moe069^2 + b25095_moe070^2 + b25095_moe071^2),
OO_50_moe = sqrt(b25095_moe009^2 + b25095_moe018^2 + b25095_moe027^2 + b25095_moe036^2 +
b25095_moe045^2 + b25095_moe054^2 + b25095_moe063^2 + b25095_moe072^2),
PROP_OO_3049_moe = OO_3049_moe/tot_OO,
PCT_OO_3049_moe = percent(PROP_OO_3049_moe*100),
PROP_OO_50_moe = OO_50_moe/tot_OO,
PCT_OO_50_moe = percent(PROP_OO_50_moe*100),
Med_Rent_moe = b25064_moe001,
RT_3049_moe = sqrt(b25074_moe006^2 + b25074_moe007^2 + b25074_moe008^2 +
b25074_moe015^2 + b25074_moe016^2 + b25074_moe017^2 +
b25074_moe024^2 + b25074_moe025^2 + b25074_moe026^2 +
b25074_moe033^2 + b25074_moe034^2 + b25074_moe035^2 +
b25074_moe042^2 + b25074_moe043^2 + b25074_moe044^2 +
b25074_moe051^2 + b25074_moe052^2 + b25074_moe053^2 +
b25074_moe060^2 + b25074_moe061^2 + b25074_moe062^2),
RT_50_moe = sqrt(b25074_moe009^2 + b25074_moe018^2 + b25074_moe027^2 +
b25074_moe036^2 + b25074_moe045^2 + b25074_moe054^2 + b25074_moe063^2),
PROP_RT_3049_moe = RT_3049_moe/tot_RT,
PCT_RT_3049_moe = percent(PROP_RT_3049_moe*100),
PROP_RT_50_moe = RT_50_moe/tot_RT,
PCT_RT_50_moe = percent(PROP_RT_50_moe*100)
)
f.ACSCTY_val <- f.ACSCTY[,c(2,129,132:135,136,139:142)]
f.ACSCTY_moe <- f.ACSCTY[,c(2,143,146:149,150,153:156)]
f.ACSCTY_valL <- as.data.frame(t(f.ACSCTY_val))
names(f.ACSCTY_valL)[1] <- "CTY_VAL"
f.ACSCTY_valL <- rownames_to_column( f.ACSCTY_valL,"value")
f.ACSCTY_moeL <- as.data.frame(t(f.ACSCTY_moe))
names(f.ACSCTY_moeL)[1] <- "CTY_MOE"
f.ACSCTY_moeL <- rownames_to_column( f.ACSCTY_moeL,"value")
f.ACSCTY_moeL$value <- gsub("_moe","",f.ACSCTY_moeL$value)
f.ACSCTY_Fin <- merge(f.ACSCTY_valL,f.ACSCTY_moeL, by="value")
if(nchar(placefips) == 0) {
# Assemble State Data
tot_OO <- as.numeric(f.ACSST$b25095001)
tot_RT <- as.numeric(f.ACSST$b25074001)
f.ACSST <- f.ACSST %>% mutate(
#Values
Med_val = b25077001,
OO_3049 = b25095006 + b25095007 + b25095008 +
b25095015 + b25095016 + b25095017 +
b25095024 + b25095025 + b25095026 +
b25095033 + b25095034 + b25095035 +
b25095042 + b25095043 + b25095044 +
b25095051 + b25095052 + b25095053 +
b25095060 + b25095061 + b25095062 +
b25095069 + b25095070 + b25095071,
OO_50 = b25095009 + b25095018 + b25095027 + b25095036 + b25095045 + b25095054 + b25095063 + b25095072,
PROP_OO_3049 = OO_3049/tot_OO,
PCT_OO_3049 = percent(PROP_OO_3049*100),
PROP_OO_50 = OO_50/tot_OO,
PCT_OO_50 = percent(PROP_OO_50*100),
Med_Rent = b25064001,
RT_3049 = b25074006 + b25074007 + b25074008 +
b25074015 + b25074016 + b25074017 +
b25074024 + b25074025 + b25074026 +
b25074033 + b25074034 + b25074035 +
b25074042 + b25074043 + b25074044 +
b25074051 + b25074052 + b25074053 +
b25074060 + b25074061 + b25074062,
RT_50 = b25074009 + b25074018 + b25074027 + b25074036 + b25074045 + b25074054 + b25074063,
PROP_RT_3049 = RT_3049/tot_RT,
PCT_RT_3049 = percent(PROP_RT_3049*100),
PROP_RT_50 = RT_50/tot_RT,
PCT_RT_50 = percent(PROP_RT_50*100),
#MOE
Med_val_moe = b25077_moe001,
OO_3049_moe = sqrt(b25095_moe006^2 + b25095_moe007^2 + b25095_moe008^2 +
b25095_moe015^2 + b25095_moe016^2 + b25095_moe017^2 +
b25095_moe024^2 + b25095_moe025^2 + b25095_moe026^2 +
b25095_moe033^2 + b25095_moe034^2 + b25095_moe035^2 +
b25095_moe042^2 + b25095_moe043^2 + b25095_moe044^2 +
b25095_moe051^2 + b25095_moe052^2 + b25095_moe053^2 +
b25095_moe060^2 + b25095_moe061^2 + b25095_moe062^2 +
b25095_moe069^2 + b25095_moe070^2 + b25095_moe071^2),
OO_50_moe = sqrt(b25095_moe009^2 + b25095_moe018^2 + b25095_moe027^2 + b25095_moe036^2 +
b25095_moe045^2 + b25095_moe054^2 + b25095_moe063^2 + b25095_moe072^2),
PROP_OO_3049_moe = OO_3049_moe/tot_OO,
PCT_OO_3049_moe = percent(PROP_OO_3049_moe*100),
PROP_OO_50_moe = OO_50_moe/tot_OO,
PCT_OO_50_moe = percent(PROP_OO_50_moe*100),
Med_Rent_moe = b25064_moe001,
RT_3049_moe = sqrt(b25074_moe006^2 + b25074_moe007^2 + b25074_moe008^2 +
b25074_moe015^2 + b25074_moe016^2 + b25074_moe017^2 +
b25074_moe024^2 + b25074_moe025^2 + b25074_moe026^2 +
b25074_moe033^2 + b25074_moe034^2 + b25074_moe035^2 +
b25074_moe042^2 + b25074_moe043^2 + b25074_moe044^2 +
b25074_moe051^2 + b25074_moe052^2 + b25074_moe053^2 +
b25074_moe060^2 + b25074_moe061^2 + b25074_moe062^2),
RT_50_moe = sqrt(b25074_moe009^2 + b25074_moe018^2 + b25074_moe027^2 +
b25074_moe036^2 + b25074_moe045^2 + b25074_moe054^2 + b25074_moe063^2),
PROP_RT_3049_moe = RT_3049_moe/tot_RT,
PCT_RT_3049_moe = percent(PROP_RT_3049_moe*100),
PROP_RT_50_moe = RT_50_moe/tot_RT,
PCT_RT_50_moe = percent(PROP_RT_50_moe*100)
)
f.ACSST_val <- f.ACSST[,c(2,129,132:135,136,139:142)]
f.ACSST_moe <- f.ACSST[,c(2,143,146:149,150,153:156)]
f.ACSST_valL <- as.data.frame(t(f.ACSST_val))
names(f.ACSST_valL)[1] <- "ST_VAL"
f.ACSST_valL <- rownames_to_column( f.ACSST_valL,"value")
f.ACSST_moeL <- as.data.frame(t(f.ACSST_moe))
names(f.ACSST_moeL)[1] <- "ST_MOE"
f.ACSST_moeL <- rownames_to_column( f.ACSST_moeL,"value")
f.ACSST_moeL$value <- gsub("_moe","",f.ACSST_moeL$value)
f.ACSST_Fin <- merge(f.ACSST_valL,f.ACSST_moeL, by="value")
} else {
# Assemble Place Data
#Raw values
tot_OO <- as.numeric(f.ACSPL$b25095001)
tot_RT <- as.numeric(f.ACSPL$b25074001)
f.ACSPL <- f.ACSPL %>% mutate(
#Values
Med_val = b25077001,
OO_3049 = b25095006 + b25095007 + b25095008 +
b25095015 + b25095016 + b25095017 +
b25095024 + b25095025 + b25095026 +
b25095033 + b25095034 + b25095035 +
b25095042 + b25095043 + b25095044 +
b25095051 + b25095052 + b25095053 +
b25095060 + b25095061 + b25095062 +
b25095069 + b25095070 + b25095071,
OO_50 = b25095009 + b25095018 + b25095027 + b25095036 + b25095045 + b25095054 + b25095063 + b25095072,
PROP_OO_3049 = OO_3049/tot_OO,
PCT_OO_3049 = percent(PROP_OO_3049*100),
PROP_OO_50 = OO_50/tot_OO,
PCT_OO_50 = percent(PROP_OO_50*100),
Med_Rent = b25064001,
RT_3049 = b25074006 + b25074007 + b25074008 +
b25074015 + b25074016 + b25074017 +
b25074024 + b25074025 + b25074026 +
b25074033 + b25074034 + b25074035 +
b25074042 + b25074043 + b25074044 +
b25074051 + b25074052 + b25074053 +
b25074060 + b25074061 + b25074062,
RT_50 = b25074009 + b25074018 + b25074027 + b25074036 + b25074045 + b25074054 + b25074063,
PROP_RT_3049 = RT_3049/tot_RT,
PCT_RT_3049 = percent(PROP_RT_3049*100),
PROP_RT_50 = RT_50/tot_RT,
PCT_RT_50 = percent(PROP_RT_50*100),
#MOE
Med_val_moe = b25077_moe001,
OO_3049_moe = sqrt(b25095_moe006^2 + b25095_moe007^2 + b25095_moe008^2 +
b25095_moe015^2 + b25095_moe016^2 + b25095_moe017^2 +
b25095_moe024^2 + b25095_moe025^2 + b25095_moe026^2 +
b25095_moe033^2 + b25095_moe034^2 + b25095_moe035^2 +
b25095_moe042^2 + b25095_moe043^2 + b25095_moe044^2 +
b25095_moe051^2 + b25095_moe052^2 + b25095_moe053^2 +
b25095_moe060^2 + b25095_moe061^2 + b25095_moe062^2 +
b25095_moe069^2 + b25095_moe070^2 + b25095_moe071^2),
OO_50_moe = sqrt(b25095_moe009^2 + b25095_moe018^2 + b25095_moe027^2 + b25095_moe036^2 +
b25095_moe045^2 + b25095_moe054^2 + b25095_moe063^2 + b25095_moe072^2),
PROP_OO_3049_moe = OO_3049_moe/tot_OO,
PCT_OO_3049_moe = percent(PROP_OO_3049_moe*100),
PROP_OO_50_moe = OO_50_moe/tot_OO,
PCT_OO_50_moe = percent(PROP_OO_50_moe*100),
Med_Rent_moe = b25064_moe001,
RT_3049_moe = sqrt(b25074_moe006^2 + b25074_moe007^2 + b25074_moe008^2 +
b25074_moe015^2 + b25074_moe016^2 + b25074_moe017^2 +
b25074_moe024^2 + b25074_moe025^2 + b25074_moe026^2 +
b25074_moe033^2 + b25074_moe034^2 + b25074_moe035^2 +
b25074_moe042^2 + b25074_moe043^2 + b25074_moe044^2 +
b25074_moe051^2 + b25074_moe052^2 + b25074_moe053^2 +
b25074_moe060^2 + b25074_moe061^2 + b25074_moe062^2),
RT_50_moe = sqrt(b25074_moe009^2 + b25074_moe018^2 + b25074_moe027^2 +
b25074_moe036^2 + b25074_moe045^2 + b25074_moe054^2 + b25074_moe063^2),
PROP_RT_3049_moe = RT_3049_moe/tot_RT,
PCT_RT_3049_moe = percent(PROP_RT_3049_moe*100),
PROP_RT_50_moe = RT_50_moe/tot_RT,
PCT_RT_50_moe = percent(PROP_RT_50_moe*100)
)
f.ACSPL_val <- f.ACSPL[,c(2,129,132:135,136,139:142)]
f.ACSPL_moe <- f.ACSPL[,c(2,143,146:149,150,153:156)]
f.ACSPL_valL <- as.data.frame(t(f.ACSPL_val))
names(f.ACSPL_valL)[1] <- "PL_VAL"
f.ACSPL_valL <- rownames_to_column( f.ACSPL_valL,"value")
f.ACSPL_moeL <- as.data.frame(t(f.ACSPL_moe))
names(f.ACSPL_moeL)[1] <- "PL_MOE"
f.ACSPL_moeL <- rownames_to_column( f.ACSPL_moeL,"value")
f.ACSPL_moeL$value <- gsub("_moe","",f.ACSPL_moeL$value)
f.ACSPL_Fin <- merge(f.ACSPL_valL,f.ACSPL_moeL, by="value")
}
# Joining Fles and calculating tests
if(nchar(placefips) == 0) {
f.HouseVal <- merge(f.ACSCTY_Fin,f.ACSST_Fin, by="value")
} else {
f.HouseVal <- merge(f.ACSPL_Fin,f.ACSCTY_Fin, by="value")
}
m.HouseVal <- as.matrix(f.HouseVal)
m.test <- matrix(nrow=11, ncol=2)
#calculating statistical test
for(i in c(2,3,8:11)) {
m.test[i,1] <-abs(as.numeric(m.HouseVal[i,2]) - as.numeric(m.HouseVal[i,4]))/sqrt((as.numeric(m.HouseVal[i,3])^2) + (as.numeric(m.HouseVal[i,5])^2))
m.test[i,2] <- ifelse(m.test[i,1] < 1,"No","Yes")
}
for(i in 4:7) {
m.test[i,1] = m.test[i+4,1]
m.test[i,2] = m.test[i+4,2]
}
m.HouseVal <- cbind(m.HouseVal,m.test)
# Bulding Renter and owner-occupied mataices for table
m.rental <- m.HouseVal[c(2,6,7),c(1:5,7)]
m.oocc <- m.HouseVal[c(3:5),c(1:5,7)]
m.rental[1,1] <- "Median Gross Rent of Rental Households (Current Dollars)"
m.rental[2,1] <- "Percentage of Rental Households paying 30-49% of income on housing"
m.rental[3,1] <- "Percentage of Rental Households paying 50% or more of income on housing"
m.oocc[1.1] <- "Median Value of Owner-Occupied Households (Current Dollars)"
m.oocc[2,1] <- "Percentage of Owner-Occupied Households paying 30-49% of income on housing"
m.oocc[3,1] <- "Percentage of Owner-Occupied Households paying 50% or more of income on housing"
#formatting values
m.rental[1,2:5] <- paste0("$",formatC(as.numeric(m.rental[1,2:5]),format="f",digits=0,big.mark=","))
m.oocc[1,2:5] <- paste0("$",formatC(as.numeric(m.oocc[1,2:5]),format="f",digits=0,big.mark=","))
#Creating output data set
m.FinTab <- rbind(m.oocc,m.rental)
f.HouseVal_Fin <- as.data.frame(m.FinTab)
if(nchar(placefips) == 0) {
names(f.HouseVal_Fin) <- c("Variable",paste0("Value: ",ctyname), paste0("MOE: ",ctyname),
paste0("Value: Colorado"), paste0("MOE: Colorado"), "Siginficant Difference?")
} else {
names(f.HouseVal_Fin) <- c("Variable",paste0("Value: ",placename), paste0("MOE: ",placename),
paste0("Value: ",ctyname), paste0("MOE: ",ctyname), "Siginficant Difference?")
}
# Setting up table
#Column Names
names_spaced <- c("Variable","Value","MOE","Value","MOE","Sig. Diff.?")
#Span Header
# create vector with colspan
if(nchar(placefips) == 0) {
tblHead1 <- c(" " = 1, ctyname = 2, "Colorado" = 2, " " = 1)
# set vector names
names(tblHead1) <- c(" ", ctyname, "Colorado", " ")
} else {
tblHead1 <- c(" " = 1, placename = 2, ctyname = 2, " " = 1)
# set vector names
names(tblHead1) <- c(" ", placename, ctyname, " ")
}
if(oType == "html") {
Housing_tab1 <- m.oocc %>%
kable(format='html', table.attr='class="cleanTable"',
row.names=FALSE,
align='lrrrrr',
caption="Comparative Owner-Occupied Housing Values",
col.names = names_spaced,
escape = FALSE) %>%
kable_styling(bootstrap_options = "condensed",full_width = F,font_size=11) %>%
row_spec(0, align = "c") %>%
column_spec(1, width = "3in") %>%
column_spec(column=2:6, width = "0.33in") %>%
add_header_above(header=tblHead1) %>%
footnote(captionSrc("ACS",ACS))
Housing_tab2 <- m.rental %>%
kable(format='html', table.attr='class="cleanTable"',
row.names=FALSE,
align='lrrrrr',
caption="Comparative Rental Housing Values",
col.names = names_spaced,
escape = FALSE) %>%
kable_styling(bootstrap_options = "condensed",full_width = F,font_size=11) %>%
row_spec(0, align = "c") %>%
column_spec(1, width = "3in") %>%
column_spec(2:6, width = "0.33in") %>%
add_header_above(header=tblHead1) %>%
footnote(captionSrc("ACS",ACS))
#Building FlexTables
# Owner-Occupied
OO <- data.frame(m.oocc)
names(OO) <- c("V1","V2","V3","V4","V5","V6")
FTOO <- regulartable(OO)
FTOO <- set_header_labels(FTOO, V1 = "Variable",
V2="Value", V3="Margin of Error",
V4="Value",V5="Margin of Error",
V6= "Significant Difference?"
)
if(nchar(placefips) == 0) {
FTOO <- add_header(FTOO,V2=ctyname,V4="Colorado",top=TRUE)
} else {
FTOO <- add_header(FTOO,V2=placename,V4=ctyname,top=TRUE)
}
FTOO <- add_header(FTOO,V1="Comparative Owner-Occupied Housing Values",top=TRUE)
FTOO <- add_footer(FTOO,V1=captionSrc("ACS",ACS))
FTOO <- merge_at(FTOO,i=1, j = 1:6, part = "header")
FTOO <- merge_at(FTOO,i=2,j=2:3,part="header")
FTOO <- merge_at(FTOO,i=2,j=4:5,part="header")
FTOO <- merge_at(FTOO,i=1, j = 1:6, part = "footer")
FTOO <- align(FTOO,i=1, j=1, align="left",part="header")
FTOO <- align(FTOO,i=2:3, j=2:6, align="center",part="header")
FTOO <- align(FTOO,i=1, align="left",part="footer")
FTOO <- align(FTOO, j=1, align="left", part="body")
FTOO <- autofit(FTOO)
FTOO <- width(FTOO,j=1, width=3)
FTOO <- width(FTOO,j=2:6, width=1)
# Rental
RT <- data.frame(m.rental)
names(RT) <- c("V1","V2","V3","V4","V5","V6")
FTRT <- regulartable(RT)
FTRT <- set_header_labels(FTRT, V1 = "Variable",
V2="Value", V3="Margin of Error",
V4="Value",V5="Margin of Error",
V6= "Significant Difference?"
)
if(nchar(placefips) == 0) {
FTRT <- add_header(FTRT,V2=ctyname,V4="Colorado",top=TRUE)
} else {
FTRT <- add_header(FTRT,V2=placename,V4=ctyname,top=TRUE)
}
FTRT <- add_header(FTRT,V1="Comparative Rental Housing Values",top=TRUE)
FTRT <- add_footer(FTRT,V1=captionSrc("ACS",ACS))
FTRT <- merge_at(FTRT,i=1, j = 1:6, part = "header")
FTRT <- merge_at(FTRT,i=2,j=2:3,part="header")
FTRT <- merge_at(FTRT,i=2,j=4:5,part="header")
FTRT <- merge_at(FTRT,i=1, j = 1:6, part = "footer")
FTRT <- align(FTRT,i=1, j=1, align="left",part="header")
FTRT <- align(FTRT,i=2:3, j=2:6, align="center",part="header")
FTRT <- align(FTRT,i=1, align="left",part="footer")
FTRT <- align(FTRT, j=1, align="left", part="body")
FTRT <- autofit(FTRT)
FTRT <- width(FTRT,j=1, width=3)
FTRT <- width(FTRT,j=2:6, width=1)
outList <- list("OOTab" = Housing_tab1, "RTTab" = Housing_tab2, "data" = f.HouseVal_Fin,
"FlexTableOO" = FTOO,"FlexTableRT" = FTRT)
return(outList)
}
if(oType == "latex") {
Housing_tab1 <- kable(m.oocc,
col.names = names_spaced,
align=c("lrrrrr"),
caption="Comparison of Owner-Occupied Housing Values", row.names=FALSE,
format="latex", booktabs=TRUE) %>%
kable_styling(latex_options=c("scale_down","HOLD_position"),font_size=10) %>%
row_spec(0, align = "c") %>%
column_spec(1,width="3in") %>%
add_header_above(header=tblHead1) %>%
footnote(captionSrc("ACS",ACS))
Housing_tab2 <- kable(m.rental,
col.names = names_spaced,
align=c("lrrrrr"),
caption="Comparison of Rental Housing Values", row.names=FALSE,
format="latex", booktabs=TRUE) %>%
kable_styling(latex_options=c("scale_down","HOLD_position"),font_size=10) %>%
row_spec(0, align = "c") %>%
column_spec(1,width="3in") %>%
add_header_above(header=tblHead1) %>%
footnote(captionSrc("ACS",ACS))
outList <- list("OOTab" = Housing_tab1, "RTTab" = Housing_tab2)
return(outList)
}
}
R Package Documentation
Browse R Packages
We want your feedback!
Note that we can't provide technical support on individual packages. You should contact the package authors for that.
Embedding an R snippet on your website
Add the following code to your website.
For more information on customizing the embed code, read Embedding Snippets.