From c9c76fcc1d9bd16d6cdc98461b4f5741c6068ee5 Mon Sep 17 00:00:00 2001
From: Paul-Corbalan <58653590+Paul-Corbalan@users.noreply.github.com>
Date: Tue, 8 Mar 2022 10:42:55 +0100
Subject: [PATCH] Local score compute with Daudin

---
 Dataset_study.rmd | 14 ++++++--------
 1 file changed, 6 insertions(+), 8 deletions(-)

diff --git a/Dataset_study.rmd b/Dataset_study.rmd
index 8952e4c..174bb11 100644
--- a/Dataset_study.rmd
+++ b/Dataset_study.rmd
@@ -259,17 +259,15 @@ sst
 Local score
 ```{r}
 library("localScore")
-
-X = floor(log(dexp(tbe1, rate = lambda1)/dexp(tbe1, rate = lambda0)))
+E = 10
+X = floor(E*log(dexp(tbe1, rate = lambda1)/dexp(tbe1, rate = lambda0)))
 
 max_X = max(X)
 min_X = min(X)
-P_X = table(X)/length(X)
-hist(X)
+P_X = table(factor(X, levels = min_X:max_X))/length(X)
 
-localScoreC_double(X)
+LS=localScoreC(X)$localScore[1]
 
-l = list()
-l[["X"]] = X
-result = automatic_analysis(l, "iid")
+result = daudin(localScore = LS, score_probabilities = P_X, sequence_length = length(x), sequence_min = min_X, sequence_max = max_X)
+result
 ```