Simulation

[Gauthier Dietrich]

Bonjour,
je voudrais faire la simulation suivante sur n pas :

"Faire un saut de +1 ou -1 selon une loi de Bernoulli."
si c'est +1 générer un nombre selon une loi uniforme sur [0, 1] et l'ajouter à L.
si c'est -1 supprimer le inf de la liste L et le placer dans une liste R.


x = sample(c(-1, 1), size=1, replace = T, prob = c(4/9, 5/9))
L = c(0)
u = runif(1, 0, 1)
if (x == 1) {
L = c(L, u)
} else {L = L[L!= min(L)]}

Comment avec n pas compléter les listes L et R ?

[jean lobry]

Bonjour,

je ne suis pas certain d'avoir compris, mais voici une piste :

Code : Tout sélectionner

set.seed(2021)
n <- 10
R <- L <- numeric(0)
for(i in seq_len(n)){
  if(sample(c(-1, 1), 1) == 1){
    L <- c(L, runif(1))
  } else {
    imin <- which.min(L)
    R <- c(R, L[imin])
    L <- L[-imin]
  }
}

Amicalement,
jean lobry