From de2ada97d5bac9950efce07b62627460483e4bd0 Mon Sep 17 00:00:00 2001
From: elisaduz <elisa.duz1807@gmail.com>
Date: Tue, 19 Apr 2022 08:09:37 +0200
Subject: [PATCH 1/2] Update Dataset_study.rmd

---
 Dataset_study.rmd | 1 +
 1 file changed, 1 insertion(+)

diff --git a/Dataset_study.rmd b/Dataset_study.rmd
index b16529d..cf52f85 100644
--- a/Dataset_study.rmd
+++ b/Dataset_study.rmd
@@ -505,6 +505,7 @@ x.verif=seq(range(x)[1],range(x)[2],1)
 #dist.theo.scores=lambda0*exp(-lambda0*(A*x.verif-B))
 #dist.theo.scores
 dist.emp.scores
+barplot(dist.emp.scores)
 
 ```
 

From f5a35a44ce5a1a71f770757594179a55e8e2cb07 Mon Sep 17 00:00:00 2001
From: elisaduz <elisa.duz1807@gmail.com>
Date: Tue, 19 Apr 2022 11:24:55 +0200
Subject: [PATCH 2/2] Comparison_Lambda_Values

---
 Comparaison_of_methods.rmd | 27 +++++++++++++++++++++++----
 1 file changed, 23 insertions(+), 4 deletions(-)

diff --git a/Comparaison_of_methods.rmd b/Comparaison_of_methods.rmd
index b0f0a00..e89f2bd 100644
--- a/Comparaison_of_methods.rmd
+++ b/Comparaison_of_methods.rmd
@@ -263,14 +263,35 @@ ScoreDistribEmpiric <- function(lambda0, lambda1, n_sample, T){
 }
 ```
 
+```{r}
+
+lambda0=5
+lambda1=7
+distrib_mc=ScoreDistribEmpiric(lambda0,lambda1,10000,T)
+score_moyen=mean(distrib_mc[,1])
+print(score_moyen)
+score_max=max(distrib_mc[,1])
+print(score_max)
+score_min=min(distrib_mc[,1])
+print(score_min)
+amplitude=abs(score_max-score_min)
+print(amplitude)
+E=ComputeE(lambda0, lambda1)
+print(E)
+barplot(distrib_mc[,2])
+```
+
+
 ```{r}
 ScoreDistribElisa <- function(lambda0, lambda1, T){
     E = ComputeE(lambda0, lambda1)
 
     score_max = floor(E*log(lambda1/lambda0))
+    
 
     ## score_min compute
     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))
@@ -289,11 +310,9 @@ distrib_score_mc=ScoreDistribEmpiric(2,3,10000,T)
 
 distrib_score_theo=ScoreDistribElisa(2,3,T)
 
+distrib_score_mc
+distrib_score_theo
 
-length(distrib_score_mc[,2])
-length(distrib_score_theo[,2])
-
-#diff_distrib_score=abs(distrib_score_mc[,2]-distrib_score_theo[,2])
 
 #par(mfrow = c(1,2))
 barplot(distrib_score_mc[,2],col="blue",axes=F)