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

Update_new_function_score

Browse Source
This commit is contained in:
elisaduz 2022-04-11 17:03:10 +02:00
parent 8b192f36eb
commit b0dd092c34
1 changed files with 9 additions and 3 deletions
Show Stats Download Patch File Download Diff File Expand all files Collapse all files

12
Dataset_study.rmd
View File

@ -361,7 +361,7 @@ max
#0)>=max)<10^-9 #0)>=max)<10^-9
t=seq(0,max,0.1) t=seq(0,max,0.1)
probaseuilmin = 10^-9 probaseuilmin = 10^-9
x=floor(E*log(lambda1/lambda0)+(lambda0-lambda1)*t) x=floor(E*log(lambda1/lambda0)+E*(lambda0-lambda1)*t)
head(x) head(x)
range(x) range(x)
@ -371,8 +371,14 @@ s
##tronquer la queue des x (x négatifs) ##tronquer la queue des x (x négatifs)
A=1/(lambda0-lambda1) A=1/(lambda0-lambda1)
B=A*log(lambda1/lambda0) B=A*log(lambda1/lambda0)
proba.s=lambda0*exp(-lambda0*(A*s-B)) proba.s=(lambda1/lambda0)*(exp(lambda0/(lambda0-lambda1)*(1-t/E))*exp(-lambda0/(E*lambda0-lambda1))-1)
proba.s = proba.s/sum(proba.s)
score_min_c=floor(E*log(lambda1/lambda0)+E*(lambda0-lambda1)*max)
somme=(lambda1/lambda0)*(exp(-lambda0/(E*(lambda0-lambda1))-1))* exp(lambda0/(lambda0-lambda1))*exp(lambda0/(E*(lambda0-lambda1))*(1-score_min_c))/(1-exp(lambda0/(E*(lambda0-lambda1))))
#proba.s = proba.s/sum(proba.s)
proba.s= proba.s/(1-somme)
proba.s proba.s
barplot(proba.s) barplot(proba.s)