add SampleSizes output to pairwise_test_bin

Author: lemireg

R/pairwise_comparisons.R

@@ -610,6 +610,7 @@ pairwise_test_bin <- function(x,
                           rep(j_group, nrow(data_here)))
         groups_here <- droplevels(factor(groups_here, levels = levels_here))
       } else {
+
         i_vals <- x[group == i_group]
         j_vals <- x[group == j_group]
         vals_here <-  x[group  %in% c(i_group, j_group)]
@@ -640,6 +641,8 @@ pairwise_test_bin <- function(x,
 
       stats_by_group <- data.frame(Group1 = i_group,
                                    Group2 = j_group,
+                                   Group1_n = sum(!is.na(i_vals)),
+                                   Group2_n = sum(!is.na(j_vals)),
                                    Group1_rr = response_info_here_by_group[[1]],
                                    Group2_rr = response_info_here_by_group[[2]],
                                    stringsAsFactors = FALSE
@@ -688,8 +691,9 @@ pairwise_test_bin <- function(x,
 
 
   data.frame(Comparison = pasted_results$Comparison,
+             SampleSizes = pasted_results$n_comparison,
              ResponseStats = pasted_results$rr_comparison,
-             ResponseTest = results$ResponseTest ,
+             ResponseTest = results$ResponseTest,
              PerfectSeparation = results$PerfectSeparation,
              stringsAsFactors = FALSE)
 

man/get_full_name.Rd

@@ -597,19 +597,20 @@ test_that("pairwise_comparisons_bin testing two groups", {
     group_by(group) %>%
     mutate(num_pos = sum(x), n = n()) %>%
     group_by(group,num_pos, n) %>%
-    group_modify( ~ wilson_ci(.$x, .95)) %>% ungroup() %>%
+    group_modify( ~ wilson_ci(.$x, .95),
+                  sum(.$x)) %>% ungroup() %>%
     mutate(rr = paste0(num_pos, '/', n, ' = ',
                        stat_paste(mean * 100, lower * 100, upper * 100, digits = 1, suffix = '%')))
 
   testing_stats <- bind_cols(
-    testing_stats_pre %>% filter(group == 'a') %>% select(Group1 = group, Group1_rr = rr),
-    testing_stats_pre %>% filter(group == 'b') %>% select(Group2 = group, Group2_rr = rr))
+    testing_stats_pre %>% filter(group == 'a') %>% select(Group1 = group, Group1_n = n, Group1_rr = rr),
+    testing_stats_pre %>% filter(group == 'b') %>% select(Group2 = group, Group2_n = n, Group2_rr = rr))
 
 
   # testing multiple methods
   test_pasting <- paste_tbl_grp(data = testing_stats)
 
-  names(test_pasting) <- c('Comparison', 'ResponseStats')
+  names(test_pasting) <- c('Comparison', 'SampleSizes', 'ResponseStats')
 
   lapply(c("barnard", "fisher", "chi.sq"),
               function(method_in){
@@ -629,18 +630,20 @@ test_that("pairwise_comparisons_bin testing two groups", {
     group_by(group) %>%
     mutate(num_pos = sum(x), n = n()) %>%
     group_by(group,num_pos, n) %>%
-    group_modify( ~ wilson_ci(.$x, .95)) %>% ungroup() %>%
+    group_modify( ~ wilson_ci(.$x, .95),
+                  sum(.$x)) %>%
+    ungroup() %>%
     mutate(rr = paste0(num_pos, '/', n, ' = ',
                        stat_paste(mean * 100, lower * 100, upper * 100, digits = 3, suffix = '%')))
 
   testing_stats_3digits <- bind_cols(
     testing_stats_pre_3digits %>%
     filter(group == 'a') %>%
-    select(Group1 = group, Group1_rr = rr), testing_stats_pre_3digits %>%
-    filter(group == 'b') %>% select(Group2 = group, Group2_rr = rr))
+    select(Group1 = group, Group1_n = n, Group1_rr = rr), testing_stats_pre_3digits %>%
+    filter(group == 'b') %>% select(Group2 = group, Group2_n = n, Group2_rr = rr))
 
   test_pasting <- paste_tbl_grp(data = testing_stats_3digits)
-  names(test_pasting) <- c('Comparison', 'ResponseStats')
+  names(test_pasting) <- c('Comparison', 'SampleSizes', 'ResponseStats')
   testing_results <- data.frame(test_pasting,
     ResponseTest  = two_samp_bin_test(x = x, y = group),
     PerfectSeparation = ifelse(diff(testing_stats_pre_3digits$mean) == 1,
@@ -652,7 +655,7 @@ test_that("pairwise_comparisons_bin testing two groups", {
   # One-sided test comparison
   test_pasting <- paste_tbl_grp(data = testing_stats, alternative = 'less')
 
-  names(test_pasting) <- c('Comparison', 'ResponseStats')
+  names(test_pasting) <- c('Comparison', 'SampleSizes', 'ResponseStats')
 
   testing_results <- data.frame(test_pasting,
     # Need reverse testing direction
@@ -667,7 +670,7 @@ test_that("pairwise_comparisons_bin testing two groups", {
   # sorted group greater than comparison
   test_pasting <- paste_tbl_grp(data = testing_stats, alternative = 'greater',
                                 first_name = 'Group2', second_name = 'Group1')
-  names(test_pasting) <- c('Comparison', 'ResponseStats')
+  names(test_pasting) <- c('Comparison', 'SampleSizes', 'ResponseStats')
   testing_results <- data.frame(
     test_pasting,
     # Need reverse testing direction
@@ -682,7 +685,7 @@ test_that("pairwise_comparisons_bin testing two groups", {
 
   # High number needed for testing
   test_pasting <- paste_tbl_grp(data = testing_stats)
-  names(test_pasting) <- c('Comparison', 'ResponseStats')
+  names(test_pasting) <- c('Comparison', 'SampleSizes', 'ResponseStats')
   testing_results <- data.frame(
     test_pasting,
     ResponseTest  = NA_integer_,
@@ -701,7 +704,8 @@ test_that("pairwise_comparisons_bin testing two groups", {
     group_by(group) %>%
     mutate(num_pos = sum(x), n = n()) %>%
     group_by(group,num_pos, n) %>%
-    group_modify( ~ wilson_ci(.$x, .95)) %>% ungroup() %>%
+    group_modify( ~ wilson_ci(.$x, .95),
+                  sum(.$x)) %>% ungroup() %>%
     mutate(rr = paste0(num_pos, '/', n, ' = ', stat_paste(mean * 100,
                                                           lower * 100,
                                                           upper * 100,
@@ -710,9 +714,10 @@ test_that("pairwise_comparisons_bin testing two groups", {
 
   paired_stats <- bind_cols(paired_stats_pre %>%
                               filter(group == 'a') %>%
-                              select(Group1 = group, Group1_rr = rr), paired_stats_pre %>%
+                              select(Group1 = group, Group1_n = n, Group1_rr = rr),
+                            paired_stats_pre %>%
                               filter(group == 'b') %>%
-                              select(Group2 = group, Group2_rr = rr))
+                              select(Group2 = group, Group2_n = n, Group2_rr = rr))
 
 
   test_pasting <- paste_tbl_grp(data = paired_stats)
@@ -722,8 +727,8 @@ test_that("pairwise_comparisons_bin testing two groups", {
                                     second_name = "Group1")
   expect_false(identical(test_pasting, test_pasting_rev))
 
-  names(test_pasting) <- c('Comparison', 'ResponseStats')
-  names(test_pasting_rev) <- c('Comparison', 'ResponseStats')
+  names(test_pasting) <- c('Comparison', 'SampleSizes', 'ResponseStats')
+  names(test_pasting_rev) <- c('Comparison', 'SampleSizes', 'ResponseStats')
 
   testing_results <- data.frame(test_pasting,
     ResponseTest  = two_samp_bin_test(x = x, y = group, method = 'mcnemar'),
@@ -759,11 +764,12 @@ test_that("pairwise_test_bin testing 3+ groups", {
               ci = wilson_ci(response),
               r1 = sum(response),
               r0 = abs(sum(response - 1)),
+              n = n(),
               .groups = "keep")  %>%
     pivot_wider(id_cols = c(antigen, visitno),
                 names_from = group,
                 names_prefix = "grp",
-                values_from = c(rfraction, ci, r0, r1)) %>%
+                values_from = c(rfraction, ci, r0, r1, n)) %>%
     mutate(pval = ifelse(((r0_grp1 == 0 & r0_grp2 == 0) | (r1_grp1 == 0 & r1_grp2 == 0)),
                               1,
                               as.double(Exact::exact.test(matrix(c(r1_grp1,
@@ -776,6 +782,7 @@ test_that("pairwise_test_bin testing 3+ groups", {
                                               alternative = "two.sided")$p.value)),
            PerfectSeparation = FALSE,
            Comparison = "1 vs. 2",
+           SampleSizes = paste0(n_grp1, " vs. ", n_grp2),
            ResponseStats = paste0(rfraction_grp1,
                                   " = ",
                                   round_away_0(ci_grp1$mean*100, 1, trailing_zeros = TRUE),
@@ -792,7 +799,7 @@ test_that("pairwise_test_bin testing 3+ groups", {
                                   "%, ",
                                   round_away_0(ci_grp2$upper*100, 1, trailing_zeros = TRUE),
                                   "%)"))  %>%
-    select(antigen, visitno, Comparison, ResponseStats, ResponseTest = pval, PerfectSeparation)
+    select(antigen, visitno, Comparison, SampleSizes, ResponseStats, ResponseTest = pval, PerfectSeparation)
 
   function_obj <- exampleData_BAMA %>%
     group_by(antigen, visitno) %>%