From 99cb19a9d941a83bec7242494d2e1dece7b4f4fd Mon Sep 17 00:00:00 2001
From: Paul-Corbalan <58653590+Paul-Corbalan@users.noreply.github.com>
Date: Tue, 12 Apr 2022 10:25:21 +0200
Subject: [PATCH] Add Elisa code

---
 Comparaison_of_methods.rmd | 16 +++++++++++++---
 1 file changed, 13 insertions(+), 3 deletions(-)

diff --git a/Comparaison_of_methods.rmd b/Comparaison_of_methods.rmd
index fe903e5..6091118 100644
--- a/Comparaison_of_methods.rmd
+++ b/Comparaison_of_methods.rmd
@@ -295,9 +295,19 @@ LocaScoreMC <- function(lambda0, lambda1, NbSeq, tbe0, T){
   pvalue = c()
   X = c()
   
-  Score = ScoreDistrib(lambda0, lambda1, NbSeq, T)
-  xp = Score$X
-  P_X = Score$P_X
+score_max=floor(E*log(lambda1/lambda0))
+
+## Calcul score_min 
+score_min_c=floor(E*log(lambda1/lambda0)+E*(lambda0-lambda1)*T)
+
+  l=seq(score_min_c,score_max,1)
+  borne_inf=(l-E*log(lambda1/lambda0))/(E*(lambda0-lambda1))
+  borne_sup=(l+1-E*log(lambda1/lambda0))/(E*(lambda0-lambda1))
+  proba.l=pexp(rate=lambda0,borne_inf)-pexp(rate=lambda0,borne_sup)
+  S=sum(proba.l)
+  new.proba.s=proba.l/S
+  xp = l
+  P_X = proba.l/S
   
   min_X = min(xp)
   max_X = max(xp)