I'm trying to come out a series of normal distribution and tested by Shapiro-Wilks test. The code is as below: DATA RAND; DO J=0 TO 9; CALL STREAMINIT(123); DO i=1 to 50; X=Rand("normal",1.5+0.001*j,0.3); OUTPUT; END; END; RUN; /************************************************* *input = dataset to check; *vars = Specify variable(s) for which you want to check normality; *output = dataset to output with normality status; **************************************************/ %macro normal(input=, vars=, output=); ods output TestsForNormality = Normal; proc univariate data = &input normal; var &vars; run; ods output close; data &output; set Normal ( where = (Test = 'Shapiro-Wilk')); if pValue > 0.05 then Status ="Normal"; else Status = "Non-normal"; drop TestLab Stat pType pSign; run; %mend; %normal(input=RAND, vars=X, output=Normality); But the output data set is only one item by one Varname X from 500 data, not for j between 0 and 9 and from 50 data. Thank you if you can support. ... View more

Dear all, I'm writing a code by sas university edition (9.4) to come out a dataset with column i, j, X, Y. The code is as below: data work.Rand; %macro Normal_Simulation; call streaminit(567); do i=1 to 50; X=Rand("normal",0.02*j-1,1.0); IF X>3 OR X<-3 THEN DO Y=X; X=0; i=i-1; END; ELSE Y=0; /*u=Rand("uniform");*/ output; end; %mend; do j=0 to 99; %Normal_Simulation; end; run; PROC SQL; CREATE TABLE work.query AS SELECT j , i , X , Y FROM work.rand; /*WHERE J=&j AND Y=0;*/ RUN; QUIT; Can I use any sql expression in Where statement here to come out like only the first 10 percentile dataset of X of each j in this table? Thank you.   MKW   ... View more

Dear all, I'm trying to generate a set of random data and export then Numpy in Python environment can read and work on. So I would export the data without var label and learned below and code as attached with label x='00'x, but x is still there, see attached csv data MYSECOND.BIRANDOM; call streaminit(1234); do i=1 to 100; x=rand("normal", 20, 3); /*y=0;*/ output; end; do i=101 to 200; x=rand("normal", 30, 3); output; /*y=0;*/ end; Drop i; run; proc print data=mysecond.birandom; var x; label x='00'x; run; /* Option group 5 (GRAPH SIZE) parameters. */ /*--Set output size--*/ ods graphics / reset imagemap; /*--SGPLOT proc statement--*/ proc sgplot data=mysecond.BIRANDOM; /*--Histogram settings--*/ histogram x / scale=COUNT SHOWBINS; /*--Vertical or Response Axis--*/ yaxis grid; run; /*ods graphics / reset; ods noproctitle; ods graphics / imagemap=on; proc princomp data=mysecond.BIRANDOM plots(only)=(scree); var x y; run;*/ /* Stream a CSV representation of SASHELP.CARS directly to the user's browser. */ proc export data=mysecond.birandom outfile="/folders/myfolders/MySimulation/birandom.csv" dbms=csv replace; label; run; %let _DATAOUT_MIME_TYPE=text/csv; %let _DATAOUT_NAME=birandom.csv; http://support.sas.com/documentation/cdl/en/proc/65145/HTML/default/viewer.htm#n17dcq1elcvpvkn1pkecj41cva6j.htm   Can someone advise where went wrong and how to fix it, thank you! ... View more

Dear all, Following my last trial to output a series of histogram through SQL with different parameters and SAS code as attached. %MACRO SPCDIST; /*Macro file is needed for repetitive SQL statement*/ %do j=1 %to 10; PROC SQL; CREATE TABLE work.query AS SELECT j , i , X , Y FROM MYSECOND.SPC; WHERE J=&j AND Y=0; RUN; QUIT; ods graphics / reset imagemap; /*--SGPLOT proc statement--*/ proc sgplot data=WORK.QUERY; /*--Histogram settings--*/ histogram X /; /*--Vertical or Response Axis--*/ yaxis grid; run; /*end; run;*/ ods graphics / reset; %end; %MEND; %SPCDIST; The result was output 10 identifical plots (like j=10). I wonder where went wrong and please advise.   Thank you! ... View more