/******************************************************************************
 * Programme : Analyse non-paramétrique des revenus des candidats
 * Reference : ANALYS3B63
 * Source    : https://www.wearecas.eu/fr/sampleCode/ANALYS3B63
 ******************************************************************************/

/* --- BLOC 1 --- */
%web_drop_table(WORK.repincome);

FILENAME REFFILE '/home/myFolder/STAT1/RepIncome..csv';

PROC IMPORT DATAFILE=REFFILE
	DBMS=CSV
	OUT=WORK.repincome;
	GETNAMES=YES;
RUN;

PROC CONTENTS DATA=WORK.repincome; RUN;

%web_open_table(WORK.repincome);

/* --- BLOC 2 --- */
proc sgplot data=work.repincome;
vbox income / category=candidate;
run;

proc sgplot data=work.repincome;
scatter y=income x=candidate;
run;

/* --- BLOC 3 --- */
*Trump Carson comparison;
DATA TrumpCarson; SET work.repincome; 
  IF code NE 3;
RUN;

proc npar1way data=TrumpCarson wilcoxon alpha=.05;
var income;
class candidate;
exact wilcoxon HL;
title 'Trump Carson';
run;
title;

/* --- BLOC 4 --- */
*Trump Cruz Comparison;
DATA TrumpCruz; SET work.repincome; 
  IF code NE 2;
RUN;

proc npar1way data=TrumpCruz wilcoxon alpha=.05;
var income;
class candidate;
exact wilcoxon HL;
title 'Trump Cruz';
run;
title;

/* --- BLOC 5 --- */
*Cruz Carson comparison;
DATA CruzCarson; SET work.repincome; 
  IF code NE 1;
RUN;

proc npar1way data=CruzCarson wilcoxon alpha=.05;
var income;
class candidate;
exact wilcoxon HL;
title 'Cruz Carson';
run;
title;