#' Change Names of Columns in a Model Selection Table from SAS GLIMMIX Model Output
#'
#' @param y
#'
#' @export
processColumnNames <- function(y) {
# change column names
if ("Columns in X" %in% names(y)) {
setnames(y, "Columns in X", "Number of Fixed Effects Parameters")
}
if (length(grep("Columns in Z", names(y), fixed=T)) > 0) {
setnames(
y,
names(y)[grep("Columns in Z", names(y), fixed=T)],
"ColumnsZ"
)
}
# plant size
if ("C_t" %in% names(y)) {setnames(y, "C_t", "$C_t$")}
# insects
if ("CA" %in% names(y)) {setnames(y, "CA", "Invasive Moth")}
if ("CH" %in% names(y)) {setnames(y, "CH", "Native Bug")}
if ("DA" %in% names(y)) {setnames(y, "DA", "Native Scale")}
if ("ME" %in% names(y)) {setnames(y, "ME", "Native Moth")}
if ("MEhistory" %in% names(y)) {setnames(y, "MEhistory", "$\\nu_h$")}
if ("CAhistory" %in% names(y)) {setnames(y, "CAhistory", "$\\mu_h$")}
if ("Old_Moth_Evidence_t_" %in% names(y)) {setnames(y, "Old_Moth_Evidence_t_", "Moth Damage Evidence")}
if ("OldMothPlantPres" %in% names(y)) {setnames(y, "OldMothPlantPres", "Moth Damage Presence During the Study")}
if ("OldMothhistory" %in% names(y)) {setnames(y, "OldMothhistory", "$\\delta_h$")}
if ("CA_t_1" %in% names(y)) {setnames(y, "CA_t_1", "Invasive Moth")}#"$\\mu_t$")}
if ("CAyr_t_1" %in% names(y)) {setnames(y, "CAyr_t_1", "Invasive Moth")}
if ("CH_t_1" %in% names(y)) {setnames(y, "CH_t_1", "Native Bug")}
if ("CHyr_t_1" %in% names(y)) {setnames(y, "CHyr_t_1", "Native Bug")}
if ("DA_t_1" %in% names(y)) {setnames(y, "DA_t_1", "Native Scale")}
if ("DAyr_t_1" %in% names(y)) {setnames(y, "DAyr_t_1", "Native Scale")}
if ("ME_t_1" %in% names(y)) {setnames(y, "ME_t_1", "Native Moth")} # "$\\nu_t$")}
if ("MEyr_t_1" %in% names(y)) {setnames(y, "MEyr_t_1", "Native Moth")}
if ("NatInsect_t_1" %in% names(y)) {
setnames(y, "NatInsect_t_1", "Native Insects")
}
if ("NatInsectyr_t_" %in% names(y)) {
setnames(y, "NatInsectyr_t_", "Native Insects")
}
# plant size
if ("Ln_Size_t_1_st" %in% names(y)) {setnames(y, "Ln_Size_t_1_st", "C_t")}
if ("Ln_Cone_t_1_st" %in% names(y)) {
setnames(
y,
"Ln_Cone_t_1_st",
"ln(Cone Volume), Stand."
)
}
if (length(grep("Ln_Cylinder", names(y), fixed=T)) > 0) {
setnames(
y,
names(y)[grep("Ln_Cylinder", names(y), fixed=T)],
"ln(Cylinder Volume), Stand."
)
}
# interactions
if ("CA_t_1*CH_t_1" %in% names(y)) {
setnames(y, "CA_t_1*CH_t_1", "Invasive Moth x Native Bug")
}
if (length(grep("CA_t_1*NatInse", names(y), fixed=T)) > 0) {
setnames(
y,
names(y)[grep("CA_t_1*NatInse", names(y), fixed=T)],
"Invasive Moth x Native Insects"
)
}
if ("P1*CA_t_1" %in% names(y)) {
setnames(y, "P1*CA_t_1", "Invasive Moth x Precipitation")
}
if ("T1*CA_t_1" %in% names(y)) {
setnames(y, "T1*CA_t_1", "Invasive Moth x Temperature")
}
if ("P1*CH_t_1" %in% names(y)) {
setnames(y, "P1*CH_t_1", "Native Bug x Precipitation")
}
if ("T1*CH_t_1" %in% names(y)) {
setnames(y, "T1*CH_t_1", "Native Bug x Temperature")
}
if ("P1*DA_t_1" %in% names(y)) {
setnames(y, "P1*DA_t_1", "Native Scale x Precipitation")
}
if ("T1*DA_t_1" %in% names(y)) {
setnames(y, "T1*DA_t_1", "Native Scale x Temperature")
}
if ("P1*ME_t_1" %in% names(y)) {
setnames(y, "P1*ME_t_1", "Native Moth x Precipitation")
}
if ("T1*ME_t_1" %in% names(y)) {
setnames(y, "T1*ME_t_1", "Native Moth x Temperature")
}
return(y)
}
#' Change Names of Precipitation and Temperature Columns
#'
#' @param y
#'
#' @export
processT1P1ColumnNames <- function(y) {
if ("T1" %in% names(y)) {setnames(y, "T1", "T")}
if ("P1" %in% names(y)) {setnames(y, "P1", "P")}
}
#' Change Names of Model Selection Columns
#'
#' @param y
#'
#' @export
processModSelectionColumnNames <- function(y) {
if ("Intercept" %in% names(y)) {setnames(y, "Intercept", "$\\beta_{0}$")}
if ("Value" %in% names(y)) {setnames(y, "Value", "AIC")}
if ("modelVars" %in% names(y)) {setnames(y, "modelVars", "Model")}
if ("delta AIC" %in% names(y)) {setnames(y, "delta AIC", "$\\Delta$AIC")}
if ("Model Lik" %in% names(y)) {setnames(y, "Model Lik", "Model $\\Lagr$")}
if ("Prob(Model)" %in% names(y)) {setnames(y, "Prob(Model)", "$P$(Model)")}
return(y)
}
#' Change Names of Statistics Columns
#'
#' @param y
#'
#' @export
processStatisticsColumnNames <- function(y) {
if ("Pearson Chi-Square/DF" %in% names(y)) {setnames(y, "Pearson Chi-Square/DF", "Pearson $\\chi^2$/DF")}
if ("ProbChiSq" %in% names(y)) {setnames(y, "ProbChiSq", "$P$($\\chi^2$)")}
if ("ChiSq" %in% names(y)) {setnames(y, "ChiSq", "$\\chi^2$")}
if ("StdErr" %in% names(y)) {setnames(y, "StdErr", "Standard Error")}
if ("LowerWaldCL" %in% names(y)) {setnames(y, "LowerWaldCL", "Lower Wald CL")}
if ("UpperWaldCL" %in% names(y)) {setnames(y, "UpperWaldCL", "Upper Wald CL")}
return(y)
}
processModSelectionColumnNames <- function(y) {
if ("Intercept" %in% names(y)) {setnames(y, "Intercept", "$\\beta_{0}$")}
if ("Value" %in% names(y)) {setnames(y, "Value", "AIC")}
if ("modelVars" %in% names(y)) {setnames(y, "modelVars", "Model")}
if ("delta AIC" %in% names(y)) {setnames(y, "delta AIC", "$\\Delta$AIC")}
if ("Model Lik" %in% names(y)) {setnames(y, "Model Lik", "Model $\\Lagr$")}
if ("Prob(Model)" %in% names(y)) {setnames(y, "Prob(Model)", "$P$(Model)")}
return(y)
}
#' Abbreviate Insect Column Names
#'
#' @param y
#'
#' @export
abbrInsectNames <- function(y) {
if ("Native Bug" %in% names(y)) {setnames(y, "Native Bug", "Bug")}
if ("Native Moth" %in% names(y)) {setnames(y, "Native Moth", "Moth")}
if ("Native Scale" %in% names(y)) {setnames(y, "Native Scale", "Scale")}
return(y)
}
Add the following code to your website.
For more information on customizing the embed code, read Embedding Snippets.