#' Create table of survival model results
#'
#' @param convergence.status
#' @param parameter.estimates
#' @param fit.statistics
#' @param replace_list
#' @param select_list
#'
#' @export
survival_model_results_function_stricta <- function(
convergence.status,
parameter.estimates,
fit.statistics
)
{
# filter out models that didn't converge
y = merge(convergence.status, parameter.estimates) %>%
filter(Reason=="Algorithm converged.") %>%
# change effects to columns
dcast(modelVars~Parameter, value.var="Estimate", fun=sum)
# fit statistics
z = fit.statistics %>%
filter(Criterion=="AIC (smaller is better)" | Criterion=="Scaled Deviance") %>%
dcast(modelVars~Criterion, value.var="Value") %>%
setnames("AIC (smaller is better)", "AIC")
y %<>% merge(z, by="modelVars") %>%
arrange(modelVars)
y %<>%
mutate(
Intercept = paste(
Intercept %>% round(2),
" [",
parameter.estimates[which(parameter.estimates$Parameter=="Intercept"), ]$LowerLRCL[1] %>% round(2),
", ",
parameter.estimates[which(parameter.estimates$Parameter=="Intercept"), ]$UpperLRCL[1] %>% round(2),
"]",
sep=""
),
`C_t` = paste(
Ln_Size_t_1_st %>% round(2),
" [",
parameter.estimates[which(parameter.estimates$Parameter=="Ln_Size_t_1_st"), ]$LowerLRCL[1] %>% round(2),
", ",
parameter.estimates[which(parameter.estimates$Parameter=="Ln_Size_t_1_st"), ]$UpperLRCL[1] %>% round(2),
"]",
sep=""
)
) %>%
dplyr::select(
Intercept,
C_t,
`Scaled Deviance`
)
return(y)
}
Add the following code to your website.
For more information on customizing the embed code, read Embedding Snippets.