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 TRUE.

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]