MACRO pour générer échantillon bootstrap

[glaconutella]

Bonjour,

Voilà une macro permettant d'effectuer un échantillon bootstrap ("out") à partir d'un jeu de données ("in").
J'ai trouvé cette macro (qui fonctionne très bien) sur :
http://www.math.univ-toulouse.fr/~be...s/tp6ssbvc.pdf

Code :

1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
 
/* Macro générant des échantillons bootstrap */
 
%macro bootstrp(IN,out);
DATA seed (keep=seed) &out (DROP=seed i);
IF _n_=1 then SET seed;
do i=1 TO ntot;
  call ranuni(seed,num);
  num=floor(ntot*num+1);
  SET &IN point=num nobs=ntot;
  output &out;
end;
output seed;
stop;
run;
%mend;
 
DATA seed;
retain seed 13;
run;

Malheureusement, étant novice en la matière , j'ai quelques difficultés à comprendre :
- à quoi correspond "seed" qu'on appelle la graine.
-pourquoi la graine est égale à 13 et et quelle valeur choisir?
data seed;
retain seed 13;
run;
- d'où vient la variable point dans "set &in point=num nobs=ntot;"
- pourquoi faire un "stop;"

Merci d'avance pour vos réponses.

[la2002]

Salut Glaconutella,

Le seed permet de générer tes nombre (pseudo) aléatoires.
En choisissant un seed fixe, tu auras les mêmes résultats entre 2 exécutions de ton programme.
je te conseille de prendre une autre valeur (Ex le nombre de secondes écoulés depuis le 1/1/2010) pour éviter cet effet. Ainsi 2 exécutions de ton programme ne donnera pas les mêmes résultats.


Le point et le stop vont ensemble. Le point permet l'accession directe à la num ème observation (dixit ton programme)

Stop permet de quitter l'étape data. sinon, il boucle sans cesse sur la num ème observation.

Ces instructions sont plus performantes que un where=(_n_=num)

[datametric]

vas voir la proc surveyselect qui fait du bootstrap toute seule avec REPS=

[RemiBousquet]

Tu peux trouver plein d'infos et de programmes sur le site d'Olivier DECOURT
http://www.od-datamining.com/

[datametric]

J'en avais parlé avec Olivier cet été et nous avions conclu que REPS= était la meilleur stratégie pour le bootstrap