Discussion :

[wsquare]

Bonjour,
Je cherche depuis longtemps.
je voudrais qu'une variable prenne les valeurs 0 ou 1 sachant que 1 par exemple 12% de chance d'apparaître.

Bonjour,
je cherche depuis longtemps à créer une variable prenant les valeurs 0 ou 1, sachant que 1 a 12% de chance d'apparaître.

je m'explique:

Pour chaque ligne (âge), j'ai une probabilité connue d'occurrence d'un évènement admettons 12% pour la ligne 1 25% pour la ligne 2 etc..
Je veux créer une variable aléatoire qui pour chaque ligne sur l'ensemble de l'échantillon prenne affiche un 0 ou un 1 en fonction de la probabilité renseignée.
Je sais qu'il faut que j'utilise soit un call logistic ou un ranuni, mais comment par exemple le ranuni m'affiche un 1 et sur l'ensemble de l'échantillon l'espérance soit de 12% de rencontrer un 1.


merci

[xav2229]

Tu peux utiliser une loi UNIFORME

X

[wsquare]

oui, je sais, mais comment faire pour que l’occurrence du 1 respect une certaine probabilité?

voici une partie de mon code :

data simperp;
set simulation2;
if finna = 0 then x1=int(ranuni(0)<tauxperp);
run;

[wsquare]

par exemple
j'ai une colonne avec l'âge de l'individu et les probabilités :
à 35 ans un individu a 0.12 de chance que un 1 lui soit affecté.


âge y x1
35 0.12 ?
36 0.5 ?
37 0.06 ?

[xav2229]

comme ceci?

Code :

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

1
2
3
4
5
 
data WORK.ALEAS_2;
   set WORK.ALEAS;
   x = ifn(ranuni(123) <= p, 1, 0);
run;

X

[RemiBousquet]

Bonjour

Tu verras que ce genre de simulation est assez simple avec ranuni()

Dans le code d'exemple, je crée mon "top" de réalisation de l'évènement si mon ranuni() est inférieur à ma proba.
La proc means permet de vérifier le résultat

Code :

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

1
2
3
4
5
6
7
8
9
10
11
12
13
14
data WORK.TEST ;
    do loop = 1 to 1000 ;
        PROBA_EVT = 0.12 ; output ;
    end ; 
run ; 
 
data WORK.SIMULATION;
    set WORK.TEST ;
    TP_EVT = (ranuni(1) < PROBA_EVT ) ;
run ; 
 
proc means data = WORK.SIMULATION ;
    var TP_EVT ;
run ;

A noter que le paramètre du ranuni() est la "racine" du pseudo aléa.
Si tu souhaites pouvoir resimuler de façon identique le même échantillon (trace d'audit), il faut alors choisir une valeur spécifique > 0.