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

Improve printing for LocaScoreMC experience plan

Browse Source
This commit is contained in:
Paul-Corbalan 2022-04-12 18:16:06 +02:00
parent 18dc50fa3f
commit 4d5bd81071
1 changed files with 3 additions and 2 deletions
Show Stats Download Patch File Download Diff File Expand all files Collapse all files

5
Comparaison_of_methods.rmd
View File

@ -301,7 +301,7 @@ T = 10
for (lambda0 in (2:5)){ for (lambda0 in (2:5)){
for (lambda1 in c(2,4,6)){ for (lambda1 in c(2,4,6)){
if (lambda0 < lambda1){ if (lambda0 < lambda1){
cat("Nb = ", NbSeq, ", lambda0 = ", lambda0, ", lambda1 = ", lambda1, "\n", sep = "") cat("For T = ", T, ", Nb = ", NbSeq, "lambda0 = ", lambda0, "and lambda1 = ", lambda1, ":\n", sep = "")
tbe0=vector("list",length=NbSeq) tbe0=vector("list",length=NbSeq)
for (i in (1:NbSeq)) { for (i in (1:NbSeq)) {
@ -311,6 +311,7 @@ for (lambda0 in (2:5)){
tbe0[[i]]=tbei tbe0[[i]]=tbei
} }
cat("- Empiric version:\n")
Score = ScoreDistribEmpiric(lambda0, lambda1, NbSeq, T) Score = ScoreDistribEmpiric(lambda0, lambda1, NbSeq, T)
X_seq = Score$Score_X X_seq = Score$Score_X
P_X = Score$P_X P_X = Score$P_X
@ -318,8 +319,8 @@ for (lambda0 in (2:5)){
LS_H0 = LocaScoreMC(lambda0, lambda1, NbSeq, T, X_seq, P_X, tbe0) LS_H0 = LocaScoreMC(lambda0, lambda1, NbSeq, T, X_seq, P_X, tbe0)
print(summary(LS_H0)) print(summary(LS_H0))
cat("-\n")
cat("- Elisa version:\n")
Score = ScoreDistribElisa(lambda0, lambda1, T) Score = ScoreDistribElisa(lambda0, lambda1, T)
X_seq = Score$Score_X X_seq = Score$Score_X
P_X = Score$P_X P_X = Score$P_X