Add Elisa code

This commit is contained in:
Paul-Corbalan 2022-04-19 11:29:43 +02:00
parent 6e3480d478
commit 6a66069040
1 changed files with 17 additions and 0 deletions

View File

@ -272,6 +272,23 @@ 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} ```{r}
ScoreDistribElisa <- function(lambda0, lambda1, T){ ScoreDistribElisa <- function(lambda0, lambda1, T){
E = ComputeE(lambda0, lambda1) E = ComputeE(lambda0, lambda1)