Check if the sum of transition probabilities equal to one. — check_sum_of_transition

check_sum_of_transition_array checks if each of the rows of the transition matrices sum to one.

check_sum_of_transition_array(
  a_P,
  n_states = dim(a_P)[1],
  n_t = dim(a_P)[3],
  err_stop = FALSE,
  verbose = FALSE,
  tolerance = 1e-08
)

Arguments

a_P: A transition probability array.
n_states: Number of health states. Defaults to the first dimension of a_P.
n_t: Number of cycles. Defaults to the third dimension of a_P.
err_stop: Logical variable to stop model run if set up as TRUE. Default = FALSE.
verbose: Logical variable to indicate print out of messages. Default = FALSE
tolerance: Numeric tolerance for the comparison of each row sum with one. Default = 1e-08.

Value

Invisibly, a data frame with one row per row of a_P whose transition probabilities do not sum to one, empty when the array is valid. With err_stop = TRUE the function stops instead, and with verbose = TRUE it warns, listing the offending rows and cycles.

Examples

l_params_all <- load_all_params()
a_P <- decision_model(l_params_all)$a_P
check_sum_of_transition_array(a_P, verbose = TRUE)