Discussion :

[celia89]

Bonjour,

je cherche comment réaliser un tirage aléatoire de nombre selon une loi uniforme sans remise. Pour faire simple, disons que j'ai une table avec x lignes et pour chaque ligne une valeur min et une valeur max. ex :

table =

id min max
1 1 10
2 1 15
3 1 8 ...

J'aimerais pour chaque ligne tirer 3 chiffres DIFFERENTS compris entre le min et max de la ligne. Avec la fonction ranuni j'y arrive mais il peut me tirer plusieurs fois le meme chiffre pour une ligne donnée. J'utilise le code ci-dessous :

Code :

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

1
2
3
4
5
6
7
data tirage (drop=i);
set table;
   do i=1 to 3;
     tirage=int(ranuni(45123)*max)+1;
     output;
   end;
 run;

qui marche très bien, cela me duplique mes lignes 3 fois avec un chiffre compris entre 1 et le max mais il n'est pas sans remise.

Je n'ai rien trouver sur internet pour faire cela...

Merci de votre aide !

[Manoutz]

Une solution en remodelant legerement le dataset en sortie:

Code :

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

1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
data table;
input id min max 3.;
cards;
1 1 10
2 1 15
3 1 8
;
run;
 
 
 
data tirage ;
set table;
tirage1=-999;
tirage2=-999;
tirage3=-999;
do until (tirage1 ne tirage2 and  tirage1 ne tirage3 and  tirage2 ne tirage3) ;
tirage1=int(ranuni(-3)*(max-min))+min;
tirage2=int(ranuni(-3)*(max-min))+min;
tirage3=int(ranuni(-3)*(max-min))+min;
end;
run;

[olivier.decourt]

Je proposerais une solution en 2 temps ; d'abord la démultiplication de la table (une ligne = une combinaison id x nombre) puis un tirage aléatoire sans remise avec la proc Surveyselect (avec un STRATA id et un SAMPSIZE=3).
La démultiplication :

Code :

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

1
2
3
4
5
6
DATA mult ;
  SET table ;
  DO i=min TO max ;
    OUTPUT ;
  END ;
RUN ;

Le tirage

Code :

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

1
2
3
PROC SURVEYSELECT DATA=mult OUT=ech SAMPSIZE=3 SEED=1 NOPRINT ;
  STRATA id ;
RUN ;