DATA A;
   INFILE 'd:\data\sitkin\study2\sitkin2.csv' DSD;
   INPUT con$ q11-q16 q21-q26 q31-q36 q41-q46 q51-q56 q61-q66 q71-q76 q81-q86;
   *****
   The database is according to questionnaire A
   qij - i -the prospect number, j- the question number. Thus q12 is the answer to the second question regarding 
   	the first prospect.
   The values of the outcomes of the 6 prospect are in data b. g1 and g2 are the gains. g3 is the loss
   The answers were made on a 7 point Likert scale. The first two are the simple measures. The last four are
   	Sitkin and Wingart scales
   1 -risk, 2 -return, 3- threat 4- potentail for loss, 5 - negative investment, 6, low success probatility;
 sub = _n_;
 	***Reversing items 4 5 and 8;
 q14=8-q14; q15=8-q15; q16=8-q16;
 q24=8-q24; q25=8-q25; q26=8-q26;
 q34=8-q34; q35=8-q35; q36=8-q36;
 q44=8-q44; q45=8-q45; q46=8-q46;
 q54=8-q54; q55=8-q55; q56=8-q56;
 q64=8-q64; q65=8-q65; q66=8-q66;
 q74=8-q74; q75=8-q75; q76=8-q76;
 q84=8-q84; q85=8-q85; q86=8-q86;

 data a1; set a; keep con sub q11-q16 qq1-qq6 stim;stim=1;rename q11=qq1 q12=qq2 q13=qq3 q14=qq4 q15=qq5 q16=qq6;
 data a2; set a; keep con sub q21-q26 qq1-qq6 stim;stim=2;rename q21=qq1 q22=qq2 q23=qq3 q24=qq4 q25=qq5 q26=qq6;
 data a3; set a; keep con sub q31-q36 qq1-qq6 stim;stim=3;rename q31=qq1 q32=qq2 q33=qq3 q34=qq4 q35=qq5 q36=qq6;
 data a4; set a; keep con sub q41-q46 qq1-qq6 stim;stim=4;rename q41=qq1 q42=qq2 q43=qq3 q44=qq4 q45=qq5 q46=qq6;
 data a5; set a; keep con sub q51-q56 qq1-qq6 stim;stim=5;rename q51=qq1 q52=qq2 q53=qq3 q54=qq4 q55=qq5 q56=qq6;
 data a6; set a; keep con sub q61-q66 qq1-qq6 stim;stim=6;rename q61=qq1 q62=qq2 q63=qq3 q64=qq4 q65=qq5 q66=qq6;
 data a7; set a; keep con sub q71-q76 qq1-qq6 stim;stim=7;rename q71=qq1 q72=qq2 q73=qq3 q74=qq4 q75=qq5 q76=qq6;
 data a8; set a; keep con sub q81-q86 qq1-qq6 stim;stim=8;rename q81=qq1 q82=qq2 q83=qq3 q84=qq4 q85=qq5 q86=qq6;

data b; set a1 a2 a3 a4 a5 a6 a7 a8;
if stim=1 then do;g1=100;g2=190;g3=-60 ;end;
if stim=2 then do;g1=150;g2=190;g3=-50 ;end;
if stim=3 then do;g1=130;g2=240;g3=-80 ;end;
if stim=4 then do;g1=70;g2=110;g3=-50 ;end;
if stim=5 then do;g1=50;g2=160;g3=-80 ;end;
if stim=6 then do;g1=40;g2=60;g3=-30 ;end;
if stim=7 then do;g1=120;g2=140;g3=-30 ;end;
if stim=8 then do;g1=20;g2=110;g3=-60 ;end;
sit=mean(of qq3-qq6);
sd=std(of g1-g3); ev=mean(of g1-g3);

proc sort; by sub stim;
proc corr outp=c; var qq1 qq2 sit sd ev g3; by sub; run;
data c1; set c; if _name_='sit'; keep sub qq1 qq2; rename qq1=risit qq2=resit;
data c2; set c;if _name_='qq2'; keep sub qq1 rire; rename qq1= rire;
data c3; set c;if _name_='qq1';keep sub sd ev g3 risd riev rilos; rename sd=risd ev=riev g3=rilos;
data c4; set c;if _name_='qq2';keep sub sd ev g3 resd reev relos; rename sd=resd ev=reev g3=relos;
data c5; set c;if _name_='sit';keep sub sd ev g3 sitsd sitev sitlos; rename sd=sitsd ev=sitev g3=sitlos;
data d; merge c1 c2 c3 c4 c5; by sub;
if sub=8 then delete; * risit=0;
if sub=10 then delete; if sub=12 thend delete; if sub=14 then delete;if sub=16 then delete; if sub=17 then delete;
rilos=-rilos; relos=-relos;sitlos=-sitlos;
risit_resit = risit+resit;
sitlos_rilos = sitlos-rilos;
sitev_sitlos=sitev+sitlos;
resit_rire=resit-rire;
riev_rilos=riev+rilos;
proc means mean t prt stderr n; run;