plot(0:100,ylab="Qol ratings",xlab="",typ="n",xaxt="n")
mtext(c("No context","Context","No context","Context"),at=c(15,35,65,85),side=1,padj=1)
mtext(c("Healthy","Patients"),at=c(25,75),side=1,padj=3)
segments(c(15,15,65,65),c(54.89,39.80,70.27,41.58),c(35,35,85,85),c(64,18.85,78.72,11.83),
         lty=c(1,2,1,2))
points(c(15,15,65,65,35,35,85,85),c(54.89,39.80,70.27,41.58,64,18.85,78.72,11.83))
legend(x=10,y=90,legend=c("Mild","Severe"),lty=c(1,2),cex=.9)
dev.copy2eps(file="fig2.eps")
dev.off()

#,"Healthy",,"Patients",
#,"No Context","Context","No Context","Context"
#"Mild",54.89,64,70.27,78.72
#"Severe",39.08,18.85,41.58,11.83