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

Cahnge lambda numbers for hypothesis

This commit is contained in:
Paul-Corbalan 2022-02-15 08:37:29 +01:00
parent 175160d7e9
commit 70e798cf33
1 changed files with 8 additions and 8 deletions
Show Stats Download Patch File Download Diff File Expand all files Collapse all files

16
Dataset_study.rmd
View File

@ -140,14 +140,14 @@ PoissonProcess <- function(lambda,T) {
lambda1=2 lambda0=2
lambda2=3 lambda1=3
Ti=10 Ti=10
pp1=PoissonProcess(lambda1,Ti) pp1=PoissonProcess(lambda0,Ti)
print(pp1) print(pp1)
plot(c(0,pp1),0:length(pp1),type="s",xlab="time t",ylab="number of events by time t") plot(c(0,pp1),0:length(pp1),type="s",xlab="time t",ylab="number of events by time t")
pp2=PoissonProcess(lambda2,Ti) pp2=PoissonProcess(lambda1,Ti)
print(pp2) print(pp2)
plot(c(0,pp2),0:length(pp2),type="s",xlab="time t",ylab="number of events by time t") plot(c(0,pp2),0:length(pp2),type="s",xlab="time t",ylab="number of events by time t")
@ -160,9 +160,9 @@ n2=length(pp2)
tbe2=pp2[2:n2]-pp2[1:n2-1] tbe2=pp2[2:n2]-pp2[1:n2-1]
tbe2 tbe2
ks.test(tbe1,pexp,lambda1, alternative="two.sided") ks.test(tbe1,pexp,lambda0, alternative="two.sided")
ks.test(tbe2,pexp,lambda2, alternative="two.sided") ks.test(tbe2,pexp,lambda1, alternative="two.sided")
@ -177,13 +177,13 @@ Je reprends votre code pour faire un data set :
```{r} ```{r}
# Etape 1 : simu Poisson process sous H0 # Etape 1 : simu Poisson process sous H0
ppH0=PoissonProcess(lambda1,Ti) ppH0=PoissonProcess(lambda0,Ti)
ppH0 ppH0
length(ppH0) length(ppH0)
# Etape 2 : creation d'un segment sous H1 # Etape 2 : creation d'un segment sous H1
tau= 2.5 # longeur de l'intervalle modifie, a fortiori tau < Ti tau= 2.5 # longeur de l'intervalle modifie, a fortiori tau < Ti
ppH1.segt=PoissonProcess(lambda2,tau) ppH1.segt=PoissonProcess(lambda1,tau)
ppH1.segt ppH1.segt
length(ppH1.segt) length(ppH1.segt)