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

Add Elisa code

This commit is contained in:
Paul-Corbalan 2022-04-12 10:25:21 +02:00
parent cfc4097a1a
commit 99cb19a9d9
1 changed files with 13 additions and 3 deletions
Show Stats Download Patch File Download Diff File Expand all files Collapse all files

16
Comparaison_of_methods.rmd
View File

@ -295,9 +295,19 @@ LocaScoreMC <- function(lambda0, lambda1, NbSeq, tbe0, T){
pvalue = c() pvalue = c()
X = c() X = c()
Score = ScoreDistrib(lambda0, lambda1, NbSeq, T) score_max=floor(E*log(lambda1/lambda0))
xp = Score$X
P_X = Score$P_X ## Calcul score_min
score_min_c=floor(E*log(lambda1/lambda0)+E*(lambda0-lambda1)*T)
l=seq(score_min_c,score_max,1)
borne_inf=(l-E*log(lambda1/lambda0))/(E*(lambda0-lambda1))
borne_sup=(l+1-E*log(lambda1/lambda0))/(E*(lambda0-lambda1))
proba.l=pexp(rate=lambda0,borne_inf)-pexp(rate=lambda0,borne_sup)
S=sum(proba.l)
new.proba.s=proba.l/S
xp = l
P_X = proba.l/S
min_X = min(xp) min_X = min(xp)
max_X = max(xp) max_X = max(xp)