@@ -5,8 +5,6 @@ subtitle: Event Rate Summary for Recurrent Events
5
5
6
6
------------------------------------------------------------------------
7
7
8
- {{< include ../../test-utils/envir_hook.qmd >}}
9
-
10
8
This example will focus on data from a single visit "WEEK 1 DAY 8" and a single endpoint "FKSI-FWB".
11
9
Only evaluable populations will be analyzed.
12
10
@@ -31,65 +29,9 @@ anl <- df_explicit_na(anl)
31
29
## Customized Table
32
30
33
31
``` {r variant1, test = list(result_v1 = "result")}
34
- lyt <- basic_table(show_colcounts = TRUE) %>%
35
- split_cols_by("ARM", ref_group = "B: Placebo", split_fun = ref_group_position("first")) %>%
36
- analyze_vars(
37
- "AVAL_f",
38
- var_labels = "Number of exacerbations per patient",
39
- .stats = c("count_fraction"),
40
- .formats = c("count_fraction" = "xx (xx.xx%)"),
41
- .label = c("Number of exacerbations per patient")
42
- ) %>%
43
- summarize_glm_count(
44
- vars = "AVAL",
45
- variables = list(arm = "ARM", offset = "lgTMATRSK", covariates = NULL),
46
- conf_level = 0.95,
47
- distribution = "poisson",
48
- rate_mean_method = "emmeans",
49
- var_labels = "Unadjusted exacerbation rate (per year)",
50
- table_names = "unadj",
51
- .stats = c("rate"),
52
- .labels = c(rate = "Rate")
53
- ) %>%
54
- summarize_glm_count(
55
- vars = "AVAL",
56
- variables = list(arm = "ARM", offset = "lgTMATRSK", covariates = c("REGION1")),
57
- conf_level = 0.95,
58
- distribution = "quasipoisson",
59
- rate_mean_method = "ppmeans",
60
- var_labels = "Adjusted (QP) exacerbation rate (per year)",
61
- table_names = "adj-qp",
62
- .stats = c("rate", "rate_ci", "rate_ratio", "rate_ratio_ci", "pval"),
63
- .labels = c(
64
- rate = "Rate", rate_ci = "Rate CI", rate_ratio = "Rate Ratio",
65
- rate_ratio_ci = "Rate Ratio CI", pval = "p value"
66
- )
67
- ) %>%
68
- summarize_glm_count(
69
- vars = "AVAL",
70
- variables = list(arm = "ARM", offset = "lgTMATRSK", covariates = c("REGION1")),
71
- conf_level = 0.95,
72
- distribution = "negbin",
73
- rate_mean_method = "emmeans",
74
- var_labels = "Adjusted (NB) exacerbation rate (per year)",
75
- table_names = "adj-nb",
76
- .stats = c("rate", "rate_ci", "rate_ratio", "rate_ratio_ci", "pval"),
77
- .labels = c(
78
- rate = "Rate", rate_ci = "Rate CI", rate_ratio = "Rate Ratio",
79
- rate_ratio_ci = "Rate Ratio CI", pval = "p value"
80
- )
81
- )
82
-
83
- result <- build_table(
84
- lyt = lyt,
85
- df = anl
86
- )
87
-
88
- result
32
+ # Under development
89
33
```
90
34
91
- {{< include ../../test-utils/save_results.qmd >}}
92
-
93
35
## ` teal ` App
94
36
95
37
``` {r teal, message=FALSE, opts.label=c('skip_if_testing', 'app')}
0 commit comments