Description Usage Arguments Value Examples
inter_frailty
returns a list with the estimated parameters par
and their asymptotic covariance matrix mcov
. The list also contains a
dummy variable stop_c
assuming 0 if algorithm converged and 1 if a
stop criteria ended the process.
1 2 3 | inter_frailty(dataset, left, right, delta, cov_theta, cov_beta, M, b = 0.001,
tol = 0.001, max_n = 100, par_cl = NULL, burn_in = 30,
output_files = FALSE)
|
dataset |
Dataset used to fit the model. |
left |
Vector containing the last check times before event. |
right |
Vector containing the first check times after event. |
delta |
Flag vector indicating failure inside interval. |
cov_theta |
String vector containing the column names to be used on the cure rate predictor. |
cov_beta |
String vector containing the column names to be used on the predictor associated with the hazard function. |
M |
Number of replicates generated by each iteration on the ANDA (Asymptotic Normal Data Augmentation) algorithm. |
b |
Parameter for initial theta and beta variances. |
tol |
Numeric for tolerance of convergence. |
max_n |
Maximum number of algorithm's iterations without the burn in. |
par_cl |
Registered SOCK cluster for parallel process. If NULL (default) the program loops are executed sequentially. |
burn_in |
Number of burn in iterations. |
output_files |
Boolean indicating if text outputs for the estimates and variances should be generated. |
The inter_frailty
function returns an list containing the
following outputs:
|
estimates of theta and beta parameters. |
|
estimates for the covariance matrix of theta and beta parameters. |
|
stop criteria indicator assuming 1 when process is stopped for a non-convergence criteria. Assumes 0 when convergence is reached. |
1 2 3 4 5 6 7 8 9 10 11 12 13 14 | ## few iterations just to check how to use the function
set.seed(3)
sample_set <- sim_frailty(80)
inter_frailty(sample_set, sample_set$L, sample_set$R, sample_set$delta,
c("xi1","xi2"), c("xi1","xi2"), M = 10, max_n = 3, burn_in = 0)
## precise estimate (computationally intensive)
## Not run:
inter_frailty(sample_set, sample_set$L, sample_set$R, sample_set$delta,
c("xi1"), c("xi2"), M = 50, max_n = 50, burn_in = 10)
## End(Not run)
|
Add the following code to your website.
For more information on customizing the embed code, read Embedding Snippets.