---
title: RATET01
subtitle: Event Rate Summary for Recurrent Events
---
------------------------------------------------------------------------
{{< include ../../_utils/envir_hook.qmd >}}
```{r setup, echo = FALSE, warning = FALSE, message = FALSE}
library (tern)
library (dplyr)
anl <- random.cdisc.data:: cadtte %>%
filter (PARAMCD == "TNE" )
# Ensure number of exacerbation is a factor and NAs are explicit missing levels.
anl$ AVAL_f <- as.factor (anl$ AVAL)
anl <- df_explicit_na (anl)
```
```{r include = FALSE}
webr_code_labels <- c ("setup" )
```
{{< include ../../_utils/webr_no_include.qmd >}}
## Output
This example will focus on data from a single visit "WEEK 1 DAY 8" and a single endpoint "FKSI-FWB". Only evaluable populations will be analyzed.
:::: panel-tabset
## Customized Table
::: {.panel-tabset .nav-justified group="webr"}
## {{< fa regular file-lines sm fw >}} Preview
```{r variant1, test = list(result_v1 = "result")}
lyt <- basic_table (show_colcounts = TRUE ) %>%
split_cols_by ("ARM" , ref_group = "B: Placebo" , split_fun = ref_group_position ("first" )) %>%
analyze_vars (
"AVAL_f" ,
var_labels = "Number of exacerbations per patient" ,
.stats = c ("count_fraction" ),
.formats = c ("count_fraction" = "xx (xx.xx%)" ),
.labels = c ("Number of exacerbations per patient" )
) %>%
summarize_glm_count (
vars = "AVAL" ,
variables = list (arm = "ARM" , offset = "lgTMATRSK" , covariates = NULL ),
conf_level = 0.95 ,
distribution = "poisson" ,
rate_mean_method = "emmeans" ,
var_labels = "Unadjusted exacerbation rate (per year)" ,
table_names = "unadj" ,
.stats = c ("rate" ),
.labels = c (rate = "Rate" )
) %>%
summarize_glm_count (
vars = "AVAL" ,
variables = list (arm = "ARM" , offset = "lgTMATRSK" , covariates = c ("REGION1" )),
conf_level = 0.95 ,
distribution = "quasipoisson" ,
rate_mean_method = "ppmeans" ,
var_labels = "Adjusted (QP) exacerbation rate (per year)" ,
table_names = "adj-qp" ,
.stats = c ("rate" , "rate_ci" , "rate_ratio" , "rate_ratio_ci" , "pval" ),
.labels = c (
rate = "Rate" , rate_ci = "Rate CI" , rate_ratio = "Rate Ratio" ,
rate_ratio_ci = "Rate Ratio CI" , pval = "p-value"
)
) %>%
summarize_glm_count (
vars = "AVAL" ,
variables = list (arm = "ARM" , offset = "lgTMATRSK" , covariates = c ("REGION1" )),
conf_level = 0.95 ,
distribution = "negbin" ,
rate_mean_method = "emmeans" ,
var_labels = "Adjusted (NB) exacerbation rate (per year)" ,
table_names = "adj-nb" ,
.stats = c ("rate" , "rate_ci" , "rate_ratio" , "rate_ratio_ci" , "pval" ),
.labels = c (
rate = "Rate" , rate_ci = "Rate CI" , rate_ratio = "Rate Ratio" ,
rate_ratio_ci = "Rate Ratio CI" , pval = "p-value"
)
)
result <- build_table (
lyt = lyt,
df = anl
)
result
```
```{r include = FALSE}
webr_code_labels <- c ("variant1" )
```
{{< include ../../_utils/webr.qmd >}}
:::
## Data Setup
```{r setup}
#| code-fold: show
```
::::
{{< include ../../_utils/save_results.qmd >}}
{{< include ../../repro.qmd >}}