From b0dd092c3470f5f802ccb4c0b7f504300cec4cd3 Mon Sep 17 00:00:00 2001 From: elisaduz Date: Mon, 11 Apr 2022 17:03:10 +0200 Subject: [PATCH] Update_new_function_score --- Dataset_study.rmd | 12 +++++++++--- 1 file changed, 9 insertions(+), 3 deletions(-) diff --git a/Dataset_study.rmd b/Dataset_study.rmd index 32c5455..243a9fc 100644 --- a/Dataset_study.rmd +++ b/Dataset_study.rmd @@ -361,7 +361,7 @@ max #0)>=max)<10^-9 t=seq(0,max,0.1) probaseuilmin = 10^-9 -x=floor(E*log(lambda1/lambda0)+(lambda0-lambda1)*t) +x=floor(E*log(lambda1/lambda0)+E*(lambda0-lambda1)*t) head(x) range(x) @@ -371,8 +371,14 @@ s ##tronquer la queue des x (x négatifs) A=1/(lambda0-lambda1) B=A*log(lambda1/lambda0) -proba.s=lambda0*exp(-lambda0*(A*s-B)) -proba.s = proba.s/sum(proba.s) +proba.s=(lambda1/lambda0)*(exp(lambda0/(lambda0-lambda1)*(1-t/E))*exp(-lambda0/(E*lambda0-lambda1))-1) + +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 barplot(proba.s)