Répondre à la discussion
Affichage des résultats 1 à 7 sur 7

résolution d'un exercice de probabilité sous le logiciel R



  1. #1
    Lza

    résolution d'un exercice de probabilité sous le logiciel R


    ------

    bonjour,
    je travaille actuellement sur un projet à traiter sous R.
    j'aimerai savoir si certains travaillent dessus pour avoir de l'aide.
    Je dois illustrer sur des exemples le théorème central limite:
    Si X1,X2,…Xn sont des variables indépendantes de même loi, d’espérance m et d’écart-type s, on démontre que quand n tend vers l’infini, la loi de

    [(X1+…+Xn)-n*m]/s*sqrt(n) tend vers la loi normale centrée réduite.

    On choisira diverses lois pour les Xi, et on considèrera en particulier la convergence de la loi
    Binomiale vers la loi normale, de la loi de Poisson vers la loi normale.

    Merci d'avance.

    -----

  2. Publicité
  3. #2
    HigginsVincent

    Re : résolution d'un exercice de probabilité sous le logiciel R

    Bonjour,

    Pour faire des tirages d'une variable aléatoire suivant une loi de Poisson ou Binomiale, tu peux essayer les fonctions rpois et rbinom...

    Je n'ai pas compris trop en quoi précisément tu voulais de l'aide !

    Bon courage,
    V.

  4. #3
    Lza

    Re : résolution d'un exercice de probabilité sous le logiciel R

    merci d'avoir répondu. en fait je me creuse la tête pour savoir exactement ce qui est demandé c'est un travail surtout de maths et non d'informatique enfin il ne restera plus qu'à faire les graphes . Ce que j'ai compris de ce sujet est qu'il va falloir sans doute comparer les courbes afin de montrer que la courbe de la loi de binomiale avec des n grands tend vers la courbe de la loi normale il en va de même pour la loi de poisson. Mais je ne sais comment m'y prendre
    enfin je vais continuer à réfléchir...
    merci ^^

  5. #4
    HigginsVincent

    Re : résolution d'un exercice de probabilité sous le logiciel R

    Je ne comprends pas ce qui est compliqué :

    1 - déclarer (l'entier du théorème)
    2 - générer réalisations d'une variable aléatoire X "de Poisson", par exemple, avec le paramètre que tu veux
    3 - faire la somme de toutes ces réalisations
    4 - répéter les étapes 2 et 3 fois : tu obtiens réalisations d'une variable que tu peux appeler Y, par exemple
    5 - tracer l'histogramme des réalisations de Y

    Mais peut-être que je n'ai pas bien compris le problème ?

    Cordialement,
    V.

  6. A voir en vidéo sur Futura
  7. #5
    invite986312212
    Invité

    Re : résolution d'un exercice de probabilité sous le logiciel R

    Citation Envoyé par Lza Voir le message
    merci d'avoir répondu. en fait je me creuse la tête pour savoir exactement ce qui est demandé c'est un travail surtout de maths et non d'informatique

    quand j'étais en math sup, le prof de physique m'a demandé de dessiner un triangle rectangle sur le tableau, puis à l'aide d'une règle graduée longue d'1m il m'a fait mesurer les côtés, enfin il m'a demandé de vérifier expérimentalement le théorème de Pythagore. Je lui ai répondu que cela n'avait pas lieu d'être et je pense que j'ai dû avoir une très mauvaise note...

    autrement, tu peux peut-être, une fois que tu auras un échantillon de moyennes, utiliser le test de Kolomogorov-Smirnov (une distribuion empririque contre une théorique). C'est un test peu puissant alors tu n'as pas besoin d'un n énorme (n=10 devrait suffire largement)

    sous R: ?ks.test

  8. #6
    HigginsVincent

    Re : résolution d'un exercice de probabilité sous le logiciel R

    Il s'agit de redémontrer le Théorème ?

    Dans ce cas tu peux te servir de la propriété suivante :
    La loi d'une somme de variables aléatoires indépendantes est le produit de convolution de toutes les lois de ces variables...


    Citation Envoyé par ambrosio
    C'est un test peu puissant alors tu n'as pas besoin d'un n énorme (n=10 devrait suffire largement)
    Je ne comprends pas pourquoi le fait d'avoir un test peu puissant permet de se limiter à n=10

  9. Publicité
  10. #7
    invite986312212
    Invité

    Re : résolution d'un exercice de probabilité sous le logiciel R

    et ben j'avais tort, avec n=10 ça passe pas:

    > ks.test((apply(matrix(rbinom(2 0000,20,0.5),ncol=10),1,mean)-10)/sqrt(20*0.25/10),"pnorm")

    One-sample Kolmogorov-Smirnov test

    data: (apply(matrix(rbinom(20000, 20, 0.5), ncol = 10), 1, mean) - 10)/sqrt(20 * 0.25/10)
    D = 0.0403, p-value = 0.003005
    alternative hypothesis: two-sided


    et avec n=20 ça marche:

    > ks.test((apply(matrix(rbinom(2 0000,20,0.5),ncol=20),1,mean)-10)/sqrt(20*0.25/20),"pnorm")

    One-sample Kolmogorov-Smirnov test

    data: (apply(matrix(rbinom(20000, 20, 0.5), ncol = 20), 1, mean) - 10)/sqrt(20 * 0.25/20)
    D = 0.0369, p-value = 0.1322
    alternative hypothesis: two-sided

Discussions similaires

  1. resolution d'un system sous MATLAB
    Par ABN84 dans le forum Mathématiques du supérieur
    Réponses: 1
    Dernier message: 16/12/2007, 18h15
  2. Résolution algébrique d'un système d'équations sous Matlab
    Par jbollagnier dans le forum Mathématiques du supérieur
    Réponses: 6
    Dernier message: 22/10/2007, 10h45
  3. Un exercice de probabilité
    Par Bleyblue dans le forum Mathématiques du supérieur
    Réponses: 2
    Dernier message: 28/03/2007, 10h03
  4. Résolution d'un exercice de G.A
    Par GUYSTAVE dans le forum Mathématiques du supérieur
    Réponses: 3
    Dernier message: 16/03/2007, 12h40
  5. Exercice de probabilité
    Par R is R dans le forum Mathématiques du supérieur
    Réponses: 4
    Dernier message: 28/09/2005, 14h53