Discussion :

[selver057]

Bonjour à tous,

J'ai une variable x qui prend les valeurs suivantes : -10, 20, 100. Chaque valeur est associée à la probabilité : 0.7, 0.2, 0.1.

Je souhaite faire une simulation à partir de plusieurs échantillons de taille 100 (par exemple une cinquantaine d'échantillons) pour voir comment évolue la somme des valeurs de chaque échantillon (un échantillon = 1 somme, deux échantillons = somme de l'échantillon 1 + somme de l'échantillon 2 etc). je souhaite représenter les différents résultats obtenus (résultat = somme) sur un histogramme. Voici comment j'ai fait, j'ai créé:

- un vecteur x pour les valeurs de la variable,
- un vecteur y pour les probabilités,
- un vecteur z1 qui représente un échantillon,
- un vecteur z2 qui recueille la somme de chaque échantillon créé

Pour le moment mon code se limite aux commandes suivantes :

Code :

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

1
2
3
4
5
x<-c(-10,20,100)
y<-c(0.7,0.2,0.1)
z1<-sample(x,100,replace=T,prob=y)
z2=sumz1
hist(z2, freq=FALSE, col="dark green")

Vous constaterez qu'il ne me permet pas d'arriver au résultat que je cherche. J'imagine qu'il y a moyen de l'améliorer en introduisant une boucle qui me permettrait de récupérer la somme de(s) échantillon(s) dans un vecteur puis de créer mon histogramme à partir de ce vecteur. Pouvez-vous m'indiquer comment faire ? D'avance je vous remercie

[hoccha]

Bonsoir,
Voici une proposition:

Code :

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

1
2
3
4
5
6
x  <- c(-10 , 20 , 100)
y  <- c(0.7 , 0.2, 0.1)
z1 <- matrix(sample(x,100*50,replace=T,prob=y),ncol=50)    # matrice de 100 lignes et 50 colonnes.
z2 <- apply(z1,2,sum)      # Je somme les colonnes --> resulat: vecteur de longeur 50.
res <- cumsum(z2)          # Somme cumulée du vecteur
barplot(res)

J'espère avoir pu vous aider.