View source: R/multibias_plot.R
| multibias_plot | R Documentation |
This function generates a forest plot comparing the observed effect estimate with adjusted effect estimates from sensitivity analyses. The plot includes point estimates and confidence intervals for each analysis.
multibias_plot(data_observed, multibias_result_list, log_scale = FALSE)
data_observed |
Object of class |
multibias_result_list |
A named list of sensitivity analysis results.
Each element should be a result from |
log_scale |
Boolean indicating whether to display the x-axis on the log scale. Default is FALSE. |
A ggplot object showing a forest plot with:
Point estimates (blue dots)
Confidence intervals (gray horizontal lines)
A vertical reference line at x=1 (dashed)
Appropriate labels and title
## Not run:
df_observed <- data_observed(
data = df_em,
bias = "em",
exposure = "Xstar",
outcome = "Y",
confounders = "C1"
)
bp1 <- bias_params(coef_list = list(x = c(-2.10, 1.62, 0.63, 0.35)))
bp2 <- bias_params(coef_list = list(x = c(-2.10 * 2, 1.62 * 2, 0.63 * 2, 0.35 * 2)))
result1 <- multibias_adjust(
data_observed = df_observed,
bias_params = bp1
)
result2 <- multibias_adjust(
data_observed = df_observed,
bias_params = bp2
)
multibias_plot(
data_observed = df_observed,
multibias_result_list = list(
"Adjusted with bias params" = result1,
"Adjusted with bias params doubled" = result2
)
)
## End(Not run)
Add the following code to your website.
For more information on customizing the embed code, read Embedding Snippets.