```Date: Fri, 16 Jun 2000 11:14:05 -0700 Reply-To: barrere bendia Sender: "SAS(r) Discussion" From: barrere bendia Subject: Re: random number?Thank you all Comments: To: Bill_Knowlton@baxter.com Content-Type: text/plain Thank you all. The SAS-List are good friends. I coldn't resist to ask Bill Knowlton, where can I find this usuful information about RANGAM function. And I would ask about PROC BOX PLOT: I wante plot juste one BOXPLOT for a variable, but more BOXPLOT. Excactly I have data set like this: data Relief1; input Level Code Acupun Relief @@; datalines; 1 1 1 0.0 1 2 1 0.5 1 1 2 0.6 1 2 2 1.2 2 1 1 0.3 2 2 1 0.6 2 1 2 0.7 2 2 2 1.3 3 1 1 0.4 3 2 1 0.8 3 1 2 0.8 3 2 2 1.6 4 1 1 0.4 4 2 1 0.7 4 1 2 0.9 4 2 2 1.5 5 1 1 0.6 5 2 1 1.0 5 1 2 1.5 5 2 2 1.9 6 1 1 0.9 6 2 1 1.4 6 1 2 1.6 6 2 2 2.3 7 1 1 1.0 7 2 1 1.8 7 1 2 1.7 7 2 2 2.1 8 1 1 1.2 8 2 1 1.7 8 1 2 1.6 8 2 2 2.4 ; I wante get juste one BOXPLOT for the variable Relief . I tried some combinaison with PROC BOXPLOT, but not succes. Here some thing that tried to program: proc sort data=Relief1 out=Relief2; by level Acupun ; run; symbol color = salmon h = .8; goptions ftext=swiss; axis1 minor=none color=black label=(angle=90 rotate=0); title 'Here I have 8 Box Plot. How to I have juste one Box Plot for Relief or other? '; proc boxplot data=Relief2; plot Relief*level/ cframe = vligb cboxes = dagr cboxfill = ywh vaxis = axis1; run; The result is 8 Box Plot,. And 16 8 Box Plot. I mean Relief by level and by Acupun . Thanks a other time. > > > Here's some info from the SAS v8 CD Manual on the > > 1. RANGAM function > > and > > 2. RANGAM call routine. > > 'Hope this helps. > > -Bill > > > 1. RANGAM > > Syntax > > > > > RANGAM(seed,a) > > > > > > > Arguments > > > seed > is an integer. For more information on seeds, see Seed Values. > > Range: seed < 2(superscript: 31)-1 > > Note: If seed = 0, the time of day is used to > initialize the seed stream. > > > > > > a > is a numeric shape parameter. > > Range: a > > 0. > > > > > > > Details > > > > > > > The RANGAM function returns a variate that is generated from a gamma > distribution with parameter a. For a > 1, an acceptance-rejection method > due to Cheng (1977) (See References) is used. For a = 1, an > acceptance-rejection method due to Fishman is used (1978, Algorithm G2) > (See References). > > > A gamma variate X with shape parameter ALPHA and scale BETA can be > generated: > x=beta*rangam(seed,alpha); > > > If 2*ALPHA is an integer, a chi-square variate X with 2*ALPHA degrees of > freedom can be generated: > x=2*rangam(seed,alpha); > > > If N is a positive integer, an Erlang variate X can be generated: > x=beta*rangam(seed,N); > > > It has the distribution of the sum of N independent exponential variates > whose means are BETA. > > > And finally, a beta variate X with parameters ALPHA and BETA can be > generated: > y1=rangam(seed,alpha); > y2=rangam(seed,beta); > x=y1/(y1+y2); > > > > > > > > > > > The CALL RANGAM routine, an alternative to the RANGAM function, gives > greater control of the seed and random number streams. > > > 2. CALL RANGAM > > > > Returns a random variate from a gamma distribution > > Category: Random > Number > > > > > > > > Syntax > > > > > CALL RANGAM(seed,a,x > ); > > > > > > Arguments > > > seed > is the seed value. For more information about seeds, see Seed Values. A new > value for seed is returned each time CALL RANGAM is executed. > > Range: seed < 2(superscript: 31) - 1 > > Note: If seed =0, the time of day is used to > initialize the seed stream. > > > > > > a > is a numeric shape parameter. > > Range: a > > 0 > > > > > > x > is a numeric variable. A new value for the random variate x is > returned each time CALL RANGAM is executed. > > > > Details > > > > > > > The CALL RANGAM routine updates seed and returns a variate x that is > generated from a gamma distribution with parameter a. > > > By adjusting the seeds, you can force streams of variates to agree or > disagree for some or all of the observations in the same, or in subsequent, > DATA steps. > > > For a>1, an acceptance-rejection method (Cheng 1977) (See References) is > used. For a =1, an acceptance-rejection method due to (Fishman 1978) (See > References) is used. > > > > Comparisons > > > > > > > The CALL RANGAM routine gives greater control of the seed and random number > streams than does the RANGAM function. > > > > Examples > > > > > > > This example uses the CALL RANGAM routine: > options nodate pageno=1 linesize=80 pagesize=60; > > data case; > retain Seed_1 Seed_2 Seed_3 45; > a=2; > do i=1 to 10; > call rangam(Seed_1,a,X1); > call rangam(Seed_2,a,X2); > X3=rangam(Seed_3,a); > if i=5 then > do; > Seed_2=18; > Seed_3=18; > end; > output; > end; > run; > > proc print; > id i; > var Seed_1-Seed_3 X1-X3; > run; > > > > > > > > > > barrere bendia > VISTA.COM> cc: > Sent by: "SAS(r) Subject: random number? > Discussion" > GA.EDU> > > > 06/15/00 01:40 PM > Please respond to > barrere bendia > > > > > > Hi > I wante to generate number according gamma distribution with mean=mu and > variance= var. > > I tried the functin: RANGAM(seed,a) but what mens the parameter a. The > GAMMA distribution has two parameters. > > Can I do like this: > do i=1 to 100; > gama= mu + sqrt(var)*RANGA(seed,mu); > output; > end; > drop i; > ? > Thank you all. > > > _______________________________________________________________________ > > Why pay when you don't have to? Get AltaVista Free Internet Access now! > http://jump.altavista.com/freeaccess4.go > > _______________________________________________________________________ _______________________________________________________________________ Why pay when you don't have to? Get AltaVista Free Internet Access now! http://jump.altavista.com/freeaccess4.go _______________________________________________________________________ ```

Back to: Top of message | Previous page | Main SAS-L page