paul-corbalan
/
Scan-Statistics-Project-4Y-INSA
1
0
Fork 0
Code Issues Pull Requests Packages Projects Releases Wiki Activity

Passing to function plot_graph_distrib_score

This commit is contained in:
Paul-Corbalan 2022-04-19 10:30:49 +02:00
parent ff0f6cf253
commit b7e079a55e
1 changed files with 15 additions and 13 deletions
Show Stats Download Patch File Download Diff File Expand all files Collapse all files

28
Comparaison_of_methods.rmd
View File

@ -285,23 +285,25 @@ ScoreDistribElisa <- function(lambda0, lambda1, T){
```
```{r}
distrib_score_mc=ScoreDistribEmpiric(2,3,10000,T)
distrib_score_mc = ScoreDistribEmpiric(2,3,10000,T)
distrib_score_theo = ScoreDistribElisa(2,3,T)
distrib_score_theo=ScoreDistribElisa(2,3,T)
plot_graph_distrib_score <- function(distrib_score_theo, distrib_score_mc){
# length(distrib_score_mc[,2])
# length(distrib_score_theo[,2])
#diff_distrib_score=abs(distrib_score_mc[,2]-distrib_score_theo[,2])
length(distrib_score_mc[,2])
length(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")
axis(2, ylim=c(0,10))
par(new = T)
barplot(distrib_score_theo[,2],col="red",axes=F)
mtext("Distribution of scores using the theoretical method",side=1,line=4,col="red")
}
#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 des scores via Monte_Carlo",side=1,line=2.5,col="blue")
axis(2, ylim=c(0,10))
par(new = T)
barplot(distrib_score_theo[,2],col="red",axes=F)
mtext("Distribution des scores via la méthode théorique",side=1,line=4,col="red")
plot_graph_distrib_score(distrib_score_theo, distrib_score_mc)
```