From 922f42cf76b8b123cfc5b35c229b2c7eb9765ef9 Mon Sep 17 00:00:00 2001
From: corentinchoisy <corentin.choisy@proton.me>
Date: Thu, 16 May 2024 11:16:45 +0200
Subject: [PATCH] Added lrt proportion to ROSALI aggregation

---
 RProject/Scripts/aggregation.R | 4 +++-
 1 file changed, 3 insertions(+), 1 deletion(-)

diff --git a/RProject/Scripts/aggregation.R b/RProject/Scripts/aggregation.R
index 4bb16b2..65184ae 100644
--- a/RProject/Scripts/aggregation.R
+++ b/RProject/Scripts/aggregation.R
@@ -621,6 +621,7 @@ compile_simulation2_rosali <- function(scenario) {
     thay.fpr <- mean(sapply(1:1000,function(x) ifelse(sum(!is.na(s[x,paste0("dif_detect_",1:7)]))!=0,sum(!(unlist(s[x,paste0("dif_detect_",1:7)])%in%c(s[x,"real_dif_1"],s[x,"real_dif_2"],s[x,"real_dif_3"])),na.rm=F)/7,0)  ))
     
   }
+  lrt.pass <- mean(s$lrt_passed)
   z <- data.frame(m.beta=mean(s$beta),
                   se.empirical.beta=sd(s$beta),
                   se.analytical.beta=mean(s$se_beta),
@@ -636,7 +637,8 @@ compile_simulation2_rosali <- function(scenario) {
                   moreflexible.detect=moreflexible.detect,
                   any.detect=any.detect,
                   thay.tpr=thay.tpr,
-                  thay.fpr=thay.fpr
+                  thay.fpr=thay.fpr,
+                  lrt.pass=lrt.pass
   )
   d <- cbind(b,a,z)
   d$prop.