*****************************************************;
* Traitement 3 - Modele Méthode Van Westendorp;
* Version 2.0 ;
* Date : 5/5/2009 ;
* remarques :   ;
*****************************************************;


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;


data in; 
  set in0 ;
array TC{160} TC001-TC160 ;
array C{160}  C001-C160 ;
array E{160}  E001-E160 ;
array TE{160} TE001-TE160 ;

array M{160}  M001-M160 ;

do i=1 to TC0-1  ; TC{i}=1 ; end ; 
do i=TC0 to 160  ; TC{i}=0 ; end ; 

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

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

do i=1 to TE0-1  ; TE{i}=0 ; end ; 
do i=TE0 to 160  ; TE{i}=1 ; end ; 

Do i=1 to 160 ;
if i<TC0+.01 then M{i}=1; 
 else if i<C0+.01 then M{i}=2; 
  else if i<E0+.01 then M{i}=3;
   else if i<TE0+.01 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;
 keep quest prix log_prix TC0 C0 E0 TE0  ;

 array TC{160} TC001-TC160 ;
 array C{160}  C001-C160 ;
 array E{160}  E001-E160 ;
 array TE{160} TE001-TE160 ;

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

proc sort data=out2;
 	by prix;
Proc logistic data=out2 ;
 	model TC0 (ref="0") = log_prix   / link=logit ;
 	output out=out_log1 p=prevu1;
Proc logistic data=out2 ;
 	model C0 (ref="0") = log_prix   / link=logit ;
 	output out=out_log2 p=prevu2;
 Proc logistic data=out2 ;
 model E0 (ref="0") = log_prix   / link=logit ;
 output out=out_log3 p=prevu3;
 Proc logistic data=out2 ;
 model TE0 (ref="0") = log_prix   / link=logit ;
 output out=out_log4 p=prevu4;

data tout ;
 merge out_log1 out_log2 out_log3 out_log4;
 by prix ;

 BR= TC0-C0;
 OK= C0-E0;
 PR= E0-TE0 ;
 Prix_AcceptR  = TC0+ (1-TE0)-1;
 Prix_AcceptT  = prevu1+ (1-prevu4)-1;

 bargain       = prevu1-prevu2 ;
 acceptable    = prevu2-prevu3 ;
 premium       = prevu3-prevu4;

 run;
proc sort data=tout;
 by prix;
title2 "Data out4";
proc means data=tout noprint ;
 by prix ;
 var  TC0 C0 E0 TE0 BR OK PR  prevu1 prevu2 prevu3 prevu4 bargain acceptable premium Prix_AcceptR Prix_AcceptT;
 output out= out4 mean= TC C E TE BR OK PR  prevu1 prevu2 prevu3 prevu4 bargain acceptable premium Prix_AcceptR Prix_AcceptT;

 data out4; 
	set out4 ;
 QTC=1-TC ; QC=1-C ; QE=1-E ; QTE=1-TE ;
 BR= TC-C;
 OK= C-E;
 PR= E-TE ;
title2 "Fig 1 : Courbes TC C E TE reelles ";
proc gplot data=out4 ;
 plot (TC C E TE )*prix/overlay legend ;
 run;
title2 "Fig 2 : QTC C QE TE";
proc gplot data=out4 ;
 plot (QTC C QE TE)*prix/overlay ;
 run;

legend1 value=("too cheap" "bargain" "acceptable" "premium" "too expensive");

 title2 "Fig 3 : Courbes TC C E TE  prevues";
proc gplot data=out4 ;
 plot (prevu1 prevu2 prevu3 prevu4)*prix/overlay  legend ;
 run;
 title2 "Fig 4 : Courbes TC C E TE  prevues et reelles";
proc gplot data=out4 ;
 plot (TC C E TE prevu1 prevu2 prevu3 prevu4)*prix/overlay legend ;
 run;

 title2 "Fig 5 : Bargain Acceptable Premium ";
 proc gplot data=out4 ;
 plot (BR OK PR  bargain acceptable premium)*prix/overlay  legend ;
 run;
 title2 "Fig 5-1 : Bargain (entre TC < < C ";
 proc gplot data=out4 ;
 plot (BR  bargain )*prix/overlay  legend ;
 run;
 title2 "Fig 5-2 : Acceptable entre C< < E";
 proc gplot data=out4 ;
 plot ( OK  acceptable )*prix/overlay  legend ;
 run;
 title2 "Fig 5-3 :  Premium entre E< < TE";
  proc gplot data=out4 ;
 plot (PR  premium)*prix/overlay  legend ;
 run;

title2 "Fig 6 : Prix Acceptable théorique et réel entre TC < < TE";
 proc gplot data=out4 ;
 plot (Prix_AcceptR Prix_AcceptT)*prix/overlay  legend ;
 run;