This commit is contained in:
elisaduz 2022-05-08 17:09:20 +02:00
commit 00642de65a
1 changed files with 4 additions and 4 deletions

View File

@ -292,7 +292,7 @@ barplot(distrib_mc[,2])
```{r} ```{r}
ScoreDistribElisa <- function(lambda0, lambda1, T){ ScoreDistribTheo <- function(lambda0, lambda1, T){
E = ComputeE(lambda0, lambda1) E = ComputeE(lambda0, lambda1)
score_max = floor(E*log(lambda1/lambda0)) score_max = floor(E*log(lambda1/lambda0))
@ -326,7 +326,7 @@ distrib_score_theo
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 = ScoreDistribTheo(2,3,T)
plot_graph_distrib_score <- function(distrib_score_theo, distrib_score_mc){ plot_graph_distrib_score <- function(distrib_score_theo, distrib_score_mc){
# length(distrib_score_mc[,2]) # length(distrib_score_mc[,2])
@ -434,7 +434,7 @@ for (lambda0 in (2)){
#print(confusionMatrix(factor(LS_H0$class), factor(SS_H0$class))) #print(confusionMatrix(factor(LS_H0$class), factor(SS_H0$class)))
#cat("- Elisa version:\n") #cat("- Elisa version:\n")
Score = ScoreDistribElisa(lambda0, lambda1, T) Score = ScoreDistribTheo(lambda0, lambda1, T)
Emp = EmpDistrib(lambda0,n_sample,T,tau) Emp = EmpDistrib(lambda0,n_sample,T,tau)
X_seq = Score$Score_X X_seq = Score$Score_X
@ -521,7 +521,7 @@ for (i in (1:NbSeq)) {
#print(confusionMatrix(factor(LS_H0$class), factor(SS_H0$class))) #print(confusionMatrix(factor(LS_H0$class), factor(SS_H0$class)))
#cat("- Elisa version:\n") #cat("- Elisa version:\n")
Score = ScoreDistribElisa(lambda0, lambda1, T) Score = ScoreDistribTheo(lambda0, lambda1, T)
Emp = EmpDistrib(lambda0,n_sample,T,tau) Emp = EmpDistrib(lambda0,n_sample,T,tau)
X_seq = Score$Score_X X_seq = Score$Score_X