diff --git a/Comparaison_of_methods.rmd b/Comparaison_of_methods.rmd index aa43869..c5a5e8f 100644 --- a/Comparaison_of_methods.rmd +++ b/Comparaison_of_methods.rmd @@ -316,6 +316,7 @@ ScoreDistribElisa <- function(lambda0, lambda1, T){ ```{r} +T=10 distrib_score_mc=ScoreDistribEmpiric(2,3,10000,T) distrib_score_theo=ScoreDistribElisa(2,3,T) @@ -333,6 +334,7 @@ plot_graph_distrib_score <- function(distrib_score_theo, distrib_score_mc){ #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) mtext("Distribution of scores via Monte Carlo",side=1,line=2.5,col="blue") @@ -346,7 +348,6 @@ plot_graph_distrib_score <- function(distrib_score_theo, distrib_score_mc){ plot_graph_distrib_score(distrib_score_theo, distrib_score_mc) ``` - ### 3.2. Local score calculation ```{r} LocalScoreMC <- function(lambda0, lambda1, NbSeq, T, X_seq, P_X, tbe0){