plot(0:100,ty="n",axes=F,xlab="",ylab="",main="Non-patient")
arrows(c(50,0),c(0,50),c(50,100),c(100,50),lwd=2)
curve(.5*sqrt(50)*sqrt(x-50)+50,from=50,to=100,add=T)
curve(-sqrt(50)*sqrt(50-x)+50,from=0,to=50,add=T)
text(42,90,labels="HRQol")
text(88,53,labels="Health")
segments(c(20,35),c(50,50),c(20,35),c(50-sqrt(50)*sqrt(30),50-sqrt(50)*sqrt(15)),lty=2)
text(c(70,70,70),c(35,20,2),labels=c("no lung disease",
                               "mild lung disease","severe lung disease"))
arrows(c(59,57,55),c(37,20,2),c(51,37,21),c(49,23,11))
dev.copy2eps(file="fig1b.eps")
dev.off()

plot(0:100,ty="n",axes=F,xlab="",ylab="",main="Lung disease patient")
arrows(c(50,0),c(0,50),c(50,100),c(100,50),lwd=2)
curve(.5*sqrt(50)*sqrt(x-50)+50,from=50,to=100,add=T)
curve(-sqrt(50)*sqrt(50-x)+50,from=0,to=50,add=T)
text(42,90,labels="HRQol")
text(88,53,labels="Health")
segments(c(35,65),c(50,50),c(35,65),c(50-sqrt(50)*sqrt(15),.5*sqrt(50)*sqrt(15)+50),lty=2)
text(c(80,80,20),c(35,20,65),labels=c("mild lung disease",
                               "moderate lung disease","severe lung disease"))
arrows(c(65,69,30),c(22,37,63),c(51,65,35),c(49,49,51))
dev.copy2eps(file="fig1a.eps")
dev.off()