This commit is contained in:
elisaduz 2022-05-08 17:08:57 +02:00
commit f9850545b6
1 changed files with 2 additions and 1 deletions

View File

@ -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){