Discussion :

[Jennn]

Bonjour ,

j'ai un petit soucis au niveau d'un macro programme que j'ai et j'ai besoin de votre aide ...

Dans une base de données j'ai créé une variable réunissant deux variables de type numérique, donc mes conso A et conso B sont réunies en "conso".
Grâce à cette macro (voir ci dessous) je tire aléatoirement sans remise nb individus de ma variable conso :

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
%macro tasr(libref=,entree=,sortie=,nb=,nbsim=);
%do i=1 %to &nbsim;
data t&libref. (drop= i j count);
	count=0;
	array obsnum(&nb) _temporary_;
   do i=1 to &nb;
   redo:
     select=ceil(ranuni(0)*n);
	 set &libref..&entree point=select nobs=n;
	 do j=1 to count;
	 if obsnum(j)=select then goto redo;
	 end;
	 position=select;
	 count=count+1;
	 obsnum(count)=select;
     output;
   end;
   stop;
   set &libref..&entree;
%tasr(libref=work,entree=reech,sortie=donnees2,nb=36,nbsim=10);

Je précise que j'ai trouvé cette macro sur le net. Bon alors il me sélectionne bien 36 données SR. Simplement maintenant j'aimerai qu'ils me les stockent dans une nouvelle variable et que les 36 autres valeurs non sélectionnées soient stockées dans une autre variable.. J'espere être clair

Pouvez vous m'aider s'il vous plait ? Je tiens également a préciser que je suis débutante en SAS

En vous remerciant ,

[Fabien G.]

Bonjour,
ta macro n'est pas complète visiblement.

En partant depuis le début, tu voudrais à partir d'une variable de n observations, faire x tirages aléatoire sans remise pour obtenir :
1 - une table avec tes x observations
2 - une table avec tes n-x observations

C'est bien ça ?

[Fabien G.]

Si c'est ça je te propose la macro suivante, c'est franchement pas top mais c'est un début :

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
%macro tasr(lib,table,nbtasr);
 
* Nb observations de ta table;
data _null_; set &lib..&table.; call symput('nb',compress(_n_)); run;
 
* Teste que le nombre d observations est strictement supérieur au nombre de tirages;
%if &nbtasr.>=&nb. %then %put La table contient &nb. observations, vous devez faire un nombre de tirages strictement inférieur.;
%else %do;
 
* Crée les tables "reste" et "tasr" dans la work;
data reste; set &lib..&table.; run; data tasr; set _null_; run;
 
* tire une observation aléatoire par étape jusqu a nbtasr;
%do i=&nb.-1 %to &nb.-&nbtasr. %by -1;
data _null_; call symput('v',int(ranuni(0)*&i.)+1); run;
data reste tasr_prov; set reste;if _n_=&v. then output tasr_prov; else output reste; run;
data tasr; set tasr tasr_prov; run;
%end;
 
%end;
%mend;
%tasr(sashelp,class,18);

[datametric]

Ca ne serait pas plus simple d'utiliser la proc SURVEYSELECT ?

[Jennn]

Merci pour vos réponses !!!
la PROC SURVEYSELECT j'y ai pensé mais est ce réellement un tirage aléatoire sans remise ? je n'en étais pas certaine ...

[datametric]

c'est LA proc pour faire des tirages avec (with replacement) ou sans remise (without replacement)

http://support.sas.com/documentation...ct_sect007.htm

C'est une procédure a privilégier car elle gère correctement le tirage aléatoire. Lorsque l'on fait un tirage avec une étape DATA, je me souviens de cas contestables.
Je pourrais retrouver de très très vieilles discussions sur ce sujet mais je ne jouerai pas au vieux con ce soir.

[olivier.decourt]

J'ajoute que dans SURVEYSELECT METHOD=URS (tirage avec remise) on doit utiliser de préférence l'option OUTHITS pour avoir 2 observations si le même individu est tiré 2 fois ; et que l'option OUTALL permet d'avoir une variable indiquant la sélection ou non d'un individu, ce qui permet ensuite de scinder facilement le résultat en une table de sélectionnés et une table de non-sélectionnés.

[Jennn]

Ok merci pour vos réponses, je vais me rabattre sur la PROC SURVEYSELECT mais le problème reste le même : un tirage aléatoire sans remise va être fait mais comment stocker les valeurs restantes dans une autre variables ?

[olivier.decourt]

un tirage aléatoire sans remise va être fait mais comment stocker les valeurs restantes dans une autre variables ?

l'option OUTALL permet d'avoir une variable indiquant la sélection ou non d'un individu

Code :

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

1
2
3
4
PROC SURVEYSELECT DATA=sashelp.class OUT=work.echantillon 
                METHOD=URS SEED=1 SAMPSIZE=15
                OUTALL OUTHITS NOPRINT ;
RUN ;

Tu constateras que la table en sortie contient les individus sélectionnés ET ceux qui ne le sont pas. La variable SELECTED permet de faire le tri entre les deux.

[Jennn]

Je reviens sur cette procédure que j'ai un peu laissé tomber quelques temps...
La méthode URS est bien sauf que c'est avec remise, moi il me faut du sans remise absolument. sinon tout fonctionne très bien je vous remercie !

[datametric]

tu es allé voir le lien que j'ai mis vers la doc ?

Si URS est un tirage with replacement, tu n'as pas trouvé la méthode without replacement ?

[Jennn]

Oui oui j'ai été voir bien sur mais il n'y a pas réellement d'équivalent ?

[datametric]

d'équivalent à quoi ?
il y a le tirage avec ou sans remise simple.

[Jennn]

Oui la méthode SRS , seulement si j'utilise cette méthode je ne peux plus mettre les options telles que OUTHITS pour récupérer mes valeurs non sélectionnées...

[datametric]

SRS te permet de savoir qui est sélectionné une fois grâce à OUTALL et non pas OUTHITS.
OUTHITS te dit combien de fois la ligne est prise puisque le URS est avec remise.

[Jennn]

Oui sauf que mes variables non utilisées je les sélectionne de cette façon pour les placer dans une nouvelle variable :

if NumberHits = 0 then consoa = consos;
else if NumberHits NE 0 then consob = consos;

et ça c'est bien avec l'option OUTHITS que j'obtiens cette variable NumberHits...

[datametric]

certes mais ton numberhits est le nombre de fois où l'observation est sélectionnée. Je ne vois donc pas le rapport avec un tirage aléatoire et ton IF. Qu'as tu en ligne ? des individus ou des variables ?

[Jennn]

En ligne j'ai des individus. Ceux qui ont été tirés au sort je les place dans consob, et ceux qui n'ont jamais été tirés au sort forment consoa.

et ça je le fais à chaque tirage aléatoire...