plot.ICE {gfoRmulaICE} | R Documentation |
Plot method for ICE estimator objects
Description
This function provides visualization of estimated risk for all specified interventions, estimated natural course risk, and observed risk at each time point.
Usage
## S3 method for class 'ICE'
plot(..., plot_obs = TRUE, label = 0)
Arguments
... |
ICE estimator objects. |
plot_obs |
a logical value indicating whether to plot the observed risk over time.
Default is |
label |
a number specifying which time label is used in x-axis. 0 represents using generic numerical time index, and 1 represents using the original time label from the data. Default is 0. |
Value
a plot for risks of all the interventions specified in ...
.
Examples
data <- gfoRmulaICE::compData
ice_fit1 <- ice(
data = data,
time_points = 4,
id = "id",
time_name = "t0",
censor_name = "C",
outcome_name = "Y",
compevent_name = "D",
comp_effect = 0,
outcome_model = Y ~ L1 + L2 + A1 + A2,
censor_model = C ~ L1 + L2 + A1 + A2,
ref_idx = 0,
estimator = pool(hazard = FALSE),
int_descript = "Static Intervention",
intervention1.A1 = list(static(3)),
intervention1.A2 = list(static(1))
)
ice_fit2 <- ice(
data = data,
time_points = 4,
id = "id",
time_name = "t0",
censor_name = "C",
outcome_name = "Y",
compevent_name = "D",
comp_effect = 0,
outcome_model = Y ~ L1 + L2 + A1 + A2,
censor_model = C ~ L1 + L2 + A1 + A2,
competing_model = D ~ L1 + L2 + A1 + A2,
ref_idx = 0,
estimator = pool(hazard = TRUE),
int_descript = "Static Intervention",
intervention1.A1 = list(static(3)),
intervention1.A2 = list(static(1))
)
plot(ice_fit1, ice_fit2)
[Package gfoRmulaICE version 0.1.0 Index]