Discussion :

[Puppet_Master]

Bonjour.

Voici un petit bout de code:

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
v <- c()
for( i in runif(n = 10, min = -7, max = 7) ){
    v <- c( v, rnorm(n = 30, mean = i, sd = 3))
}
 
id <- rep(seq(from = 1, to = 10, by = 1), each = 30)
 
df <- data.frame(
      cat  = as.factor(id)
    , data = v
)
 
library(ggplot2)
ggplot(df) +
    aes(x = cat, y = data) +
    geom_boxplot()

La question est simple:

Peut-on obtenir la même chose sans boucle for ?

La boucle me sert à créer un vecteur de 10 tirages de 30 valeurs d'une distribution normale, avec 10 moyennes aléatoires entre deux bornes.

[xavier-Pierre]

bonjour ,
avec le tidyverse (en particulier purrr)
tu peux creer ton dataframe

Code :

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

1
2
3
library(tidyverse)
df <- map2_df(rep(seq(from = 1, to = 10, by = 1)),runif(n = 10, min = -7, max = 7),
             ~tibble(cat =.x,data=rnorm(n= 30,mean = .y,sd = 3)))

[tototode]

Bonjour,

tu as deux possibilités pour le faire sans boucle, soit tu utilises le recyclage sur le vecteur de moyenne, soit tu répètes chaque moyenne 30 fois chacune.
Les deux te donneront le même résultat mais par contre pas dans le même ordre. Dans le premier cas les données seront dans l'ordre g1, g2, g3, ..., g10, g1, g2, ... alors que dans le deuxième cas les données d'un même groupe seront à la suite puis une fois les 30 valeurs obtenues tu passeras au groupe suivant :

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
18
19
set.seed(100)
v <- c()
for( i in runif(n = 10, min = -7, max = 7) ){
    v <- c( v, rnorm(n = 30, mean = i, sd = 3))
}
 
# recyclage
set.seed(100)
mns <- runif(10, -7, 7)
v1 <- rnorm(30*length(mns), mns, 3)
id1 <- rep(1:10, 30)
 
# repetition des moyennes
set.seed(100)
mns <- runif(10, -7, 7)
v2 <- rnorm(30*length(mns), rep(mns, each = 30), 3)
id2 <- rep(1:10, each = 30)
 
plot(v, v2)

cdlt

[Puppet_Master]

Bonjour.

Merci pour les solutions !