Update Dataset_study.rmd

Update_Feedback_List

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(ppH0<ppH1.repo[1])]
+ppH0_apres=ppH0[which(ppH0>ppH1.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 = ";")