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

Change indentation

Browse Source
This commit is contained in:
Paul-Corbalan 2022-04-05 08:14:51 +02:00
parent fb66d1307a
commit bf50e2294c
1 changed files with 13 additions and 13 deletions
Show Stats Download Patch File Download Diff File Expand all files Collapse all files

18
Comparaison_of_methods.rmd
View File

@ -172,17 +172,17 @@ E = 10
pvalue=c() pvalue=c()
for (i in 1:NbSeqH0){ for (i in 1:NbSeqH0){
X = floor(E*log(dexp(tbe0[[i]], rate = lambda1)/dexp(tbe0[[i]], rate = lambda0))) X = floor(E*log(dexp(tbe0[[i]], rate = lambda1)/dexp(tbe0[[i]], rate = lambda0)))
max_X = max(X) max_X = max(X)
min_X = min(X) min_X = min(X)
P_X = table(factor(X, levels = min_X:max_X))/length(X) P_X = table(factor(X, levels = min_X:max_X))/length(X)
LS=localScoreC(X)$localScore[1] LS=localScoreC(X)$localScore[1]
result = daudin(localScore = LS, score_probabilities = P_X, sequence_length = length(X), sequence_min = min_X, sequence_max =max_X) result = daudin(localScore = LS, score_probabilities = P_X, sequence_length = length(X), sequence_min = min_X, sequence_max =max_X)
pvalue=c(pvalue,result) pvalue=c(pvalue,result)
} }
LS_H0=data.frame(num=1:NbSeqH0, pvalue_scan=pvalue, class=(pvalue<0.05)) LS_H0=data.frame(num=1:NbSeqH0, pvalue_scan=pvalue, class=(pvalue<0.05))
LS_H0 LS_H0
@ -211,8 +211,8 @@ E = 1
maxXk = floor(E*(log(lambda1/lambda0))) maxXk = floor(E*(log(lambda1/lambda0)))
maxXk maxXk
while (maxXk < 3) { while (maxXk < 3) {
E = E+1 E = E+1
maxXk = floor(E*(log(lambda1/lambda0))) maxXk = floor(E*(log(lambda1/lambda0)))
} }
print(E) print(E)