sas_sampling

"http://cmq.esalq.usp.br/wiki/lib/exe/fetch.php?media=publico:syllabvs:lcf5759:dados_grupo01_amazonia.xls"

FILENAME website url "http://www2.census.gov/acs2005/Tables_Profiles_Subject_Tables/010Nation/United%20States.xls"
debug;
website DBMS=xls REPLACE;
SHEET="Base Tables";
GETNAMES=YES;
RUN;

PROC IMPORT OUT= Readin DATAFILE=
ods rtf file='ex3edwards.rtf';
options center nodate nonumber;

/*Importação dos dados */

PROC IMPORT OUT= WORK.A1
DATAFILE= "T:\dados_grupo01_amazonia.xls"
DBMS=EXCEL REPLACE;
SHEET="DADOS$";
GETNAMES=YES;
MIXED=NO;
SCANTEXT=YES;
USADATE=YES;
SCANTIME=YES;
RUN;

PROC PRINT DATA=A1;RUN;

/*Testar Normalidade para todas as plantas*/

PROC UNIVARIATE DATA=A1 NORMAL;
VAR ALTURA_M DAP_CM;
histogram / normal(color=red fill mu=est sigma=est);
inset n = 'Sample Size'
         normal / pos=ne  cfill=blank;
RUN;

/*Amostragem ao acaso*/

PROC SURVEYSELECT DATA=A1 METHOD=SRS N=5 OUT=A5;
PROC PRINT DATA=A5;
RUN;

PROC SURVEYSELECT DATA=A1 METHOD=SRS N=10 OUT=A10;
PROC PRINT DATA=A10;
RUN;

PROC SURVEYSELECT DATA=A1 METHOD=SRS N=15 OUT=A15;
PROC PRINT DATA=A15;
RUN;

title 'Enhancing a Histogram';
PROC UNIVARIATE DATA=A5 NORMAL;
VAR ALTURA_M DAP_CM;
histogram / normal(color=red fill mu=est sigma=est);
inset n = 'Sample Size'
         title="Amostragem aleatorio"
         title2="Midpoints"
         normal / pos=ne  cfill=blank;
RUN;

PROC UNIVARIATE DATA=A10 NORMAL;
VAR ALTURA_M DAP_CM;
histogram / normal(color=red fill mu=est sigma=est);
inset n = 'Sample Size'
         normal / pos=ne  cfill=blank;  /*normal (ksdpval)*/
RUN;

PROC UNIVARIATE DATA=A15 NORMAL;
VAR ALTURA_M DAP_CM;
histogram / normal(color=red fill mu=est sigma=est);
inset n = 'Sample Size'
         normal / pos=ne  cfill=blank;
RUN;

/*Amostragem sistemática*/

PROC SURVEYSELECT DATA=A1 METHOD=SYS N=5 OUT=S5;
PROC PRINT DATA=S5;
RUN;

PROC SURVEYSELECT DATA=A1 METHOD=SYS N=10 OUT=S10;
PROC PRINT DATA=S10;
RUN;

PROC SURVEYSELECT DATA=A1 METHOD=SYS N=15 OUT=S15;
PROC PRINT DATA=S15;
RUN;

PROC UNIVARIATE DATA=S5 NORMAL;
VAR ALTURA_M DAP_CM;
histogram / normal(color=red fill mu=est sigma=est);
inset n = 'Sample Size'
         normal / pos=ne  cfill=blank;
RUN;

PROC UNIVARIATE DATA=S10 NORMAL;
VAR ALTURA_M DAP_CM;
histogram / normal(color=red fill mu=est sigma=est);
inset n = 'Sample Size'
         normal / pos=ne  cfill=blank;
RUN;

PROC UNIVARIATE DATA=S15 NORMAL;
VAR ALTURA_M DAP_CM;
histogram / normal(color=red fill mu=est sigma=est);
inset n = 'Sample Size'
         normal / pos=ne  cfill=blank;
RUN;

ODS rtf CLOSE;