sas_reshaping wide_long

/*PROC IMPORT OUT= WORK.PF */
/*            DATAFILE= "T:\Profundidades.xlsx" */
/*            DBMS=EXCEL REPLACE;*/
/*     RANGE="Plan1$"; */
/*     GETNAMES=YES;*/
/*     MIXED=NO;*/
/*     SCANTEXT=YES;*/
/*     USEDATE=YES;*/
/*     SCANTIME=YES;*/
/*RUN;*/

DATA PF;
INPUT PONTO	P20	P150 PONTO1 P201 P1501;
DATALINES;
1	12	6	21	3	0
2	9	4	22	9	1
3	14	5	23	13	11
4	10	7	24	11	3
5	11	6	25	12	10
6	13	7	26	16	6
7	14	8	27	9	10
8	13	7	28	6	0
9	11	13	29	10	7
10	12	7	30	7	3
11	11	10	31	11	11
12	11	5	32	14	6
13	13	10	33	10	11
14	13	8	34	13	14
15	11	5	35	15	17
16	8	11	36	3	0
17	7	7	37	13	10
18	11	3	38	10	0
19	13	4	39	2	0
20	7	0	40	14	9
;
PROC PRINT DATA=PF; RUN;

DATA A;
SET PF;
DROP Ponto1 P201 P1501; OUTPUT;
RUN;
PROC PRINT DATA=A; RUN;

DATA B;
SET PF;
DROP Ponto P20 P150; OUTPUT;
RENAME PONTO1=PONTO;
RENAME P201=P20;
RENAME P1501=P150;
RUN;
PROC PRINT DATA=B; RUN;

DATA C ;
SET A B;
DIF=P20-P150;
PROC PRINT DATA=C; RUN;

PROC UNIVARIATE DATA=C NORMAL ;
VAR DIF ;
RUN ;
/*goptions reset=all border cback=white htitle=12pt htext=10pt;*/

DATA C;
do PONTO=1 to 40;
P20=round(ranuni (20);
P150=round(ranuni(20);
output;  end;
run;

symbol1 interpol=join value=squarefilled   color=vibg height=1;
symbol2 interpol=join value=trianglefilled color=depk height=1;
legend1 label=none frame;
axis1 label=("Pontos") minor=none offset=(1,1);
axis2 label=(angle=90 "No de organismos por litro de água")
      order=(0 to 20 by 2) minor=(n=1);

proc gplot data=C;
   plot (P20 P150)*PONTO / overlay legend=legend1
                                   haxis=axis1 vaxis=axis2;
run;
quit;