PROC UNIVARIATE : Modifier l'ordre d'affichage des valeurs pour tracer une courbe théorique sur le graph

[avalen]

Bonjour,

J'ai une distribution(voir pièces jointes) qui commence à 0 heures.
J'ai créé une autre distribution (voir pièces jointes) avec une proc gchart qui me convient et qui commence à 5h et termine à 4h.

J'aimerais appliquer sur cette deuxième distribution une proc univariate pour y faire apparaître les courbes théoriques (même si à vue d'oeil c'est une loi normale).

Mon pb : proc univariate n'accepte pas que les heures soient classées par ordre autre que croissant. Donc MIDPOINTS ne marche pas, ORDER non plus...

J'espère que quelqu'un a déjà eu ce problème et pourra donc remédier au mien

Merci d'avance.

PS: Est-il possible d'ajouter des titres aux graphiques dans le cadre d'une univariate??

[MEGAMIND2]

Code :

ODS GRAPHICS ON

te permet d'avoir les graphiques en V9.2.
ça ne coute rien d'assayer!!

[oualina]

Bonjour,
Utilise cette macro SAS, permet de tracer l'histogramme, la courbe du noyau et le QQplot.

Code :

1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
29
30
31
32
33
34
35
36
37
38
39
40
41
42
43
44
45
46
47
48
49
50
51
52
53
54
55
56
57
58
59
60
61
62
63
64
65
66
67
68
69
70
71
72
73
74
75
76
77
%macro normalityplot(DATA=, var=, varlabel=); 
quit;
*ods SELECT NONE;
ods output Moments=moment TestsforNormality=normal;
proc univariate DATA=&DATA normal;
  var &var;
run;
 
DATA ChoixStat; SET normal;
  IF _n_ = 1;
run;
 
DATA size; SET moment;
  WHERE Label1 = 'N';
  n = cValue1; keep n;
run;
 
DATA skew; SET moment;
  WHERE Label1 = 'Skewness';
  skewness = cValue1; keep skewness;
run;
 
DATA kurt; SET moment;
  WHERE Label2 = 'Kurtosis';
  kurtosis = cValue2; keep kurtosis;
run;
 
DATA affiche;
  merge size skew kurt ChoixStat;
  DROP VarName pSign;
run;
 
DATA anno;
  length text $20.4;
  retain xsys '3' ysys '3' FUNCTION 'label' when 'a' y_pct 85
         hsys '4' size .75;
  SET affiche;
  IF _n_ = 1 then do;
    x=23; y=85; text=put(Test, 20.4);   output;
    x=18; y=82; text=put(TestLab, 2.);  output;
    x=28; y=82; text=put(pType, 6.);    output;
  end;
    y_pct=y_pct-3;
   x=18; y=79; text=put(Stat, 6.4);   output;
    x=28; y=79; text=put(pValue, 6.4); output;
 
   y_pct=y_pct-8;
   x=18; y=72; text='Skewness'; output;
   x=28; y=72; text='Kurtosis'; output;
   x=18; y=69; text=put(skewness, 8.4); output;
    x=28; y=69; text=put(kurtosis, 8.4); output;
 
run;
 
 
*** Probability Plots AND Histograms***;                                                                                                              
goptions ftext=SWISS ctext=BLACK htext=1 cells;                                                                                         
symbol v=star c=BLUE h=0.5 cells;                                                                                                       
axis1 label=("&varlabel");
 
proc univariate DATA=&DATA noprint;                                                                                               
   var &var; 
   title "Verification de l'hypothese de normalite &varlabel";
   probplot / caxes=BLACK cframe=CXF7E1C2 waxis= 1                                                                                      
              hminor=0 vminor=0 vaxislabel="&varlabel"                                                                                       
              normal( mu=est sigma=est color=BLUE l=1 w=1);                                                                                                             
   inset normal;  
   histogram / normal (noprint color=red)
               caxes=BLACK cframe=CXF7E1C2 waxis= 1                                                                                     
               cbarline=BLACK cfill=BLUE pfill=SOLID                                                                                    
               vscale=percent hminor=0 vminor=0 anno=anno nmidpoints=11;                               
run;                                                                                                                                    
   title;
ods SELECT ALL;
%mend;
 
ods SELECT ALL;

[avalen]

Mais qu'est ce que c'est que cela mon dieu!!!
J'ai testé, je suppose que je n'avais qu'à remplir les espaces vides.
J'ai plus ou moins compris le programme mais je ne suis jamais passée par les macros donc :S

Par contre le fichier fait 0KB et le message d'erreur suivant apparaît :

Aucune action par défaut pour le type de données Macro

[olivier.decourt]

Je proposerais de créer tout simplement une nouvelle variable.
Si actuellement heure >= 5 alors nouvelle_heure = heure-5
Sinon nouvelle_heure = heure + 19
Ou un truc du genre.
Olivier