Discussion :

[Nifaste]

Bonjour à tous !

je dois construire un échantillon de taille N=1000 issu d’une loi exponentielle de paramètre 2 et calculer la moyenne empirique m et la variance empirique s² de E.
Mon problème est que je ne connais pas la fonction R qui permet de creer une loi exponentielle de paramètre 2.
La moyenne se calcule avec mean et la variance avec sd, c'est bien ca ?
J'ai regardé, il y a rexp(n,rate) mais je ne sais pas où mettre mon paramètre ...
Quelqu'un pourrait m'eclairer ?
Merci d'avance!

[Ayest]

bonjour,

La première chose à faire pour avoir de l'aide sur une fonction :

Code :

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

1
2
 
R> ?rexp

Tu y verras que 'n' représente le nombre de tirage aléatoire que tu veux obtenir, et 'rate' est le paramètre de la loi exponenetielle.

Pour compléter cette réponse : dans l'aide, certains paramètres sont écrit avec '=_valeur_' à coté de leur nom, ce la signifie que ces paramètres ont une valeur par défaut.
Ex :

Code :

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

1
2
 
rexp(n, rate = 1)

signifie que si tu n'indique pas de valeur pour le paramètre 'rate', il vaudra 1, c'est à dire :

Code :

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

1
2
 
rexp(1000)

va générer 1000 réalisation d'une loi exponentielle de paramètre 1.

[Nifaste]

merci beaucoup !! je ne savais pas pour "rate" !

Une autre question :
je dois effectuer 10 tirages aléatoires de taille 1000 avec remise dans E ( = rexp(100,2), j'ai fait :

Code :

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

1
2
E=rexp(100,2) 
replicate(10,sample(E,1000,replace='TRUE'))

cela vous semble t'il juste ?

Merci d'avance

[Ayest]

C'est juste.
Si tu veux vraiment t'en convaincre, et voir ce que ca donne, fais le avec des chiffres plus petits que 100 et 1000 (genre 5 et 7), ca permet de voir mieux le résultats.

Merci de mettre ton message comme "Résolu" en cliquant sur le bouton... "Résolu".

[Nifaste]

Un autre question :
Je dois decouper mon vecteur E en 10 parts egales : je dois utiliser cut () ou seq () ?
Car ensuite je dois calculer la moyenne des 10 parts, et cut() me sort qq chose du type "factor" donc je n'y arrive pas ...

[Ayest]

'cut' est utilisé pour faire des intervalles à partir d'un vecteur numérique.

Mais vraiment, il faut lire la documentation des fonctions avant de poster sur le forum.

Tu peux utiliser les fonctions, 'rep' et 'tapply' pour calculer tes moyennes :

Code :

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

1
2
 
R> tapply(E,rep(1:10,each=10),mean)

Une autre solution consiste à transformer ton vecteur de taille 100 en une matrice de taille 10x10, et utiliser la fonction 'apply'

Code :

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

1
2
3
 
R> dim(E)=c(10,10)
R> apply(E,1,mean)