Nothing
tests/testthat/_snaps/summarize_glm_count.md
h_glm_poisson glm-fit works with healthy input
Code
res
Output
Estimate SE z_value Pr coefs
1 2.10487843 0.07035975 29.9159442 1.220818e-196 (Intercept)
2 0.10946932 0.09613709 1.1386793 2.548369e-01 ARMB: Placebo
3 -0.04371457 0.10265030 -0.4258591 6.702105e-01 ARMC: Combination
h_glm_poisson emmeans-fit works with healthy input no offset
Code
res
Output
ARMCD rate std.error df null statistic p.value
1 ARM A 2.927536 0.2059807 Inf 1 15.26670 1.274641e-52
2 ARM B 3.191781 0.2091005 Inf 1 17.71547 3.185540e-70
3 ARM C 3.086207 0.2306739 Inf 1 15.07747 2.278432e-51
h_glm_poisson glm-fit works with healthy input with covariates
Code
res
Output
Estimate SE z_value Pr coefs
1 2.01065582 0.18541942 10.8438255 2.133586e-27 (Intercept)
2 0.07631174 0.17896220 0.4264126 6.698072e-01 REGION1Asia
3 0.64425750 0.22389462 2.8775033 4.008358e-03 REGION1Eurasia
4 2.13096720 0.36521976 5.8347533 5.387022e-09 REGION1Europe
5 -0.07449500 0.20314837 -0.3667024 7.138410e-01 REGION1North America
6 0.38101695 0.21554753 1.7676703 7.711605e-02 REGION1South America
7 0.11047866 0.09872549 1.1190490 2.631192e-01 ARMB: Placebo
8 -0.17694419 0.10873176 -1.6273459 1.036637e-01 ARMC: Combination
h_glm_poisson emmeans-fit works with healthy input with covariates
Code
res
Output
ARMCD rate std.error df null statistic p.value
1 ARM A 12.64167 1.2378669 Inf 1 25.90902 5.270655e-148
2 ARM B 14.11838 1.2848735 Inf 1 29.09088 4.682722e-186
3 ARM C 10.59153 0.9708089 Inf 1 25.74821 3.375733e-146
h_glm_quasipoisson glm-fit works with healthy input
Code
res
Output
Estimate SE z_value Pr coefs
1 2.01065582 0.7781805 2.5837912 0.01051488 (Intercept)
2 0.07631174 0.7510804 0.1016026 0.91917812 REGION1Asia
3 0.64425750 0.9396557 0.6856314 0.49377257 REGION1Eurasia
4 2.13096720 1.5327784 1.3902643 0.16605816 REGION1Europe
5 -0.07449500 0.8525865 -0.0873753 0.93046426 REGION1North America
6 0.38101695 0.9046241 0.4211881 0.67408881 REGION1South America
7 0.11047866 0.4143377 0.2666392 0.79003312 ARMB: Placebo
8 -0.17694419 0.4563327 -0.3877526 0.69862863 ARMC: Combination
h_glm_quasipoisson emmeans-fit works with healthy input
Code
res
Output
ARM rate std.error df null statistic p.value
1 A: Drug X 12.64167 5.195162 Inf 1 6.173420 6.682841e-10
2 B: Placebo 14.11838 5.392442 Inf 1 6.931571 4.161914e-12
3 C: Combination 10.59153 4.074355 Inf 1 6.135104 8.510352e-10
h_glm_negbin glm-fit works with healthy input with offset
Code
res
Output
Estimate SE z_value Pr coefs
1 2.56211826 0.3969418 6.4546439 1.084737e-10 (Intercept)
2 0.12359121 0.2085554 0.5926060 5.534448e-01 ARMB: Placebo
3 0.05569230 0.2247874 0.2477554 8.043237e-01 ARMC: Combination
4 0.05684153 0.3845557 0.1478109 8.824920e-01 REGION1Asia
5 1.84806287 0.4849801 3.8105953 1.386326e-04 REGION1Eurasia
6 2.04520468 0.9029032 2.2651428 2.350392e-02 REGION1Europe
7 0.44720215 0.4417735 1.0122884 3.114002e-01 REGION1North America
8 1.21153294 0.4690914 2.5827223 9.802419e-03 REGION1South America
h_glm_negbin emmeans-fit works with healthy input no offset
Code
res
Output
ARM response std.error df null statistic p.value
1 A: Drug X 3.322579 0.3367532 Inf 1 11.84712 2.227054e-32
2 B: Placebo 3.720373 0.3782682 Inf 1 12.92183 3.390023e-38
3 C: Combination 3.412887 0.3424577 Inf 1 12.23369 2.054037e-34
h_glm_count glm-fit works with healthy input
Code
res
Output
Estimate SE z_value Pr coefs
1 2.10487843 0.07035975 29.9159442 1.220818e-196 (Intercept)
2 0.10946932 0.09613709 1.1386793 2.548369e-01 ARMCDARM B
3 -0.04371457 0.10265030 -0.4258591 6.702105e-01 ARMCDARM C
h_glm_count emmeans-fit works with healthy input
Code
res
Output
ARMCD rate std.error df null statistic p.value
1 ARM A 8.206105 0.5773795 Inf 1 29.91594 1.220818e-196
2 ARM B 9.155436 0.5997925 Inf 1 33.80055 1.935734e-250
3 ARM C 7.855107 0.5871181 Inf 1 27.57650 2.129731e-167
h_ppmeans works with healthy input
Code
fits2
Output
$glm_fit
Call: stats::glm(formula = formula, family = stats::quasipoisson(link = "log"),
data = .df_row, offset = offset)
Coefficients:
(Intercept) REGION1Asia REGION1Eurasia
2.01066 0.07631 0.64426
REGION1Europe REGION1North America REGION1South America
2.13097 -0.07450 0.38102
ARMCDARM B ARMCDARM C
0.11048 -0.17694
Degrees of Freedom: 199 Total (i.e. Null); 192 Residual
Null Deviance: 983.8
Residual Deviance: 939 AIC: NA
$emmeans_fit
ARMCD rate SE df asymp.LCL asymp.UCL
ARM A 12.6 5.20 Inf 5.65 28.3
ARM B 14.1 5.39 Inf 6.68 29.8
ARM C 10.6 4.07 Inf 4.98 22.5
Results are averaged over the levels of: REGION1
Confidence level used: 0.95
Intervals are back-transformed from the log scale
s_glm_count works with healthy input
Code
res
Output
$n
[1] 73
$rate
[1] 14.11838
attr(,"label")
[1] "Adjusted Rate"
$rate_ci
[1] 11.81189 16.87525
attr(,"label")
[1] "95% CI"
$rate_ratio
character(0)
attr(,"label")
[1] "Adjusted Rate Ratio"
$rate_ratio_ci
character(0)
attr(,"label")
[1] "95% CI"
$pval
character(0)
attr(,"label")
[1] "p-value"
s_glm_count (negative binomial) works with healthy input
Code
res
Output
$n
[1] 73
$rate
[1] 37.35687
attr(,"label")
[1] "Adjusted Rate"
$rate_ci
[1] 24.59065 56.75067
attr(,"label")
[1] "95% CI"
$rate_ratio
character(0)
attr(,"label")
[1] "Adjusted Rate Ratio"
$rate_ratio_ci
character(0)
attr(,"label")
[1] "95% CI"
$pval
character(0)
attr(,"label")
[1] "p-value"
s_glm_count works with no reference group selected.
Code
res
Output
$n
[1] 73
$rate
[1] 14.11838
attr(,"label")
[1] "Adjusted Rate"
$rate_ci
[1] 11.81189 16.87525
attr(,"label")
[1] "95% CI"
$rate_ratio
[1] 0.8954054 0.7501944
attr(,"label")
[1] "Adjusted Rate Ratio"
$rate_ratio_ci
[1] 0.7378778 0.6062152 1.0865633 0.9283695
attr(,"label")
[1] "95% CI"
$pval
[1] 0.263119218 0.008203621
attr(,"label")
[1] "p-value"
s_glm_count (negative binomial) works with no reference group selected.
Code
res
Output
$n
[1] 73
$rate
[1] 37.35687
attr(,"label")
[1] "Adjusted Rate"
$rate_ci
[1] 24.59065 56.75067
attr(,"label")
[1] "95% CI"
$rate_ratio
[1] 0.8837410 0.9343549
attr(,"label")
[1] "Adjusted Rate Ratio"
$rate_ratio_ci
[1] 0.5872220 0.6062612 1.3299880 1.4400049
attr(,"label")
[1] "95% CI"
$pval
[1] 0.5534448 0.7583367
attr(,"label")
[1] "p-value"
summarize_glm_count works with healthy inputs
Code
res
Output
B: Placebo A: Drug X C: Combination
(N=73) (N=69) (N=58)
————————————————————————————————————————————————————————————————————————————————————
Number of exacerbations per patient
0 8 (10.96%) 3 (4.35%) 6 (10.34%)
1 9 (12.33%) 11 (15.94%) 6 (10.34%)
2 15 (20.55%) 18 (26.09%) 9 (15.52%)
3 11 (15.07%) 14 (20.29%) 15 (25.86%)
4 9 (12.33%) 10 (14.49%) 9 (15.52%)
5 9 (12.33%) 7 (10.14%) 8 (13.79%)
6 4 (5.48%) 4 (5.80%) 4 (6.90%)
7 8 (10.96%) 2 (2.90%) 0 (0.00%)
10 0 (0.00%) 0 (0.00%) 1 (1.72%)
Unadjusted exacerbation rate (per year)
Rate 9.1554 8.2061 7.8551
summarize_glm_count (negative binomial) works with healthy inputs
Code
res
Output
B: Placebo A: Drug X C: Combination
(N=73) (N=69) (N=58)
————————————————————————————————————————————————————————————————————————————————————
Number of exacerbations per patient
0 8 (10.96%) 3 (4.35%) 6 (10.34%)
1 9 (12.33%) 11 (15.94%) 6 (10.34%)
2 15 (20.55%) 18 (26.09%) 9 (15.52%)
3 11 (15.07%) 14 (20.29%) 15 (25.86%)
4 9 (12.33%) 10 (14.49%) 9 (15.52%)
5 9 (12.33%) 7 (10.14%) 8 (13.79%)
6 4 (5.48%) 4 (5.80%) 4 (6.90%)
7 8 (10.96%) 2 (2.90%) 0 (0.00%)
10 0 (0.00%) 0 (0.00%) 1 (1.72%)
Unadjusted exacerbation rate (per year)
Rate 18.6539 21.8760 22.6006
Try the tern package in your browser
Any scripts or data that you put into this service are public.
