[Scilab] Simulation de la loi géométrique

[melbachir]

Bonjour à tous,

Je veux simuler la loi géométrique à l'aide de Scilab. Donc j'écris :

Code : Tout sélectionner

//- Loi géométrique de paramètre p
function X = geo(p)
  X = 1;
  while %t
  //- Si échec +1, si Succès sort
    if(rand() < p) then break; else, X = X + 1;  end 
  end  
endfunction

Cette fonction si je l'ai bien codée, nous donne le rang du premier succès. Mon problème est celui là :
Comment puis-je représenter P(X=k) en fonction de k.

Merci d'avance,

[guiguiche]

Il faut répéter ta fonction un grand nombre de fois et faire compter la fréquence d'apparition de k (plus le nombre est grand plus sa probabilité d'être proche de la P(X=k) est importante).

[melbachir]

Merci pour ta réponse,

Alors ce qui me vient à l'esprit :

Code : Tout sélectionner

//- Taille de l'échantillon
  N = 10000;
//- Le paramètre p
  p = 0.3;
//- Créer un vecteur colonne de composantes nulles
  Z = zeros(N,1);
//- Remplir ce vecteur
  for k=1:N
     Z(k) = geo(p);
  end

Mais après je ne sais pas quoi faire pour obtenir P(X=k) en fonction de k :(

[guiguiche]

J'aurai mis un truc du genre (je ne connais pas la syntaxe de scilab) :

Code : Tout sélectionner

k=geo(p) ;
Z(k)=Z(k)+1 ;

[melbachir]

J'aurai mis un truc du genre (je ne connais pas la syntaxe de scilab) :

Code : Tout sélectionner

k=geo(p) ;
Z(k)=Z(k)+1 ;

Mon gros problème, c'est que j'obtiens sur l'axe des ordonnées des nombres > 1 alors que je m'attends à voir une probabilité...

[guiguiche]

Quelle est la définition d'une fréquence ?