Performs checks to gather a list of guesses and iterations

Gather_Guesses_CPP called from within R, uses a list of options and the model definition to generate a list of parameters and iterations that do not produce errors

Usage

Gather_Guesses_CPP(
  df,
  dfc,
  names,
  term_n,
  tform,
  keep_constant,
  a_n,
  x_all,
  a_n_default,
  modelform,
  fir,
  control,
  guesses_control,
  model_control = list()
)

Arguments

df: a data.table containing the columns of interest
dfc: vector matching subterm number to matrix column
names: columns for elements of the model, used to identify data columns
term_n: term numbers for each element of the model
tform: list of string function identifiers, used for linear/step
keep_constant: binary values to denote which parameters to change
a_n: list of initial parameter values, used to determine number of parameters. May be either a list of vectors or a single vector.
x_all: covariate matrix
a_n_default: center of parameter distribution guessing scope
modelform: string specifying the model type: M, ME, A, PA, PAE, GMIX, GMIX-R, GMIX-E
fir: term number for the initial term, used for models of the form T0*f(Ti) in which the order matters
control: list of parameters controlling the convergence, see Def_Control() for options or vignette("Control_Options")
guesses_control: list of parameters to control how the guessing works, see Def_Control_Guess() for options or vignette("Control_Options")
model_control: controls which alternative model options are used, see Def_model_control() for options and vignette("Control_Options") for further details

Value

returns a list of the final results

Examples

library(data.table)
a <- c(0, 1, 2, 3, 4, 5, 6)
b <- c(1, 2, 3, 4, 5, 6, 7)
c <- c(0, 1, 0, 0, 0, 1, 0)
d <- c(3, 4, 5, 6, 7, 8, 9)
df <- data.table::data.table("a" = a, "b" = b, "c" = c, "d" = d)
time1 <- "a"
time2 <- "b"
event <- "c"
names <- c("d")
term_n <- c(0)
tform <- c("loglin")
keep_constant <- c(0)
a_n <- c(-0.1)
a_n_default <- a_n
modelform <- "M"
fir <- 0
der_iden <- 0
control <- list(
  "ncores" = 2, "lr" = 0.75, "maxiter" = -1,
  "halfmax" = 5, "epsilon" = 1e-9,
  "deriv_epsilon" = 1e-9, "abs_max" = 1.0, "change_all" = TRUE,
  "dose_abs_max" = 100.0, "verbose" = FALSE, "ties" = "breslow",
  "double_step" = 1
)
guesses_control <- list()
model_control <- list()
all_names <- unique(names(df))
dfc <- match(names, all_names)
term_tot <- max(term_n) + 1
x_all <- as.matrix(df[, all_names, with = FALSE])
control <- Def_Control(control)
guesses_control <- Def_Control_Guess(guesses_control, a_n)
model_control <- Def_model_control(model_control)
options(warn = -1)
Gather_Guesses_CPP(
  df, dfc, names, term_n, tform, keep_constant,
  a_n, x_all, a_n_default,
  modelform, fir, control, guesses_control
)
#> $a_ns
#>  [1] -0.100000000  0.613150009  0.753092897  0.004755189  0.054131389
#>  [6]  0.274895491 -0.610563541 -1.010765683  0.719709113 -0.958637562
#> 
#> $maxiters
#>  [1] 5 5 5 5 5 5 5 5 5 5
#>