Discussion :

[vicjohns]

Bonjour à tous,

Je viens solliciter votre aide sur cette problématique, je dispose d'un échantillon de taille 40 et j'essaye de créer 1000 échantillons bootstrap de taille 40 chacunes, mais je ne sais pas trop comment le faire.

Merci d'avance

[tototode]

Bonjour,

La fonction sample est parfaite pour faire des tirages aléatoires dans un vecteur ou pour faire un tirage aléatoire sur le numéro de tes lignes d'un data.frame par exemple. Par contre il faut bien spécifier l'argument replace = TRUE pour que la même valeur puisse être choisie plusieurs fois.

Code :

Sélectionner tout - Visualiser dans une fenêtre à part

1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
x <- rnorm(40)
 
sims <- replicate(1000, sample(x, 40, replace = T))
# chaque colonne correspond à un tirage
 
# ou encore
sims <- sapply(1:1000, function(u) sample(x, 40, replace = TRUE))
 
# ou encore
sims <- matrix(numeric(), 40, 1000)
for (i in 1:1000) {
  sims[, i] <- sample(x, 40, replace = TRUE)
  } 
 
# ou encore
sims <- sample(x, 40 * 1000,  replace = TRUE)
dim(sims) <- c(40, 1000)

cdlt