Interpret_Output uses the list output from a regression, prints off a table of results and summarizes the score and convergence.
See also
Other Output and Information Functions:
System_Version()
Examples
library(data.table)
## basic example code reproduced from the starting-description vignette
df <- data.table::data.table(
"UserID" = c(112, 114, 213, 214, 115, 116, 117),
"Starting_Age" = c(18, 20, 18, 19, 21, 20, 18),
"Ending_Age" = c(30, 45, 57, 47, 36, 60, 55),
"Cancer_Status" = c(0, 0, 1, 0, 1, 0, 0),
"a" = c(0, 1, 1, 0, 1, 0, 1),
"b" = c(1, 1.1, 2.1, 2, 0.1, 1, 0.2),
"c" = c(10, 11, 10, 11, 12, 9, 11),
"d" = c(0, 0, 0, 1, 1, 1, 1),
"e" = c(0, 0, 1, 0, 0, 0, 1)
)
# For the interval case
time1 <- "Starting_Age"
time2 <- "Ending_Age"
event <- "Cancer_Status"
names <- c("a", "b", "c", "d")
a_n <- list(c(1.1, -0.1, 0.2, 0.5), c(1.6, -0.12, 0.3, 0.4))
# used to test at a specific point
term_n <- c(0, 1, 1, 2)
tform <- c("loglin", "lin", "lin", "plin")
modelform <- "M"
keep_constant <- c(0, 0, 0, 0)
control <- list(
"ncores" = 2, "lr" = 0.75, "maxiters" = c(5, 5, 5),
"halfmax" = 5, "epsilon" = 1e-3, "deriv_epsilon" = 1e-3,
"abs_max" = 1.0, "dose_abs_max" = 100.0,
"verbose" = FALSE,
"ties" = "breslow", "double_step" = 1, "guesses" = 2
)
e <- RunCoxRegression_Omnibus(df, time1, time2, event,
names, term_n, tform, keep_constant,
a_n, modelform, control,
model_control = list(
"single" = FALSE,
"basic" = FALSE, "cr" = FALSE, "null" = FALSE
)
)
Interpret_Output(e)
#> |-------------------------------------------------------------------|
#> Final Results
#> Covariate Subterm Term Number Central Estimate Standard Error 2-tail p-value
#> <char> <char> <int> <num> <num> <num>
#> 1: a loglin 0 7.6000000 21.75009 0.7267711
#> 2: b lin 1 0.1541029 333.48693 0.9996313
#> 3: c lin 1 3.3758662 3963.75979 0.9993205
#> 4: d plin 2 6.4000000 57.28494 0.9110437
#>
#> Cox Model Used
#> -2*Log-Likelihood: 1.7, AIC: 9.7
#> Iterations run: 5
#> maximum step size: 1.00e+00, maximum first derivative: 2.12e-02
#> Analysis did not converge, check convergence criteria or run further
#> Run finished in 0.05 seconds
#> |-------------------------------------------------------------------|