data in0 ;

input quest TC0 C0 E0 TE0 ;

cards ;
1 10 25 35 40
2 11 26 36 45
3 13 30 38 50
4 14 32 40 55
5 18 34 41 56
6 20 35 42 58
7 23 38 45 60
8 25 40 48 65
9 26 41 49 70
10 27 42 50 75
11 28 44 55 80
12 29 45 60 85
;
symbol1 i=join L=1;
symbol2 i=join L=1;
symbol3 i=join L=1;
symbol4 i=join L=1;
symbol5 i=join L=1;
proc print data=in0 ;
run;
proc tabulate data=in0;
var TC0 C0 E0 TE0;
table (TC0 C0 E0 TE0), mean;
run ;

 data in; 
 	set in0 ;
array TC{100} TC001-TC100 ;
array C{100}  C001-C100 ;
array E{100}  E001-E100 ;
array TE{100} TE001-TE100 ;

array M{100}  M001-M100 ;

do i=1 to TC0  ; TC{i}=1 ; end ; 
do i=TC0+1 to 100  ; TC{i}=0 ; end ; 

do i=1 to C0-1  ; C{i}=0 ; end ; 
do i=C0 to 100  ; C{i}=1 ; end ; 

do i=1 to E0-1  ; E{i}=0 ; end ; 
do i=E0 to 100  ; E{i}=1 ; end ; 

do i=1 to TE0  ; TE{i}=0 ; end ; 
do i=TE0+1 to 100  ; TE{i}=1 ; end ; 

Do i=1 to 100 ;
if i<TC0+1 then M{i}=1; 
	else if i<C0+1 then M{i}=2; 
		else if i<E0+1 then M{i}=3;
			else if i<TE0+1 then M{i}=4;
			else M{i}=5;
end;

output ;

*******************************************************;
* DATA OUT2 ;
* sortie des variables par individu  x 100 lignes ;
*******************************************************;
title2 "OUT2 variables par individu  x 100 lignes " ;
data out2 ;
	set in;
	array TC{100} TC001-TC100 ;
	array C{100}  C001-C100 ;
	array E{100}  E001-E100 ;
	array TE{100} TE001-TE100 ;
	array QE{100} QE001-QE100 ;
	array QTC{100} QTC001-QTC100 ;

	do i=1 to 100 ;
	prix=i ;
	QTC{i}=1-TC{i} ;
	QE{i}=1-E{i} ;
		output ;
	end;
	run;



data out2 ;
	set in;
	keep quest prix log_prix TC0 C0 E0 TE0 ACC BR OK PR QTC QE;

	array TC{100} TC001-TC100 ;
	array C{100}  C001-C100 ;
	array E{100}  E001-E100 ;
	array TE{100} TE001-TE100 ;

	do i=1 to 100 ;
	prix=i ;
	log_prix=log(i) ;
	TC0= 1- TC{i};C0=C{i};E0=E{i};TE0=TE{i};

	BR= TC0-C0;
	OK= C0-E0;
	PR= E0-TE0 ;

	ACC=TC0-TE0 ;

	QTC=1-TC0 ;
	QE=1-E0;

output ;

	end;
proc print data=out2 ;
	where quest<2 ;
	var quest prix TC0 C0 E0 TE0 ;

proc sort data=out2;
	by prix;
title2 "Data out4";
proc means data=out2 noprint ;
	by prix ;
	var  TC0 C0 E0 TE0 BR OK PR ACC  QTC QE;
	output out= out4 mean= TC C E TE BR OK PR ACC QTC QE;
proc gplot data=out4 ;
	plot (TC C E TE )*prix/overlay legend ;
	run;
	proc gplot data=out4 ;
	plot (QTC C QE TE )*prix/overlay legend ;
	run;
title2 "Fig 5 :  bargain correct premium";
proc gplot data=out4 ;
	plot ( BR OK PR  )*prix/overlay  legend ;
	run;
	title2 "Fig 6 : Volumes (quantités) à partir de l'acceptabilité (Hyp QA/NA homogène)";
proc gplot data=out4 ;
	plot (ACC )*prix/overlay  legend ;
	run;
data out4; set out4 ;
	CA= prix*acc;
	title2 "Fig 7 : CA à partir de l'acceptabilité ";
proc gplot data=out4 ;
	plot (CA )*prix/overlay  legend ;
	run;

	/* pour détermination automatique des points d'intersection / A FINIR
			if i>1 then do ;
			if (QTC{i}<C{i} and QTC{i-1}>C{i-1}) or (QTC{i}>C{i} and QTC{i-1}<C{i-1}) then MGP= i;
			if (TE{i}<QE{i} and TE{i-1}>QE{i-1}) or (TE{i}>QE{i} and TE{i-1}>QE{i-1})then MDP= i;
			if (C{i}<QE{i} and C{i-1}>QE{i-1}) or (C{i}>QE{i} and C{i-1}<QE{i-1})  then IGP= i;
			if (QTC{i}<TE{i} and QTC{i-1}>TE{i-1}) or (QTC{i}>TE{i} and QTC{i-1}<TE{i-1}) then OPS= i;
		end ;
	proc print ;
	var MGP MDP IGP OPS ;
	run;
	*/