data dem;
   infile 'Z:\Vicari-Lab\Esercizi\tassi_dem.txt'  pad; *obs=10;
   input nat mort minf paese $22.;
run;
proc print data=dem noobs;
var paese nat; 
run;
data sample;
set dem nobs=N; /* N = variabile temporanea*/
x=nat*N;
x1= ranuni(1); /* (1) specificazione del valore da cui deve partire; random di uniforme*/
run;
proc print data=sample;
run;
data sample;
retain K 4 N; /* creo una variabile k di valore 4 e una variabile N*/
set dem nobs=N;
x=ranuni(1);
y=K/N;
run;
proc print data=sample;
run;

data sample1; /* campione casuale senza ripetizioni*/
set dem;
x=ranuni(0);
if x<=0.20; /* e sottinteso then output*/
run;
proc print data=sample1;run;

data sample2;/* selezione di un campione casuale senza ripetizioni contenente esattamente 10 ripetizione*/
set dem;
x=ranuni(0);
run;
proc sort data=sample2;
by x;
run;
data sample2;
set sample2;
if _N_<=10;/* _ _ rappresenta il contatore delle osservazioni*/
run;
proc print data=sample2;
run;

data sample3;
set dem nobs=N;
y=ranuni(1)*N;
x=int(y)+1; /* int() funzione interna che prende la parte intera del valore indicato ;
            generazione di un numero da 1 a N*/
run;
proc sort data=sample3;by x; run;
proc print data=sample3;
run;

/* selezione di un campione casuale con ripetizione contenente esattamente 50 osservazioni*/
data sample5;
x= int(ranuni(1)*N)+1;
set dem point=x nobs=N; /* point=puntatore= accesso diretto a una certa osservazione*/
i+1;/* i=contatore; valore default e 0;*/
if i > 50 then stop;
drop i;
run;
proc print data=sample5;run;

/* verifica di ipotesi; si puo fare se il campione e con ripetizione*/
proc means data=dem mean;
var mort;
run;
/* normalita*/
proc univariate data=dem plot normal; /* con normal vediamo se l'andamento della mortalita puo essere considerato normale*/
var mort;
run;

proc univariate data=sample5 mu0=10.83 loccount plot; /* loccount considera i quantili*/
var mort;
run;

proc univariate data=sample5 normal plot mu0=10.83 ;
 var mort ;
 histogram mort;
 histogram mort/ midpercents endpoints=70 to 130 by 5;
 run;
proc univariate data=sample5  noprint;
histogram mort/ normal (mu=10.83, sigma=4.647);
 inset normal(ksd ksdpva1);
 run;