From 70e798cf336d89c1198114b69f20ed0b56c39110 Mon Sep 17 00:00:00 2001 From: Paul-Corbalan <58653590+Paul-Corbalan@users.noreply.github.com> Date: Tue, 15 Feb 2022 08:37:29 +0100 Subject: [PATCH] Cahnge lambda numbers for hypothesis --- Dataset_study.rmd | 16 ++++++++-------- 1 file changed, 8 insertions(+), 8 deletions(-) diff --git a/Dataset_study.rmd b/Dataset_study.rmd index 262b344..207a413 100644 --- a/Dataset_study.rmd +++ b/Dataset_study.rmd @@ -140,14 +140,14 @@ PoissonProcess <- function(lambda,T) { -lambda1=2 -lambda2=3 +lambda0=2 +lambda1=3 Ti=10 -pp1=PoissonProcess(lambda1,Ti) +pp1=PoissonProcess(lambda0,Ti) print(pp1) 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) 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 -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} # Etape 1 : simu Poisson process sous H0 -ppH0=PoissonProcess(lambda1,Ti) +ppH0=PoissonProcess(lambda0,Ti) ppH0 length(ppH0) # Etape 2 : creation d'un segment sous H1 tau= 2.5 # longeur de l'intervalle modifie, a fortiori tau < Ti -ppH1.segt=PoissonProcess(lambda2,tau) +ppH1.segt=PoissonProcess(lambda1,tau) ppH1.segt length(ppH1.segt)