Creando distribuciones en SAS

* Curvas de funciones normales de densidad con diferentes medias;

data normal;

do  n = 1 to 101;                 * loop 101 veces;

    x = -5+ (n-1)/10;             * x = -5 to 5;

  n01 = pdf('normal',x,0,1);      * N(0,1);

  n11 = pdf('normal',x,1,1);      * N(1,1);

  n21 = pdf('normal',x,2,1);      * N(2,1);

  n04 = pdf('normal',x,0,2);      * N(0,4);

  n09 = pdf('normal',x,0,3);      * N(0,9);

  output;                         * observation to data set;

end;                              * finaliza loop;

symbol1 interpol=join value=none color=blue;

symbol2 interpol=join value=dot color=green  height=.5;

symbol3 interpol=join value=square color=red height=.5;

axis1 label=("x");

axis2 label=("pdf");

legend1 label=(justify=c "Distribucion" justify=c "Parametros")

value=(tick=1 justify=c "N(0,1)" tick=2 justify=c "N(1,1)" tick=3 justify=c "N(2,1)");

proc gplot data=normal;

       plot n01*x=1n11*x=2n21*x=3/ overlay

       haxis=axis1 vaxis=axis2

       legend=legend1;

       title1 "Funciones normales de densidad: cambiando medias";

run;quit;

run;

* Curvas de funciones normales de densidad con diferentes varianzas;

data normal;

do   n = 1 to 101;                * loop 101 veces;

     x = -5+ (n-1)/10;            * x = -5 to 5;

   n01 = pdf('normal',x,0,1);     * N(0,1);

   n11 = pdf('normal',x,1,1);     * N(1,1);

   n21 = pdf('normal',x,2,1);     * N(2,1);

   n04 = pdf('normal',x,0,2);     * N(0,4);

   n09 = pdf('normal',x,0,3);     * N(0,9);

   output;                              

end;                              * finaliza loop;

symbol1 interpol=join value=none color=blue;

symbol2 interpol=join value=dot color=green  height=.5;

symbol3 interpol=join value=square color=red height=.5;

axis1 label=("x");

axis2 label=("pdf");

legend2 label=(justify=c "Distribucion" justify=c "Parametros")

value=(tick=1 justify=c "N(0,1)" tick=2 justify=c "N(0,4)" tick=3 justify=c "N(0,9)");

proc gplot data=normal;

       plot n01*x=1n04*x=2n09*x=3/ overlay

       haxis=axis1 vaxis=axis2

       legend=legend2;

       title1 "Funciones normales de densidad: cambiando varianzas";

run; quit;