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

Add NbSeq.NonNulles to LocalScoreMC

Browse Source
This commit is contained in:
Paul-Corbalan 2022-05-10 18:45:12 +02:00
parent 29efa02528
commit 0054d077be
1 changed files with 22 additions and 18 deletions
Show Stats Download Patch File Download Diff File Expand all files Collapse all files

40
Comparaison_of_methods.rmd
View File

@ -356,24 +356,28 @@ plot_graph_distrib_score(distrib_score_theo, distrib_score_mc)
### 3.2. Local score calculation ### 3.2. Local score calculation
```{r} ```{r}
LocalScoreMC <- function(lambda0, lambda1, NbSeq, T, X_seq, P_X, tbe0){ LocalScoreMC <- function(lambda0, lambda1, NbSeq, T, X_seq, P_X, tbe0){
E = ComputeE(lambda0, lambda1) E = ComputeE(lambda0, lambda1)
pvalue = c() pvalue = c()
X = c()
min_X = min(X_seq)
min_X = min(X_seq) max_X = max(X_seq)
max_X = max(X_seq)
NbSeq.NonNulles = 0
for (i in 1:NbSeq){
x = floor(E*log(dexp(tbe0[[i]], rate = lambda1)/dexp(tbe0[[i]], rate = lambda0))) for (i in 1:NbSeq){
LS = localScoreC(x)$localScore[1] x = floor(E*log(dexp(tbe0[[i]], rate = lambda1)/dexp(tbe0[[i]], rate = lambda0)))
if (length(x)!=0){
daudin_result = daudin(localScore = LS, score_probabilities = P_X, sequence_length = length(x), sequence_min = min_X, sequence_max = max_X) LS = localScoreC(x)$localScore[1]
options(warn = -1) # Disable warnings print
daudin_result = daudin(localScore = LS, score_probabilities = P_X, sequence_length = length(x), sequence_min = min_X, sequence_max = max_X)
pvalue = c(pvalue, daudin_result) options(warn = -1) # Disable warnings print
pvalue = c(pvalue, daudin_result)
NbSeq.NonNulles = NbSeq.NonNulles + 1
}
} }
LS_H0=data.frame(num=1:NbSeq, pvalue_scan=pvalue, class=as.numeric(pvalue<0.05)) LS_H0=data.frame(num=1:NbSeq.NonNulles, pvalue_scan=pvalue, class=as.numeric(pvalue<0.05))
return(LS_H0) return(LS_H0)
} }
``` ```
@ -451,7 +455,7 @@ CompareMethods <- function(lambda0, lambda1, NbSeq, T, tau_scan, tau_H1){
``` ```
```{r} ```{r}
NbSeq = 10**2 NbSeq = 10**4
T = 10 T = 10
tau = 2 tau = 2