Beni feedback

data/preprocess-rpcb-data.R

@@ -56,6 +56,11 @@ datClean %>%
               successes = sum(resulto == "Null" &
                               resultr %in% c("Null-positive", "Null-negative",
                                              "Null")))
+## this should give the same
+datClean %>%
+    summarise(nulls = sum(resulto == "Null"),
+              successes = sum(resulto == "Null" &
+                              pr > 0.05))
 
 ## should give 112 original positive effects and 44 successful replications
 ## (see positive results in Table 1)
@@ -64,6 +69,12 @@ datClean %>%
               successes = sum(resulto == "Positive" &
                               resultr == "Positive" &
                               sign(smdo) == sign(smdr)))
+## this should give the same
+datClean %>%
+    summarise(positives = sum(resulto == "Positive"),
+              successes = sum(resulto == "Positive" &
+                              pr < 0.05 &
+                              sign(smdo) == sign(smdr)))
 
 ## save
 write.csv(datClean, "rpcb-outcome-level.csv", row.names = FALSE)
@@ -114,6 +125,11 @@ datClean2 %>%
               successes = sum(resulto == "Null" &
                               resultr %in% c("Null-positive", "Null-negative",
                                              "Null")))
+## this should give the same
+datClean2 %>%
+    summarise(nulls = sum(resulto == "Null"),
+              successes = sum(resulto == "Null" &
+                              pr > 0.05))
 
 ## should give 97 original positive effects and 42 successful replications
 ## (see positive results in Table 1)
@@ -122,6 +138,12 @@ datClean2 %>%
               successes = sum(resulto == "Positive" &
                               resultr == "Positive" &
                               sign(smdo) == sign(smdr)))
+## this should give the same
+datClean2 %>%
+    summarise(positives = sum(resulto == "Positive"),
+              successes = sum(resulto == "Positive" &
+                              pr < 0.05 &
+                              sign(smdo) == sign(smdr)))
 
 ## replicate Figure 1
 ## https://iiif.elifesciences.org/lax/71601%2Felife-71601-fig2-v3.tif/full/1500,/0/default.jpg