Renorme probas s

This commit is contained in:
Paul-Corbalan 2022-03-15 10:48:14 +01:00
parent 88a52200a8
commit fb7958adae
1 changed files with 2 additions and 1 deletions

View File

@ -315,10 +315,11 @@ for (lambda0 in c(5,50,500,5000,50000) ) {
B=A*log(lambda1/lambda0)
probaS=floor(E*lambda0*exp(-lambda0*(A*x[j]-B)))
print(S)
# print(s)
probaS_vect[j]=probaS
}
probaS_vect = probaS_vect/sum(probaS_vect)
print(probaS_vect)
p <- barplot(probaS_vect, main=str_c("lambda0=", as.character(lambda0), "et", "lambda1=", as.character(lambda1), sep=" " ))
p