## 1a ## Objective Probability Diagonal C/X ## 0.05 0.05 0.1 ## 0.1 0.1 0.13 ## 0.15 0.15 0.18 ## 0.25 0.25 0.22 ## 0.3 0.3 0.28 ## 0.4 0.4 0.29 ## 0.6 0.6 0.43 ## 0.7 0.7 0.49 ## 0.75 0.75 0.48 ## 0.85 0.85 0.62 ## 0.9 0.9 0.59 ## 0.95 0.95 0.76 ## 1b ## Diagonal Circle Triangle ## 0.01 0.01 0.05 0.03 ## 0.05 0.05 0.14 NA ## 0.1 0.1 0.18 0.1 ## 0.25 0.25 0.25 NA ## 0.5 0.5 0.42 0.38 ## 0.5 0.5 0.36 NA ## 0.75 0.75 0.52 NA ## 0.9 0.9 0.74 0.66 ## 0.95 0.95 0.78 NA ## 0.99 0.99 0.94 0.94 ## 2a ## Outcomes Events ## 3 vs. 0,32 65 73 ## 3 vs. 0 (incomplete) 9 0 ## 3 vs. 0 (complete) 62 66 ## 2b ## Outcomes Events ## -3 vs. 0,-32 70 69 ## -3 vs. 0 (incomplete) 0 13 ## -3 vs. 0 (complete) 60 46 ## 2c ## Outcomes Events ## 0,3 vs. 0,32 93 86 ## 0,3 vs. 0 (incomplete) 14 5 ## 0,3 vs. 0 (complete) 100 89 ## 3 ## A: (3,1) > (4,.8) B: (3,.25) > (4,.2) ## A then B 73 59 ## B then A 29 42 f1a <- matrix(c(0.05,0.1 ,0.15,0.25,0.3 ,0.4 ,0.6 ,0.7 ,0.75,0.85,0.9 ,0.95, 0.1 ,0.13,0.18,0.22,0.28,0.29,0.43,0.49,0.48,0.62,0.59,0.76),,2) postscript(file="fig1a.eps",width=8,height=8, horiz=F,onefile=F,pointsize=18,paper="special") plot(f1a[,1:2],xlab="Objective probability (p)", ylab="Median normalized certainty equivalent (c/x)", xlim=c(0,1),ylim=c(0,1)) lines(c(0,1),c(0,1),lty=2) text(c(.2,.7),c(.7,.2),labels=c("Risk-seeking","Risk-aversion"), col=c("darkgreen","red")) dev.off() system("bbox fig1a.eps") f1b <- matrix(c(0.01,0.05,0.1,0.25,0.5,0.5,0.75,0.9,0.95,0.99, 0.05,0.14,0.18,0.25,0.42,0.36,0.52,0.74,0.78,0.94, 0.03,NA,0.1,NA,0.38,NA,NA,0.66,NA,0.94),,3) postscript(file="fig1b.eps",width=8,height=8, horiz=F,onefile=F,pointsize=18,paper="special") plot(f1b[,1:2],xlab="Objective probability (p)", ylab="Median normalized certainty equivalent (c/x)", xlim=c(0,1),ylim=c(0,1),col="blue") points(f1b[,c(1,3)],pch=2,col="brown") lines(c(0,1),c(0,1),lty=2) # curve(1.239*x^4 - .4141*x^3 - 1.0697*x^2 + 1.1265*x + .0624,add=T,col="blue") dev.off() system("bbox fig1b.eps") f2b <- matrix(c(65,9,62,73,0,66),3,2) rownames(f2b)=c("3 vs. 0,32\n","3 vs. 0\nincomplete","3 vs. 0\ncomplete") postscript(file="fig2b.eps",width=8,height=6, horiz=F,onefile=F,pointsize=18,paper="special") barplot(t(f2b),beside=T,ylim=c(0,100)) legend("top",c("Outcomes","Events"),fill=c("gray40","gray80")) dev.off() system("bbox fig2b.eps") f2a <- matrix(c(70,0,60,69,13,46),3,2) rownames(f2a) <- c("-3 vs. 0,-32\n","-3 vs. 0\nincomplete","-3 vs. 0\ncomplete") postscript(file="fig2a.eps",width=8,height=6, horiz=F,onefile=F,pointsize=18,paper="special") barplot(t(f2a),beside=T,ylim=c(0,100)) legend("top",c("Outcomes","Events"),fill=c("gray40","gray80")) dev.off() system("bbox fig2a.eps") f2c <- matrix(c(93,13,89,86,5,89),3,2) rownames(f2c) <- c("3 vs. 0,32\n","3 vs. 0\nincomplete","3 vs. 0\ncomplete") postscript(file="fig2c.eps",width=8,height=6, horiz=F,onefile=F,pointsize=18,paper="special") barplot(t(f2c),beside=T,ylim=c(0,100)) legend("top",c("Outcomes","Events"),fill=c("gray40","gray80")) dev.off() system("bbox fig2c.eps") f3 <- matrix(c(73,29,59,42),2,2) # rownames(f3) <- c("A: (3,1) > (4,.8)","B: (3,.25) > (4,.2)") postscript(file="fig3.eps",width=8,height=6, horiz=F,onefile=F,pointsize=18,paper="special") barplot(f3,beside=T,ylim=c(0,100)) legend("top",c("A then B","B then A"),fill=c("gray40","gray80")) dev.off() system("bbox fig3.eps")