diff --git a/Dataset_study.rmd b/Dataset_study.rmd index 44c550d..262b344 100644 --- a/Dataset_study.rmd +++ b/Dataset_study.rmd @@ -170,6 +170,58 @@ ks.test(tbe2,pexp,lambda2, alternative="two.sided") ``` The Kolmogorov-Smirnov test rejects the hypothesis that the time between events sequence is following an exponential distribution. +\section{Proposition de simulation sous H1} + +Je reprends votre code pour faire un data set : + + +```{r} +# Etape 1 : simu Poisson process sous H0 +ppH0=PoissonProcess(lambda1,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 +length(ppH1.segt) + +# Etape 3 : insertion du segment dans la sequence H0 +dbt=runif(1,0,Ti-tau) # choix de l'indice de temps ou va commencer le segment modifie +dbt +ppH1.repo=dbt+ppH1.segt # repositionnement des observations dans le temps +ppH1.repo +ppH0_avant=ppH0[which(ppH0ppH1.repo[length(ppH1.repo)])] +ppH1=c(ppH0_avant,ppH1.repo,ppH0_apres) +ppH1 +length(ppH1) + + +#time between events +n1=length(ppH1) +tbe1=ppH1[2:n1]-ppH1[1:n1-1] + +n0=length(ppH0) +tbe0=ppH0[2:n0]-ppH0[1:n0-1] + +tbe1=c(0,tbe1) +tbe1 + +list1=data.frame(ProcessusPoissonH1=ppH1, + TimeBetweenEventH1=tbe1) +list1 + +tbe0=c(0,tbe0) +tbe0 + +list0=data.frame(ProcessusPoissonH0=ppH0, + TimeBetweenEventH0=tbe0) +list0 + +``` + Import data of rainfall in France every 3 hours. ```{r} Rain_Dataset = read.csv("data/synop.202202.csv", sep = ";")