Browse code

prevented NA F statistics from being considered for FDR estimation

Nils Kurzawa authored on 19/03/2020 15:29:04
Showing2 changed files

... ...
@@ -1,6 +1,6 @@
1 1
 Package: TPP2D
2 2
 Title: FDR-controlled analysis of 2D-TPP experiments
3
-Version: 1.3.9
3
+Version: 1.3.10
4 4
 Authors@R: c(
5 5
     person("Nils", "Kurzawa", email = "nils.kurzawa@embl.de", 
6 6
         role = c("aut", "cre")),
... ...
@@ -36,6 +36,7 @@ getFDR <- function(df_out, df_null){
36 36
                         df_null) %>%
37 37
         mutate(nObsRound = round(nObs/10)*10) %>%
38 38
         group_by(nObsRound) %>%
39
+        filter(!is.na(F_statistic)) %>% 
39 40
         arrange(desc(F_statistic)) %>%
40 41
         mutate(max_rank = n(),
41 42
                rank = dense_rank(desc(F_statistic)),