R/processColumnNames.R

#' 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)
}
ksauby/modresproc documentation built on May 20, 2019, 7:02 p.m.