pdf("roc1.pdf",width=10,height=10,paper="special")#
par(cex.axis=1.7,cex.lab=1.7,pty="s")#
plot(pca.perf.r50,col="black",lty=1,lwd=1,cex=1.7)#
plot(pca.perf.r100,add=TRUE,col="black",lty=2,lwd=1)#
plot(pca.perf.r200,add=TRUE,col="black",lty=3,lwd=1)#
plot(rasch.perf.r50,add=TRUE,col="black",lty=1,lwd=3)#
plot(rasch.perf.r100,add=TRUE,col="black",lty=2,lwd=3)#
plot(rasch.perf.r200,add=TRUE,col="black",lty=3,lwd=3)#
legend("bottomright",c(expression(RM~(paste(p,"*",sep="")==50)),expression(RM~(paste(p,"*",sep="")==100)),expression(RM~(paste(p,"*",sep="")==200)),expression(PCA~(paste(p,"*",sep="")==50)),expression(PCA~(paste(p,"*",sep="")==100)),expression(PCA~(paste(p,"*",sep="")==200))),lty=c(1,2,3,1,2,3),lwd=c(2,2,2,1,1,1),bty="n",cex=1.7)#
dev.off()
pdf("roc1.pdf",width=10,height=10,paper="special")#
par(cex.axis=1.7,cex.lab=1.7,pty="s")#
plot(pca.perf.r50,col="black",lty=1,lwd=2,cex=1.7)#
plot(pca.perf.r100,add=TRUE,col="black",lty=2,lwd=2)#
plot(pca.perf.r200,add=TRUE,col="black",lty=3,lwd=2)#
plot(rasch.perf.r50,add=TRUE,col="red",lty=1,lwd=2)#
plot(rasch.perf.r100,add=TRUE,col="blue",lty=2,lwd=2)#
plot(rasch.perf.r200,add=TRUE,col="green",lty=3,lwd=2)#
legend("bottomright",c(expression(RM~(paste(p,"*",sep="")==50)),expression(RM~(paste(p,"*",sep="")==100)),expression(RM~(paste(p,"*",sep="")==200)),expression(PCA~(paste(p,"*",sep="")==50)),expression(PCA~(paste(p,"*",sep="")==100)),expression(PCA~(paste(p,"*",sep="")==200))),lty=c(1,2,3,1,2,3),lwd=c(2,2,2,1,1,1),bty="n",cex=1.7)#
dev.off()
pdf("roc1.pdf",width=10,height=10,paper="special")#
par(cex.axis=1.7,cex.lab=1.7,pty="s")#
plot(pca.perf.r50,col="black",lty=1,lwd=1,cex=1.7)#
plot(pca.perf.r100,add=TRUE,col="black",lty=2,lwd=1)#
plot(pca.perf.r200,add=TRUE,col="black",lty=3,lwd=1)#
plot(rasch.perf.r50,add=TRUE,col="black",lty=1,lwd=3)#
plot(rasch.perf.r100,add=TRUE,col="black",lty=2,lwd=3)#
plot(rasch.perf.r200,add=TRUE,col="black",lty=3,lwd=3)#
legend("bottomright",c(expression(RM~(paste(p,"*",sep="")==50)),expression(RM~(paste(p,"*",sep="")==100)),expression(RM~(paste(p,"*",sep="")==200)),expression(PCA~(paste(p,"*",sep="")==50)),expression(PCA~(paste(p,"*",sep="")==100)),expression(PCA~(paste(p,"*",sep="")==200))),lty=c(1,2,3,1,2,3),lwd=c(3,3,3,1,1,1),bty="n",cex=1.7)#
dev.off()
pdf("roc1.pdf",width=10,height=10,paper="special")#
par(cex.axis=1.7,cex.lab=1.7,pty="s")#
plot(pca.perf.r50,col="bray",lty=1,lwd=3,cex=1.7)#
plot(pca.perf.r100,add=TRUE,col="gray",lty=2,lwd=3)#
plot(pca.perf.r200,add=TRUE,col="gray",lty=3,lwd=3)#
plot(rasch.perf.r50,add=TRUE,col="black",lty=1,lwd=3)#
plot(rasch.perf.r100,add=TRUE,col="black",lty=2,lwd=3)#
plot(rasch.perf.r200,add=TRUE,col="black",lty=3,lwd=3)#
legend("bottomright",c(expression(RM~(paste(p,"*",sep="")==50)),expression(RM~(paste(p,"*",sep="")==100)),expression(RM~(paste(p,"*",sep="")==200)),expression(PCA~(paste(p,"*",sep="")==50)),expression(PCA~(paste(p,"*",sep="")==100)),expression(PCA~(paste(p,"*",sep="")==200))),lty=c(1,2,3,1,2,3),lwd=c(3,3,3,1,1,1),col=c("black","black","black","gray","gray","gray"),bty="n",cex=1.7)#
dev.off()
pdf("roc1.pdf",width=10,height=10,paper="special")#
par(cex.axis=1.7,cex.lab=1.7,pty="s")#
plot(pca.perf.r50,col="bray",lty=1,lwd=3,cex=1.7)#
plot(pca.perf.r100,add=TRUE,col="gray",lty=2,lwd=3)#
plot(pca.perf.r200,add=TRUE,col="gray",lty=3,lwd=3)#
plot(rasch.perf.r50,add=TRUE,col="black",lty=1,lwd=3)#
plot(rasch.perf.r100,add=TRUE,col="black",lty=2,lwd=3)#
plot(rasch.perf.r200,add=TRUE,col="black",lty=3,lwd=3)#
legend("bottomright",c(expression(RM~(paste(p,"*",sep="")==50)),expression(RM~(paste(p,"*",sep="")==100)),expression(RM~(paste(p,"*",sep="")==200)),expression(PCA~(paste(p,"*",sep="")==50)),expression(PCA~(paste(p,"*",sep="")==100)),expression(PCA~(paste(p,"*",sep="")==200))),lty=c(1,2,3,1,2,3),lwd=c(3,3,3,3,3,3),col=c("black","black","black","gray","gray","gray"),bty="n",cex=1.7)#
dev.off()
pdf("roc1.pdf",width=10,height=10,paper="special")#
par(cex.axis=1.7,cex.lab=1.7,pty="s")#
plot(pca.perf.r50,col="bray",lty=1,lwd=3,cex=1.7)#
plot(pca.perf.r100,add=TRUE,col="gray",lty=2,lwd=3)#
plot(pca.perf.r200,add=TRUE,col="gray",lty=3,lwd=3)#
plot(rasch.perf.r50,add=TRUE,col="black",lty=1,lwd=3)#
plot(rasch.perf.r100,add=TRUE,col="black",lty=2,lwd=3)#
plot(rasch.perf.r200,add=TRUE,col="black",lty=3,lwd=3)#
legend("bottomright",c(expression(RM~(paste(p,"*",sep="")==50)),expression(RM~(paste(p,"*",sep="")==100)),expression(RM~(paste(p,"*",sep="")==200)),expression(PCA~(paste(p,"*",sep="")==50)),expression(PCA~(paste(p,"*",sep="")==100)),expression(PCA~(paste(p,"*",sep="")==200))),lty=c(1,2,3,1,2,3),lwd=c(3,3,3,3,3,3),bty="n",cex=1.7)#
dev.off()
pdf("roc1.pdf",width=10,height=10,paper="special")#
par(cex.axis=1.7,cex.lab=1.7,pty="s")#
plot(pca.perf.r50,col="gray",lty=1,lwd=3,cex=1.7)#
plot(pca.perf.r100,add=TRUE,col="gray",lty=2,lwd=3)#
plot(pca.perf.r200,add=TRUE,col="gray",lty=3,lwd=3)#
plot(rasch.perf.r50,add=TRUE,col="black",lty=1,lwd=3)#
plot(rasch.perf.r100,add=TRUE,col="black",lty=2,lwd=3)#
plot(rasch.perf.r200,add=TRUE,col="black",lty=3,lwd=3)#
legend("bottomright",c(expression(RM~(paste(p,"*",sep="")==50)),expression(RM~(paste(p,"*",sep="")==100)),expression(RM~(paste(p,"*",sep="")==200)),expression(PCA~(paste(p,"*",sep="")==50)),expression(PCA~(paste(p,"*",sep="")==100)),expression(PCA~(paste(p,"*",sep="")==200))),lty=c(1,2,3,1,2,3),lwd=c(3,3,3,3,3,3),col=c("black","black","black","gray","gray","gray"),bty="n",cex=1.7)#
dev.off()
pdf("roc1.pdf",width=10,height=10,paper="special")#
par(cex.axis=1.7,cex.lab=1.7,pty="s")#
plot(pca.perf.r50,col="gray",lty=1,lwd=3,cex=1.7)#
plot(pca.perf.r100,add=TRUE,col="gray",lty=2,lwd=3)#
plot(pca.perf.r200,add=TRUE,col="gray",lty=3,lwd=3)#
plot(rasch.perf.r50,add=TRUE,col="black",lty=1,lwd=3)#
plot(rasch.perf.r100,add=TRUE,col="black",lty=2,lwd=3)#
plot(rasch.perf.r200,add=TRUE,col="black",lty=3,lwd=3)#
legend("bottomright",c(expression(RM\t(paste(p,"*",sep="")==50)),expression(RM~(paste(p,"*",sep="")==100)),expression(RM~(paste(p,"*",sep="")==200)),expression(PCA~(paste(p,"*",sep="")==50)),expression(PCA~(paste(p,"*",sep="")==100)),expression(PCA~(paste(p,"*",sep="")==200))),lty=c(1,2,3,1,2,3),lwd=c(3,3,3,3,3,3),col=c("black","black","black","gray","gray","gray"),bty="n",cex=1.7)#
dev.off()
pdf("roc1.pdf",width=10,height=10,paper="special")#
par(cex.axis=1.7,cex.lab=1.7,pty="s")#
plot(pca.perf.r50,col="gray",lty=1,lwd=3,cex=1.7)#
plot(pca.perf.r100,add=TRUE,col="gray",lty=2,lwd=3)#
plot(pca.perf.r200,add=TRUE,col="gray",lty=3,lwd=3)#
plot(rasch.perf.r50,add=TRUE,col="black",lty=1,lwd=3)#
plot(rasch.perf.r100,add=TRUE,col="black",lty=2,lwd=3)#
plot(rasch.perf.r200,add=TRUE,col="black",lty=3,lwd=3)#
legend("bottomright",c(expression(RM"\t"(paste(p,"*",sep="")==50)),expression(RM~(paste(p,"*",sep="")==100)),expression(RM~(paste(p,"*",sep="")==200)),expression(PCA~(paste(p,"*",sep="")==50)),expression(PCA~(paste(p,"*",sep="")==100)),expression(PCA~(paste(p,"*",sep="")==200))),lty=c(1,2,3,1,2,3),lwd=c(3,3,3,3,3,3),col=c("black","black","black","gray","gray","gray"),bty="n",cex=1.7)#
dev.off()
pdf("roc1.pdf",width=10,height=10,paper="special")#
par(cex.axis=1.7,cex.lab=1.7,pty="s")#
plot(pca.perf.r50,col="gray",lty=1,lwd=3,cex=1.7)#
plot(pca.perf.r100,add=TRUE,col="gray",lty=2,lwd=3)#
plot(pca.perf.r200,add=TRUE,col="gray",lty=3,lwd=3)#
plot(rasch.perf.r50,add=TRUE,col="black",lty=1,lwd=3)#
plot(rasch.perf.r100,add=TRUE,col="black",lty=2,lwd=3)#
plot(rasch.perf.r200,add=TRUE,col="black",lty=3,lwd=3)#
legend("bottomright",c(expression(RM~(paste(p,"*",sep="")==50)),expression(RM~(paste(p[RM],"*",sep="")==100)),expression(RM~(paste(p,"*",sep="")==200)),expression(PCA~(paste(p,"*",sep="")==50)),expression(PCA~(paste(p,"*",sep="")==100)),expression(PCA~(paste(p,"*",sep="")==200))),lty=c(1,2,3,1,2,3),lwd=c(3,3,3,3,3,3),col=c("black","black","black","gray","gray","gray"),bty="n",cex=1.7)#
dev.off()
pdf("roc1.pdf",width=10,height=10,paper="special")#
par(cex.axis=1.7,cex.lab=1.7,pty="s")#
plot(pca.perf.r50,col="gray",lty=1,lwd=3,cex=1.7)#
plot(pca.perf.r100,add=TRUE,col="gray",lty=2,lwd=3)#
plot(pca.perf.r200,add=TRUE,col="gray",lty=3,lwd=3)#
plot(rasch.perf.r50,add=TRUE,col="black",lty=1,lwd=3)#
plot(rasch.perf.r100,add=TRUE,col="black",lty=2,lwd=3)#
plot(rasch.perf.r200,add=TRUE,col="black",lty=3,lwd=3)#
legend("bottomright",c(expression(RM~(paste(p,"*",sep="")==50)),expression(RM~~(paste(p,"*",sep="")==100)),expression(RM~(paste(p,"*",sep="")==200)),expression(PCA~(paste(p,"*",sep="")==50)),expression(PCA~(paste(p,"*",sep="")==100)),expression(PCA~(paste(p,"*",sep="")==200))),lty=c(1,2,3,1,2,3),lwd=c(3,3,3,3,3,3),col=c("black","black","black","gray","gray","gray"),bty="n",cex=1.7)#
dev.off()
pdf("roc1.pdf",width=10,height=10,paper="special")#
par(cex.axis=1.7,cex.lab=1.7,pty="s")#
plot(pca.perf.r50,col="gray",lty=1,lwd=3,cex=1.7)#
plot(pca.perf.r100,add=TRUE,col="gray",lty=2,lwd=3)#
plot(pca.perf.r200,add=TRUE,col="gray",lty=3,lwd=3)#
plot(rasch.perf.r50,add=TRUE,col="black",lty=1,lwd=3)#
plot(rasch.perf.r100,add=TRUE,col="black",lty=2,lwd=3)#
plot(rasch.perf.r200,add=TRUE,col="black",lty=3,lwd=3)#
legend("bottomright",c(expression(RM~(paste(p,"*",sep="")==50)),expression(RM~[paste(p,"*",sep="")==100]),expression(RM~(paste(p,"*",sep="")==200)),expression(PCA~(paste(p,"*",sep="")==50)),expression(PCA~(paste(p,"*",sep="")==100)),expression(PCA~(paste(p,"*",sep="")==200))),lty=c(1,2,3,1,2,3),lwd=c(3,3,3,3,3,3),col=c("black","black","black","gray","gray","gray"),bty="n",cex=1.7)#
dev.off()
pdf("roc1.pdf",width=10,height=10,paper="special")#
par(cex.axis=1.7,cex.lab=1.7,pty="s")#
plot(pca.perf.r50,col="gray",lty=1,lwd=3,cex=1.7)#
plot(pca.perf.r100,add=TRUE,col="gray",lty=2,lwd=3)#
plot(pca.perf.r200,add=TRUE,col="gray",lty=3,lwd=3)#
plot(rasch.perf.r50,add=TRUE,col="black",lty=1,lwd=3)#
plot(rasch.perf.r100,add=TRUE,col="black",lty=2,lwd=3)#
plot(rasch.perf.r200,add=TRUE,col="black",lty=3,lwd=3)#
legend("bottomright",c(expression(RM~(paste(p,"*",sep="")==50)),expression(RM--(paste(p,"*",sep="")==100)),expression(RM~(paste(p,"*",sep="")==200)),expression(PCA~(paste(p,"*",sep="")==50)),expression(PCA~(paste(p,"*",sep="")==100)),expression(PCA~(paste(p,"*",sep="")==200))),lty=c(1,2,3,1,2,3),lwd=c(3,3,3,3,3,3),col=c("black","black","black","gray","gray","gray"),bty="n",cex=1.7)#
dev.off()
pdf("roc1.pdf",width=10,height=10,paper="special")#
par(cex.axis=1.7,cex.lab=1.7,pty="s")#
plot(pca.perf.r50,col="gray",lty=1,lwd=3,cex=1.7)#
plot(pca.perf.r100,add=TRUE,col="gray",lty=2,lwd=3)#
plot(pca.perf.r200,add=TRUE,col="gray",lty=3,lwd=3)#
plot(rasch.perf.r50,add=TRUE,col="black",lty=1,lwd=3)#
plot(rasch.perf.r100,add=TRUE,col="black",lty=2,lwd=3)#
plot(rasch.perf.r200,add=TRUE,col="black",lty=3,lwd=3)#
legend("bottomright",c(expression(RM~(paste(p,"*",sep="")==50)),expression(RM-paste(p,"*",sep="")==100),expression(RM~(paste(p,"*",sep="")==200)),expression(PCA~(paste(p,"*",sep="")==50)),expression(PCA~(paste(p,"*",sep="")==100)),expression(PCA~(paste(p,"*",sep="")==200))),lty=c(1,2,3,1,2,3),lwd=c(3,3,3,3,3,3),col=c("black","black","black","gray","gray","gray"),bty="n",cex=1.7)#
dev.off()
pdf("roc1.pdf",width=10,height=10,paper="special")#
par(cex.axis=1.7,cex.lab=1.7,pty="s")#
plot(pca.perf.r50,col="gray",lty=1,lwd=3,cex=1.7)#
plot(pca.perf.r100,add=TRUE,col="gray",lty=2,lwd=3)#
plot(pca.perf.r200,add=TRUE,col="gray",lty=3,lwd=3)#
plot(rasch.perf.r50,add=TRUE,col="black",lty=1,lwd=3)#
plot(rasch.perf.r100,add=TRUE,col="black",lty=2,lwd=3)#
plot(rasch.perf.r200,add=TRUE,col="black",lty=3,lwd=3)#
legend("bottomright",c(expression(RM~(paste(p,"*",sep="")==50)),expression(RM-(paste(p,"*",sep="")==100)),expression(RM~(paste(p,"*",sep="")==200)),expression(PCA~(paste(p,"*",sep="")==50)),expression(PCA~(paste(p,"*",sep="")==100)),expression(PCA~(paste(p,"*",sep="")==200))),lty=c(1,2,3,1,2,3),lwd=c(3,3,3,3,3,3),col=c("black","black","black","gray","gray","gray"),bty="n",cex=1.7)#
dev.off()
pdf("roc1.pdf",width=10,height=10,paper="special")#
par(cex.axis=1.7,cex.lab=1.7,pty="s")#
plot(pca.perf.r50,col="gray",lty=1,lwd=3,cex=1.7)#
plot(pca.perf.r100,add=TRUE,col="gray",lty=2,lwd=3)#
plot(pca.perf.r200,add=TRUE,col="gray",lty=3,lwd=3)#
plot(rasch.perf.r50,add=TRUE,col="black",lty=1,lwd=3)#
plot(rasch.perf.r100,add=TRUE,col="black",lty=2,lwd=3)#
plot(rasch.perf.r200,add=TRUE,col="black",lty=3,lwd=3)#
legend("bottomright",c(expression(RM~(paste(p,"*",sep="")==50)),expression(RM-(paste(p,"*",sep="")==100)),expression(RM~(paste(p,"*",sep="")==200)),expression(PCA~(paste(p,"*",sep="")==50)),expression(PCA~(paste(p,"*",sep="")==100)),expression(PCA~(paste(p,"*",sep="")==200))),lty=c(1,2,3,1,2,3),lwd=c(3,3,3,3,3,3),col=c("black","black","black","gray","gray","gray60"),bty="n",cex=1.7)#
dev.off()
pdf("roc1.pdf",width=10,height=10,paper="special")#
par(cex.axis=1.7,cex.lab=1.7,pty="s")#
plot(pca.perf.r50,col="gray",lty=1,lwd=3,cex=1.7)#
plot(pca.perf.r100,add=TRUE,col="gray",lty=2,lwd=3)#
plot(pca.perf.r200,add=TRUE,col="gray",lty=3,lwd=3)#
plot(rasch.perf.r50,add=TRUE,col="black",lty=1,lwd=3)#
plot(rasch.perf.r100,add=TRUE,col="black",lty=2,lwd=3)#
plot(rasch.perf.r200,add=TRUE,col="black",lty=3,lwd=3)#
legend("bottomright",c(expression(RM~(paste(p,"*",sep="")==50)),expression(RM-(paste(p,"*",sep="")==100)),expression(RM~(paste(p,"*",sep="")==200)),expression(PCA~(paste(p,"*",sep="")==50)),expression(PCA~(paste(p,"*",sep="")==100)),expression(PCA~(paste(p,"*",sep="")==200))),lty=c(1,2,3,1,2,3),lwd=c(3,3,3,3,3,3),col=c("black","black","black","gray50","gray50","gray50"),bty="n",cex=1.7)#
dev.off()
pdf("roc1.pdf",width=10,height=10,paper="special")#
par(cex.axis=1.7,cex.lab=1.7,pty="s")#
plot(pca.perf.r50,col="gray",lty=1,lwd=3,cex=1.7)#
plot(pca.perf.r100,add=TRUE,col="gray",lty=2,lwd=3)#
plot(pca.perf.r200,add=TRUE,col="gray",lty=3,lwd=3)#
plot(rasch.perf.r50,add=TRUE,col="black",lty=1,lwd=3)#
plot(rasch.perf.r100,add=TRUE,col="black",lty=2,lwd=3)#
plot(rasch.perf.r200,add=TRUE,col="black",lty=3,lwd=3)#
legend("bottomright",c(expression(RM~(paste(p,"*",sep="")==50)),expression(RM-(paste(p,"*",sep="")==100)),expression(RM~(paste(p,"*",sep="")==200)),expression(PCA~(paste(p,"*",sep="")==50)),expression(PCA~(paste(p,"*",sep="")==100)),expression(PCA~(paste(p,"*",sep="")==200))),lty=c(1,2,3,1,2,3),lwd=c(3,3,3,3,3,3),col=c("black","black","black","gray55","gray55","gray55"),bty="n",cex=1.7)#
dev.off()
pdf("roc1.pdf",width=10,height=10,paper="special")#
par(cex.axis=1.7,cex.lab=1.7,pty="s")#
plot(pca.perf.r50,col="gray",lty=1,lwd=3,cex=1.7)#
plot(pca.perf.r100,add=TRUE,col="gray",lty=2,lwd=3)#
plot(pca.perf.r200,add=TRUE,col="gray",lty=3,lwd=3)#
plot(rasch.perf.r50,add=TRUE,col="black",lty=1,lwd=3)#
plot(rasch.perf.r100,add=TRUE,col="black",lty=2,lwd=3)#
plot(rasch.perf.r200,add=TRUE,col="black",lty=3,lwd=3)#
legend("bottomright",c(expression(RM~(paste(p,"*",sep="")==50)),expression(RM~(paste(p,"*",sep="")==100)),expression(RM~(paste(p,"*",sep="")==200)),expression(PCA~(paste(p,"*",sep="")==50)),expression(PCA~(paste(p,"*",sep="")==100)),expression(PCA~(paste(p,"*",sep="")==200))),lty=c(1,2,3,1,2,3),lwd=c(3,3,3,3,3,3),col=c("black","black","black","gray55","gray55","gray55"),bty="n",cex=1.7)#
dev.off()
pdf("roc1.pdf",width=10,height=10,paper="special")#
par(cex.axis=1.7,cex.lab=1.7,pty="s")#
plot(pca.perf.r50,col="gray",lty=1,lwd=3,cex=1.7)#
plot(pca.perf.r100,add=TRUE,col="gray",lty=2,lwd=3)#
plot(pca.perf.r200,add=TRUE,col="gray",lty=3,lwd=3)#
plot(rasch.perf.r50,add=TRUE,col="black",lty=1,lwd=3)#
plot(rasch.perf.r100,add=TRUE,col="black",lty=2,lwd=3)#
plot(rasch.perf.r200,add=TRUE,col="black",lty=3,lwd=3)#
legend("bottomright",c(expression(RM[{paste(p,"*",sep="")==50}]),expression(RM~(paste(p,"*",sep="")==100)),expression(RM~(paste(p,"*",sep="")==200)),expression(PCA~(paste(p,"*",sep="")==50)),expression(PCA~(paste(p,"*",sep="")==100)),expression(PCA~(paste(p,"*",sep="")==200))),lty=c(1,2,3,1,2,3),lwd=c(3,3,3,3,3,3),col=c("black","black","black","gray55","gray55","gray55"),bty="n",cex=1.7)#
dev.off()
pdf("roc1.pdf",width=10,height=10,paper="special")#
par(cex.axis=1.7,cex.lab=1.7,pty="s")#
plot(pca.perf.r50,col="gray",lty=1,lwd=3,cex=1.7)#
plot(pca.perf.r100,add=TRUE,col="gray",lty=2,lwd=3)#
plot(pca.perf.r200,add=TRUE,col="gray",lty=3,lwd=3)#
plot(rasch.perf.r50,add=TRUE,col="black",lty=1,lwd=3)#
plot(rasch.perf.r100,add=TRUE,col="black",lty=2,lwd=3)#
plot(rasch.perf.r200,add=TRUE,col="black",lty=3,lwd=3)#
legend("bottomright",c(expression(RM[{paste(p,"*",sep="")~==~50}]),expression(RM~(paste(p,"*",sep="")==100)),expression(RM~(paste(p,"*",sep="")==200)),expression(PCA~(paste(p,"*",sep="")==50)),expression(PCA~(paste(p,"*",sep="")==100)),expression(PCA~(paste(p,"*",sep="")==200))),lty=c(1,2,3,1,2,3),lwd=c(3,3,3,3,3,3),col=c("black","black","black","gray55","gray55","gray55"),bty="n",cex=1.7)#
dev.off()
pdf("roc1.pdf",width=10,height=10,paper="special")#
par(cex.axis=1.7,cex.lab=1.7,pty="s")#
plot(pca.perf.r50,col="gray",lty=1,lwd=3,cex=1.7)#
plot(pca.perf.r100,add=TRUE,col="gray",lty=2,lwd=3)#
plot(pca.perf.r200,add=TRUE,col="gray",lty=3,lwd=3)#
plot(rasch.perf.r50,add=TRUE,col="black",lty=1,lwd=3)#
plot(rasch.perf.r100,add=TRUE,col="black",lty=2,lwd=3)#
plot(rasch.perf.r200,add=TRUE,col="black",lty=3,lwd=3)#
legend("bottomright",c(expression(RM[{paste(p,"*",sep="") == 50}]),expression(RM~(paste(p,"*",sep="")==100)),expression(RM~(paste(p,"*",sep="")==200)),expression(PCA~(paste(p,"*",sep="")==50)),expression(PCA~(paste(p,"*",sep="")==100)),expression(PCA~(paste(p,"*",sep="")==200))),lty=c(1,2,3,1,2,3),lwd=c(3,3,3,3,3,3),col=c("black","black","black","gray55","gray55","gray55"),bty="n",cex=1.7)#
dev.off()
pdf("roc1.pdf",width=10,height=10,paper="special")#
par(cex.axis=1.7,cex.lab=1.7,pty="s")#
plot(pca.perf.r50,col="gray",lty=1,lwd=3,cex=1.7)#
plot(pca.perf.r100,add=TRUE,col="gray",lty=2,lwd=3)#
plot(pca.perf.r200,add=TRUE,col="gray",lty=3,lwd=3)#
plot(rasch.perf.r50,add=TRUE,col="black",lty=1,lwd=3)#
plot(rasch.perf.r100,add=TRUE,col="black",lty=2,lwd=3)#
plot(rasch.perf.r200,add=TRUE,col="black",lty=3,lwd=3)#
legend("bottomright",c(expression(RM[{paste(p,"*")==50}]),expression(RM~(paste(p,"*",sep="")==100)),expression(RM~(paste(p,"*",sep="")==200)),expression(PCA~(paste(p,"*",sep="")==50)),expression(PCA~(paste(p,"*",sep="")==100)),expression(PCA~(paste(p,"*",sep="")==200))),lty=c(1,2,3,1,2,3),lwd=c(3,3,3,3,3,3),col=c("black","black","black","gray55","gray55","gray55"),bty="n",cex=1.7)#
dev.off()
pdf("roc1.pdf",width=10,height=10,paper="special")#
par(cex.axis=1.7,cex.lab=1.7,pty="s")#
plot(pca.perf.r50,col="gray",lty=1,lwd=3,cex=1.7)#
plot(pca.perf.r100,add=TRUE,col="gray",lty=2,lwd=3)#
plot(pca.perf.r200,add=TRUE,col="gray",lty=3,lwd=3)#
plot(rasch.perf.r50,add=TRUE,col="black",lty=1,lwd=3)#
plot(rasch.perf.r100,add=TRUE,col="black",lty=2,lwd=3)#
plot(rasch.perf.r200,add=TRUE,col="black",lty=3,lwd=3)#
legend("bottomright",c(bquote(RM~(paste(p,"*",sep="")==50)),expression(RM~(paste(p,"*",sep="")==100)),expression(RM~(paste(p,"*",sep="")==200)),expression(PCA~(paste(p,"*",sep="")==50)),expression(PCA~(paste(p,"*",sep="")==100)),expression(PCA~(paste(p,"*",sep="")==200))),lty=c(1,2,3,1,2,3),lwd=c(3,3,3,3,3,3),col=c("black","black","black","gray55","gray55","gray55"),bty="n",cex=1.7)#
dev.off()
pdf("roc1.pdf",width=10,height=10,paper="special")#
par(cex.axis=1.7,cex.lab=1.7,pty="s")#
plot(pca.perf.r50,col="gray",lty=1,lwd=3,cex=1.7)#
plot(pca.perf.r100,add=TRUE,col="gray",lty=2,lwd=3)#
plot(pca.perf.r200,add=TRUE,col="gray",lty=3,lwd=3)#
plot(rasch.perf.r50,add=TRUE,col="black",lty=1,lwd=3)#
plot(rasch.perf.r100,add=TRUE,col="black",lty=2,lwd=3)#
plot(rasch.perf.r200,add=TRUE,col="black",lty=3,lwd=3)#
legend("bottomright",c(bquote(RM~[paste(p,"*",sep="")==50]),expression(RM~(paste(p,"*",sep="")==100)),expression(RM~(paste(p,"*",sep="")==200)),expression(PCA~(paste(p,"*",sep="")==50)),expression(PCA~(paste(p,"*",sep="")==100)),expression(PCA~(paste(p,"*",sep="")==200))),lty=c(1,2,3,1,2,3),lwd=c(3,3,3,3,3,3),col=c("black","black","black","gray55","gray55","gray55"),bty="n",cex=1.7)#
dev.off()
pdf("roc1.pdf",width=10,height=10,paper="special")#
par(cex.axis=1.7,cex.lab=1.7,pty="s")#
plot(pca.perf.r50,col="gray",lty=1,lwd=3,cex=1.7)#
plot(pca.perf.r100,add=TRUE,col="gray",lty=2,lwd=3)#
plot(pca.perf.r200,add=TRUE,col="gray",lty=3,lwd=3)#
plot(rasch.perf.r50,add=TRUE,col="black",lty=1,lwd=3)#
plot(rasch.perf.r100,add=TRUE,col="black",lty=2,lwd=3)#
plot(rasch.perf.r200,add=TRUE,col="black",lty=3,lwd=3)#
legend("bottomright",c(bquote(RM~["p*"==50]),expression(RM~(paste(p,"*",sep="")==100)),expression(RM~(paste(p,"*",sep="")==200)),expression(PCA~(paste(p,"*",sep="")==50)),expression(PCA~(paste(p,"*",sep="")==100)),expression(PCA~(paste(p,"*",sep="")==200))),lty=c(1,2,3,1,2,3),lwd=c(3,3,3,3,3,3),col=c("black","black","black","gray55","gray55","gray55"),bty="n",cex=1.7)#
dev.off()
pdf("roc1.pdf",width=10,height=10,paper="special")#
par(cex.axis=1.7,cex.lab=1.7,pty="s")#
plot(pca.perf.r50,col="gray",lty=1,lwd=3,cex=1.7)#
plot(pca.perf.r100,add=TRUE,col="gray",lty=2,lwd=3)#
plot(pca.perf.r200,add=TRUE,col="gray",lty=3,lwd=3)#
plot(rasch.perf.r50,add=TRUE,col="black",lty=1,lwd=3)#
plot(rasch.perf.r100,add=TRUE,col="black",lty=2,lwd=3)#
plot(rasch.perf.r200,add=TRUE,col="black",lty=3,lwd=3)#
legend("bottomright",c(bquote("RM"),expression(RM~(paste(p,"*",sep="")==100)),expression(RM~(paste(p,"*",sep="")==200)),expression(PCA~(paste(p,"*",sep="")==50)),expression(PCA~(paste(p,"*",sep="")==100)),expression(PCA~(paste(p,"*",sep="")==200))),lty=c(1,2,3,1,2,3),lwd=c(3,3,3,3,3,3),col=c("black","black","black","gray55","gray55","gray55"),bty="n",cex=1.7)#
dev.off()
?bquote
pdf("roc1.pdf",width=10,height=10,paper="special")#
par(cex.axis=1.7,cex.lab=1.7,pty="s")#
plot(pca.perf.r50,col="gray",lty=1,lwd=3,cex=1.7)#
plot(pca.perf.r100,add=TRUE,col="gray",lty=2,lwd=3)#
plot(pca.perf.r200,add=TRUE,col="gray",lty=3,lwd=3)#
plot(rasch.perf.r50,add=TRUE,col="black",lty=1,lwd=3)#
plot(rasch.perf.r100,add=TRUE,col="black",lty=2,lwd=3)#
plot(rasch.perf.r200,add=TRUE,col="black",lty=3,lwd=3)#
legend("bottomright",c(bquote("RM"["p"]),expression(RM~(paste(p,"*",sep="")==100)),expression(RM~(paste(p,"*",sep="")==200)),expression(PCA~(paste(p,"*",sep="")==50)),expression(PCA~(paste(p,"*",sep="")==100)),expression(PCA~(paste(p,"*",sep="")==200))),lty=c(1,2,3,1,2,3),lwd=c(3,3,3,3,3,3),col=c("black","black","black","gray55","gray55","gray55"),bty="n",cex=1.7)#
dev.off()
pdf("roc1.pdf",width=10,height=10,paper="special")#
par(cex.axis=1.7,cex.lab=1.7,pty="s")#
plot(pca.perf.r50,col="gray",lty=1,lwd=3,cex=1.7)#
plot(pca.perf.r100,add=TRUE,col="gray",lty=2,lwd=3)#
plot(pca.perf.r200,add=TRUE,col="gray",lty=3,lwd=3)#
plot(rasch.perf.r50,add=TRUE,col="black",lty=1,lwd=3)#
plot(rasch.perf.r100,add=TRUE,col="black",lty=2,lwd=3)#
plot(rasch.perf.r200,add=TRUE,col="black",lty=3,lwd=3)#
legend("bottomright",c(bquote("RM"["p*"]),expression(RM~(paste(p,"*",sep="")==100)),expression(RM~(paste(p,"*",sep="")==200)),expression(PCA~(paste(p,"*",sep="")==50)),expression(PCA~(paste(p,"*",sep="")==100)),expression(PCA~(paste(p,"*",sep="")==200))),lty=c(1,2,3,1,2,3),lwd=c(3,3,3,3,3,3),col=c("black","black","black","gray55","gray55","gray55"),bty="n",cex=1.7)#
dev.off()
pdf("roc1.pdf",width=10,height=10,paper="special")#
par(cex.axis=1.7,cex.lab=1.7,pty="s")#
plot(pca.perf.r50,col="gray",lty=1,lwd=3,cex=1.7)#
plot(pca.perf.r100,add=TRUE,col="gray",lty=2,lwd=3)#
plot(pca.perf.r200,add=TRUE,col="gray",lty=3,lwd=3)#
plot(rasch.perf.r50,add=TRUE,col="black",lty=1,lwd=3)#
plot(rasch.perf.r100,add=TRUE,col="black",lty=2,lwd=3)#
plot(rasch.perf.r200,add=TRUE,col="black",lty=3,lwd=3)#
legend("bottomright",c(bquote("RM"["p*"==50]),expression(RM~(paste(p,"*",sep="")==100)),expression(RM~(paste(p,"*",sep="")==200)),expression(PCA~(paste(p,"*",sep="")==50)),expression(PCA~(paste(p,"*",sep="")==100)),expression(PCA~(paste(p,"*",sep="")==200))),lty=c(1,2,3,1,2,3),lwd=c(3,3,3,3,3,3),col=c("black","black","black","gray55","gray55","gray55"),bty="n",cex=1.7)#
dev.off()
pdf("roc1.pdf",width=10,height=10,paper="special")#
par(cex.axis=1.7,cex.lab=1.7,pty="s")#
plot(pca.perf.r50,col="gray",lty=1,lwd=3,cex=1.7)#
plot(pca.perf.r100,add=TRUE,col="gray",lty=2,lwd=3)#
plot(pca.perf.r200,add=TRUE,col="gray",lty=3,lwd=3)#
plot(rasch.perf.r50,add=TRUE,col="black",lty=1,lwd=3)#
plot(rasch.perf.r100,add=TRUE,col="black",lty=2,lwd=3)#
plot(rasch.perf.r200,add=TRUE,col="black",lty=3,lwd=3)#
legend("bottomright",c(bquote("RM"["p* "==50]),expression(RM~(paste(p,"*",sep="")==100)),expression(RM~(paste(p,"*",sep="")==200)),expression(PCA~(paste(p,"*",sep="")==50)),expression(PCA~(paste(p,"*",sep="")==100)),expression(PCA~(paste(p,"*",sep="")==200))),lty=c(1,2,3,1,2,3),lwd=c(3,3,3,3,3,3),col=c("black","black","black","gray55","gray55","gray55"),bty="n",cex=1.7)#
dev.off()
pdf("roc1.pdf",width=10,height=10,paper="special")#
par(cex.axis=1.7,cex.lab=1.7,pty="s")#
plot(pca.perf.r50,col="gray",lty=1,lwd=3,cex=1.7)#
plot(pca.perf.r100,add=TRUE,col="gray",lty=2,lwd=3)#
plot(pca.perf.r200,add=TRUE,col="gray",lty=3,lwd=3)#
plot(rasch.perf.r50,add=TRUE,col="black",lty=1,lwd=3)#
plot(rasch.perf.r100,add=TRUE,col="black",lty=2,lwd=3)#
plot(rasch.perf.r200,add=TRUE,col="black",lty=3,lwd=3)#
legend("bottomright",c(bquote("RM"["p* "==.50]),expression(RM~(paste(p,"*",sep="")==100)),expression(RM~(paste(p,"*",sep="")==200)),expression(PCA~(paste(p,"*",sep="")==50)),expression(PCA~(paste(p,"*",sep="")==100)),expression(PCA~(paste(p,"*",sep="")==200))),lty=c(1,2,3,1,2,3),lwd=c(3,3,3,3,3,3),col=c("black","black","black","gray55","gray55","gray55"),bty="n",cex=1.7)#
dev.off()
pdf("roc1.pdf",width=10,height=10,paper="special")#
par(cex.axis=1.7,cex.lab=1.7,pty="s")#
plot(pca.perf.r50,col="gray",lty=1,lwd=3,cex=1.7)#
plot(pca.perf.r100,add=TRUE,col="gray",lty=2,lwd=3)#
plot(pca.perf.r200,add=TRUE,col="gray",lty=3,lwd=3)#
plot(rasch.perf.r50,add=TRUE,col="black",lty=1,lwd=3)#
plot(rasch.perf.r100,add=TRUE,col="black",lty=2,lwd=3)#
plot(rasch.perf.r200,add=TRUE,col="black",lty=3,lwd=3)#
legend("bottomright",c(bquote("RM"["p* "==~50]),expression(RM~(paste(p,"*",sep="")==100)),expression(RM~(paste(p,"*",sep="")==200)),expression(PCA~(paste(p,"*",sep="")==50)),expression(PCA~(paste(p,"*",sep="")==100)),expression(PCA~(paste(p,"*",sep="")==200))),lty=c(1,2,3,1,2,3),lwd=c(3,3,3,3,3,3),col=c("black","black","black","gray55","gray55","gray55"),bty="n",cex=1.7)#
dev.off()
pdf("roc1.pdf",width=10,height=10,paper="special")#
par(cex.axis=1.7,cex.lab=1.7,pty="s")#
plot(pca.perf.r50,col="gray",lty=1,lwd=3,cex=1.7)#
plot(pca.perf.r100,add=TRUE,col="gray",lty=2,lwd=3)#
plot(pca.perf.r200,add=TRUE,col="gray",lty=3,lwd=3)#
plot(rasch.perf.r50,add=TRUE,col="black",lty=1,lwd=3)#
plot(rasch.perf.r100,add=TRUE,col="black",lty=2,lwd=3)#
plot(rasch.perf.r200,add=TRUE,col="black",lty=3,lwd=3)#
legend("bottomright",c(bquote("RM"["p*"~==~50]),expression(RM~(paste(p,"*",sep="")==100)),expression(RM~(paste(p,"*",sep="")==200)),expression(PCA~(paste(p,"*",sep="")==50)),expression(PCA~(paste(p,"*",sep="")==100)),expression(PCA~(paste(p,"*",sep="")==200))),lty=c(1,2,3,1,2,3),lwd=c(3,3,3,3,3,3),col=c("black","black","black","gray55","gray55","gray55"),bty="n",cex=1.7)#
dev.off()
pdf("roc1.pdf",width=10,height=10,paper="special")#
par(cex.axis=1.7,cex.lab=1.7,pty="s")#
plot(pca.perf.r50,col="gray",lty=1,lwd=3,cex=1.7)#
plot(pca.perf.r100,add=TRUE,col="gray",lty=2,lwd=3)#
plot(pca.perf.r200,add=TRUE,col="gray",lty=3,lwd=3)#
plot(rasch.perf.r50,add=TRUE,col="black",lty=1,lwd=3)#
plot(rasch.perf.r100,add=TRUE,col="black",lty=2,lwd=3)#
plot(rasch.perf.r200,add=TRUE,col="black",lty=3,lwd=3)#
legend("bottomright",c(bquote("RM"["p* "==" 50"]),expression(RM~(paste(p,"*",sep="")==100)),expression(RM~(paste(p,"*",sep="")==200)),expression(PCA~(paste(p,"*",sep="")==50)),expression(PCA~(paste(p,"*",sep="")==100)),expression(PCA~(paste(p,"*",sep="")==200))),lty=c(1,2,3,1,2,3),lwd=c(3,3,3,3,3,3),col=c("black","black","black","gray55","gray55","gray55"),bty="n",cex=1.7)#
dev.off()
pdf("roc1.pdf",width=10,height=10,paper="special")#
par(cex.axis=1.7,cex.lab=1.7,pty="s")#
plot(pca.perf.r50,col="gray",lty=1,lwd=3,cex=1.7)#
plot(pca.perf.r100,add=TRUE,col="gray",lty=2,lwd=3)#
plot(pca.perf.r200,add=TRUE,col="gray",lty=3,lwd=3)#
plot(rasch.perf.r50,add=TRUE,col="black",lty=1,lwd=3)#
plot(rasch.perf.r100,add=TRUE,col="black",lty=2,lwd=3)#
plot(rasch.perf.r200,add=TRUE,col="black",lty=3,lwd=3)#
legend("bottomright",c(bquote("RM"["p*~"=="~50"]),expression(RM~(paste(p,"*",sep="")==100)),expression(RM~(paste(p,"*",sep="")==200)),expression(PCA~(paste(p,"*",sep="")==50)),expression(PCA~(paste(p,"*",sep="")==100)),expression(PCA~(paste(p,"*",sep="")==200))),lty=c(1,2,3,1,2,3),lwd=c(3,3,3,3,3,3),col=c("black","black","black","gray55","gray55","gray55"),bty="n",cex=1.7)#
dev.off()
pdf("roc1.pdf",width=10,height=10,paper="special")#
par(cex.axis=1.7,cex.lab=1.7,pty="s")#
plot(pca.perf.r50,col="gray",lty=1,lwd=3,cex=1.7)#
plot(pca.perf.r100,add=TRUE,col="gray",lty=2,lwd=3)#
plot(pca.perf.r200,add=TRUE,col="gray",lty=3,lwd=3)#
plot(rasch.perf.r50,add=TRUE,col="black",lty=1,lwd=3)#
plot(rasch.perf.r100,add=TRUE,col="black",lty=2,lwd=3)#
plot(rasch.perf.r200,add=TRUE,col="black",lty=3,lwd=3)#
legend("bottomright",c(bquote("RM"["p*\."=="~50"]),expression(RM~(paste(p,"*",sep="")==100)),expression(RM~(paste(p,"*",sep="")==200)),expression(PCA~(paste(p,"*",sep="")==50)),expression(PCA~(paste(p,"*",sep="")==100)),expression(PCA~(paste(p,"*",sep="")==200))),lty=c(1,2,3,1,2,3),lwd=c(3,3,3,3,3,3),col=c("black","black","black","gray55","gray55","gray55"),bty="n",cex=1.7)#
dev.off()
pdf("roc1.pdf",width=10,height=10,paper="special")#
par(cex.axis=1.7,cex.lab=1.7,pty="s")#
plot(pca.perf.r50,col="gray",lty=1,lwd=3,cex=1.7)#
plot(pca.perf.r100,add=TRUE,col="gray",lty=2,lwd=3)#
plot(pca.perf.r200,add=TRUE,col="gray",lty=3,lwd=3)#
plot(rasch.perf.r50,add=TRUE,col="black",lty=1,lwd=3)#
plot(rasch.perf.r100,add=TRUE,col="black",lty=2,lwd=3)#
plot(rasch.perf.r200,add=TRUE,col="black",lty=3,lwd=3)#
legend("bottomright",c(bquote("RM"["p*"\.=="~50"]),expression(RM~(paste(p,"*",sep="")==100)),expression(RM~(paste(p,"*",sep="")==200)),expression(PCA~(paste(p,"*",sep="")==50)),expression(PCA~(paste(p,"*",sep="")==100)),expression(PCA~(paste(p,"*",sep="")==200))),lty=c(1,2,3,1,2,3),lwd=c(3,3,3,3,3,3),col=c("black","black","black","gray55","gray55","gray55"),bty="n",cex=1.7)#
dev.off()
pdf("roc1.pdf",width=10,height=10,paper="special")#
par(cex.axis=1.7,cex.lab=1.7,pty="s")#
plot(pca.perf.r50,col="gray",lty=1,lwd=3,cex=1.7)#
plot(pca.perf.r100,add=TRUE,col="gray",lty=2,lwd=3)#
plot(pca.perf.r200,add=TRUE,col="gray",lty=3,lwd=3)#
plot(rasch.perf.r50,add=TRUE,col="black",lty=1,lwd=3)#
plot(rasch.perf.r100,add=TRUE,col="black",lty=2,lwd=3)#
plot(rasch.perf.r200,add=TRUE,col="black",lty=3,lwd=3)#
legend("bottomright",c(bquote("RM"["p*"=="50"]),expression(RM~(paste(p,"*",sep="")==100)),expression(RM~(paste(p,"*",sep="")==200)),expression(PCA~(paste(p,"*",sep="")==50)),expression(PCA~(paste(p,"*",sep="")==100)),expression(PCA~(paste(p,"*",sep="")==200))),lty=c(1,2,3,1,2,3),lwd=c(3,3,3,3,3,3),col=c("black","black","black","gray55","gray55","gray55"),bty="n",cex=1.7)#
dev.off()
pdf("roc1.pdf",width=10,height=10,paper="special")#
par(cex.axis=1.7,cex.lab=1.7,pty="s")#
plot(pca.perf.r50,col="gray",lty=1,lwd=3,cex=1.7)#
plot(pca.perf.r100,add=TRUE,col="gray",lty=2,lwd=3)#
plot(pca.perf.r200,add=TRUE,col="gray",lty=3,lwd=3)#
plot(rasch.perf.r50,add=TRUE,col="black",lty=1,lwd=3)#
plot(rasch.perf.r100,add=TRUE,col="black",lty=2,lwd=3)#
plot(rasch.perf.r200,add=TRUE,col="black",lty=3,lwd=3)#
legend("bottomright",c(bquote("RM"["p*"="50"]),expression(RM~(paste(p,"*",sep="")==100)),expression(RM~(paste(p,"*",sep="")==200)),expression(PCA~(paste(p,"*",sep="")==50)),expression(PCA~(paste(p,"*",sep="")==100)),expression(PCA~(paste(p,"*",sep="")==200))),lty=c(1,2,3,1,2,3),lwd=c(3,3,3,3,3,3),col=c("black","black","black","gray55","gray55","gray55"),bty="n",cex=1.7)#
dev.off()
pdf("roc1.pdf",width=10,height=10,paper="special")#
par(cex.axis=1.7,cex.lab=1.7,pty="s")#
plot(pca.perf.r50,col="gray",lty=1,lwd=3,cex=1.7)#
plot(pca.perf.r100,add=TRUE,col="gray",lty=2,lwd=3)#
plot(pca.perf.r200,add=TRUE,col="gray",lty=3,lwd=3)#
plot(rasch.perf.r50,add=TRUE,col="black",lty=1,lwd=3)#
plot(rasch.perf.r100,add=TRUE,col="black",lty=2,lwd=3)#
plot(rasch.perf.r200,add=TRUE,col="black",lty=3,lwd=3)#
legend("bottomright",c(bquote("RM"["p*".=="50"]),expression(RM~(paste(p,"*",sep="")==100)),expression(RM~(paste(p,"*",sep="")==200)),expression(PCA~(paste(p,"*",sep="")==50)),expression(PCA~(paste(p,"*",sep="")==100)),expression(PCA~(paste(p,"*",sep="")==200))),lty=c(1,2,3,1,2,3),lwd=c(3,3,3,3,3,3),col=c("black","black","black","gray55","gray55","gray55"),bty="n",cex=1.7)#
dev.off()
pdf("roc1.pdf",width=10,height=10,paper="special")#
par(cex.axis=1.7,cex.lab=1.7,pty="s")#
plot(pca.perf.r50,col="gray",lty=1,lwd=3,cex=1.7)#
plot(pca.perf.r100,add=TRUE,col="gray",lty=2,lwd=3)#
plot(pca.perf.r200,add=TRUE,col="gray",lty=3,lwd=3)#
plot(rasch.perf.r50,add=TRUE,col="black",lty=1,lwd=3)#
plot(rasch.perf.r100,add=TRUE,col="black",lty=2,lwd=3)#
plot(rasch.perf.r200,add=TRUE,col="black",lty=3,lwd=3)#
legend("bottomright",c(bquote("RM"["p*"~=="50"]),expression(RM~(paste(p,"*",sep="")==100)),expression(RM~(paste(p,"*",sep="")==200)),expression(PCA~(paste(p,"*",sep="")==50)),expression(PCA~(paste(p,"*",sep="")==100)),expression(PCA~(paste(p,"*",sep="")==200))),lty=c(1,2,3,1,2,3),lwd=c(3,3,3,3,3,3),col=c("black","black","black","gray55","gray55","gray55"),bty="n",cex=1.7)#
dev.off()
pdf("roc1.pdf",width=10,height=10,paper="special")#
par(cex.axis=1.7,cex.lab=1.7,pty="s")#
plot(pca.perf.r50,col="gray",lty=1,lwd=3,cex=1.7)#
plot(pca.perf.r100,add=TRUE,col="gray",lty=2,lwd=3)#
plot(pca.perf.r200,add=TRUE,col="gray",lty=3,lwd=3)#
plot(rasch.perf.r50,add=TRUE,col="black",lty=1,lwd=3)#
plot(rasch.perf.r100,add=TRUE,col="black",lty=2,lwd=3)#
plot(rasch.perf.r200,add=TRUE,col="black",lty=3,lwd=3)#
legend("bottomright",c(bquote("RM"["p*"~~=="50"]),expression(RM~(paste(p,"*",sep="")==100)),expression(RM~(paste(p,"*",sep="")==200)),expression(PCA~(paste(p,"*",sep="")==50)),expression(PCA~(paste(p,"*",sep="")==100)),expression(PCA~(paste(p,"*",sep="")==200))),lty=c(1,2,3,1,2,3),lwd=c(3,3,3,3,3,3),col=c("black","black","black","gray55","gray55","gray55"),bty="n",cex=1.7)#
dev.off()
pdf("roc1.pdf",width=10,height=10,paper="special")#
par(cex.axis=1.7,cex.lab=1.7,pty="s")#
plot(pca.perf.r50,col="gray",lty=1,lwd=3,cex=1.7)#
plot(pca.perf.r100,add=TRUE,col="gray",lty=2,lwd=3)#
plot(pca.perf.r200,add=TRUE,col="gray",lty=3,lwd=3)#
plot(rasch.perf.r50,add=TRUE,col="black",lty=1,lwd=3)#
plot(rasch.perf.r100,add=TRUE,col="black",lty=2,lwd=3)#
plot(rasch.perf.r200,add=TRUE,col="black",lty=3,lwd=3)#
legend("bottomright",c(bquote("RM"["p* "==" 50"]),expression(RM~(paste(p,"*",sep="")==100)),expression(RM~(paste(p,"*",sep="")==200)),expression(PCA~(paste(p,"*",sep="")==50)),expression(PCA~(paste(p,"*",sep="")==100)),expression(PCA~(paste(p,"*",sep="")==200))),lty=c(1,2,3,1,2,3),lwd=c(3,3,3,3,3,3),col=c("black","black","black","gray55","gray55","gray55"),bty="n",cex=1.7)#
dev.off()
pdf("roc1.pdf",width=10,height=10,paper="special")#
par(cex.axis=1.7,cex.lab=1.7,pty="s")#
plot(pca.perf.r50,col="gray",lty=1,lwd=3,cex=1.7)#
plot(pca.perf.r100,add=TRUE,col="gray",lty=2,lwd=3)#
plot(pca.perf.r200,add=TRUE,col="gray",lty=3,lwd=3)#
plot(rasch.perf.r50,add=TRUE,col="black",lty=1,lwd=3)#
plot(rasch.perf.r100,add=TRUE,col="black",lty=2,lwd=3)#
plot(rasch.perf.r200,add=TRUE,col="black",lty=3,lwd=3)#
legend("bottomright",c(bquote("RM"["p*==50"]),expression(RM~(paste(p,"*",sep="")==100)),expression(RM~(paste(p,"*",sep="")==200)),expression(PCA~(paste(p,"*",sep="")==50)),expression(PCA~(paste(p,"*",sep="")==100)),expression(PCA~(paste(p,"*",sep="")==200))),lty=c(1,2,3,1,2,3),lwd=c(3,3,3,3,3,3),col=c("black","black","black","gray55","gray55","gray55"),bty="n",cex=1.7)#
dev.off()
pdf("roc1.pdf",width=10,height=10,paper="special")#
par(cex.axis=1.7,cex.lab=1.7,pty="s")#
plot(pca.perf.r50,col="gray",lty=1,lwd=3,cex=1.7)#
plot(pca.perf.r100,add=TRUE,col="gray",lty=2,lwd=3)#
plot(pca.perf.r200,add=TRUE,col="gray",lty=3,lwd=3)#
plot(rasch.perf.r50,add=TRUE,col="black",lty=1,lwd=3)#
plot(rasch.perf.r100,add=TRUE,col="black",lty=2,lwd=3)#
plot(rasch.perf.r200,add=TRUE,col="black",lty=3,lwd=3)#
legend("bottomright",c(bquote("RM"["p*=50"]),expression(RM~(paste(p,"*",sep="")==100)),expression(RM~(paste(p,"*",sep="")==200)),expression(PCA~(paste(p,"*",sep="")==50)),expression(PCA~(paste(p,"*",sep="")==100)),expression(PCA~(paste(p,"*",sep="")==200))),lty=c(1,2,3,1,2,3),lwd=c(3,3,3,3,3,3),col=c("black","black","black","gray55","gray55","gray55"),bty="n",cex=1.7)#
dev.off()
pdf("roc1.pdf",width=10,height=10,paper="special")#
par(cex.axis=1.7,cex.lab=1.7,pty="s")#
plot(pca.perf.r50,col="gray",lty=1,lwd=3,cex=1.7)#
plot(pca.perf.r100,add=TRUE,col="gray",lty=2,lwd=3)#
plot(pca.perf.r200,add=TRUE,col="gray",lty=3,lwd=3)#
plot(rasch.perf.r50,add=TRUE,col="black",lty=1,lwd=3)#
plot(rasch.perf.r100,add=TRUE,col="black",lty=2,lwd=3)#
plot(rasch.perf.r200,add=TRUE,col="black",lty=3,lwd=3)#
legend("bottomright",c(bquote("RM"["p* = 50"]),expression(RM~(paste(p,"*",sep="")==100)),expression(RM~(paste(p,"*",sep="")==200)),expression(PCA~(paste(p,"*",sep="")==50)),expression(PCA~(paste(p,"*",sep="")==100)),expression(PCA~(paste(p,"*",sep="")==200))),lty=c(1,2,3,1,2,3),lwd=c(3,3,3,3,3,3),col=c("black","black","black","gray55","gray55","gray55"),bty="n",cex=1.7)#
dev.off()
pdf("roc1.pdf",width=10,height=10,paper="special")#
par(cex.axis=1.7,cex.lab=1.7,pty="s")#
plot(pca.perf.r50,col="gray",lty=1,lwd=3,cex=1.7)#
plot(pca.perf.r100,add=TRUE,col="gray",lty=2,lwd=3)#
plot(pca.perf.r200,add=TRUE,col="gray",lty=3,lwd=3)#
plot(rasch.perf.r50,add=TRUE,col="black",lty=1,lwd=3)#
plot(rasch.perf.r100,add=TRUE,col="black",lty=2,lwd=3)#
plot(rasch.perf.r200,add=TRUE,col="black",lty=3,lwd=3)#
legend("bottomright",c(expression(RM~(paste(p,"*",sep="")==50)),expression(RM~(paste(p,"*",sep="")==100)),expression(RM~(paste(p,"*",sep="")==200)),expression(PCA~(paste(p,"*",sep="")==50)),expression(PCA~(paste(p,"*",sep="")==100)),expression(PCA~(paste(p,"*",sep="")==200))),lty=c(1,2,3,1,2,3),lwd=c(3,3,3,3,3,3),col=c("black","black","black","gray55","gray55","gray55"),bty="n",cex=1.7)#
dev.off()
pdf("roc1.pdf",width=20,height=10,paper="special")#
par(cex.axis=1.7,cex.lab=1.7,pty="s",mfrow=c(2,1))#
plot(pca.perf.r50,col="gray",lty=1,lwd=3,cex=1.7)#
plot(pca.perf.r100,add=TRUE,col="gray",lty=2,lwd=3)#
plot(pca.perf.r200,add=TRUE,col="gray",lty=3,lwd=3)#
plot(rasch.perf.r50,add=TRUE,col="black",lty=1,lwd=3)#
plot(rasch.perf.r100,add=TRUE,col="black",lty=2,lwd=3)#
plot(rasch.perf.r200,add=TRUE,col="black",lty=3,lwd=3)#
legend("bottomright",c(expression(RM~(paste(p,"*",sep="")==50)),expression(RM~(paste(p,"*",sep="")==100)),expression(RM~(paste(p,"*",sep="")==200)),expression(PCA~(paste(p,"*",sep="")==50)),expression(PCA~(paste(p,"*",sep="")==100)),expression(PCA~(paste(p,"*",sep="")==200))),lty=c(1,2,3,1,2,3),lwd=c(3,3,3,3,3,3),col=c("black","black","black","gray55","gray55","gray55"),bty="n",cex=1.7)#
#dev.off()#
#
plot(pca.perf.p50,col="black",lty=1,lwd=1)#
plot(pca.perf.p100,add=TRUE,col="black",lty=2,lwd=1)#
plot(pca.perf.p200,add=TRUE,col="black",lty=3,lwd=1)#
plot(rasch.perf.p50,add=TRUE,col="black",lty=1,lwd=3)#
plot(rasch.perf.p100,add=TRUE,col="black",lty=2,lwd=3)#
plot(rasch.perf.p200,add=TRUE,col="black",lty=3,lwd=3)#
legend("bottomright",c(expression(RM~(paste(p,"*",sep="")==50)),expression(RM~(paste(p,"*",sep="")==100)),expression(RM~(paste(p,"*",sep="")==200)),expression(PCA~(paste(p,"*",sep="")==50)),expression(PCA~(paste(p,"*",sep="")==100)),expression(PCA~(paste(p,"*",sep="")==200))),lty=c(1,2,3,1,2,3),lwd=c(2,2,2,1,1,1),bty="n")#
dev.off()
pdf("roc1.pdf",width=20,height=10,paper="special")#
par(cex.axis=1.7,cex.lab=1.7,pty="s",mfrow=c(1,2))#
plot(pca.perf.r50,col="gray",lty=1,lwd=3,cex=1.7)#
plot(pca.perf.r100,add=TRUE,col="gray",lty=2,lwd=3)#
plot(pca.perf.r200,add=TRUE,col="gray",lty=3,lwd=3)#
plot(rasch.perf.r50,add=TRUE,col="black",lty=1,lwd=3)#
plot(rasch.perf.r100,add=TRUE,col="black",lty=2,lwd=3)#
plot(rasch.perf.r200,add=TRUE,col="black",lty=3,lwd=3)#
legend("bottomright",c(expression(RM~(paste(p,"*",sep="")==50)),expression(RM~(paste(p,"*",sep="")==100)),expression(RM~(paste(p,"*",sep="")==200)),expression(PCA~(paste(p,"*",sep="")==50)),expression(PCA~(paste(p,"*",sep="")==100)),expression(PCA~(paste(p,"*",sep="")==200))),lty=c(1,2,3,1,2,3),lwd=c(3,3,3,3,3,3),col=c("black","black","black","gray55","gray55","gray55"),bty="n",cex=1.7)#
#dev.off()#
#
plot(pca.perf.p50,col="black",lty=1,lwd=1)#
plot(pca.perf.p100,add=TRUE,col="black",lty=2,lwd=1)#
plot(pca.perf.p200,add=TRUE,col="black",lty=3,lwd=1)#
plot(rasch.perf.p50,add=TRUE,col="black",lty=1,lwd=3)#
plot(rasch.perf.p100,add=TRUE,col="black",lty=2,lwd=3)#
plot(rasch.perf.p200,add=TRUE,col="black",lty=3,lwd=3)#
legend("bottomright",c(expression(RM~(paste(p,"*",sep="")==50)),expression(RM~(paste(p,"*",sep="")==100)),expression(RM~(paste(p,"*",sep="")==200)),expression(PCA~(paste(p,"*",sep="")==50)),expression(PCA~(paste(p,"*",sep="")==100)),expression(PCA~(paste(p,"*",sep="")==200))),lty=c(1,2,3,1,2,3),lwd=c(2,2,2,1,1,1),bty="n")#
dev.off()
pdf("roc1.pdf",width=20,height=10,paper="special")#
par(cex.axis=1.7,cex.lab=1.7,pty="s",mfrow=c(1,2))#
plot(pca.perf.r50,col="gray",lty=1,lwd=3,cex=1.7)#
plot(pca.perf.r100,add=TRUE,col="gray",lty=2,lwd=3)#
plot(pca.perf.r200,add=TRUE,col="gray",lty=3,lwd=3)#
plot(rasch.perf.r50,add=TRUE,col="black",lty=1,lwd=3)#
plot(rasch.perf.r100,add=TRUE,col="black",lty=2,lwd=3)#
plot(rasch.perf.r200,add=TRUE,col="black",lty=3,lwd=3)#
legend("bottomright",c(expression(RM~(paste(p,"*",sep="")==50)),expression(RM~(paste(p,"*",sep="")==100)),expression(RM~(paste(p,"*",sep="")==200)),expression(PCA~(paste(p,"*",sep="")==50)),expression(PCA~(paste(p,"*",sep="")==100)),expression(PCA~(paste(p,"*",sep="")==200))),lty=c(1,2,3,1,2,3),lwd=c(3,3,3,3,3,3),col=c("black","black","black","gray55","gray55","gray55"),bty="n",cex=1.7)#
#
plot(pca.perf.p50,col="black",lty=1,lwd=1)#
plot(pca.perf.p100,add=TRUE,col="black",lty=2,lwd=1)#
plot(pca.perf.p200,add=TRUE,col="black",lty=3,lwd=1)#
plot(rasch.perf.p50,add=TRUE,col="black",lty=1,lwd=3)#
plot(rasch.perf.p100,add=TRUE,col="black",lty=2,lwd=3)#
plot(rasch.perf.p200,add=TRUE,col="black",lty=3,lwd=3)#
legend("bottomright",c(expression(RM~(paste(p,"*",sep="")==50)),expression(RM~(paste(p,"*",sep="")==100)),expression(RM~(paste(p,"*",sep="")==200)),expression(PCA~(paste(p,"*",sep="")==50)),expression(PCA~(paste(p,"*",sep="")==100)),expression(PCA~(paste(p,"*",sep="")==200))),lty=c(1,2,3,1,2,3),lwd=c(3,3,3,3,3,3),col=c("black","black","black","gray55","gray55","gray55"),bty="n",cex=1.7)#
dev.off()
pdf("roc1.pdf",width=20,height=10,paper="special")#
par(cex.axis=1.7,cex.lab=1.7,pty="s",mfrow=c(1,2))#
plot(pca.perf.r50,col="gray",lty=1,lwd=3,cex=1.7)#
plot(pca.perf.r100,add=TRUE,col="gray",lty=2,lwd=3)#
plot(pca.perf.r200,add=TRUE,col="gray",lty=3,lwd=3)#
plot(rasch.perf.r50,add=TRUE,col="black",lty=1,lwd=3)#
plot(rasch.perf.r100,add=TRUE,col="black",lty=2,lwd=3)#
plot(rasch.perf.r200,add=TRUE,col="black",lty=3,lwd=3)#
legend("bottomright",c(expression(RM~(paste(p,"*",sep="")==50)),expression(RM~(paste(p,"*",sep="")==100)),expression(RM~(paste(p,"*",sep="")==200)),expression(PCA~(paste(p,"*",sep="")==50)),expression(PCA~(paste(p,"*",sep="")==100)),expression(PCA~(paste(p,"*",sep="")==200))),lty=c(1,2,3,1,2,3),lwd=c(3,3,3,3,3,3),col=c("black","black","black","gray55","gray55","gray55"),bty="n",cex=1.7)#
#
plot(pca.perf.p50,col="black",lty=1,lwd=3)#
plot(pca.perf.p100,add=TRUE,col="black",lty=2,lwd=3)#
plot(pca.perf.p200,add=TRUE,col="black",lty=3,lwd=3)#
plot(rasch.perf.p50,add=TRUE,col="black",lty=1,lwd=3)#
plot(rasch.perf.p100,add=TRUE,col="black",lty=2,lwd=3)#
plot(rasch.perf.p200,add=TRUE,col="black",lty=3,lwd=3)#
legend("bottomright",c(expression(RM~(paste(p,"*",sep="")==50)),expression(RM~(paste(p,"*",sep="")==100)),expression(RM~(paste(p,"*",sep="")==200)),expression(PCA~(paste(p,"*",sep="")==50)),expression(PCA~(paste(p,"*",sep="")==100)),expression(PCA~(paste(p,"*",sep="")==200))),lty=c(1,2,3,1,2,3),lwd=c(3,3,3,3,3,3),col=c("black","black","black","gray55","gray55","gray55"),bty="n",cex=1.7)#
dev.off()
pdf("roc1_cl.pdf")#
par(cex.axis=1.7,cex.lab=1.7,pty="s",mfrow=c(1,2))#
plot(pca.perf.r50,col="blue",lty=1,lwd=3,cex=1.7)#
plot(pca.perf.r100,add=TRUE,col="blue",lty=2,lwd=3)#
plot(pca.perf.r200,add=TRUE,col="blue",lty=3,lwd=3)#
plot(rasch.perf.r50,add=TRUE,col="red",lty=1,lwd=3)#
plot(rasch.perf.r100,add=TRUE,col="red",lty=2,lwd=3)#
plot(rasch.perf.r200,add=TRUE,col="red",lty=3,lwd=3)#
legend("bottomright",c(expression(RM~(paste(p,"*",sep="")==50)),expression(RM~(paste(p,"*",sep="")==100)),expression(RM~(paste(p,"*",sep="")==200)),expression(PCA~(paste(p,"*",sep="")==50)),expression(PCA~(paste(p,"*",sep="")==100)),expression(PCA~(paste(p,"*",sep="")==200))),lty=c(1,2,3,1,2,3),lwd=c(3,3,3,3,3,3),col=c("red","red","red","blue","blue","blue"),bty="n",cex=1.7)#
plot(pca.perf.p50,col="blue",lty=1,lwd=3)#
plot(pca.perf.p100,add=TRUE,col="blue",lty=2,lwd=3)#
plot(pca.perf.p200,add=TRUE,col="blue",lty=3,lwd=3)#
plot(rasch.perf.p50,add=TRUE,col="red",lty=1,lwd=3)#
plot(rasch.perf.p100,add=TRUE,col="red",lty=2,lwd=3)#
plot(rasch.perf.p200,add=TRUE,col="red",lty=3,lwd=3)#
legend("bottomright",c(expression(RM~(paste(p,"*",sep="")==50)),expression(RM~(paste(p,"*",sep="")==100)),expression(RM~(paste(p,"*",sep="")==200)),expression(PCA~(paste(p,"*",sep="")==50)),expression(PCA~(paste(p,"*",sep="")==100)),expression(PCA~(paste(p,"*",sep="")==200))),lty=c(1,2,3,1,2,3),lwd=c(3,3,3,3,3,3),col=c("red","red","red","blue","blue","blue"),bty="n",cex=1.7)#
dev.off()
?pdf
pdf("roc1_cl.pdf",width=20,height=10,paper="special")#
par(cex.axis=1.7,cex.lab=1.7,pty="s",mfrow=c(1,2))#
plot(pca.perf.r50,col="blue",lty=1,lwd=3,cex=1.7)#
plot(pca.perf.r100,add=TRUE,col="blue",lty=2,lwd=3)#
plot(pca.perf.r200,add=TRUE,col="blue",lty=3,lwd=3)#
plot(rasch.perf.r50,add=TRUE,col="red",lty=1,lwd=3)#
plot(rasch.perf.r100,add=TRUE,col="red",lty=2,lwd=3)#
plot(rasch.perf.r200,add=TRUE,col="red",lty=3,lwd=3)#
legend("bottomright",c(expression(RM~(paste(p,"*",sep="")==50)),expression(RM~(paste(p,"*",sep="")==100)),expression(RM~(paste(p,"*",sep="")==200)),expression(PCA~(paste(p,"*",sep="")==50)),expression(PCA~(paste(p,"*",sep="")==100)),expression(PCA~(paste(p,"*",sep="")==200))),lty=c(1,2,3,1,2,3),lwd=c(3,3,3,3,3,3),col=c("red","red","red","blue","blue","blue"),bty="n",cex=1.7)#
plot(pca.perf.p50,col="blue",lty=1,lwd=3)#
plot(pca.perf.p100,add=TRUE,col="blue",lty=2,lwd=3)#
plot(pca.perf.p200,add=TRUE,col="blue",lty=3,lwd=3)#
plot(rasch.perf.p50,add=TRUE,col="red",lty=1,lwd=3)#
plot(rasch.perf.p100,add=TRUE,col="red",lty=2,lwd=3)#
plot(rasch.perf.p200,add=TRUE,col="red",lty=3,lwd=3)#
legend("bottomright",c(expression(RM~(paste(p,"*",sep="")==50)),expression(RM~(paste(p,"*",sep="")==100)),expression(RM~(paste(p,"*",sep="")==200)),expression(PCA~(paste(p,"*",sep="")==50)),expression(PCA~(paste(p,"*",sep="")==100)),expression(PCA~(paste(p,"*",sep="")==200))),lty=c(1,2,3,1,2,3),lwd=c(3,3,3,3,3,3),col=c("red","red","red","blue","blue","blue"),bty="n",cex=1.7)#
dev.off()
rasch.perf.p200 <- performance(rasch.pred.p200,"tpr","fpr")
pdf("roc1_cl.pdf",width=20,height=10,paper="special")#
par(cex.axis=1.7,cex.lab=1.7,pty="s",mfrow=c(1,2))#
plot(pca.perf.r50,col="blue",lty=1,lwd=3,cex=1.7)#
plot(pca.perf.r100,add=TRUE,col="blue",lty=2,lwd=3)#
plot(pca.perf.r200,add=TRUE,col="blue",lty=3,lwd=3)#
plot(rasch.perf.r50,add=TRUE,col="red",lty=1,lwd=3)#
plot(rasch.perf.r100,add=TRUE,col="red",lty=2,lwd=3)#
plot(rasch.perf.r200,add=TRUE,col="red",lty=3,lwd=3)#
legend("bottomright",c(expression(RM~(paste(p,"*",sep="")==50)),expression(RM~(paste(p,"*",sep="")==100)),expression(RM~(paste(p,"*",sep="")==200)),expression(PCA~(paste(p,"*",sep="")==50)),expression(PCA~(paste(p,"*",sep="")==100)),expression(PCA~(paste(p,"*",sep="")==200))),lty=c(1,2,3,1,2,3),lwd=c(3,3,3,3,3,3),col=c("red","red","red","blue","blue","blue"),bty="n",cex=1.7)#
plot(pca.perf.p50,col="blue",lty=1,lwd=3)#
plot(pca.perf.p100,add=TRUE,col="blue",lty=2,lwd=3)#
plot(pca.perf.p200,add=TRUE,col="blue",lty=3,lwd=3)#
plot(rasch.perf.p50,add=TRUE,col="red",lty=1,lwd=3)#
plot(rasch.perf.p100,add=TRUE,col="red",lty=2,lwd=3)#
plot(rasch.perf.p200,add=TRUE,col="red",lty=3,lwd=3)#
legend("bottomright",c(expression(RM~(paste(p,"*",sep="")==50)),expression(RM~(paste(p,"*",sep="")==100)),expression(RM~(paste(p,"*",sep="")==200)),expression(PCA~(paste(p,"*",sep="")==50)),expression(PCA~(paste(p,"*",sep="")==100)),expression(PCA~(paste(p,"*",sep="")==200))),lty=c(1,2,3,1,2,3),lwd=c(3,3,3,3,3,3),col=c("red","red","red","blue","blue","blue"),bty="n",cex=1.7)#
dev.off()
rasch.perf.p200 <- performance(rasch.pred.p200,"tpr","fpr")#
#
#
pdf("roc1_bw.pdf",width=20,height=10,paper="special")#
par(cex.axis=1.7,cex.lab=1.7,pty="s",mfrow=c(1,2))#
plot(pca.perf.r50,col="gray",lty=1,lwd=3,cex=1.7)#
plot(pca.perf.r100,add=TRUE,col="gray",lty=2,lwd=3)#
plot(pca.perf.r200,add=TRUE,col="gray",lty=3,lwd=3)#
plot(rasch.perf.r50,add=TRUE,col="black",lty=1,lwd=3)#
plot(rasch.perf.r100,add=TRUE,col="black",lty=2,lwd=3)#
plot(rasch.perf.r200,add=TRUE,col="black",lty=3,lwd=3)#
legend("bottomright",c(expression(RM~(paste(p,"*",sep="")==50)),expression(RM~(paste(p,"*",sep="")==100)),expression(RM~(paste(p,"*",sep="")==200)),expression(PCA~(paste(p,"*",sep="")==50)),expression(PCA~(paste(p,"*",sep="")==100)),expression(PCA~(paste(p,"*",sep="")==200))),lty=c(1,2,3,1,2,3),lwd=c(3,3,3,3,3,3),col=c("black","black","black","gray55","gray55","gray55"),bty="n",cex=1.7)#
plot(pca.perf.p50,col="black",lty=1,lwd=3)#
plot(pca.perf.p100,add=TRUE,col="black",lty=2,lwd=3)#
plot(pca.perf.p200,add=TRUE,col="black",lty=3,lwd=3)#
plot(rasch.perf.p50,add=TRUE,col="black",lty=1,lwd=3)#
plot(rasch.perf.p100,add=TRUE,col="black",lty=2,lwd=3)#
plot(rasch.perf.p200,add=TRUE,col="black",lty=3,lwd=3)#
legend("bottomright",c(expression(RM~(paste(p,"*",sep="")==50)),expression(RM~(paste(p,"*",sep="")==100)),expression(RM~(paste(p,"*",sep="")==200)),expression(PCA~(paste(p,"*",sep="")==50)),expression(PCA~(paste(p,"*",sep="")==100)),expression(PCA~(paste(p,"*",sep="")==200))),lty=c(1,2,3,1,2,3),lwd=c(3,3,3,3,3,3),col=c("black","black","black","gray55","gray55","gray55"),bty="n",cex=1.7)#
dev.off()
str(rasch.perf.p200)
unlist(slot(rasch.perf.p200,x.values))
unlist(slot(rasch.perf.p200,"x.values"))
pdf("roc1_bw.pdf",width=20,height=10,paper="special")#
par(cex.axis=1.7,cex.lab=1.7,pty="s",mfrow=c(1,2))#
plot(pca.perf.r50,col="gray",lty=1,lwd=3,cex=1.7)#
plot(pca.perf.r100,add=TRUE,col="gray",lty=2,lwd=3)#
plot(pca.perf.r200,add=TRUE,col="gray",lty=3,lwd=3)#
plot(rasch.perf.r50,add=TRUE,col="black",lty=1,lwd=3)#
plot(rasch.perf.r100,add=TRUE,col="black",lty=2,lwd=3)#
plot(rasch.perf.r200,add=TRUE,col="black",lty=3,lwd=3)#
legend("bottomright",c(expression(RM~(paste(p,"*",sep="")==50)),expression(RM~(paste(p,"*",sep="")==100)),expression(RM~(paste(p,"*",sep="")==200)),expression(PCA~(paste(p,"*",sep="")==50)),expression(PCA~(paste(p,"*",sep="")==100)),expression(PCA~(paste(p,"*",sep="")==200))),lty=c(1,2,3,1,2,3),lwd=c(3,3,3,3,3,3),col=c("black","black","black","gray55","gray55","gray55"),bty="n",cex=1.7)#
plot(pca.perf.p50,col="black",lty=1,lwd=3)#
plot(pca.perf.p100,add=TRUE,col="black",lty=2,lwd=3)#
plot(pca.perf.p200,add=TRUE,col="black",lty=3,lwd=3)#
plot(rasch.perf.p50,add=TRUE,col="black",lty=1,lwd=3)#
plot(rasch.perf.p100,add=TRUE,col="black",lty=2,lwd=3)#
plot(unlist(slot(rasch.perf.p200,"x.values")),unlist(slot(rasch.perf.p200,"y.values")),add=TRUE,col="black",lty=3,lwd=3)#
legend("bottomright",c(expression(RM~(paste(p,"*",sep="")==50)),expression(RM~(paste(p,"*",sep="")==100)),expression(RM~(paste(p,"*",sep="")==200)),expression(PCA~(paste(p,"*",sep="")==50)),expression(PCA~(paste(p,"*",sep="")==100)),expression(PCA~(paste(p,"*",sep="")==200))),lty=c(1,2,3,1,2,3),lwd=c(3,3,3,3,3,3),col=c("black","black","black","gray55","gray55","gray55"),bty="n",cex=1.7)#
dev.off()
plot(unlist(slot(rasch.perf.p200,"x.values")),unlist(slot(rasch.perf.p200,"y.values")),col="black",lty=3,lwd=3)
str(rasch.perf.p200)
plot(1-unlist(slot(rasch.perf.p200,"x.values")),unlist(slot(rasch.perf.p200,"y.values")),col="black",lty=3,lwd=3)
plot(unlist(slot(rasch.perf.p200,"x.values")),1-unlist(slot(rasch.perf.p200,"y.values")),col="black",lty=3,lwd=3)
plot(unlist(slot(rasch.perf.p200,"x.values")),unlist(slot(rasch.perf.p200,"y.values")),col="black",lty=3,lwd=3)
rasch.perf.p200 <- performance(rasch.pred.p200,"sens","spec")#
plot(unlist(slot(rasch.perf.p200,"x.values")),unlist(slot(rasch.perf.p200,"y.values")),col="black",lty=3,lwd=3)
rasch.perf.p200 <- performance(rasch.pred.p200,"spec","sens")#
plot(unlist(slot(rasch.perf.p200,"x.values")),unlist(slot(rasch.perf.p200,"y.values")),col="black",lty=3,lwd=3)
rasch.perf.p200 <- performance(rasch.pred.p200,"spec","sens")#
plot(unlist(slot(rasch.perf.p200,"x.values")),1-unlist(slot(rasch.perf.p200,"y.values")),col="black",lty=3,lwd=3)
rasch.perf.p200 <- performance(rasch.pred.p200,"spec","sens")
str(rasch.perf.p200)
plot(unlist(slot(rasch.perf.p200,"x.values")),1-unlist(slot(rasch.perf.p200,"y.values")),col="black",lty=3,lwd=3)
plot(1-unlist(slot(rasch.perf.p200,"y.values")),unlist(slot(rasch.perf.p200,"x.values")),col="black",lty=3,lwd=3)
plot(unlist(slot(rasch.perf.p200,"y.values")),unlist(slot(rasch.perf.p200,"x.values")),col="black",lty=3,lwd=3)
rasch.perf.p200 <- performance(rasch.pred.p200,"tpr","fpr")#
plot(unlist(slot(rasch.perf.p200,"y.values")),unlist(slot(rasch.perf.p200,"x.values")),col="black",lty=3,lwd=3)
str(rasch.perf.p200)
rasch.perf.p200 <- performance(rasch.pred.p200,"fpr","tpr")
pdf("roc1_bw.pdf",width=20,height=10,paper="special")#
par(cex.axis=1.7,cex.lab=1.7,pty="s",mfrow=c(1,2))#
plot(pca.perf.r50,col="gray",lty=1,lwd=3,cex=1.7,xlab="False positive rate",ylab="True positive rate",main="(a)")#
plot(pca.perf.r100,add=TRUE,col="gray",lty=2,lwd=3)#
plot(pca.perf.r200,add=TRUE,col="gray",lty=3,lwd=3)#
plot(rasch.perf.r50,add=TRUE,col="black",lty=1,lwd=3)#
plot(rasch.perf.r100,add=TRUE,col="black",lty=2,lwd=3)#
plot(rasch.perf.r200,add=TRUE,col="black",lty=3,lwd=3)#
legend("bottomright",c(expression(RM~(paste(p,"*",sep="")==50)),expression(RM~(paste(p,"*",sep="")==100)),expression(RM~(paste(p,"*",sep="")==200)),expression(PCA~(paste(p,"*",sep="")==50)),expression(PCA~(paste(p,"*",sep="")==100)),expression(PCA~(paste(p,"*",sep="")==200))),lty=c(1,2,3,1,2,3),lwd=c(3,3,3,3,3,3),col=c("black","black","black","gray55","gray55","gray55"),bty="n",cex=1.7)#
plot(pca.perf.p50,col="black",lty=1,lwd=3)#
plot(pca.perf.p100,add=TRUE,col="black",lty=2,lwd=3)#
plot(pca.perf.p200,add=TRUE,col="black",lty=3,lwd=3)#
plot(rasch.perf.p50,add=TRUE,col="black",lty=1,lwd=3)#
plot(rasch.perf.p100,add=TRUE,col="black",lty=2,lwd=3)#
plot(unlist(slot(rasch.perf.p200,"x.values")),unlist(slot(rasch.perf.p200,"y.values")),add=TRUE,col="black",lty=3,lwd=3)#
legend("bottomright",c(expression(RM~(paste(p,"*",sep="")==50)),expression(RM~(paste(p,"*",sep="")==100)),expression(RM~(paste(p,"*",sep="")==200)),expression(PCA~(paste(p,"*",sep="")==50)),expression(PCA~(paste(p,"*",sep="")==100)),expression(PCA~(paste(p,"*",sep="")==200))),lty=c(1,2,3,1,2,3),lwd=c(3,3,3,3,3,3),col=c("black","black","black","gray55","gray55","gray55"),bty="n",cex=1.7)#
dev.off()
pdf("roc1_bw.pdf",width=20,height=10,paper="special")#
par(cex.axis=1.7,cex.lab=1.7,pty="s",mfrow=c(1,2))#
plot(pca.perf.r50,col="gray",lty=1,lwd=3,cex=1.7,xlab="False positive rate",ylab="True positive rate",main="(a)")#
plot(pca.perf.r100,add=TRUE,col="gray",lty=2,lwd=3)#
plot(pca.perf.r200,add=TRUE,col="gray",lty=3,lwd=3)#
plot(rasch.perf.r50,add=TRUE,col="black",lty=1,lwd=3)#
plot(rasch.perf.r100,add=TRUE,col="black",lty=2,lwd=3)#
plot(rasch.perf.r200,add=TRUE,col="black",lty=3,lwd=3)#
legend("bottomright",c(expression(RM~(paste(p,"*",sep="")==50)),expression(RM~(paste(p,"*",sep="")==100)),expression(RM~(paste(p,"*",sep="")==200)),expression(PCA~(paste(p,"*",sep="")==50)),expression(PCA~(paste(p,"*",sep="")==100)),expression(PCA~(paste(p,"*",sep="")==200))),lty=c(1,2,3,1,2,3),lwd=c(3,3,3,3,3,3),col=c("black","black","black","gray55","gray55","gray55"),bty="n",cex=1.7)#
plot(pca.perf.p50,col="black",lty=1,lwd=3,xlab="False positive rate",ylab="True positive rate",main="(b)")#
plot(pca.perf.p100,add=TRUE,col="black",lty=2,lwd=3)#
plot(pca.perf.p200,add=TRUE,col="black",lty=3,lwd=3)#
plot(rasch.perf.p50,add=TRUE,col="black",lty=1,lwd=3)#
plot(rasch.perf.p100,add=TRUE,col="black",lty=2,lwd=3)#
plot(rasch.perf.p200,add=TRUE,col="black",lty=3,lwd=3)#
legend("bottomright",c(expression(RM~(paste(p,"*",sep="")==50)),expression(RM~(paste(p,"*",sep="")==100)),expression(RM~(paste(p,"*",sep="")==200)),expression(PCA~(paste(p,"*",sep="")==50)),expression(PCA~(paste(p,"*",sep="")==100)),expression(PCA~(paste(p,"*",sep="")==200))),lty=c(1,2,3,1,2,3),lwd=c(3,3,3,3,3,3),col=c("black","black","black","gray55","gray55","gray55"),bty="n",cex=1.7)#
dev.off()
pdf("roc1_bw.pdf",width=20,height=10,paper="special")#
par(cex.axis=1.7,cex.lab=1.7,cex.main=1.7,pty="s",mfrow=c(1,2))#
plot(pca.perf.r50,col="gray",lty=1,lwd=3,cex=1.7,xlab="False positive rate",ylab="True positive rate",main="(a)")#
plot(pca.perf.r100,add=TRUE,col="gray",lty=2,lwd=3)#
plot(pca.perf.r200,add=TRUE,col="gray",lty=3,lwd=3)#
plot(rasch.perf.r50,add=TRUE,col="black",lty=1,lwd=3)#
plot(rasch.perf.r100,add=TRUE,col="black",lty=2,lwd=3)#
plot(rasch.perf.r200,add=TRUE,col="black",lty=3,lwd=3)#
legend("bottomright",c(expression(RM~(paste(p,"*",sep="")==50)),expression(RM~(paste(p,"*",sep="")==100)),expression(RM~(paste(p,"*",sep="")==200)),expression(PCA~(paste(p,"*",sep="")==50)),expression(PCA~(paste(p,"*",sep="")==100)),expression(PCA~(paste(p,"*",sep="")==200))),lty=c(1,2,3,1,2,3),lwd=c(3,3,3,3,3,3),col=c("black","black","black","gray55","gray55","gray55"),bty="n",cex=1.7)#
plot(pca.perf.p50,col="black",lty=1,lwd=3,xlab="False positive rate",ylab="True positive rate",main="(b)")#
plot(pca.perf.p100,add=TRUE,col="black",lty=2,lwd=3)#
plot(pca.perf.p200,add=TRUE,col="black",lty=3,lwd=3)#
plot(rasch.perf.p50,add=TRUE,col="black",lty=1,lwd=3)#
plot(rasch.perf.p100,add=TRUE,col="black",lty=2,lwd=3)#
plot(rasch.perf.p200,add=TRUE,col="black",lty=3,lwd=3)#
legend("bottomright",c(expression(RM~(paste(p,"*",sep="")==50)),expression(RM~(paste(p,"*",sep="")==100)),expression(RM~(paste(p,"*",sep="")==200)),expression(PCA~(paste(p,"*",sep="")==50)),expression(PCA~(paste(p,"*",sep="")==100)),expression(PCA~(paste(p,"*",sep="")==200))),lty=c(1,2,3,1,2,3),lwd=c(3,3,3,3,3,3),col=c("black","black","black","gray55","gray55","gray55"),bty="n",cex=1.7)#
dev.off()
pdf("roc1_bw.pdf",width=20,height=10,paper="special")#
par(cex.axis=1.7,cex.lab=1.7,cex.main=1.7,pty="s",mfrow=c(1,2))#
plot(pca.perf.r50,col="gray",lty=1,lwd=3,cex=1.7,main="(a)",xlab="False positive rate",ylab="True positive rate")#
plot(pca.perf.r100,add=TRUE,col="gray",lty=2,lwd=3)#
plot(pca.perf.r200,add=TRUE,col="gray",lty=3,lwd=3)#
plot(rasch.perf.r50,add=TRUE,col="black",lty=1,lwd=3)#
plot(rasch.perf.r100,add=TRUE,col="black",lty=2,lwd=3)#
plot(rasch.perf.r200,add=TRUE,col="black",lty=3,lwd=3)#
legend("bottomright",c(expression(RM~(paste(p,"*",sep="")==50)),expression(RM~(paste(p,"*",sep="")==100)),expression(RM~(paste(p,"*",sep="")==200)),expression(PCA~(paste(p,"*",sep="")==50)),expression(PCA~(paste(p,"*",sep="")==100)),expression(PCA~(paste(p,"*",sep="")==200))),lty=c(1,2,3,1,2,3),lwd=c(3,3,3,3,3,3),col=c("black","black","black","gray55","gray55","gray55"),bty="n",cex=1.7)#
plot(pca.perf.p50,col="black",lty=1,lwd=3,xlab="False positive rate",ylab="True positive rate",main="(b)")#
plot(pca.perf.p100,add=TRUE,col="black",lty=2,lwd=3)#
plot(pca.perf.p200,add=TRUE,col="black",lty=3,lwd=3)#
plot(rasch.perf.p50,add=TRUE,col="black",lty=1,lwd=3)#
plot(rasch.perf.p100,add=TRUE,col="black",lty=2,lwd=3)#
plot(rasch.perf.p200,add=TRUE,col="black",lty=3,lwd=3)#
legend("bottomright",c(expression(RM~(paste(p,"*",sep="")==50)),expression(RM~(paste(p,"*",sep="")==100)),expression(RM~(paste(p,"*",sep="")==200)),expression(PCA~(paste(p,"*",sep="")==50)),expression(PCA~(paste(p,"*",sep="")==100)),expression(PCA~(paste(p,"*",sep="")==200))),lty=c(1,2,3,1,2,3),lwd=c(3,3,3,3,3,3),col=c("black","black","black","gray55","gray55","gray55"),bty="n",cex=1.7)#
dev.off()#
#
pdf("roc1_cl.pdf",width=20,height=10,paper="special")#
par(cex.axis=1.7,cex.lab=1.7,cex.main=1.7,pty="s",mfrow=c(1,2))#
plot(pca.perf.r50,col="blue",lty=1,lwd=3,cex=1.7,main="(a)",xlab="False positive rate",ylab="True positive rate")#
plot(pca.perf.r100,add=TRUE,col="blue",lty=2,lwd=3)#
plot(pca.perf.r200,add=TRUE,col="blue",lty=3,lwd=3)#
plot(rasch.perf.r50,add=TRUE,col="red",lty=1,lwd=3)#
plot(rasch.perf.r100,add=TRUE,col="red",lty=2,lwd=3)#
plot(rasch.perf.r200,add=TRUE,col="red",lty=3,lwd=3)#
legend("bottomright",c(expression(RM~(paste(p,"*",sep="")==50)),expression(RM~(paste(p,"*",sep="")==100)),expression(RM~(paste(p,"*",sep="")==200)),expression(PCA~(paste(p,"*",sep="")==50)),expression(PCA~(paste(p,"*",sep="")==100)),expression(PCA~(paste(p,"*",sep="")==200))),lty=c(1,2,3,1,2,3),lwd=c(3,3,3,3,3,3),col=c("red","red","red","blue","blue","blue"),bty="n",cex=1.7)#
plot(pca.perf.p50,col="blue",lty=1,lwd=3,,main="(a)",xlab="False positive rate",ylab="True positive rate")#
plot(pca.perf.p100,add=TRUE,col="blue",lty=2,lwd=3)#
plot(pca.perf.p200,add=TRUE,col="blue",lty=3,lwd=3)#
plot(rasch.perf.p50,add=TRUE,col="red",lty=1,lwd=3)#
plot(rasch.perf.p100,add=TRUE,col="red",lty=2,lwd=3)#
plot(rasch.perf.p200,add=TRUE,col="red",lty=3,lwd=3)#
legend("bottomright",c(expression(RM~(paste(p,"*",sep="")==50)),expression(RM~(paste(p,"*",sep="")==100)),expression(RM~(paste(p,"*",sep="")==200)),expression(PCA~(paste(p,"*",sep="")==50)),expression(PCA~(paste(p,"*",sep="")==100)),expression(PCA~(paste(p,"*",sep="")==200))),lty=c(1,2,3,1,2,3),lwd=c(3,3,3,3,3,3),col=c("red","red","red","blue","blue","blue"),bty="n",cex=1.7)#
dev.off()
pdf("roc1_bw_temp.pdf",width=20,height=10,paper="special")#
par(cex.axis=1.7,cex.lab=1.7,cex.main=1.7,pty="s",mfrow=c(1,2))#
plot(pca.perf.r50,col="gray",lty=1,lwd=3,cex=1.7,main="(a)",xlab="False positive rate",ylab="True positive rate")#
plot(pca.perf.r100,add=TRUE,col="gray",lty=2,lwd=3)#
plot(pca.perf.r200,add=TRUE,col="gray",lty=3,lwd=3)#
plot(rasch.perf.r50,add=TRUE,col="black",lty=1,lwd=3)#
plot(rasch.perf.r100,add=TRUE,col="black",lty=2,lwd=3)#
plot(rasch.perf.r200,add=TRUE,col="black",lty=3,lwd=3)#
legend("bottomright",c(expression(RM~(paste(p,"*",sep="")==50)),expression(RM~(paste(p,"*",sep="")==100)),expression(RM~(paste(p,"*",sep="")==200)),expression(PCA~(paste(p,"*",sep="")==50)),expression(PCA~(paste(p,"*",sep="")==100)),expression(PCA~(paste(p,"*",sep="")==200))),lty=c(1,2,3,1,2,3),lwd=c(3,3,3,3,3,3),col=c("black","black","black","gray55","gray55","gray55"),bty="n",cex=1.7)#
plot(pca.perf.p50,col="black",lty=1,lwd=3,xlab="False positive rate",ylab="True positive rate",main="(b)")#
plot(pca.perf.p100,add=TRUE,col="black",lty=2,lwd=3)#
plot(pca.perf.p200,add=TRUE,col="black",lty=3,lwd=3)#
plot(rasch.perf.p50,add=TRUE,col="black",lty=1,lwd=3)#
plot(rasch.perf.p100,add=TRUE,col="black",lty=2,lwd=3)#
plot(rasch.perf.p200,add=TRUE,col="black",lty=3,lwd=3)#
legend("bottomright",c(expression(RM~(paste(p,"*",sep="")==50)),expression(RM~(paste(p,"*",sep="")==100)),expression(RM~(paste(p,"*",sep="")==200)),expression(PCA~(paste(p,"*",sep="")==50)),expression(PCA~(paste(p,"*",sep="")==100)),expression(PCA~(paste(p,"*",sep="")==200))),lty=c(1,2,3,1,2,3),lwd=c(3,3,3,3,3,3),col=c("black","black","black","gray55","gray55","gray55"),bty="n",cex=1.7)#
dev.off()#
#
pdf("roc1_cl_temp.pdf",width=20,height=10,paper="special")#
par(cex.axis=1.7,cex.lab=1.7,cex.main=1.7,pty="s",mfrow=c(1,2))#
plot(pca.perf.r50,col="blue",lty=1,lwd=3,cex=1.7,main="(a)",xlab="False positive rate",ylab="True positive rate")#
plot(pca.perf.r100,add=TRUE,col="blue",lty=2,lwd=3)#
plot(pca.perf.r200,add=TRUE,col="blue",lty=3,lwd=3)#
plot(rasch.perf.r50,add=TRUE,col="red",lty=1,lwd=3)#
plot(rasch.perf.r100,add=TRUE,col="red",lty=2,lwd=3)#
plot(rasch.perf.r200,add=TRUE,col="red",lty=3,lwd=3)#
legend("bottomright",c(expression(RM~(paste(p,"*",sep="")==50)),expression(RM~(paste(p,"*",sep="")==100)),expression(RM~(paste(p,"*",sep="")==200)),expression(PCA~(paste(p,"*",sep="")==50)),expression(PCA~(paste(p,"*",sep="")==100)),expression(PCA~(paste(p,"*",sep="")==200))),lty=c(1,2,3,1,2,3),lwd=c(3,3,3,3,3,3),col=c("red","red","red","blue","blue","blue"),bty="n",cex=1.7)#
plot(pca.perf.p50,col="blue",lty=1,lwd=3,,main="(a)",xlab="False positive rate",ylab="True positive rate")#
plot(pca.perf.p100,add=TRUE,col="blue",lty=2,lwd=3)#
plot(pca.perf.p200,add=TRUE,col="blue",lty=3,lwd=3)#
plot(rasch.perf.p50,add=TRUE,col="red",lty=1,lwd=3)#
plot(rasch.perf.p100,add=TRUE,col="red",lty=2,lwd=3)#
plot(rasch.perf.p200,add=TRUE,col="red",lty=3,lwd=3)#
legend("bottomright",c(expression(RM~(paste(p,"*",sep="")==50)),expression(RM~(paste(p,"*",sep="")==100)),expression(RM~(paste(p,"*",sep="")==200)),expression(PCA~(paste(p,"*",sep="")==50)),expression(PCA~(paste(p,"*",sep="")==100)),expression(PCA~(paste(p,"*",sep="")==200))),lty=c(1,2,3,1,2,3),lwd=c(3,3,3,3,3,3),col=c("red","red","red","blue","blue","blue"),bty="n",cex=1.7)#
dev.off()
pdf("roc1_bw_temp.pdf",width=20,height=10,paper="special")#
par(cex.axis=1.7,cex.lab=1.7,cex.main=1.7,pty="s",mfrow=c(1,2))#
plot(pca.perf.r50,col="gray",lty=1,lwd=3,cex=1.7,main="(a)",xlab="False positive rate",ylab="True positive rate")#
plot(pca.perf.r100,add=TRUE,col="gray",lty=2,lwd=3)#
plot(pca.perf.r200,add=TRUE,col="gray",lty=3,lwd=3)#
plot(rasch.perf.r50,add=TRUE,col="black",lty=1,lwd=3)#
plot(rasch.perf.r100,add=TRUE,col="black",lty=2,lwd=3)#
plot(rasch.perf.r200,add=TRUE,col="black",lty=3,lwd=3)#
legend("bottomright",c(expression(RM~(paste(p,"*",sep="")==50)),expression(RM~(paste(p,"*",sep="")==100)),expression(RM~(paste(p,"*",sep="")==200)),expression(PCA~(paste(p,"*",sep="")==50)),expression(PCA~(paste(p,"*",sep="")==100)),expression(PCA~(paste(p,"*",sep="")==200))),lty=c(1,2,3,1,2,3),lwd=c(3,3,3,3,3,3),col=c("black","black","black","gray55","gray55","gray55"),bty="n",cex=1.7)#
plot(pca.perf.p50,col="black",lty=1,lwd=3,main="(b)",xlab="False positive rate",ylab="True positive rate")#
plot(pca.perf.p100,add=TRUE,col="black",lty=2,lwd=3)#
plot(pca.perf.p200,add=TRUE,col="black",lty=3,lwd=3)#
plot(rasch.perf.p50,add=TRUE,col="black",lty=1,lwd=3)#
plot(rasch.perf.p100,add=TRUE,col="black",lty=2,lwd=3)#
plot(rasch.perf.p200,add=TRUE,col="black",lty=3,lwd=3)#
legend("bottomright",c(expression(RM~(paste(p,"*",sep="")==50)),expression(RM~(paste(p,"*",sep="")==100)),expression(RM~(paste(p,"*",sep="")==200)),expression(PCA~(paste(p,"*",sep="")==50)),expression(PCA~(paste(p,"*",sep="")==100)),expression(PCA~(paste(p,"*",sep="")==200))),lty=c(1,2,3,1,2,3),lwd=c(3,3,3,3,3,3),col=c("black","black","black","gray55","gray55","gray55"),bty="n",cex=1.7)#
dev.off()#
#
pdf("roc1_cl_temp.pdf",width=20,height=10,paper="special")#
par(cex.axis=1.7,cex.lab=1.7,cex.main=1.7,pty="s",mfrow=c(1,2))#
plot(pca.perf.r50,col="blue",lty=1,lwd=3,cex=1.7,main="(a)",xlab="False positive rate",ylab="True positive rate")#
plot(pca.perf.r100,add=TRUE,col="blue",lty=2,lwd=3)#
plot(pca.perf.r200,add=TRUE,col="blue",lty=3,lwd=3)#
plot(rasch.perf.r50,add=TRUE,col="red",lty=1,lwd=3)#
plot(rasch.perf.r100,add=TRUE,col="red",lty=2,lwd=3)#
plot(rasch.perf.r200,add=TRUE,col="red",lty=3,lwd=3)#
legend("bottomright",c(expression(RM~(paste(p,"*",sep="")==50)),expression(RM~(paste(p,"*",sep="")==100)),expression(RM~(paste(p,"*",sep="")==200)),expression(PCA~(paste(p,"*",sep="")==50)),expression(PCA~(paste(p,"*",sep="")==100)),expression(PCA~(paste(p,"*",sep="")==200))),lty=c(1,2,3,1,2,3),lwd=c(3,3,3,3,3,3),col=c("red","red","red","blue","blue","blue"),bty="n",cex=1.7)#
plot(pca.perf.p50,col="blue",lty=1,lwd=3,main="(b)",xlab="False positive rate",ylab="True positive rate")#
plot(pca.perf.p100,add=TRUE,col="blue",lty=2,lwd=3)#
plot(pca.perf.p200,add=TRUE,col="blue",lty=3,lwd=3)#
plot(rasch.perf.p50,add=TRUE,col="red",lty=1,lwd=3)#
plot(rasch.perf.p100,add=TRUE,col="red",lty=2,lwd=3)#
plot(rasch.perf.p200,add=TRUE,col="red",lty=3,lwd=3)#
legend("bottomright",c(expression(RM~(paste(p,"*",sep="")==50)),expression(RM~(paste(p,"*",sep="")==100)),expression(RM~(paste(p,"*",sep="")==200)),expression(PCA~(paste(p,"*",sep="")==50)),expression(PCA~(paste(p,"*",sep="")==100)),expression(PCA~(paste(p,"*",sep="")==200))),lty=c(1,2,3,1,2,3),lwd=c(3,3,3,3,3,3),col=c("red","red","red","blue","blue","blue"),bty="n",cex=1.7)#
dev.off()
?pdf
library(ltm)#
data(LSAT)#
fitRasch <- rasch(LSAT)
pdf("icc_cl_temp.pdf",width=10,height=8,paper="special")#
par(cex.axis=1.7,cex.lab=1.7,mar=c(5,5,4,2))#
plot(fitRasch,items=c(1,2,3),annot=F,legend=T,cx="bottomright",main=NULL,xlab="Latent factor",ylab="Probability",lwd=3,col=c("black","red","blue"),cex=1.7)#
dev.off()
pdf("icc_cl_temp.pdf",width=10,height=8,paper="special")#
par(cex.axis=1.7,cex.lab=1.7)#
plot(fitRasch,items=c(1,2,3),annot=F,legend=T,cx="bottomright",main=NULL,xlab="Latent factor",ylab="Probability",lwd=3,col=c("black","red","blue"),cex=1.7)#
dev.off()
pdf("icc_bw_temp.pdf",width=10,height=8,paper="special")#
par(cex.axis=1.7,cex.lab=1.7,mar=c(5,5,4,2))#
plot(fitRasch,items=c(1,2,3),annot=F,legend=T,cx="bottomright",main=NULL,xlab="Latent factor",ylab="Probability",lwd=3,lty=c(1,2,4),col="black",cex=1.7)#
dev.off()#
#
pdf("icc_cl_temp.pdf",width=10,height=8,paper="special")#
par(cex.axis=1.7,cex.lab=1.7,mar=c(5,5,4,2))#
plot(fitRasch,items=c(1,2,3),annot=F,legend=T,cx="bottomright",main=NULL,xlab="Latent factor",ylab="Probability",lwd=3,col=c("black","red","blue"),cex=1.7)#
dev.off()
library(MASS)#
library(ltm)#
library(biclust)#
source("./r/pca.lda.r")#
source("./r/rasch.lda.r")#
source("./r/my.binarize.r")#
#
load("./data/prostate.rda")
getwd()
r50 <- as.matrix(read.table("./subsets/r50.txt",h=F))#
pca.components.r50 <- read.table("./components/prostate_pca_components_r50.txt",h=F)$V1#
rasch.components.r50 <- read.table("./components/prostate_rasch_components_r50.txt",h=F)$V1#
r100 <- as.matrix(read.table("./subsets/r100.txt",h=F))#
pca.components.r100 <- read.table("./components/prostate_pca_components_r100.txt",h=F)$V1#
rasch.components.r100 <- read.table("./components/prostate_rasch_components_r100.txt",h=F)$V1#
r200 <- as.matrix(read.table("./subsets/r200.txt",h=F))#
pca.components.r200 <- read.table("./components/prostate_pca_components_r200.txt",h=F)$V1#
rasch.components.r200 <- read.table("./components/prostate_rasch_components_r200.txt",h=F)$V1
pca.posterior.r50 <- NULL#
pca.class.r50 <- NULL#
for(i in 1:nrow(my.train)) {#
    predict <- pca.lda(Xtrain=prostate.x[my.train[i,],r50[i,]],Ytrain=prostate.y[my.train[i,]],Xtest=prostate.x[-(my.train)[i,],r50[i,]],ncomp=pca.components.r50[i])$posterior[,1]#
    class <- prostate.y[-(my.train)[i,]]#
    pca.posterior.r50 <- c(pca.posterior.r50,predict)#
    pca.class.r50 <- c(pca.class.r50,class)#
}#
pca.pred.r50 <- prediction(pca.posterior.r50,pca.class.r50)#
pca.perf.r50 <- performance(pca.pred.r50,"fpr","tpr")#
#
pca.posterior.r100 <- NULL#
pca.class.r100 <- NULL#
for(i in 1:nrow(my.train)) {#
    predict <- pca.lda(Xtrain=prostate.x[my.train[i,],r100[i,]],Ytrain=prostate.y[my.train[i,]],Xtest=prostate.x[-(my.train)[i,],r100[i,]],ncomp=pca.components.r100[i])$posterior[,1]#
    class <- prostate.y[-(my.train)[i,]]#
    pca.posterior.r100 <- c(pca.posterior.r100,predict)#
    pca.class.r100 <- c(pca.class.r100,class)#
}#
pca.pred.r100 <- prediction(pca.posterior.r100,pca.class.r100)#
pca.perf.r100 <- performance(pca.pred.r100,"fpr","tpr")#
#
pca.posterior.r200 <- NULL#
pca.class.r200 <- NULL#
for(i in 1:nrow(my.train)) {#
    predict <- pca.lda(Xtrain=prostate.x[my.train[i,],r200[i,]],Ytrain=prostate.y[my.train[i,]],Xtest=prostate.x[-(my.train)[i,],r200[i,]],ncomp=pca.components.r200[i])$posterior[,1]#
    class <- prostate.y[-(my.train)[i,]]#
    pca.posterior.r200 <- c(pca.posterior.r200,predict)#
    pca.class.r200 <- c(pca.class.r200,class)   #
}#
pca.pred.r200 <- prediction(pca.posterior.r200,pca.class.r200)#
pca.perf.r200 <- performance(pca.pred.r200,"fpr","tpr")
rasch.posterior.r50 <- NULL#
rasch.class.r50 <- NULL#
for(i in 1:nrow(my.train)) {#
    predict <- rasch.lda.predict(Xtrain=prostate.x[my.train[i,],r50[i,]],Ytrain=prostate.y[my.train[i,]],Xtest=prostate.x[-(my.train)[i,],r50[i,]],ncomp=rasch.components.r50[i])$posterior[,1]#
    class <- prostate.y[-(my.train)[i,]]#
    rasch.posterior.r50 <- c(rasch.posterior.r50,predict)#
    rasch.class.r50 <- c(rasch.class.r50,class)#
    cat("Complete: ",i,"\n")#
}#
rasch.pred.r50 <- prediction(rasch.posterior.r50,rasch.class.r50)#
rasch.perf.r50 <- performance(rasch.pred.r50,"fpr","tpr")
rasch.posterior.r100 <- NULL#
rasch.class.r100 <- NULL#
for(i in 1:nrow(my.train)) {#
    predict <- rasch.lda.predict(Xtrain=prostate.x[my.train[i,],r100[i,]],Ytrain=prostate.y[my.train[i,]],Xtest=prostate.x[-(my.train)[i,],r100[i,]],ncomp=rasch.components.r100[i])$posterior[,1]#
    class <- prostate.y[-(my.train)[i,]]#
    rasch.posterior.r100 <- c(rasch.posterior.r100,predict)#
    rasch.class.r100 <- c(rasch.class.r100,class)#
    cat("Complete: ",i,"\n")#
}
?rm
rm(list = ls())
library(MASS)#
library(ltm)#
library(biclust)#
source("./r/pca.lda.r")#
source("./r/rasch.lda.r")#
source("./r/my.binarize.r")#
#
load("./data/prostate.rda")#
#
my.train <- as.matrix(read.table("./subsets/my.train.txt",h=F))#
#
r50 <- as.matrix(read.table("./subsets/r50.txt",h=F))#
pca.components.r50 <- read.table("./components/prostate_pca_components_r50.txt",h=F)$V1#
rasch.components.r50 <- read.table("./components/prostate_rasch_components_r50.txt",h=F)$V1#
r100 <- as.matrix(read.table("./subsets/r100.txt",h=F))#
pca.components.r100 <- read.table("./components/prostate_pca_components_r100.txt",h=F)$V1#
rasch.components.r100 <- read.table("./components/prostate_rasch_components_r100.txt",h=F)$V1#
r200 <- as.matrix(read.table("./subsets/r200.txt",h=F))#
pca.components.r200 <- read.table("./components/prostate_pca_components_r200.txt",h=F)$V1#
rasch.components.r200 <- read.table("./components/prostate_rasch_components_r200.txt",h=F)$V1
pca.posterior.r50 <- NULL#
pca.class.r50 <- NULL#
for(i in 1:nrow(my.train)) {#
    predict <- pca.lda(Xtrain=prostate.x[my.train[i,],r50[i,]],Ytrain=prostate.y[my.train[i,]],Xtest=prostate.x[-(my.train)[i,],r50[i,]],ncomp=pca.components.r50[i])$posterior[,1]#
    class <- prostate.y[-(my.train)[i,]]#
    pca.posterior.r50 <- c(pca.posterior.r50,predict)#
    pca.class.r50 <- c(pca.class.r50,class)#
}#
pca.pred.r50 <- prediction(pca.posterior.r50,pca.class.r50)#
pca.perf.r50 <- performance(pca.pred.r50,"fpr","tpr")#
#
pca.posterior.r100 <- NULL#
pca.class.r100 <- NULL#
for(i in 1:nrow(my.train)) {#
    predict <- pca.lda(Xtrain=prostate.x[my.train[i,],r100[i,]],Ytrain=prostate.y[my.train[i,]],Xtest=prostate.x[-(my.train)[i,],r100[i,]],ncomp=pca.components.r100[i])$posterior[,1]#
    class <- prostate.y[-(my.train)[i,]]#
    pca.posterior.r100 <- c(pca.posterior.r100,predict)#
    pca.class.r100 <- c(pca.class.r100,class)#
}#
pca.pred.r100 <- prediction(pca.posterior.r100,pca.class.r100)#
pca.perf.r100 <- performance(pca.pred.r100,"fpr","tpr")#
#
pca.posterior.r200 <- NULL#
pca.class.r200 <- NULL#
for(i in 1:nrow(my.train)) {#
    predict <- pca.lda(Xtrain=prostate.x[my.train[i,],r200[i,]],Ytrain=prostate.y[my.train[i,]],Xtest=prostate.x[-(my.train)[i,],r200[i,]],ncomp=pca.components.r200[i])$posterior[,1]#
    class <- prostate.y[-(my.train)[i,]]#
    pca.posterior.r200 <- c(pca.posterior.r200,predict)#
    pca.class.r200 <- c(pca.class.r200,class)   #
}#
pca.pred.r200 <- prediction(pca.posterior.r200,pca.class.r200)#
pca.perf.r200 <- performance(pca.pred.r200,"fpr","tpr")
rasch.posterior.r50 <- NULL#
rasch.class.r50 <- NULL#
for(i in 1:nrow(my.train)) {#
    predict <- rasch.lda.predict(Xtrain=prostate.x[my.train[i,],r50[i,]],Ytrain=prostate.y[my.train[i,]],Xtest=prostate.x[-(my.train)[i,],r50[i,]],ncomp=rasch.components.r50[i])$posterior[,1]#
    class <- prostate.y[-(my.train)[i,]]#
    rasch.posterior.r50 <- c(rasch.posterior.r50,predict)#
    rasch.class.r50 <- c(rasch.class.r50,class)#
    cat("Complete: ",i,"\n")#
}#
rasch.pred.r50 <- prediction(rasch.posterior.r50,rasch.class.r50)#
rasch.perf.r50 <- performance(rasch.pred.r50,"fpr","tpr")
rasch.posterior.r100 <- NULL#
rasch.class.r100 <- NULL#
for(i in 1:nrow(my.train)) {#
    predict <- rasch.lda.predict(Xtrain=prostate.x[my.train[i,],r100[i,]],Ytrain=prostate.y[my.train[i,]],Xtest=prostate.x[-(my.train)[i,],r100[i,]],ncomp=rasch.components.r100[i])$posterior[,1]#
    class <- prostate.y[-(my.train)[i,]]#
    rasch.posterior.r100 <- c(rasch.posterior.r100,predict)#
    rasch.class.r100 <- c(rasch.class.r100,class)#
    cat("Complete: ",i,"\n")#
}#
rasch.pred.r100 <- prediction(rasch.posterior.r100,rasch.class.r100)#
rasch.perf.r100 <- performance(rasch.pred.r100,"fpr","tpr")
rasch.posterior.r200 <- NULL#
rasch.class.r200 <- NULL#
for(i in 1:nrow(my.train)) {#
    predict <- rasch.lda.predict(Xtrain=prostate.x[my.train[i,],r200[i,]],Ytrain=prostate.y[my.train[i,]],Xtest=prostate.x[-(my.train)[i,],r200[i,]],ncomp=rasch.components.r200[i])$posterior[,1]#
    class <- prostate.y[-(my.train)[i,]]#
    rasch.posterior.r200 <- c(rasch.posterior.r200,predict)#
    rasch.class.r200 <- c(rasch.class.r200,class)#
    cat("Complete: ",i,"\n")#
}#
rasch.pred.r200 <- prediction(rasch.posterior.r200,rasch.class.r200)#
rasch.perf.r200 <- performance(rasch.pred.r200,"fpr","tpr")
p50 <- as.matrix(read.table("./subsets/p50.txt",h=F))#
pca.components.p50 <- read.table("./components/prostate_pca_components_p50.txt",h=F)$V1#
rasch.components.p50 <- read.table("./components/prostate_rasch_components_p50.txt",h=F)$V1#
p100 <- as.matrix(read.table("./subsets/p100.txt",h=F))#
pca.components.p100 <- read.table("./components/prostate_pca_components_p100.txt",h=F)$V1#
rasch.components.p100 <- read.table("./components/prostate_rasch_components_p100.txt",h=F)$V1#
p200 <- as.matrix(read.table("./subsets/p200.txt",h=F))#
pca.components.p200 <- read.table("./components/prostate_pca_components_p200.txt",h=F)$V1#
rasch.components.p200 <- read.table("./components/prostate_rasch_components_p200.txt",h=F)$V1
pca.posterior.p50 <- NULL#
pca.class.p50 <- NULL#
for(i in 1:nrow(my.train)) {#
    predict <- pca.lda(Xtrain=prostate.x[my.train[i,],p50[i,]],Ytrain=prostate.y[my.train[i,]],Xtest=prostate.x[-(my.train)[i,],p50[i,]],ncomp=pca.components.p50[i])$posterior[,1]#
    class <- prostate.y[-(my.train)[i,]]#
    pca.posterior.p50 <- c(pca.posterior.p50,predict)#
    pca.class.p50 <- c(pca.class.p50,class)#
}#
pca.pred.p50 <- prediction(pca.posterior.p50,pca.class.p50)#
pca.perf.p50 <- performance(pca.pred.p50,"fpr","tpr")#
#
pca.posterior.p100 <- NULL#
pca.class.p100 <- NULL#
for(i in 1:nrow(my.train)) {#
    predict <- pca.lda(Xtrain=prostate.x[my.train[i,],p100[i,]],Ytrain=prostate.y[my.train[i,]],Xtest=prostate.x[-(my.train)[i,],p100[i,]],ncomp=pca.components.p100[i])$posterior[,1]#
    class <- prostate.y[-(my.train)[i,]]#
    pca.posterior.p100 <- c(pca.posterior.p100,predict)#
    pca.class.p100 <- c(pca.class.p100,class)#
}#
pca.pred.p100 <- prediction(pca.posterior.p100,pca.class.p100)#
pca.perf.p100 <- performance(pca.pred.p100,"fpr","tpr")#
#
pca.posterior.p200 <- NULL#
pca.class.p200 <- NULL#
for(i in 1:nrow(my.train)) {#
    predict <- pca.lda(Xtrain=prostate.x[my.train[i,],p200[i,]],Ytrain=prostate.y[my.train[i,]],Xtest=prostate.x[-(my.train)[i,],p200[i,]],ncomp=pca.components.p200[i])$posterior[,1]#
    class <- prostate.y[-(my.train)[i,]]#
    pca.posterior.p200 <- c(pca.posterior.p200,predict)#
    pca.class.p200 <- c(pca.class.p200,class)#
}#
pca.pred.p200 <- prediction(pca.posterior.p200,pca.class.p200)#
pca.perf.p200 <- performance(pca.pred.p200,"fpr","tpr")
rasch.posterior.p50 <- NULL#
rasch.class.p50 <- NULL#
for(i in 1:nrow(my.train)) {#
    predict <- rasch.lda.predict(Xtrain=prostate.x[my.train[i,],p50[i,]],Ytrain=prostate.y[my.train[i,]],Xtest=prostate.x[-(my.train)[i,],p50[i,]],ncomp=rasch.components.p50[i])$posterior[,1]#
    class <- prostate.y[-(my.train)[i,]]#
    rasch.posterior.p50 <- c(rasch.posterior.p50,predict)#
    rasch.class.p50 <- c(rasch.class.p50,class)#
    cat("Complete: ",i,"\n")#
}#
rasch.pred.p50 <- prediction(rasch.posterior.p50,rasch.class.p50)#
rasch.perf.p50 <- performance(rasch.pred.p50,"fpr","tpr")
rasch.posterior.p100 <- NULL#
rasch.class.p100 <- NULL#
for(i in 1:nrow(my.train)) {#
    predict <- rasch.lda.predict(Xtrain=prostate.x[my.train[i,],p100[i,]],Ytrain=prostate.y[my.train[i,]],Xtest=prostate.x[-(my.train)[i,],p100[i,]],ncomp=rasch.components.p100[i])$posterior[,1]#
    class <- prostate.y[-(my.train)[i,]]#
    rasch.posterior.p100 <- c(rasch.posterior.p100,predict)#
    rasch.class.p100 <- c(rasch.class.p100,class)#
    cat("Complete: ",i,"\n")#
}#
rasch.pred.p100 <- prediction(rasch.posterior.p100,rasch.class.p100)#
rasch.perf.p100 <- performance(rasch.pred.p100,"fpr","tpr")
rasch.posterior.p200 <- NULL#
rasch.class.p200 <- NULL#
for(i in 1:nrow(my.train)) {#
    predict <- rasch.lda.predict(Xtrain=prostate.x[my.train[i,],p200[i,]],Ytrain=prostate.y[my.train[i,]],Xtest=prostate.x[-(my.train)[i,],p200[i,]],ncomp=rasch.components.p200[i])$posterior[,1]#
    class <- prostate.y[-(my.train)[i,]]#
    rasch.posterior.p200 <- c(rasch.posterior.p200,predict)#
    rasch.class.p200 <- c(rasch.class.p200,class)#
    cat("Complete: ",i,"\n")#
}#
rasch.pred.p200 <- prediction(rasch.posterior.p200,rasch.class.p200)#
rasch.perf.p200 <- performance(rasch.pred.p200,"fpr","tpr")
rasch.class.p200
plot(rasch.perf.p200)
rasch.pred.p200 <- prediction(rasch.posterior.p200,as.factor(rasch.class.p200))#
rasch.perf.p200 <- performance(rasch.pred.p200,"fpr","tpr")
plot(rasch.perf.p200)
rasch.pred.p200 <- prediction(rasch.posterior.p200,as.factor(rasch.class.p200))#
rasch.perf.p200 <- performance(rasch.pred.p200,"tpr","fpr")
plot(rasch.perf.p200)
data(ROCR.simple)
pred <- prediction(ROCR.simple$predictions, ROCR.simple$labels)
perf <- performance(pred,"tpr","fpr")
plot(perf)
perf
ROCR.simple$predictions
ROCR.simple$labels
rasch.posterior.p200
as.vector(rasch.posterior.p200)
rasch.pred.p200 <- prediction(as.vector(rasch.posterior.p200),as.factor(rasch.class.p200))#
rasch.perf.p200 <- performance(rasch.pred.p200,"tpr","fpr")
plot(rasch.perf.p200)
plot(pca.perf.p200)
rasch.class.p200
pred <- prediction(as.vector(pca.posterior.p200),as.factor(pca.class.p200))
perf <- performance(pca.pred.p200,"tpr","fpr")
plot(peft)
plot(perf)
perf <- performance(pca.pred.p200,"fpr","tpr")
pred <- prediction(as.vector(pca.posterior.p200),as.factor(pca.class.p200))
perf <- performance(pca.pred.p200,"fpr","tpr")
plot(perf)
library(caret)
library(caret)
roc(as.vector(pca.posterior.p200),as.factor(pca.class.p200))
roc<-roc(as.vector(pca.posterior.p200),as.factor(pca.class.p200))
plot(roc[,3],roc[,2])
plot(1-roc[,3],roc[,2])
plot(roc$specificity,roc$sensitivity)
names(roc)
str(roc)
roc[,3]
str(roc)
plot(roc[,3],roc[,2])
plot(1-roc[,3],roc[,2])
?plot
plot(1-roc[,3],roc[,2],type="l")
pca.roc.r50 <- roc(as.vector(pca.posterior.r50),as.factor(pca.class.r50))
pca.roc.r100 <- roc(as.vector(pca.posterior.r100),as.factor(pca.class.r100))
pca.roc.r200 <- roc(as.vector(pca.posterior.r200),as.factor(pca.class.r200))
rasch.roc.r50 <- roc(as.vector(rasch.posterior.r50),as.factor(rasch.class.r50))
rasch.roc.r100 <- roc(as.vector(rasch.posterior.r100),as.factor(rasch.class.r100))
rasch.roc.r200 <- roc(as.vector(rasch.posterior.r200),as.factor(rasch.class.r200))
pca.roc.p50 <- roc(as.vector(pca.posterior.p50),as.factor(pca.class.p50))
pca.roc.p100 <- roc(as.vector(pca.posterior.p100),as.factor(pca.class.p100))
pca.roc.p200 <- roc(as.vector(pca.posterior.p200),as.factor(pca.class.p200))
rasch.roc.p50 <- roc(as.vector(rasch.posterior.p50),as.factor(rasch.class.p50))
rasch.roc.p100 <- roc(as.vector(rasch.posterior.p100),as.factor(rasch.class.p100))
rasch.roc.p200 <- roc(as.vector(rasch.posterior.p200),as.factor(rasch.class.p200))
str(pca.roc.r50)
plot(1-pca.roc.r50[,3],pca.roc.r50[,2],type="l",xlab="1-Specificity",ylab="Sensitivity")
curve(1-pca.roc.r100[,3],pca.roc.r100[,2],type="l")
?lines
points(1-pca.roc.r100[,3],pca.roc.r100[,2],type="l")
lines(1-pca.roc.r100[,3],pca.roc.r100[,2],type="l")
plot(1-pca.roc.r50[,3],pca.roc.r50[,2],type="l",xlab="1-Specificity",ylab="Sensitivity")
lines(1-pca.roc.r100[,3],pca.roc.r100[,2],type="l")
plot(1-pca.roc.r50[,3],pca.roc.r50[,2],type="l",xlab="1-Specificity",ylab="Sensitivity")#
curve(1-pca.roc.r100[,3],pca.roc.r100[,2],type="l")#
curve(1-pca.roc.r200[,3],pca.roc.r200[,2],type="l")
plot(1-pca.roc.r50[,3],pca.roc.r50[,2],type="l",xlab="1-Specificity",ylab="Sensitivity")#
lines(1-pca.roc.r100[,3],pca.roc.r100[,2],type="l")#
lines(1-pca.roc.r200[,3],pca.roc.r200[,2],type="l")
plot(1-pca.roc.r50[,3],pca.roc.r50[,2],type="l",xlab="1-Specificity",ylab="Sensitivity",col="blue",lty=1,lwd=3)#
lines(1-pca.roc.r100[,3],pca.roc.r100[,2],type="l",col="blue",lty=2,lwd=3)#
lines(1-pca.roc.r200[,3],pca.roc.r200[,2],type="l",col="blue",lty=3,lwd=3)
plot(1-pca.roc.r50[,3],pca.roc.r50[,2],type="l",xlab="1-Specificity",ylab="Sensitivity",col="blue",lty=1,lwd=3)#
lines(1-pca.roc.r100[,3],pca.roc.r100[,2],type="l",col="blue",lty=2,lwd=3)#
lines(1-pca.roc.r200[,3],pca.roc.r200[,2],type="l",col="blue",lty=3,lwd=3)
?lines
plot(1-pca.roc.r50[,3],pca.roc.r50[,2],type="l",xlab="1-Specificity",ylab="Sensitivity",col="blue",lty=1,lwd=3)#
lines(1-pca.roc.r100[,3],pca.roc.r100[,2],col="blue",lty=2,lwd=3)#
lines(1-pca.roc.r200[,3],pca.roc.r200[,2],col="blue",lty=3,lwd=3)
plot(1-pca.roc.r50[,3],pca.roc.r50[,2],type="l",xlab="1-Specificity",ylab="Sensitivity",col="blue",lty=1,lwd=3)#
lines(1-pca.roc.r100[,3],pca.roc.r100[,2],col="red",lty=2,lwd=3)#
lines(1-pca.roc.r200[,3],pca.roc.r200[,2],col="green",lty=3,lwd=3)
plot(1-pca.roc.r50[,3],pca.roc.r50[,2],type="l",xlab="1-Specificity",ylab="Sensitivity",col="blue",lty=1,lwd=3)#
lines(1-pca.roc.r100[,3],pca.roc.r100[,2],col="red",lty=4,lwd=3)#
lines(1-pca.roc.r200[,3],pca.roc.r200[,2],col="green",lty=3,lwd=3)
plot(1-pca.roc.r50[,3],pca.roc.r50[,2],xlab="1-Specificity",ylab="Sensitivity",col="blue",lty=1,lwd=3)#
lines(1-pca.roc.r100[,3],pca.roc.r100[,2],col="red",lty=4,lwd=3)#
lines(1-pca.roc.r200[,3],pca.roc.r200[,2],col="green",lty=3,lwd=3)
plot(1-pca.roc.r50[,3],pca.roc.r50[,2],type="l",xlab="1-Specificity",ylab="Sensitivity",col="blue",lty=1,lwd=1)#
lines(1-pca.roc.r100[,3],pca.roc.r100[,2],col="red",lty=4,lwd=1)#
lines(1-pca.roc.r200[,3],pca.roc.r200[,2],col="green",lty=3,lwd=1)
plot(1-pca.roc.r50[,3],pca.roc.r50[,2],type="l",xlab="1-Specificity",ylab="Sensitivity",col="blue",lty=1,lwd=1)#
lines(1-pca.roc.r100[,3],pca.roc.r100[,2],col="red",lty=2,lwd=1)#
lines(1-pca.roc.r200[,3],pca.roc.r200[,2],col="green",lty=3,lwd=1)
plot(1-pca.roc.r50[,3],pca.roc.r50[,2],type="l",xlab="1-Specificity",ylab="Sensitivity",col="blue",lty=1,lwd=2)#
lines(1-pca.roc.r100[,3],pca.roc.r100[,2],col="red",lty=2,lwd=2)#
lines(1-pca.roc.r200[,3],pca.roc.r200[,2],col="green",lty=3,lwd=2)
plot(1-pca.roc.r50[,3],pca.roc.r50[,2],type="l",xlab="1-Specificity",ylab="Sensitivity",col="blue",lty=1,lwd=3)#
lines(1-pca.roc.r100[,3],pca.roc.r100[,2],col="red",lty=2,lwd=3)#
lines(1-pca.roc.r200[,3],pca.roc.r200[,2],col="green",lty=3,lwd=3)
plot(1-pca.roc.r50[,3],pca.roc.r50[,2],type="l",xlab="1-Specificity",ylab="Sensitivity",col="blue",lty=1,lwd=3)#
lines(1-pca.roc.r100[,3],pca.roc.r100[,2],col="red",lty=2,lwd=3)#
lines(1-pca.roc.r200[,3],pca.roc.r200[,2],col="green",lty=3,lwd=3)
pdf("roc2_cl_temp.pdf",width=20,height=10,paper="special")#
par(cex.axis=1.7,cex.lab=1.7,cex.main=1.7,pty="s",mfrow=c(1,2))#
plot(1-pca.roc.r50[,3],pca.roc.r50[,2],type="l",cex=1.7,main="(a)",xlab="1-Specificity",ylab="Sensitivity",col="blue",lty=1,lwd=3)#
lines(1-pca.roc.r100[,3],pca.roc.r100[,2],col="blue",lty=2,lwd=3)#
lines(1-pca.roc.r200[,3],pca.roc.r200[,2],col="blue",lty=3,lwd=3)#
lines(1-rasch.roc.r50[,3],rasch.roc.r50[,2],col="red",lty=1,lwd=3)#
lines(1-rasch.roc.r100[,3],rasch.roc.r100[,2],col="red",lty=2,lwd=3)#
lines(1-rasch.roc.r200[,3],rasch.roc.r200[,2],col="red",lty=3,lwd=3)#
legend("bottomright",c(expression(RM~(paste(p,"*",sep="")==50)),expression(RM~(paste(p,"*",sep="")==100)),expression(RM~(paste(p,"*",sep="")==200)),expression(PCA~(paste(p,"*",sep="")==50)),expression(PCA~(paste(p,"*",sep="")==100)),expression(PCA~(paste(p,"*",sep="")==200))),lty=c(1,2,3,1,2,3),lwd=c(3,3,3,3,3,3),col=c("red","red","red","blue","blue","blue"),bty="n",cex=1.7)#
plot(1-pca.roc.p50[,3],pca.roc.p50[,2],type="l",cex=1.7,main="(b)",xlab="1-Specificity",ylab="Sensitivity",col="blue",lty=1,lwd=3)#
lines(1-pca.roc.p100[,3],pca.roc.p100[,2],col="blue",lty=2,lwd=3)#
lines(1-pca.roc.p200[,3],pca.roc.p200[,2],col="blue",lty=3,lwd=3)#
lines(1-rasch.roc.p50[,3],rasch.roc.p50[,2],col="red",lty=1,lwd=3)#
lines(1-rasch.roc.p100[,3],rasch.roc.p100[,2],col="red",lty=2,lwd=3)#
lines(1-rasch.roc.p200[,3],rasch.roc.p200[,2],col="red",lty=3,lwd=3)#
legend("bottomright",c(expression(RM~(paste(p,"*",sep="")==50)),expression(RM~(paste(p,"*",sep="")==100)),expression(RM~(paste(p,"*",sep="")==200)),expression(PCA~(paste(p,"*",sep="")==50)),expression(PCA~(paste(p,"*",sep="")==100)),expression(PCA~(paste(p,"*",sep="")==200))),lty=c(1,2,3,1,2,3),lwd=c(3,3,3,3,3,3),col=c("red","red","red","blue","blue","blue"),bty="n",cex=1.7)#
dev.off()
pdf("roc2_bw_temp.pdf",width=20,height=10,paper="special")#
par(cex.axis=1.7,cex.lab=1.7,cex.main=1.7,pty="s",mfrow=c(1,2))#
plot(1-pca.roc.r50[,3],pca.roc.r50[,2],type="l",cex=1.7,main="(a)",xlab="1-Specificity",ylab="Sensitivity",col="gray55",lty=1,lwd=3)#
lines(1-pca.roc.r100[,3],pca.roc.r100[,2],col="gray55",lty=2,lwd=3)#
lines(1-pca.roc.r200[,3],pca.roc.r200[,2],col="gray55",lty=3,lwd=3)#
lines(1-rasch.roc.r50[,3],rasch.roc.r50[,2],col="red",lty=1,lwd=3)#
lines(1-rasch.roc.r100[,3],rasch.roc.r100[,2],col="red",lty=2,lwd=3)#
lines(1-rasch.roc.r200[,3],rasch.roc.r200[,2],col="red",lty=3,lwd=3)#
legend("bottomright",c(expression(RM~(paste(p,"*",sep="")==50)),expression(RM~(paste(p,"*",sep="")==100)),expression(RM~(paste(p,"*",sep="")==200)),expression(PCA~(paste(p,"*",sep="")==50)),expression(PCA~(paste(p,"*",sep="")==100)),expression(PCA~(paste(p,"*",sep="")==200))),lty=c(1,2,3,1,2,3),lwd=c(3,3,3,3,3,3),col=c("black","black","black","gray55","gray55","gray55"),bty="n",cex=1.7)#
plot(1-pca.roc.p50[,3],pca.roc.p50[,2],type="l",cex=1.7,main="(b)",xlab="1-Specificity",ylab="Sensitivity",col="gray55",lty=1,lwd=3)#
lines(1-pca.roc.p100[,3],pca.roc.p100[,2],col="gray55",lty=2,lwd=3)#
lines(1-pca.roc.p200[,3],pca.roc.p200[,2],col="gray55",lty=3,lwd=3)#
lines(1-rasch.roc.p50[,3],rasch.roc.p50[,2],col="black",lty=1,lwd=3)#
lines(1-rasch.roc.p100[,3],rasch.roc.p100[,2],col="black",lty=2,lwd=3)#
lines(1-rasch.roc.p200[,3],rasch.roc.p200[,2],col="black",lty=3,lwd=3)#
legend("bottomright",c(expression(RM~(paste(p,"*",sep="")==50)),expression(RM~(paste(p,"*",sep="")==100)),expression(RM~(paste(p,"*",sep="")==200)),expression(PCA~(paste(p,"*",sep="")==50)),expression(PCA~(paste(p,"*",sep="")==100)),expression(PCA~(paste(p,"*",sep="")==200))),lty=c(1,2,3,1,2,3),lwd=c(3,3,3,3,3,3),col=c("black","black","black","gray55","gray55","gray55"),bty="n",cex=1.7)#
dev.off()#
#
pdf("roc2_cl_temp.pdf",width=20,height=10,paper="special")#
par(cex.axis=1.7,cex.lab=1.7,cex.main=1.7,pty="s",mfrow=c(1,2))#
plot(1-pca.roc.r50[,3],pca.roc.r50[,2],type="l",cex=1.7,main="(a)",xlab="1-Specificity",ylab="Sensitivity",col="blue",lty=1,lwd=3)#
lines(1-pca.roc.r100[,3],pca.roc.r100[,2],col="blue",lty=2,lwd=3)#
lines(1-pca.roc.r200[,3],pca.roc.r200[,2],col="blue",lty=3,lwd=3)#
lines(1-rasch.roc.r50[,3],rasch.roc.r50[,2],col="red",lty=1,lwd=3)#
lines(1-rasch.roc.r100[,3],rasch.roc.r100[,2],col="red",lty=2,lwd=3)#
lines(1-rasch.roc.r200[,3],rasch.roc.r200[,2],col="red",lty=3,lwd=3)#
legend("bottomright",c(expression(RM~(paste(p,"*",sep="")==50)),expression(RM~(paste(p,"*",sep="")==100)),expression(RM~(paste(p,"*",sep="")==200)),expression(PCA~(paste(p,"*",sep="")==50)),expression(PCA~(paste(p,"*",sep="")==100)),expression(PCA~(paste(p,"*",sep="")==200))),lty=c(1,2,3,1,2,3),lwd=c(3,3,3,3,3,3),col=c("red","red","red","blue","blue","blue"),bty="n",cex=1.7)#
plot(1-pca.roc.p50[,3],pca.roc.p50[,2],type="l",cex=1.7,main="(b)",xlab="1-Specificity",ylab="Sensitivity",col="blue",lty=1,lwd=3)#
lines(1-pca.roc.p100[,3],pca.roc.p100[,2],col="blue",lty=2,lwd=3)#
lines(1-pca.roc.p200[,3],pca.roc.p200[,2],col="blue",lty=3,lwd=3)#
lines(1-rasch.roc.p50[,3],rasch.roc.p50[,2],col="red",lty=1,lwd=3)#
lines(1-rasch.roc.p100[,3],rasch.roc.p100[,2],col="red",lty=2,lwd=3)#
lines(1-rasch.roc.p200[,3],rasch.roc.p200[,2],col="red",lty=3,lwd=3)#
legend("bottomright",c(expression(RM~(paste(p,"*",sep="")==50)),expression(RM~(paste(p,"*",sep="")==100)),expression(RM~(paste(p,"*",sep="")==200)),expression(PCA~(paste(p,"*",sep="")==50)),expression(PCA~(paste(p,"*",sep="")==100)),expression(PCA~(paste(p,"*",sep="")==200))),lty=c(1,2,3,1,2,3),lwd=c(3,3,3,3,3,3),col=c("red","red","red","blue","blue","blue"),bty="n",cex=1.7)#
dev.off()
pdf("roc2_bw_temp.pdf",width=20,height=10,paper="special")#
par(cex.axis=1.7,cex.lab=1.7,cex.main=1.7,pty="s",mfrow=c(1,2))#
plot(1-pca.roc.r50[,3],pca.roc.r50[,2],type="l",cex=1.7,main="(a)",xlab="1-Specificity",ylab="Sensitivity",col="gray55",lty=1,lwd=3)#
lines(1-pca.roc.r100[,3],pca.roc.r100[,2],col="gray55",lty=2,lwd=3)#
lines(1-pca.roc.r200[,3],pca.roc.r200[,2],col="gray55",lty=3,lwd=3)#
lines(1-rasch.roc.r50[,3],rasch.roc.r50[,2],col="black",lty=1,lwd=3)#
lines(1-rasch.roc.r100[,3],rasch.roc.r100[,2],col="black",lty=2,lwd=3)#
lines(1-rasch.roc.r200[,3],rasch.roc.r200[,2],col="black",lty=3,lwd=3)#
legend("bottomright",c(expression(RM~(paste(p,"*",sep="")==50)),expression(RM~(paste(p,"*",sep="")==100)),expression(RM~(paste(p,"*",sep="")==200)),expression(PCA~(paste(p,"*",sep="")==50)),expression(PCA~(paste(p,"*",sep="")==100)),expression(PCA~(paste(p,"*",sep="")==200))),lty=c(1,2,3,1,2,3),lwd=c(3,3,3,3,3,3),col=c("black","black","black","gray55","gray55","gray55"),bty="n",cex=1.7)#
plot(1-pca.roc.p50[,3],pca.roc.p50[,2],type="l",cex=1.7,main="(b)",xlab="1-Specificity",ylab="Sensitivity",col="gray55",lty=1,lwd=3)#
lines(1-pca.roc.p100[,3],pca.roc.p100[,2],col="gray55",lty=2,lwd=3)#
lines(1-pca.roc.p200[,3],pca.roc.p200[,2],col="gray55",lty=3,lwd=3)#
lines(1-rasch.roc.p50[,3],rasch.roc.p50[,2],col="black",lty=1,lwd=3)#
lines(1-rasch.roc.p100[,3],rasch.roc.p100[,2],col="black",lty=2,lwd=3)#
lines(1-rasch.roc.p200[,3],rasch.roc.p200[,2],col="black",lty=3,lwd=3)#
legend("bottomright",c(expression(RM~(paste(p,"*",sep="")==50)),expression(RM~(paste(p,"*",sep="")==100)),expression(RM~(paste(p,"*",sep="")==200)),expression(PCA~(paste(p,"*",sep="")==50)),expression(PCA~(paste(p,"*",sep="")==100)),expression(PCA~(paste(p,"*",sep="")==200))),lty=c(1,2,3,1,2,3),lwd=c(3,3,3,3,3,3),col=c("black","black","black","gray55","gray55","gray55"),bty="n",cex=1.7)#
dev.off()#
#
pdf("roc2_cl_temp.pdf",width=20,height=10,paper="special")#
par(cex.axis=1.7,cex.lab=1.7,cex.main=1.7,pty="s",mfrow=c(1,2))#
plot(1-pca.roc.r50[,3],pca.roc.r50[,2],type="l",cex=1.7,main="(a)",xlab="1-Specificity",ylab="Sensitivity",col="blue",lty=1,lwd=3)#
lines(1-pca.roc.r100[,3],pca.roc.r100[,2],col="blue",lty=2,lwd=3)#
lines(1-pca.roc.r200[,3],pca.roc.r200[,2],col="blue",lty=3,lwd=3)#
lines(1-rasch.roc.r50[,3],rasch.roc.r50[,2],col="red",lty=1,lwd=3)#
lines(1-rasch.roc.r100[,3],rasch.roc.r100[,2],col="red",lty=2,lwd=3)#
lines(1-rasch.roc.r200[,3],rasch.roc.r200[,2],col="red",lty=3,lwd=3)#
legend("bottomright",c(expression(RM~(paste(p,"*",sep="")==50)),expression(RM~(paste(p,"*",sep="")==100)),expression(RM~(paste(p,"*",sep="")==200)),expression(PCA~(paste(p,"*",sep="")==50)),expression(PCA~(paste(p,"*",sep="")==100)),expression(PCA~(paste(p,"*",sep="")==200))),lty=c(1,2,3,1,2,3),lwd=c(3,3,3,3,3,3),col=c("red","red","red","blue","blue","blue"),bty="n",cex=1.7)#
plot(1-pca.roc.p50[,3],pca.roc.p50[,2],type="l",cex=1.7,main="(b)",xlab="1-Specificity",ylab="Sensitivity",col="blue",lty=1,lwd=3)#
lines(1-pca.roc.p100[,3],pca.roc.p100[,2],col="blue",lty=2,lwd=3)#
lines(1-pca.roc.p200[,3],pca.roc.p200[,2],col="blue",lty=3,lwd=3)#
lines(1-rasch.roc.p50[,3],rasch.roc.p50[,2],col="red",lty=1,lwd=3)#
lines(1-rasch.roc.p100[,3],rasch.roc.p100[,2],col="red",lty=2,lwd=3)#
lines(1-rasch.roc.p200[,3],rasch.roc.p200[,2],col="red",lty=3,lwd=3)#
legend("bottomright",c(expression(RM~(paste(p,"*",sep="")==50)),expression(RM~(paste(p,"*",sep="")==100)),expression(RM~(paste(p,"*",sep="")==200)),expression(PCA~(paste(p,"*",sep="")==50)),expression(PCA~(paste(p,"*",sep="")==100)),expression(PCA~(paste(p,"*",sep="")==200))),lty=c(1,2,3,1,2,3),lwd=c(3,3,3,3,3,3),col=c("red","red","red","blue","blue","blue"),bty="n",cex=1.7)#
dev.off()
pdf("roc2_bw_temp.pdf",width=20,height=10,paper="special")#
par(cex.axis=1.7,cex.lab=1.7,cex.main=1.7,pty="s",mfrow=c(1,2))#
plot(1-pca.roc.r50[,3],pca.roc.r50[,2],type="l",cex=1.7,main="(a)",xlab="1-Specificity",ylab="Sensitivity",col="gray55",lty=1,lwd=3)#
lines(1-pca.roc.r100[,3],pca.roc.r100[,2],col="gray55",lty=2,lwd=3)#
lines(1-pca.roc.r200[,3],pca.roc.r200[,2],col="gray55",lty=3,lwd=3)#
lines(1-rasch.roc.r50[,3],rasch.roc.r50[,2],col="black",lty=1,lwd=3)#
lines(1-rasch.roc.r100[,3],rasch.roc.r100[,2],col="black",lty=2,lwd=3)#
lines(1-rasch.roc.r200[,3],rasch.roc.r200[,2],col="black",lty=3,lwd=3)#
legend("bottomright",c(expression(RM~(paste(p,"*",sep="")==50)),expression(RM~(paste(p,"*",sep="")==100)),expression(RM~(paste(p,"*",sep="")==200)),expression(PCA~(paste(p,"*",sep="")==50)),expression(PCA~(paste(p,"*",sep="")==100)),expression(PCA~(paste(p,"*",sep="")==200))),lty=c(1,2,3,1,2,3),lwd=c(3,3,3,3,3,3),col=c("black","black","black","gray55","gray55","gray55"),bty="n",cex=1.7)#
plot(1-pca.roc.p50[,3],pca.roc.p50[,2],type="l",cex=1.7,main="(b)",xlab="1-Specificity",ylab="Sensitivity",col="gray60",lty=1,lwd=3)#
lines(1-pca.roc.p100[,3],pca.roc.p100[,2],col="gray60",lty=2,lwd=3)#
lines(1-pca.roc.p200[,3],pca.roc.p200[,2],col="gray60",lty=3,lwd=3)#
lines(1-rasch.roc.p50[,3],rasch.roc.p50[,2],col="black",lty=1,lwd=3)#
lines(1-rasch.roc.p100[,3],rasch.roc.p100[,2],col="black",lty=2,lwd=3)#
lines(1-rasch.roc.p200[,3],rasch.roc.p200[,2],col="black",lty=3,lwd=3)#
legend("bottomright",c(expression(RM~(paste(p,"*",sep="")==50)),expression(RM~(paste(p,"*",sep="")==100)),expression(RM~(paste(p,"*",sep="")==200)),expression(PCA~(paste(p,"*",sep="")==50)),expression(PCA~(paste(p,"*",sep="")==100)),expression(PCA~(paste(p,"*",sep="")==200))),lty=c(1,2,3,1,2,3),lwd=c(3,3,3,3,3,3),col=c("black","black","black","gray55","gray55","gray55"),bty="n",cex=1.7)#
dev.off()#
#
pdf("roc2_cl_temp.pdf",width=20,height=10,paper="special")#
par(cex.axis=1.7,cex.lab=1.7,cex.main=1.7,pty="s",mfrow=c(1,2))#
plot(1-pca.roc.r50[,3],pca.roc.r50[,2],type="l",cex=1.7,main="(a)",xlab="1-Specificity",ylab="Sensitivity",col="blue",lty=1,lwd=3)#
lines(1-pca.roc.r100[,3],pca.roc.r100[,2],col="blue",lty=2,lwd=3)#
lines(1-pca.roc.r200[,3],pca.roc.r200[,2],col="blue",lty=3,lwd=3)#
lines(1-rasch.roc.r50[,3],rasch.roc.r50[,2],col="red",lty=1,lwd=3)#
lines(1-rasch.roc.r100[,3],rasch.roc.r100[,2],col="red",lty=2,lwd=3)#
lines(1-rasch.roc.r200[,3],rasch.roc.r200[,2],col="red",lty=3,lwd=3)#
legend("bottomright",c(expression(RM~(paste(p,"*",sep="")==50)),expression(RM~(paste(p,"*",sep="")==100)),expression(RM~(paste(p,"*",sep="")==200)),expression(PCA~(paste(p,"*",sep="")==50)),expression(PCA~(paste(p,"*",sep="")==100)),expression(PCA~(paste(p,"*",sep="")==200))),lty=c(1,2,3,1,2,3),lwd=c(3,3,3,3,3,3),col=c("red","red","red","blue","blue","blue"),bty="n",cex=1.7)#
plot(1-pca.roc.p50[,3],pca.roc.p50[,2],type="l",cex=1.7,main="(b)",xlab="1-Specificity",ylab="Sensitivity",col="blue",lty=1,lwd=3)#
lines(1-pca.roc.p100[,3],pca.roc.p100[,2],col="blue",lty=2,lwd=3)#
lines(1-pca.roc.p200[,3],pca.roc.p200[,2],col="blue",lty=3,lwd=3)#
lines(1-rasch.roc.p50[,3],rasch.roc.p50[,2],col="red",lty=1,lwd=3)#
lines(1-rasch.roc.p100[,3],rasch.roc.p100[,2],col="red",lty=2,lwd=3)#
lines(1-rasch.roc.p200[,3],rasch.roc.p200[,2],col="red",lty=3,lwd=3)#
legend("bottomright",c(expression(RM~(paste(p,"*",sep="")==50)),expression(RM~(paste(p,"*",sep="")==100)),expression(RM~(paste(p,"*",sep="")==200)),expression(PCA~(paste(p,"*",sep="")==50)),expression(PCA~(paste(p,"*",sep="")==100)),expression(PCA~(paste(p,"*",sep="")==200))),lty=c(1,2,3,1,2,3),lwd=c(3,3,3,3,3,3),col=c("red","red","red","blue","blue","blue"),bty="n",cex=1.7)#
dev.off()
aucRoc(rasch.roc.p200)
aucRoc(rasch.roc.p100)
aucRoc(rasch.roc.p50)
aucRoc(pca.roc.p50)
aucRoc(pca.roc.p100)
aucRoc(pca.roc.p200)
aucRoc(pca.roc.r50)
aucRoc(pca.roc.r100)
aucRoc(pca.roc.r200)
aucRoc(rasch.roc.r50)
aucRoc(rasch.roc.r100)
aucRoc(rasch.roc.r200)
aucRoc(pca.roc.p50)
aucRoc(pca.roc.p100)
aucRoc(pca.roc.p200)
aucRoc(rasch.roc.p50)
aucRoc(rasch.roc.p100)
aucRoc(rasch.roc.p200)
pdf("roc2_bw_temp.pdf",width=20,height=10,paper="special")#
par(cex.axis=1.7,cex.lab=1.7,cex.main=1.7,pty="s",mfrow=c(1,2))#
plot(1-pca.roc.r50[,3],pca.roc.r50[,2],type="l",cex=1.7,main="(a)",xlab="1-Specificity",ylab="Sensitivity",col="gray55",lty=1,lwd=3)#
lines(1-pca.roc.r100[,3],pca.roc.r100[,2],col="gray55",lty=2,lwd=3)#
lines(1-pca.roc.r200[,3],pca.roc.r200[,2],col="gray55",lty=3,lwd=3)#
lines(1-rasch.roc.r50[,3],rasch.roc.r50[,2],col="black",lty=1,lwd=3)#
lines(1-rasch.roc.r100[,3],rasch.roc.r100[,2],col="black",lty=2,lwd=3)#
lines(1-rasch.roc.r200[,3],rasch.roc.r200[,2],col="black",lty=3,lwd=3)#
legend("bottomright",c(expression(paste("RM (p*",==50,")")),expression(RM~(paste(p,"*",sep="")==100)),expression(RM~(paste(p,"*",sep="")==200)),expression(PCA~(paste(p,"*",sep="")==50)),expression(PCA~(paste(p,"*",sep="")==100)),expression(PCA~(paste(p,"*",sep="")==200))),lty=c(1,2,3,1,2,3),lwd=c(3,3,3,3,3,3),col=c("black","black","black","gray55","gray55","gray55"),bty="n",cex=1.7)#
plot(1-pca.roc.p50[,3],pca.roc.p50[,2],type="l",cex=1.7,main="(b)",xlab="1-Specificity",ylab="Sensitivity",col="gray60",lty=1,lwd=3)#
lines(1-pca.roc.p100[,3],pca.roc.p100[,2],col="gray60",lty=2,lwd=3)#
lines(1-pca.roc.p200[,3],pca.roc.p200[,2],col="gray60",lty=3,lwd=3)#
lines(1-rasch.roc.p50[,3],rasch.roc.p50[,2],col="black",lty=1,lwd=3)#
lines(1-rasch.roc.p100[,3],rasch.roc.p100[,2],col="black",lty=2,lwd=3)#
lines(1-rasch.roc.p200[,3],rasch.roc.p200[,2],col="black",lty=3,lwd=3)#
legend("bottomright",c(expression(RM~(paste(p,"*",sep="")==50)),expression(RM~(paste(p,"*",sep="")==100)),expression(RM~(paste(p,"*",sep="")==200)),expression(PCA~(paste(p,"*",sep="")==50)),expression(PCA~(paste(p,"*",sep="")==100)),expression(PCA~(paste(p,"*",sep="")==200))),lty=c(1,2,3,1,2,3),lwd=c(3,3,3,3,3,3),col=c("black","black","black","gray55","gray55","gray55"),bty="n",cex=1.7)#
dev.off()
pdf("roc2_bw_temp.pdf",width=20,height=10,paper="special")#
par(cex.axis=1.7,cex.lab=1.7,cex.main=1.7,pty="s",mfrow=c(1,2))#
plot(1-pca.roc.r50[,3],pca.roc.r50[,2],type="l",cex=1.7,main="(a)",xlab="1-Specificity",ylab="Sensitivity",col="gray55",lty=1,lwd=3)#
lines(1-pca.roc.r100[,3],pca.roc.r100[,2],col="gray55",lty=2,lwd=3)#
lines(1-pca.roc.r200[,3],pca.roc.r200[,2],col="gray55",lty=3,lwd=3)#
lines(1-rasch.roc.r50[,3],rasch.roc.r50[,2],col="black",lty=1,lwd=3)#
lines(1-rasch.roc.r100[,3],rasch.roc.r100[,2],col="black",lty=2,lwd=3)#
lines(1-rasch.roc.r200[,3],rasch.roc.r200[,2],col="black",lty=3,lwd=3)#
legend("bottomright",c(expression(paste("RM",degree)),expression(RM~(paste(p,"*",sep="")==100)),expression(RM~(paste(p,"*",sep="")==200)),expression(PCA~(paste(p,"*",sep="")==50)),expression(PCA~(paste(p,"*",sep="")==100)),expression(PCA~(paste(p,"*",sep="")==200))),lty=c(1,2,3,1,2,3),lwd=c(3,3,3,3,3,3),col=c("black","black","black","gray55","gray55","gray55"),bty="n",cex=1.7)#
plot(1-pca.roc.p50[,3],pca.roc.p50[,2],type="l",cex=1.7,main="(b)",xlab="1-Specificity",ylab="Sensitivity",col="gray60",lty=1,lwd=3)#
lines(1-pca.roc.p100[,3],pca.roc.p100[,2],col="gray60",lty=2,lwd=3)#
lines(1-pca.roc.p200[,3],pca.roc.p200[,2],col="gray60",lty=3,lwd=3)#
lines(1-rasch.roc.p50[,3],rasch.roc.p50[,2],col="black",lty=1,lwd=3)#
lines(1-rasch.roc.p100[,3],rasch.roc.p100[,2],col="black",lty=2,lwd=3)#
lines(1-rasch.roc.p200[,3],rasch.roc.p200[,2],col="black",lty=3,lwd=3)#
legend("bottomright",c(expression(RM~(paste(p,"*",sep="")==50)),expression(RM~(paste(p,"*",sep="")==100)),expression(RM~(paste(p,"*",sep="")==200)),expression(PCA~(paste(p,"*",sep="")==50)),expression(PCA~(paste(p,"*",sep="")==100)),expression(PCA~(paste(p,"*",sep="")==200))),lty=c(1,2,3,1,2,3),lwd=c(3,3,3,3,3,3),col=c("black","black","black","gray55","gray55","gray55"),bty="n",cex=1.7)#
dev.off()
pdf("roc2_bw_temp.pdf",width=20,height=10,paper="special")#
par(cex.axis=1.7,cex.lab=1.7,cex.main=1.7,pty="s",mfrow=c(1,2))#
plot(1-pca.roc.r50[,3],pca.roc.r50[,2],type="l",cex=1.7,main="(a)",xlab="1-Specificity",ylab="Sensitivity",col="gray55",lty=1,lwd=3)#
lines(1-pca.roc.r100[,3],pca.roc.r100[,2],col="gray55",lty=2,lwd=3)#
lines(1-pca.roc.r200[,3],pca.roc.r200[,2],col="gray55",lty=3,lwd=3)#
lines(1-rasch.roc.r50[,3],rasch.roc.r50[,2],col="black",lty=1,lwd=3)#
lines(1-rasch.roc.r100[,3],rasch.roc.r100[,2],col="black",lty=2,lwd=3)#
lines(1-rasch.roc.r200[,3],rasch.roc.r200[,2],col="black",lty=3,lwd=3)#
legend("bottomright",c(expression(paste("RM",["p*"])),expression(RM~(paste(p,"*",sep="")==100)),expression(RM~(paste(p,"*",sep="")==200)),expression(PCA~(paste(p,"*",sep="")==50)),expression(PCA~(paste(p,"*",sep="")==100)),expression(PCA~(paste(p,"*",sep="")==200))),lty=c(1,2,3,1,2,3),lwd=c(3,3,3,3,3,3),col=c("black","black","black","gray55","gray55","gray55"),bty="n",cex=1.7)#
plot(1-pca.roc.p50[,3],pca.roc.p50[,2],type="l",cex=1.7,main="(b)",xlab="1-Specificity",ylab="Sensitivity",col="gray60",lty=1,lwd=3)#
lines(1-pca.roc.p100[,3],pca.roc.p100[,2],col="gray60",lty=2,lwd=3)#
lines(1-pca.roc.p200[,3],pca.roc.p200[,2],col="gray60",lty=3,lwd=3)#
lines(1-rasch.roc.p50[,3],rasch.roc.p50[,2],col="black",lty=1,lwd=3)#
lines(1-rasch.roc.p100[,3],rasch.roc.p100[,2],col="black",lty=2,lwd=3)#
lines(1-rasch.roc.p200[,3],rasch.roc.p200[,2],col="black",lty=3,lwd=3)#
legend("bottomright",c(expression(RM~(paste(p,"*",sep="")==50)),expression(RM~(paste(p,"*",sep="")==100)),expression(RM~(paste(p,"*",sep="")==200)),expression(PCA~(paste(p,"*",sep="")==50)),expression(PCA~(paste(p,"*",sep="")==100)),expression(PCA~(paste(p,"*",sep="")==200))),lty=c(1,2,3,1,2,3),lwd=c(3,3,3,3,3,3),col=c("black","black","black","gray55","gray55","gray55"),bty="n",cex=1.7)#
dev.off()
pdf("roc2_bw_temp.pdf",width=20,height=10,paper="special")#
par(cex.axis=1.7,cex.lab=1.7,cex.main=1.7,pty="s",mfrow=c(1,2))#
plot(1-pca.roc.r50[,3],pca.roc.r50[,2],type="l",cex=1.7,main="(a)",xlab="1-Specificity",ylab="Sensitivity",col="gray55",lty=1,lwd=3)#
lines(1-pca.roc.r100[,3],pca.roc.r100[,2],col="gray55",lty=2,lwd=3)#
lines(1-pca.roc.r200[,3],pca.roc.r200[,2],col="gray55",lty=3,lwd=3)#
lines(1-rasch.roc.r50[,3],rasch.roc.r50[,2],col="black",lty=1,lwd=3)#
lines(1-rasch.roc.r100[,3],rasch.roc.r100[,2],col="black",lty=2,lwd=3)#
lines(1-rasch.roc.r200[,3],rasch.roc.r200[,2],col="black",lty=3,lwd=3)#
legend("bottomright",c(expression(paste("RM",[{"p*"}])),expression(RM~(paste(p,"*",sep="")==100)),expression(RM~(paste(p,"*",sep="")==200)),expression(PCA~(paste(p,"*",sep="")==50)),expression(PCA~(paste(p,"*",sep="")==100)),expression(PCA~(paste(p,"*",sep="")==200))),lty=c(1,2,3,1,2,3),lwd=c(3,3,3,3,3,3),col=c("black","black","black","gray55","gray55","gray55"),bty="n",cex=1.7)#
plot(1-pca.roc.p50[,3],pca.roc.p50[,2],type="l",cex=1.7,main="(b)",xlab="1-Specificity",ylab="Sensitivity",col="gray60",lty=1,lwd=3)#
lines(1-pca.roc.p100[,3],pca.roc.p100[,2],col="gray60",lty=2,lwd=3)#
lines(1-pca.roc.p200[,3],pca.roc.p200[,2],col="gray60",lty=3,lwd=3)#
lines(1-rasch.roc.p50[,3],rasch.roc.p50[,2],col="black",lty=1,lwd=3)#
lines(1-rasch.roc.p100[,3],rasch.roc.p100[,2],col="black",lty=2,lwd=3)#
lines(1-rasch.roc.p200[,3],rasch.roc.p200[,2],col="black",lty=3,lwd=3)#
legend("bottomright",c(expression(RM~(paste(p,"*",sep="")==50)),expression(RM~(paste(p,"*",sep="")==100)),expression(RM~(paste(p,"*",sep="")==200)),expression(PCA~(paste(p,"*",sep="")==50)),expression(PCA~(paste(p,"*",sep="")==100)),expression(PCA~(paste(p,"*",sep="")==200))),lty=c(1,2,3,1,2,3),lwd=c(3,3,3,3,3,3),col=c("black","black","black","gray55","gray55","gray55"),bty="n",cex=1.7)#
dev.off()
pdf("roc2_bw_temp.pdf",width=20,height=10,paper="special")#
par(cex.axis=1.7,cex.lab=1.7,cex.main=1.7,pty="s",mfrow=c(1,2))#
plot(1-pca.roc.r50[,3],pca.roc.r50[,2],type="l",cex=1.7,main="(a)",xlab="1-Specificity",ylab="Sensitivity",col="gray55",lty=1,lwd=3)#
lines(1-pca.roc.r100[,3],pca.roc.r100[,2],col="gray55",lty=2,lwd=3)#
lines(1-pca.roc.r200[,3],pca.roc.r200[,2],col="gray55",lty=3,lwd=3)#
lines(1-rasch.roc.r50[,3],rasch.roc.r50[,2],col="black",lty=1,lwd=3)#
lines(1-rasch.roc.r100[,3],rasch.roc.r100[,2],col="black",lty=2,lwd=3)#
lines(1-rasch.roc.r200[,3],rasch.roc.r200[,2],col="black",lty=3,lwd=3)#
legend("bottomright",c(expression(paste("RM",[{,"p*",}])),expression(RM~(paste(p,"*",sep="")==100)),expression(RM~(paste(p,"*",sep="")==200)),expression(PCA~(paste(p,"*",sep="")==50)),expression(PCA~(paste(p,"*",sep="")==100)),expression(PCA~(paste(p,"*",sep="")==200))),lty=c(1,2,3,1,2,3),lwd=c(3,3,3,3,3,3),col=c("black","black","black","gray55","gray55","gray55"),bty="n",cex=1.7)#
plot(1-pca.roc.p50[,3],pca.roc.p50[,2],type="l",cex=1.7,main="(b)",xlab="1-Specificity",ylab="Sensitivity",col="gray60",lty=1,lwd=3)#
lines(1-pca.roc.p100[,3],pca.roc.p100[,2],col="gray60",lty=2,lwd=3)#
lines(1-pca.roc.p200[,3],pca.roc.p200[,2],col="gray60",lty=3,lwd=3)#
lines(1-rasch.roc.p50[,3],rasch.roc.p50[,2],col="black",lty=1,lwd=3)#
lines(1-rasch.roc.p100[,3],rasch.roc.p100[,2],col="black",lty=2,lwd=3)#
lines(1-rasch.roc.p200[,3],rasch.roc.p200[,2],col="black",lty=3,lwd=3)#
legend("bottomright",c(expression(RM~(paste(p,"*",sep="")==50)),expression(RM~(paste(p,"*",sep="")==100)),expression(RM~(paste(p,"*",sep="")==200)),expression(PCA~(paste(p,"*",sep="")==50)),expression(PCA~(paste(p,"*",sep="")==100)),expression(PCA~(paste(p,"*",sep="")==200))),lty=c(1,2,3,1,2,3),lwd=c(3,3,3,3,3,3),col=c("black","black","black","gray55","gray55","gray55"),bty="n",cex=1.7)#
dev.off()
pdf("roc2_bw_temp.pdf",width=20,height=10,paper="special")#
par(cex.axis=1.7,cex.lab=1.7,cex.main=1.7,pty="s",mfrow=c(1,2))#
plot(1-pca.roc.r50[,3],pca.roc.r50[,2],type="l",cex=1.7,main="(a)",xlab="1-Specificity",ylab="Sensitivity",col="gray55",lty=1,lwd=3)#
lines(1-pca.roc.r100[,3],pca.roc.r100[,2],col="gray55",lty=2,lwd=3)#
lines(1-pca.roc.r200[,3],pca.roc.r200[,2],col="gray55",lty=3,lwd=3)#
lines(1-rasch.roc.r50[,3],rasch.roc.r50[,2],col="black",lty=1,lwd=3)#
lines(1-rasch.roc.r100[,3],rasch.roc.r100[,2],col="black",lty=2,lwd=3)#
lines(1-rasch.roc.r200[,3],rasch.roc.r200[,2],col="black",lty=3,lwd=3)#
legend("bottomright",c(expression(paste("RM (p* ",==50,")")),expression(RM~(paste(p,"*",sep="")==100)),expression(RM~(paste(p,"*",sep="")==200)),expression(PCA~(paste(p,"*",sep="")==50)),expression(PCA~(paste(p,"*",sep="")==100)),expression(PCA~(paste(p,"*",sep="")==200))),lty=c(1,2,3,1,2,3),lwd=c(3,3,3,3,3,3),col=c("black","black","black","gray55","gray55","gray55"),bty="n",cex=1.7)#
plot(1-pca.roc.p50[,3],pca.roc.p50[,2],type="l",cex=1.7,main="(b)",xlab="1-Specificity",ylab="Sensitivity",col="gray60",lty=1,lwd=3)#
lines(1-pca.roc.p100[,3],pca.roc.p100[,2],col="gray60",lty=2,lwd=3)#
lines(1-pca.roc.p200[,3],pca.roc.p200[,2],col="gray60",lty=3,lwd=3)#
lines(1-rasch.roc.p50[,3],rasch.roc.p50[,2],col="black",lty=1,lwd=3)#
lines(1-rasch.roc.p100[,3],rasch.roc.p100[,2],col="black",lty=2,lwd=3)#
lines(1-rasch.roc.p200[,3],rasch.roc.p200[,2],col="black",lty=3,lwd=3)#
legend("bottomright",c(expression(RM~(paste(p,"*",sep="")==50)),expression(RM~(paste(p,"*",sep="")==100)),expression(RM~(paste(p,"*",sep="")==200)),expression(PCA~(paste(p,"*",sep="")==50)),expression(PCA~(paste(p,"*",sep="")==100)),expression(PCA~(paste(p,"*",sep="")==200))),lty=c(1,2,3,1,2,3),lwd=c(3,3,3,3,3,3),col=c("black","black","black","gray55","gray55","gray55"),bty="n",cex=1.7)#
dev.off()
pdf("roc2_bw_temp.pdf",width=20,height=10,paper="special")#
par(cex.axis=1.7,cex.lab=1.7,cex.main=1.7,pty="s",mfrow=c(1,2))#
plot(1-pca.roc.r50[,3],pca.roc.r50[,2],type="l",cex=1.7,main="(a)",xlab="1-Specificity",ylab="Sensitivity",col="gray55",lty=1,lwd=3)#
lines(1-pca.roc.r100[,3],pca.roc.r100[,2],col="gray55",lty=2,lwd=3)#
lines(1-pca.roc.r200[,3],pca.roc.r200[,2],col="gray55",lty=3,lwd=3)#
lines(1-rasch.roc.r50[,3],rasch.roc.r50[,2],col="black",lty=1,lwd=3)#
lines(1-rasch.roc.r100[,3],rasch.roc.r100[,2],col="black",lty=2,lwd=3)#
lines(1-rasch.roc.r200[,3],rasch.roc.r200[,2],col="black",lty=3,lwd=3)#
legend("bottomright",c(expression(paste("RM (p*=50)")),expression(RM~(paste(p,"*",sep="")==100)),expression(RM~(paste(p,"*",sep="")==200)),expression(PCA~(paste(p,"*",sep="")==50)),expression(PCA~(paste(p,"*",sep="")==100)),expression(PCA~(paste(p,"*",sep="")==200))),lty=c(1,2,3,1,2,3),lwd=c(3,3,3,3,3,3),col=c("black","black","black","gray55","gray55","gray55"),bty="n",cex=1.7)#
plot(1-pca.roc.p50[,3],pca.roc.p50[,2],type="l",cex=1.7,main="(b)",xlab="1-Specificity",ylab="Sensitivity",col="gray60",lty=1,lwd=3)#
lines(1-pca.roc.p100[,3],pca.roc.p100[,2],col="gray60",lty=2,lwd=3)#
lines(1-pca.roc.p200[,3],pca.roc.p200[,2],col="gray60",lty=3,lwd=3)#
lines(1-rasch.roc.p50[,3],rasch.roc.p50[,2],col="black",lty=1,lwd=3)#
lines(1-rasch.roc.p100[,3],rasch.roc.p100[,2],col="black",lty=2,lwd=3)#
lines(1-rasch.roc.p200[,3],rasch.roc.p200[,2],col="black",lty=3,lwd=3)#
legend("bottomright",c(expression(RM~(paste(p,"*",sep="")==50)),expression(RM~(paste(p,"*",sep="")==100)),expression(RM~(paste(p,"*",sep="")==200)),expression(PCA~(paste(p,"*",sep="")==50)),expression(PCA~(paste(p,"*",sep="")==100)),expression(PCA~(paste(p,"*",sep="")==200))),lty=c(1,2,3,1,2,3),lwd=c(3,3,3,3,3,3),col=c("black","black","black","gray55","gray55","gray55"),bty="n",cex=1.7)#
dev.off()
demo(plotmath)
pdf("roc2_bw_temp.pdf",width=20,height=10,paper="special")#
par(cex.axis=1.7,cex.lab=1.7,cex.main=1.7,pty="s",mfrow=c(1,2))#
plot(1-pca.roc.r50[,3],pca.roc.r50[,2],type="l",cex=1.7,main="(a)",xlab="1-Specificity",ylab="Sensitivity",col="gray55",lty=1,lwd=3)#
lines(1-pca.roc.r100[,3],pca.roc.r100[,2],col="gray55",lty=2,lwd=3)#
lines(1-pca.roc.r200[,3],pca.roc.r200[,2],col="gray55",lty=3,lwd=3)#
lines(1-rasch.roc.r50[,3],rasch.roc.r50[,2],col="black",lty=1,lwd=3)#
lines(1-rasch.roc.r100[,3],rasch.roc.r100[,2],col="black",lty=2,lwd=3)#
lines(1-rasch.roc.r200[,3],rasch.roc.r200[,2],col="black",lty=3,lwd=3)#
legend("bottomright",c(expression(plain(p)),expression(RM~(paste(p,"*",sep="")==100)),expression(RM~(paste(p,"*",sep="")==200)),expression(PCA~(paste(p,"*",sep="")==50)),expression(PCA~(paste(p,"*",sep="")==100)),expression(PCA~(paste(p,"*",sep="")==200))),lty=c(1,2,3,1,2,3),lwd=c(3,3,3,3,3,3),col=c("black","black","black","gray55","gray55","gray55"),bty="n",cex=1.7)#
plot(1-pca.roc.p50[,3],pca.roc.p50[,2],type="l",cex=1.7,main="(b)",xlab="1-Specificity",ylab="Sensitivity",col="gray60",lty=1,lwd=3)#
lines(1-pca.roc.p100[,3],pca.roc.p100[,2],col="gray60",lty=2,lwd=3)#
lines(1-pca.roc.p200[,3],pca.roc.p200[,2],col="gray60",lty=3,lwd=3)#
lines(1-rasch.roc.p50[,3],rasch.roc.p50[,2],col="black",lty=1,lwd=3)#
lines(1-rasch.roc.p100[,3],rasch.roc.p100[,2],col="black",lty=2,lwd=3)#
lines(1-rasch.roc.p200[,3],rasch.roc.p200[,2],col="black",lty=3,lwd=3)#
legend("bottomright",c(expression(RM~(paste(p,"*",sep="")==50)),expression(RM~(paste(p,"*",sep="")==100)),expression(RM~(paste(p,"*",sep="")==200)),expression(PCA~(paste(p,"*",sep="")==50)),expression(PCA~(paste(p,"*",sep="")==100)),expression(PCA~(paste(p,"*",sep="")==200))),lty=c(1,2,3,1,2,3),lwd=c(3,3,3,3,3,3),col=c("black","black","black","gray55","gray55","gray55"),bty="n",cex=1.7)#
dev.off()
pdf("roc2_bw_temp.pdf",width=20,height=10,paper="special")#
par(cex.axis=1.7,cex.lab=1.7,cex.main=1.7,pty="s",mfrow=c(1,2))#
plot(1-pca.roc.r50[,3],pca.roc.r50[,2],type="l",cex=1.7,main="(a)",xlab="1-Specificity",ylab="Sensitivity",col="gray55",lty=1,lwd=3)#
lines(1-pca.roc.r100[,3],pca.roc.r100[,2],col="gray55",lty=2,lwd=3)#
lines(1-pca.roc.r200[,3],pca.roc.r200[,2],col="gray55",lty=3,lwd=3)#
lines(1-rasch.roc.r50[,3],rasch.roc.r50[,2],col="black",lty=1,lwd=3)#
lines(1-rasch.roc.r100[,3],rasch.roc.r100[,2],col="black",lty=2,lwd=3)#
lines(1-rasch.roc.r200[,3],rasch.roc.r200[,2],col="black",lty=3,lwd=3)#
legend("bottomright",c(expression(plain(RM),plain(p*)),expression(RM~(paste(p,"*",sep="")==100)),expression(RM~(paste(p,"*",sep="")==200)),expression(PCA~(paste(p,"*",sep="")==50)),expression(PCA~(paste(p,"*",sep="")==100)),expression(PCA~(paste(p,"*",sep="")==200))),lty=c(1,2,3,1,2,3),lwd=c(3,3,3,3,3,3),col=c("black","black","black","gray55","gray55","gray55"),bty="n",cex=1.7)#
plot(1-pca.roc.p50[,3],pca.roc.p50[,2],type="l",cex=1.7,main="(b)",xlab="1-Specificity",ylab="Sensitivity",col="gray60",lty=1,lwd=3)#
lines(1-pca.roc.p100[,3],pca.roc.p100[,2],col="gray60",lty=2,lwd=3)#
lines(1-pca.roc.p200[,3],pca.roc.p200[,2],col="gray60",lty=3,lwd=3)#
lines(1-rasch.roc.p50[,3],rasch.roc.p50[,2],col="black",lty=1,lwd=3)#
lines(1-rasch.roc.p100[,3],rasch.roc.p100[,2],col="black",lty=2,lwd=3)#
lines(1-rasch.roc.p200[,3],rasch.roc.p200[,2],col="black",lty=3,lwd=3)#
legend("bottomright",c(expression(RM~(paste(p,"*",sep="")==50)),expression(RM~(paste(p,"*",sep="")==100)),expression(RM~(paste(p,"*",sep="")==200)),expression(PCA~(paste(p,"*",sep="")==50)),expression(PCA~(paste(p,"*",sep="")==100)),expression(PCA~(paste(p,"*",sep="")==200))),lty=c(1,2,3,1,2,3),lwd=c(3,3,3,3,3,3),col=c("black","black","black","gray55","gray55","gray55"),bty="n",cex=1.7)#
dev.off()
pdf("roc2_bw_temp.pdf",width=20,height=10,paper="special")#
par(cex.axis=1.7,cex.lab=1.7,cex.main=1.7,pty="s",mfrow=c(1,2))#
plot(1-pca.roc.r50[,3],pca.roc.r50[,2],type="l",cex=1.7,main="(a)",xlab="1-Specificity",ylab="Sensitivity",col="gray55",lty=1,lwd=3)#
lines(1-pca.roc.r100[,3],pca.roc.r100[,2],col="gray55",lty=2,lwd=3)#
lines(1-pca.roc.r200[,3],pca.roc.r200[,2],col="gray55",lty=3,lwd=3)#
lines(1-rasch.roc.r50[,3],rasch.roc.r50[,2],col="black",lty=1,lwd=3)#
lines(1-rasch.roc.r100[,3],rasch.roc.r100[,2],col="black",lty=2,lwd=3)#
lines(1-rasch.roc.r200[,3],rasch.roc.r200[,2],col="black",lty=3,lwd=3)#
legend("bottomright",c(expression(plain(p*)),expression(RM~(paste(p,"*",sep="")==100)),expression(RM~(paste(p,"*",sep="")==200)),expression(PCA~(paste(p,"*",sep="")==50)),expression(PCA~(paste(p,"*",sep="")==100)),expression(PCA~(paste(p,"*",sep="")==200))),lty=c(1,2,3,1,2,3),lwd=c(3,3,3,3,3,3),col=c("black","black","black","gray55","gray55","gray55"),bty="n",cex=1.7)#
plot(1-pca.roc.p50[,3],pca.roc.p50[,2],type="l",cex=1.7,main="(b)",xlab="1-Specificity",ylab="Sensitivity",col="gray60",lty=1,lwd=3)#
lines(1-pca.roc.p100[,3],pca.roc.p100[,2],col="gray60",lty=2,lwd=3)#
lines(1-pca.roc.p200[,3],pca.roc.p200[,2],col="gray60",lty=3,lwd=3)#
lines(1-rasch.roc.p50[,3],rasch.roc.p50[,2],col="black",lty=1,lwd=3)#
lines(1-rasch.roc.p100[,3],rasch.roc.p100[,2],col="black",lty=2,lwd=3)#
lines(1-rasch.roc.p200[,3],rasch.roc.p200[,2],col="black",lty=3,lwd=3)#
legend("bottomright",c(expression(RM~(paste(p,"*",sep="")==50)),expression(RM~(paste(p,"*",sep="")==100)),expression(RM~(paste(p,"*",sep="")==200)),expression(PCA~(paste(p,"*",sep="")==50)),expression(PCA~(paste(p,"*",sep="")==100)),expression(PCA~(paste(p,"*",sep="")==200))),lty=c(1,2,3,1,2,3),lwd=c(3,3,3,3,3,3),col=c("black","black","black","gray55","gray55","gray55"),bty="n",cex=1.7)#
dev.off()
pdf("roc2_bw_temp.pdf",width=20,height=10,paper="special")#
par(cex.axis=1.7,cex.lab=1.7,cex.main=1.7,pty="s",mfrow=c(1,2))#
plot(1-pca.roc.r50[,3],pca.roc.r50[,2],type="l",cex=1.7,main="(a)",xlab="1-Specificity",ylab="Sensitivity",col="gray55",lty=1,lwd=3)#
lines(1-pca.roc.r100[,3],pca.roc.r100[,2],col="gray55",lty=2,lwd=3)#
lines(1-pca.roc.r200[,3],pca.roc.r200[,2],col="gray55",lty=3,lwd=3)#
lines(1-rasch.roc.r50[,3],rasch.roc.r50[,2],col="black",lty=1,lwd=3)#
lines(1-rasch.roc.r100[,3],rasch.roc.r100[,2],col="black",lty=2,lwd=3)#
lines(1-rasch.roc.r200[,3],rasch.roc.r200[,2],col="black",lty=3,lwd=3)#
legend("bottomright",c(expression(plain("p*")),expression(RM~(paste(p,"*",sep="")==100)),expression(RM~(paste(p,"*",sep="")==200)),expression(PCA~(paste(p,"*",sep="")==50)),expression(PCA~(paste(p,"*",sep="")==100)),expression(PCA~(paste(p,"*",sep="")==200))),lty=c(1,2,3,1,2,3),lwd=c(3,3,3,3,3,3),col=c("black","black","black","gray55","gray55","gray55"),bty="n",cex=1.7)#
plot(1-pca.roc.p50[,3],pca.roc.p50[,2],type="l",cex=1.7,main="(b)",xlab="1-Specificity",ylab="Sensitivity",col="gray60",lty=1,lwd=3)#
lines(1-pca.roc.p100[,3],pca.roc.p100[,2],col="gray60",lty=2,lwd=3)#
lines(1-pca.roc.p200[,3],pca.roc.p200[,2],col="gray60",lty=3,lwd=3)#
lines(1-rasch.roc.p50[,3],rasch.roc.p50[,2],col="black",lty=1,lwd=3)#
lines(1-rasch.roc.p100[,3],rasch.roc.p100[,2],col="black",lty=2,lwd=3)#
lines(1-rasch.roc.p200[,3],rasch.roc.p200[,2],col="black",lty=3,lwd=3)#
legend("bottomright",c(expression(RM~(paste(p,"*",sep="")==50)),expression(RM~(paste(p,"*",sep="")==100)),expression(RM~(paste(p,"*",sep="")==200)),expression(PCA~(paste(p,"*",sep="")==50)),expression(PCA~(paste(p,"*",sep="")==100)),expression(PCA~(paste(p,"*",sep="")==200))),lty=c(1,2,3,1,2,3),lwd=c(3,3,3,3,3,3),col=c("black","black","black","gray55","gray55","gray55"),bty="n",cex=1.7)#
dev.off()
?plotmath
pdf("roc2_bw_temp.pdf",width=20,height=10,paper="special")#
par(cex.axis=1.7,cex.lab=1.7,cex.main=1.7,pty="s",mfrow=c(1,2))#
plot(1-pca.roc.r50[,3],pca.roc.r50[,2],type="l",cex=1.7,main="(a)",xlab="1-Specificity",ylab="Sensitivity",col="gray55",lty=1,lwd=3)#
lines(1-pca.roc.r100[,3],pca.roc.r100[,2],col="gray55",lty=2,lwd=3)#
lines(1-pca.roc.r200[,3],pca.roc.r200[,2],col="gray55",lty=3,lwd=3)#
lines(1-rasch.roc.r50[,3],rasch.roc.r50[,2],col="black",lty=1,lwd=3)#
lines(1-rasch.roc.r100[,3],rasch.roc.r100[,2],col="black",lty=2,lwd=3)#
lines(1-rasch.roc.r200[,3],rasch.roc.r200[,2],col="black",lty=3,lwd=3)#
legend("bottomright",c(expression(plain("p*")==50),expression(RM~(paste(p,"*",sep="")==100)),expression(RM~(paste(p,"*",sep="")==200)),expression(PCA~(paste(p,"*",sep="")==50)),expression(PCA~(paste(p,"*",sep="")==100)),expression(PCA~(paste(p,"*",sep="")==200))),lty=c(1,2,3,1,2,3),lwd=c(3,3,3,3,3,3),col=c("black","black","black","gray55","gray55","gray55"),bty="n",cex=1.7)#
plot(1-pca.roc.p50[,3],pca.roc.p50[,2],type="l",cex=1.7,main="(b)",xlab="1-Specificity",ylab="Sensitivity",col="gray60",lty=1,lwd=3)#
lines(1-pca.roc.p100[,3],pca.roc.p100[,2],col="gray60",lty=2,lwd=3)#
lines(1-pca.roc.p200[,3],pca.roc.p200[,2],col="gray60",lty=3,lwd=3)#
lines(1-rasch.roc.p50[,3],rasch.roc.p50[,2],col="black",lty=1,lwd=3)#
lines(1-rasch.roc.p100[,3],rasch.roc.p100[,2],col="black",lty=2,lwd=3)#
lines(1-rasch.roc.p200[,3],rasch.roc.p200[,2],col="black",lty=3,lwd=3)#
legend("bottomright",c(expression(RM~(paste(p,"*",sep="")==50)),expression(RM~(paste(p,"*",sep="")==100)),expression(RM~(paste(p,"*",sep="")==200)),expression(PCA~(paste(p,"*",sep="")==50)),expression(PCA~(paste(p,"*",sep="")==100)),expression(PCA~(paste(p,"*",sep="")==200))),lty=c(1,2,3,1,2,3),lwd=c(3,3,3,3,3,3),col=c("black","black","black","gray55","gray55","gray55"),bty="n",cex=1.7)#
dev.off()
pdf("roc2_bw_temp.pdf",width=20,height=10,paper="special")#
par(cex.axis=1.7,cex.lab=1.7,cex.main=1.7,pty="s",mfrow=c(1,2))#
plot(1-pca.roc.r50[,3],pca.roc.r50[,2],type="l",cex=1.7,main="(a)",xlab="1-Specificity",ylab="Sensitivity",col="gray55",lty=1,lwd=3)#
lines(1-pca.roc.r100[,3],pca.roc.r100[,2],col="gray55",lty=2,lwd=3)#
lines(1-pca.roc.r200[,3],pca.roc.r200[,2],col="gray55",lty=3,lwd=3)#
lines(1-rasch.roc.r50[,3],rasch.roc.r50[,2],col="black",lty=1,lwd=3)#
lines(1-rasch.roc.r100[,3],rasch.roc.r100[,2],col="black",lty=2,lwd=3)#
lines(1-rasch.roc.r200[,3],rasch.roc.r200[,2],col="black",lty=3,lwd=3)#
legend("bottomright",c(expression([{plain("p*")==50}]),expression(RM~(paste(p,"*",sep="")==100)),expression(RM~(paste(p,"*",sep="")==200)),expression(PCA~(paste(p,"*",sep="")==50)),expression(PCA~(paste(p,"*",sep="")==100)),expression(PCA~(paste(p,"*",sep="")==200))),lty=c(1,2,3,1,2,3),lwd=c(3,3,3,3,3,3),col=c("black","black","black","gray55","gray55","gray55"),bty="n",cex=1.7)#
plot(1-pca.roc.p50[,3],pca.roc.p50[,2],type="l",cex=1.7,main="(b)",xlab="1-Specificity",ylab="Sensitivity",col="gray60",lty=1,lwd=3)#
lines(1-pca.roc.p100[,3],pca.roc.p100[,2],col="gray60",lty=2,lwd=3)#
lines(1-pca.roc.p200[,3],pca.roc.p200[,2],col="gray60",lty=3,lwd=3)#
lines(1-rasch.roc.p50[,3],rasch.roc.p50[,2],col="black",lty=1,lwd=3)#
lines(1-rasch.roc.p100[,3],rasch.roc.p100[,2],col="black",lty=2,lwd=3)#
lines(1-rasch.roc.p200[,3],rasch.roc.p200[,2],col="black",lty=3,lwd=3)#
legend("bottomright",c(expression(RM~(paste(p,"*",sep="")==50)),expression(RM~(paste(p,"*",sep="")==100)),expression(RM~(paste(p,"*",sep="")==200)),expression(PCA~(paste(p,"*",sep="")==50)),expression(PCA~(paste(p,"*",sep="")==100)),expression(PCA~(paste(p,"*",sep="")==200))),lty=c(1,2,3,1,2,3),lwd=c(3,3,3,3,3,3),col=c("black","black","black","gray55","gray55","gray55"),bty="n",cex=1.7)#
dev.off()
pdf("roc2_bw_temp.pdf",width=20,height=10,paper="special")#
par(cex.axis=1.7,cex.lab=1.7,cex.main=1.7,pty="s",mfrow=c(1,2))#
plot(1-pca.roc.r50[,3],pca.roc.r50[,2],type="l",cex=1.7,main="(a)",xlab="1-Specificity",ylab="Sensitivity",col="gray55",lty=1,lwd=3)#
lines(1-pca.roc.r100[,3],pca.roc.r100[,2],col="gray55",lty=2,lwd=3)#
lines(1-pca.roc.r200[,3],pca.roc.r200[,2],col="gray55",lty=3,lwd=3)#
lines(1-rasch.roc.r50[,3],rasch.roc.r50[,2],col="black",lty=1,lwd=3)#
lines(1-rasch.roc.r100[,3],rasch.roc.r100[,2],col="black",lty=2,lwd=3)#
lines(1-rasch.roc.r200[,3],rasch.roc.r200[,2],col="black",lty=3,lwd=3)#
legend("bottomright",c(expression(RM plain("p*")==50),expression(RM~(paste(p,"*",sep="")==100)),expression(RM~(paste(p,"*",sep="")==200)),expression(PCA~(paste(p,"*",sep="")==50)),expression(PCA~(paste(p,"*",sep="")==100)),expression(PCA~(paste(p,"*",sep="")==200))),lty=c(1,2,3,1,2,3),lwd=c(3,3,3,3,3,3),col=c("black","black","black","gray55","gray55","gray55"),bty="n",cex=1.7)#
plot(1-pca.roc.p50[,3],pca.roc.p50[,2],type="l",cex=1.7,main="(b)",xlab="1-Specificity",ylab="Sensitivity",col="gray60",lty=1,lwd=3)#
lines(1-pca.roc.p100[,3],pca.roc.p100[,2],col="gray60",lty=2,lwd=3)#
lines(1-pca.roc.p200[,3],pca.roc.p200[,2],col="gray60",lty=3,lwd=3)#
lines(1-rasch.roc.p50[,3],rasch.roc.p50[,2],col="black",lty=1,lwd=3)#
lines(1-rasch.roc.p100[,3],rasch.roc.p100[,2],col="black",lty=2,lwd=3)#
lines(1-rasch.roc.p200[,3],rasch.roc.p200[,2],col="black",lty=3,lwd=3)#
legend("bottomright",c(expression(RM~(paste(p,"*",sep="")==50)),expression(RM~(paste(p,"*",sep="")==100)),expression(RM~(paste(p,"*",sep="")==200)),expression(PCA~(paste(p,"*",sep="")==50)),expression(PCA~(paste(p,"*",sep="")==100)),expression(PCA~(paste(p,"*",sep="")==200))),lty=c(1,2,3,1,2,3),lwd=c(3,3,3,3,3,3),col=c("black","black","black","gray55","gray55","gray55"),bty="n",cex=1.7)#
dev.off()
pdf("roc2_bw_temp.pdf",width=20,height=10,paper="special")#
par(cex.axis=1.7,cex.lab=1.7,cex.main=1.7,pty="s",mfrow=c(1,2))#
plot(1-pca.roc.r50[,3],pca.roc.r50[,2],type="l",cex=1.7,main="(a)",xlab="1-Specificity",ylab="Sensitivity",col="gray55",lty=1,lwd=3)#
lines(1-pca.roc.r100[,3],pca.roc.r100[,2],col="gray55",lty=2,lwd=3)#
lines(1-pca.roc.r200[,3],pca.roc.r200[,2],col="gray55",lty=3,lwd=3)#
lines(1-rasch.roc.r50[,3],rasch.roc.r50[,2],col="black",lty=1,lwd=3)#
lines(1-rasch.roc.r100[,3],rasch.roc.r100[,2],col="black",lty=2,lwd=3)#
lines(1-rasch.roc.r200[,3],rasch.roc.r200[,2],col="black",lty=3,lwd=3)#
legend("bottomright",c(expression(RM~(paste(p,"*",sep="")==50)),expression(RM~(paste(p,"*",sep="")==100)),expression(RM~(paste(p,"*",sep="")==200)),expression(PCA~(paste(p,"*",sep="")==50)),expression(PCA~(paste(p,"*",sep="")==100)),expression(PCA~(paste(p,"*",sep="")==200))),lty=c(1,2,3,1,2,3),lwd=c(3,3,3,3,3,3),col=c("black","black","black","gray55","gray55","gray55"),bty="n",cex=1.7)#
plot(1-pca.roc.p50[,3],pca.roc.p50[,2],type="l",cex=1.7,main="(b)",xlab="1-Specificity",ylab="Sensitivity",col="gray60",lty=1,lwd=3)#
lines(1-pca.roc.p100[,3],pca.roc.p100[,2],col="gray60",lty=2,lwd=3)#
lines(1-pca.roc.p200[,3],pca.roc.p200[,2],col="gray60",lty=3,lwd=3)#
lines(1-rasch.roc.p50[,3],rasch.roc.p50[,2],col="black",lty=1,lwd=3)#
lines(1-rasch.roc.p100[,3],rasch.roc.p100[,2],col="black",lty=2,lwd=3)#
lines(1-rasch.roc.p200[,3],rasch.roc.p200[,2],col="black",lty=3,lwd=3)#
legend("bottomright",c(expression(RM~(paste(p,"*",sep="")==50)),expression(RM~(paste(p,"*",sep="")==100)),expression(RM~(paste(p,"*",sep="")==200)),expression(PCA~(paste(p,"*",sep="")==50)),expression(PCA~(paste(p,"*",sep="")==100)),expression(PCA~(paste(p,"*",sep="")==200))),lty=c(1,2,3,1,2,3),lwd=c(3,3,3,3,3,3),col=c("black","black","black","gray55","gray55","gray55"),bty="n",cex=1.7)#
dev.off()
pdf("roc2_bw_temp.pdf",width=20,height=10,paper="special")#
par(cex.axis=1.7,cex.lab=1.7,cex.main=1.7,pty="s",mfrow=c(1,2))#
plot(1-pca.roc.r50[,3],pca.roc.r50[,2],type="l",cex=1.7,main="(a)",xlab="1-Specificity",ylab="Sensitivity",col="gray55",lty=1,lwd=3)#
lines(1-pca.roc.r100[,3],pca.roc.r100[,2],col="gray55",lty=2,lwd=3)#
lines(1-pca.roc.r200[,3],pca.roc.r200[,2],col="gray55",lty=3,lwd=3)#
lines(1-rasch.roc.r50[,3],rasch.roc.r50[,2],col="black",lty=1,lwd=3)#
lines(1-rasch.roc.r100[,3],rasch.roc.r100[,2],col="black",lty=2,lwd=3)#
lines(1-rasch.roc.r200[,3],rasch.roc.r200[,2],col="black",lty=3,lwd=3)#
legend("bottomright",c(expression(RM~(paste(p,"*",sep="")==50)),expression(RM~(paste(p,"*",sep="")==100)),expression(RM~(paste(p,"*",sep="")==200)),expression(PCA~(paste(p,"*",sep="")==50)),expression(PCA~(paste(p,"*",sep="")==100)),expression(PCA~(paste(p,"*",sep="")==200))),lty=c(1,2,3,1,2,3),lwd=c(3,3,3,3,3,3),col=c("black","black","black","gray55","gray55","gray55"),bty="n",cex=1.7)#
plot(1-pca.roc.p50[,3],pca.roc.p50[,2],type="l",cex=1.7,main="(b)",xlab="1-Specificity",ylab="Sensitivity",col="gray60",lty=1,lwd=3)#
lines(1-pca.roc.p100[,3],pca.roc.p100[,2],col="gray60",lty=2,lwd=3)#
lines(1-pca.roc.p200[,3],pca.roc.p200[,2],col="gray60",lty=3,lwd=3)#
lines(1-rasch.roc.p50[,3],rasch.roc.p50[,2],col="black",lty=1,lwd=3)#
lines(1-rasch.roc.p100[,3],rasch.roc.p100[,2],col="black",lty=2,lwd=3)#
lines(1-rasch.roc.p200[,3],rasch.roc.p200[,2],col="black",lty=3,lwd=3)#
legend("bottomright",c(expression(RM~(paste(p,"*",sep="")==50)),expression(RM~(paste(p,"*",sep="")==100)),expression(RM~(paste(p,"*",sep="")==200)),expression(PCA~(paste(p,"*",sep="")==50)),expression(PCA~(paste(p,"*",sep="")==100)),expression(PCA~(paste(p,"*",sep="")==200))),lty=c(1,2,3,1,2,3),lwd=c(3,3,3,3,3,3),col=c("black","black","black","gray55","gray55","gray55"),bty="n",cex=1.7)#
dev.off()#
#
pdf("roc2_cl_temp.pdf",width=20,height=10,paper="special")#
par(cex.axis=1.7,cex.lab=1.7,cex.main=1.7,pty="s",mfrow=c(1,2))#
plot(1-pca.roc.r50[,3],pca.roc.r50[,2],type="l",cex=1.7,main="(a)",xlab="1-Specificity",ylab="Sensitivity",col="blue",lty=1,lwd=3)#
lines(1-pca.roc.r100[,3],pca.roc.r100[,2],col="blue",lty=2,lwd=3)#
lines(1-pca.roc.r200[,3],pca.roc.r200[,2],col="blue",lty=3,lwd=3)#
lines(1-rasch.roc.r50[,3],rasch.roc.r50[,2],col="red",lty=1,lwd=3)#
lines(1-rasch.roc.r100[,3],rasch.roc.r100[,2],col="red",lty=2,lwd=3)#
lines(1-rasch.roc.r200[,3],rasch.roc.r200[,2],col="red",lty=3,lwd=3)#
legend("bottomright",c(expression(RM~(paste(p,"*",sep="")==50)),expression(RM~(paste(p,"*",sep="")==100)),expression(RM~(paste(p,"*",sep="")==200)),expression(PCA~(paste(p,"*",sep="")==50)),expression(PCA~(paste(p,"*",sep="")==100)),expression(PCA~(paste(p,"*",sep="")==200))),lty=c(1,2,3,1,2,3),lwd=c(3,3,3,3,3,3),col=c("red","red","red","blue","blue","blue"),bty="n",cex=1.7)#
plot(1-pca.roc.p50[,3],pca.roc.p50[,2],type="l",cex=1.7,main="(b)",xlab="1-Specificity",ylab="Sensitivity",col="blue",lty=1,lwd=3)#
lines(1-pca.roc.p100[,3],pca.roc.p100[,2],col="blue",lty=2,lwd=3)#
lines(1-pca.roc.p200[,3],pca.roc.p200[,2],col="blue",lty=3,lwd=3)#
lines(1-rasch.roc.p50[,3],rasch.roc.p50[,2],col="red",lty=1,lwd=3)#
lines(1-rasch.roc.p100[,3],rasch.roc.p100[,2],col="red",lty=2,lwd=3)#
lines(1-rasch.roc.p200[,3],rasch.roc.p200[,2],col="red",lty=3,lwd=3)#
legend("bottomright",c(expression(RM~(paste(p,"*",sep="")==50)),expression(RM~(paste(p,"*",sep="")==100)),expression(RM~(paste(p,"*",sep="")==200)),expression(PCA~(paste(p,"*",sep="")==50)),expression(PCA~(paste(p,"*",sep="")==100)),expression(PCA~(paste(p,"*",sep="")==200))),lty=c(1,2,3,1,2,3),lwd=c(3,3,3,3,3,3),col=c("red","red","red","blue","blue","blue"),bty="n",cex=1.7)#
dev.off()
pdf("roc2_bw_temp.pdf",width=20,height=10,paper="special")#
par(cex.axis=1.7,cex.lab=1.7,cex.main=1.7,pty="s",mfrow=c(1,2))#
plot(1-pca.roc.r50[,3],pca.roc.r50[,2],type="l",cex=1.7,main="(a)",xlab="1-Specificity",ylab="Sensitivity",col="gray60",lty=1,lwd=3)#
lines(1-pca.roc.r100[,3],pca.roc.r100[,2],col="gray60",lty=2,lwd=3)#
lines(1-pca.roc.r200[,3],pca.roc.r200[,2],col="gray60",lty=3,lwd=3)#
lines(1-rasch.roc.r50[,3],rasch.roc.r50[,2],col="black",lty=1,lwd=3)#
lines(1-rasch.roc.r100[,3],rasch.roc.r100[,2],col="black",lty=2,lwd=3)#
lines(1-rasch.roc.r200[,3],rasch.roc.r200[,2],col="black",lty=3,lwd=3)#
legend("bottomright",c(expression(RM~(paste(p,"*",sep="")==50)),expression(RM~(paste(p,"*",sep="")==100)),expression(RM~(paste(p,"*",sep="")==200)),expression(PCA~(paste(p,"*",sep="")==50)),expression(PCA~(paste(p,"*",sep="")==100)),expression(PCA~(paste(p,"*",sep="")==200))),lty=c(1,2,3,1,2,3),lwd=c(3,3,3,3,3,3),col=c("black","black","black","gray60","gray60","gray60"),bty="n",cex=1.7)#
plot(1-pca.roc.p50[,3],pca.roc.p50[,2],type="l",cex=1.7,main="(b)",xlab="1-Specificity",ylab="Sensitivity",col="gray60",lty=1,lwd=3)#
lines(1-pca.roc.p100[,3],pca.roc.p100[,2],col="gray60",lty=2,lwd=3)#
lines(1-pca.roc.p200[,3],pca.roc.p200[,2],col="gray60",lty=3,lwd=3)#
lines(1-rasch.roc.p50[,3],rasch.roc.p50[,2],col="black",lty=1,lwd=3)#
lines(1-rasch.roc.p100[,3],rasch.roc.p100[,2],col="black",lty=2,lwd=3)#
lines(1-rasch.roc.p200[,3],rasch.roc.p200[,2],col="black",lty=3,lwd=3)#
legend("bottomright",c(expression(RM~(paste(p,"*",sep="")==50)),expression(RM~(paste(p,"*",sep="")==100)),expression(RM~(paste(p,"*",sep="")==200)),expression(PCA~(paste(p,"*",sep="")==50)),expression(PCA~(paste(p,"*",sep="")==100)),expression(PCA~(paste(p,"*",sep="")==200))),lty=c(1,2,3,1,2,3),lwd=c(3,3,3,3,3,3),col=c("black","black","black","gray60","gray60","gray60"),bty="n",cex=1.7)#
dev.off()#
#
pdf("roc2_cl_temp.pdf",width=20,height=10,paper="special")#
par(cex.axis=1.7,cex.lab=1.7,cex.main=1.7,pty="s",mfrow=c(1,2))#
plot(1-pca.roc.r50[,3],pca.roc.r50[,2],type="l",cex=1.7,main="(a)",xlab="1-Specificity",ylab="Sensitivity",col="blue",lty=1,lwd=3)#
lines(1-pca.roc.r100[,3],pca.roc.r100[,2],col="blue",lty=2,lwd=3)#
lines(1-pca.roc.r200[,3],pca.roc.r200[,2],col="blue",lty=3,lwd=3)#
lines(1-rasch.roc.r50[,3],rasch.roc.r50[,2],col="red",lty=1,lwd=3)#
lines(1-rasch.roc.r100[,3],rasch.roc.r100[,2],col="red",lty=2,lwd=3)#
lines(1-rasch.roc.r200[,3],rasch.roc.r200[,2],col="red",lty=3,lwd=3)#
legend("bottomright",c(expression(RM~(paste(p,"*",sep="")==50)),expression(RM~(paste(p,"*",sep="")==100)),expression(RM~(paste(p,"*",sep="")==200)),expression(PCA~(paste(p,"*",sep="")==50)),expression(PCA~(paste(p,"*",sep="")==100)),expression(PCA~(paste(p,"*",sep="")==200))),lty=c(1,2,3,1,2,3),lwd=c(3,3,3,3,3,3),col=c("red","red","red","blue","blue","blue"),bty="n",cex=1.7)#
plot(1-pca.roc.p50[,3],pca.roc.p50[,2],type="l",cex=1.7,main="(b)",xlab="1-Specificity",ylab="Sensitivity",col="blue",lty=1,lwd=3)#
lines(1-pca.roc.p100[,3],pca.roc.p100[,2],col="blue",lty=2,lwd=3)#
lines(1-pca.roc.p200[,3],pca.roc.p200[,2],col="blue",lty=3,lwd=3)#
lines(1-rasch.roc.p50[,3],rasch.roc.p50[,2],col="red",lty=1,lwd=3)#
lines(1-rasch.roc.p100[,3],rasch.roc.p100[,2],col="red",lty=2,lwd=3)#
lines(1-rasch.roc.p200[,3],rasch.roc.p200[,2],col="red",lty=3,lwd=3)#
legend("bottomright",c(expression(RM~(paste(p,"*",sep="")==50)),expression(RM~(paste(p,"*",sep="")==100)),expression(RM~(paste(p,"*",sep="")==200)),expression(PCA~(paste(p,"*",sep="")==50)),expression(PCA~(paste(p,"*",sep="")==100)),expression(PCA~(paste(p,"*",sep="")==200))),lty=c(1,2,3,1,2,3),lwd=c(3,3,3,3,3,3),col=c("red","red","red","blue","blue","blue"),bty="n",cex=1.7)#
dev.off()
pdf("roc2_bw_temp.pdf",width=20,height=10,paper="special")#
par(cex.axis=1.7,cex.lab=1.7,cex.main=1.7,pty="s",mfrow=c(1,2))#
plot(1-pca.roc.r50[,3],pca.roc.r50[,2],type="l",cex=1.7,main="(a)",xlab=bquote(1-"Specificity"),ylab="Sensitivity",col="gray60",lty=1,lwd=3)#
lines(1-pca.roc.r100[,3],pca.roc.r100[,2],col="gray60",lty=2,lwd=3)#
lines(1-pca.roc.r200[,3],pca.roc.r200[,2],col="gray60",lty=3,lwd=3)#
lines(1-rasch.roc.r50[,3],rasch.roc.r50[,2],col="black",lty=1,lwd=3)#
lines(1-rasch.roc.r100[,3],rasch.roc.r100[,2],col="black",lty=2,lwd=3)#
lines(1-rasch.roc.r200[,3],rasch.roc.r200[,2],col="black",lty=3,lwd=3)#
legend("bottomright",c(expression(RM~(paste(p,"*",sep="")==50)),expression(RM~(paste(p,"*",sep="")==100)),expression(RM~(paste(p,"*",sep="")==200)),expression(PCA~(paste(p,"*",sep="")==50)),expression(PCA~(paste(p,"*",sep="")==100)),expression(PCA~(paste(p,"*",sep="")==200))),lty=c(1,2,3,1,2,3),lwd=c(3,3,3,3,3,3),col=c("black","black","black","gray60","gray60","gray60"),bty="n",cex=1.7)#
plot(1-pca.roc.p50[,3],pca.roc.p50[,2],type="l",cex=1.7,main="(b)",xlab="1-Specificity",ylab="Sensitivity",col="gray60",lty=1,lwd=3)#
lines(1-pca.roc.p100[,3],pca.roc.p100[,2],col="gray60",lty=2,lwd=3)#
lines(1-pca.roc.p200[,3],pca.roc.p200[,2],col="gray60",lty=3,lwd=3)#
lines(1-rasch.roc.p50[,3],rasch.roc.p50[,2],col="black",lty=1,lwd=3)#
lines(1-rasch.roc.p100[,3],rasch.roc.p100[,2],col="black",lty=2,lwd=3)#
lines(1-rasch.roc.p200[,3],rasch.roc.p200[,2],col="black",lty=3,lwd=3)#
legend("bottomright",c(expression(RM~(paste(p,"*",sep="")==50)),expression(RM~(paste(p,"*",sep="")==100)),expression(RM~(paste(p,"*",sep="")==200)),expression(PCA~(paste(p,"*",sep="")==50)),expression(PCA~(paste(p,"*",sep="")==100)),expression(PCA~(paste(p,"*",sep="")==200))),lty=c(1,2,3,1,2,3),lwd=c(3,3,3,3,3,3),col=c("black","black","black","gray60","gray60","gray60"),bty="n",cex=1.7)#
dev.off()
pdf("roc2_bw_temp.pdf",width=20,height=10,paper="special")#
par(cex.axis=1.7,cex.lab=1.7,cex.main=1.7,pty="s",mfrow=c(1,2))#
plot(1-pca.roc.r50[,3],pca.roc.r50[,2],type="l",cex=1.7,main="(a)",xlab=bquote(1-"Specificity"),ylab="Sensitivity",col="gray60",lty=1,lwd=3)#
lines(1-pca.roc.r100[,3],pca.roc.r100[,2],col="gray60",lty=2,lwd=3)#
lines(1-pca.roc.r200[,3],pca.roc.r200[,2],col="gray60",lty=3,lwd=3)#
lines(1-rasch.roc.r50[,3],rasch.roc.r50[,2],col="black",lty=1,lwd=3)#
lines(1-rasch.roc.r100[,3],rasch.roc.r100[,2],col="black",lty=2,lwd=3)#
lines(1-rasch.roc.r200[,3],rasch.roc.r200[,2],col="black",lty=3,lwd=3)#
legend("bottomright",c(expression(RM~(paste(p,"*",sep="")==50)),expression(RM~(paste(p,"*",sep="")==100)),expression(RM~(paste(p,"*",sep="")==200)),expression(PCA~(paste(p,"*",sep="")==50)),expression(PCA~(paste(p,"*",sep="")==100)),expression(PCA~(paste(p,"*",sep="")==200))),lty=c(1,2,3,1,2,3),lwd=c(3,3,3,3,3,3),col=c("black","black","black","gray60","gray60","gray60"),bty="n",cex=1.7)#
plot(1-pca.roc.p50[,3],pca.roc.p50[,2],type="l",cex=1.7,main="(b)",xlab="1-Specificity",ylab="Sensitivity",col="gray60",lty=1,lwd=3)#
lines(1-pca.roc.p100[,3],pca.roc.p100[,2],col="gray60",lty=2,lwd=3)#
lines(1-pca.roc.p200[,3],pca.roc.p200[,2],col="gray60",lty=3,lwd=3)#
lines(1-rasch.roc.p50[,3],rasch.roc.p50[,2],col="black",lty=1,lwd=3)#
lines(1-rasch.roc.p100[,3],rasch.roc.p100[,2],col="black",lty=2,lwd=3)#
lines(1-rasch.roc.p200[,3],rasch.roc.p200[,2],col="black",lty=3,lwd=3)#
legend("bottomright",c(expression(RM~(paste(p,"*",sep="")==50)),expression(RM~(paste(p,"*",sep="")==100)),expression(RM~(paste(p,"*",sep="")==200)),expression(PCA~(paste(p,"*",sep="")==50)),expression(PCA~(paste(p,"*",sep="")==100)),expression(PCA~(paste(p,"*",sep="")==200))),lty=c(1,2,3,1,2,3),lwd=c(3,3,3,3,3,3),col=c("black","black","black","gray60","gray60","gray60"),bty="n",cex=1.7)#
dev.off()#
#
pdf("roc2_cl_temp.pdf",width=20,height=10,paper="special")#
par(cex.axis=1.7,cex.lab=1.7,cex.main=1.7,pty="s",mfrow=c(1,2))#
plot(1-pca.roc.r50[,3],pca.roc.r50[,2],type="l",cex=1.7,main="(a)",xlab=bquote(1-"Specificity"),ylab="Sensitivity",col="blue",lty=1,lwd=3)#
lines(1-pca.roc.r100[,3],pca.roc.r100[,2],col="blue",lty=2,lwd=3)#
lines(1-pca.roc.r200[,3],pca.roc.r200[,2],col="blue",lty=3,lwd=3)#
lines(1-rasch.roc.r50[,3],rasch.roc.r50[,2],col="red",lty=1,lwd=3)#
lines(1-rasch.roc.r100[,3],rasch.roc.r100[,2],col="red",lty=2,lwd=3)#
lines(1-rasch.roc.r200[,3],rasch.roc.r200[,2],col="red",lty=3,lwd=3)#
legend("bottomright",c(expression(RM~(paste(p,"*",sep="")==50)),expression(RM~(paste(p,"*",sep="")==100)),expression(RM~(paste(p,"*",sep="")==200)),expression(PCA~(paste(p,"*",sep="")==50)),expression(PCA~(paste(p,"*",sep="")==100)),expression(PCA~(paste(p,"*",sep="")==200))),lty=c(1,2,3,1,2,3),lwd=c(3,3,3,3,3,3),col=c("red","red","red","blue","blue","blue"),bty="n",cex=1.7)#
plot(1-pca.roc.p50[,3],pca.roc.p50[,2],type="l",cex=1.7,main="(b)",xlab="1-Specificity",ylab="Sensitivity",col="blue",lty=1,lwd=3)#
lines(1-pca.roc.p100[,3],pca.roc.p100[,2],col="blue",lty=2,lwd=3)#
lines(1-pca.roc.p200[,3],pca.roc.p200[,2],col="blue",lty=3,lwd=3)#
lines(1-rasch.roc.p50[,3],rasch.roc.p50[,2],col="red",lty=1,lwd=3)#
lines(1-rasch.roc.p100[,3],rasch.roc.p100[,2],col="red",lty=2,lwd=3)#
lines(1-rasch.roc.p200[,3],rasch.roc.p200[,2],col="red",lty=3,lwd=3)#
legend("bottomright",c(expression(RM~(paste(p,"*",sep="")==50)),expression(RM~(paste(p,"*",sep="")==100)),expression(RM~(paste(p,"*",sep="")==200)),expression(PCA~(paste(p,"*",sep="")==50)),expression(PCA~(paste(p,"*",sep="")==100)),expression(PCA~(paste(p,"*",sep="")==200))),lty=c(1,2,3,1,2,3),lwd=c(3,3,3,3,3,3),col=c("red","red","red","blue","blue","blue"),bty="n",cex=1.7)#
dev.off()
list(ls)
?rm
rm(list=ls())
ls()
getwd()
library(MASS)#
library(ltm)#
library(biclust)#
source("./r/pca.lda.r")#
source("./r/rasch.lda.r")#
source("./r/my.binarize.r")#
#
load("./data/leukemia.rda")#
#
my.train <- as.matrix(read.table("./subsets/my.train.txt",h=F))#
#
r50 <- as.matrix(read.table("./subsets/r50.txt",h=F))#
pca.components.r50 <- read.table("./components/leukemia_pca_components_r50.txt",h=F)$V1#
rasch.components.r50 <- read.table("./components/leukemia_rasch_components_r50.txt",h=F)$V1#
r100 <- as.matrix(read.table("./subsets/r100.txt",h=F))#
pca.components.r100 <- read.table("./components/leukemia_pca_components_r100.txt",h=F)$V1#
rasch.components.r100 <- read.table("./components/leukemia_rasch_components_r100.txt",h=F)$V1#
r200 <- as.matrix(read.table("./subsets/r200.txt",h=F))#
pca.components.r200 <- read.table("./components/leukemia_pca_components_r200.txt",h=F)$V1#
rasch.components.r200 <- read.table("./components/leukemia_rasch_components_r200.txt",h=F)$V1#
#
pca.posterior.r50 <- NULL#
pca.class.r50 <- NULL#
for(i in 1:nrow(my.train)) {#
    predict <- pca.lda(Xtrain=lx.original[my.train[i,],r50[i,]],Ytrain=ly.original[my.train[i,]],Xtest=lx.original[-(my.train)[i,],r50[i,]],ncomp=pca.components.r50[i])$posterior[,1]#
    class <- ly.original[-(my.train)[i,]]#
    pca.posterior.r50 <- c(pca.posterior.r50,predict)#
    pca.class.r50 <- c(pca.class.r50,class)#
}#
pca.roc.r50 <- roc(as.vector(pca.posterior.r50),as.factor(pca.class.r50))#
#
pca.posterior.r100 <- NULL#
pca.class.r100 <- NULL#
for(i in 1:nrow(my.train)) {#
    predict <- pca.lda(Xtrain=lx.original[my.train[i,],r100[i,]],Ytrain=ly.original[my.train[i,]],Xtest=lx.original[-(my.train)[i,],r100[i,]],ncomp=pca.components.r100[i])$posterior[,1]#
    class <- ly.original[-(my.train)[i,]]#
    pca.posterior.r100 <- c(pca.posterior.r100,predict)#
    pca.class.r100 <- c(pca.class.r100,class)#
}#
pca.roc.r100 <- roc(as.vector(pca.posterior.r100),as.factor(pca.class.r100))#
#
pca.posterior.r200 <- NULL#
pca.class.r200 <- NULL#
for(i in 1:nrow(my.train)) {#
    predict <- pca.lda(Xtrain=lx.original[my.train[i,],r200[i,]],Ytrain=ly.original[my.train[i,]],Xtest=lx.original[-(my.train)[i,],r200[i,]],ncomp=pca.components.r200[i])$posterior[,1]#
    class <- ly.original[-(my.train)[i,]]#
    pca.posterior.r200 <- c(pca.posterior.r200,predict)#
    pca.class.r200 <- c(pca.class.r200,class)   #
}#
pca.roc.r200 <- roc(as.vector(pca.posterior.r200),as.factor(pca.class.r200))
rasch.posterior.r50 <- NULL#
rasch.class.r50 <- NULL#
for(i in 1:nrow(my.train)) {#
    predict <- rasch.lda.predict(Xtrain=lx.original[my.train[i,],r50[i,]],Ytrain=ly.original[my.train[i,]],Xtest=lx.original[-(my.train)[i,],r50[i,]],ncomp=rasch.components.r50[i])$posterior[,1]#
    class <- ly.original[-(my.train)[i,]]#
    rasch.posterior.r50 <- c(rasch.posterior.r50,predict)#
    rasch.class.r50 <- c(rasch.class.r50,class)#
    cat("Complete: ",i,"\n")#
}#
rasch.roc.r50 <- roc(as.vector(rasch.posterior.r50),as.factor(rasch.class.r50))#
#
#
rasch.posterior.r100 <- NULL#
rasch.class.r100 <- NULL#
for(i in 1:nrow(my.train)) {#
    predict <- rasch.lda.predict(Xtrain=lx.original[my.train[i,],r100[i,]],Ytrain=ly.original[my.train[i,]],Xtest=lx.original[-(my.train)[i,],r100[i,]],ncomp=rasch.components.r100[i])$posterior[,1]#
    class <- ly.original[-(my.train)[i,]]#
    rasch.posterior.r100 <- c(rasch.posterior.r100,predict)#
    rasch.class.r100 <- c(rasch.class.r100,class)#
    cat("Complete: ",i,"\n")#
}#
rasch.roc.r100 <- roc(as.vector(rasch.posterior.r100),as.factor(rasch.class.r100))#
#
rasch.posterior.r200 <- NULL#
rasch.class.r200 <- NULL#
for(i in 1:nrow(my.train)) {#
    predict <- rasch.lda.predict(Xtrain=lx.original[my.train[i,],r200[i,]],Ytrain=ly.original[my.train[i,]],Xtest=lx.original[-(my.train)[i,],r200[i,]],ncomp=rasch.components.r200[i])$posterior[,1]#
    class <- ly.original[-(my.train)[i,]]#
    rasch.posterior.r200 <- c(rasch.posterior.r200,predict)#
    rasch.class.r200 <- c(rasch.class.r200,class)#
    cat("Complete: ",i,"\n")#
}#
rasch.roc.r200 <- roc(as.vector(rasch.posterior.r200),as.factor(rasch.class.r200))
p50 <- as.matrix(read.table("./subsets/p50.txt",h=F))#
pca.components.p50 <- read.table("./components/leukemia_pca_components_p50.txt",h=F)$V1#
rasch.components.p50 <- read.table("./components/leukemia_rasch_components_p50.txt",h=F)$V1#
p100 <- as.matrix(read.table("./subsets/p100.txt",h=F))#
pca.components.p100 <- read.table("./components/leukemia_pca_components_p100.txt",h=F)$V1#
rasch.components.p100 <- read.table("./components/leukemia_rasch_components_p100.txt",h=F)$V1#
p200 <- as.matrix(read.table("./subsets/p200.txt",h=F))#
pca.components.p200 <- read.table("./components/leukemia_pca_components_p200.txt",h=F)$V1#
rasch.components.p200 <- read.table("./components/leukemia_rasch_components_p200.txt",h=F)$V1#
#
#
#
pca.posterior.p50 <- NULL#
pca.class.p50 <- NULL#
for(i in 1:nrow(my.train)) {#
    predict <- pca.lda(Xtrain=lx.original[my.train[i,],p50[i,]],Ytrain=ly.original[my.train[i,]],Xtest=lx.original[-(my.train)[i,],p50[i,]],ncomp=pca.components.p50[i])$posterior[,1]#
    class <- ly.original[-(my.train)[i,]]#
    pca.posterior.p50 <- c(pca.posterior.p50,predict)#
    pca.class.p50 <- c(pca.class.p50,class)#
}#
pca.roc.p50 <- roc(as.vector(pca.posterior.p50),as.factor(rasch.class.p50))#
#
pca.posterior.p100 <- NULL#
pca.class.p100 <- NULL#
for(i in 1:nrow(my.train)) {#
    predict <- pca.lda(Xtrain=lx.original[my.train[i,],p100[i,]],Ytrain=ly.original[my.train[i,]],Xtest=lx.original[-(my.train)[i,],p100[i,]],ncomp=pca.components.p100[i])$posterior[,1]#
    class <- ly.original[-(my.train)[i,]]#
    pca.posterior.p100 <- c(pca.posterior.p100,predict)#
    pca.class.p100 <- c(pca.class.p100,class)#
}#
pca.roc.p100 <- roc(as.vector(pca.posterior.p100),as.factor(rasch.class.p100))#
#
pca.posterior.p200 <- NULL#
pca.class.p200 <- NULL#
for(i in 1:nrow(my.train)) {#
    predict <- pca.lda(Xtrain=lx.original[my.train[i,],p200[i,]],Ytrain=ly.original[my.train[i,]],Xtest=lx.original[-(my.train)[i,],p200[i,]],ncomp=pca.components.p200[i])$posterior[,1]#
    class <- ly.original[-(my.train)[i,]]#
    pca.posterior.p200 <- c(pca.posterior.p200,predict)#
    pca.class.p200 <- c(pca.class.p200,class)#
}#
pca.roc.p200 <- roc(as.vector(pca.posterior.p200),as.factor(rasch.class.p200))
pca.posterior.p50 <- NULL#
pca.class.p50 <- NULL#
for(i in 1:nrow(my.train)) {#
    predict <- pca.lda(Xtrain=lx.original[my.train[i,],p50[i,]],Ytrain=ly.original[my.train[i,]],Xtest=lx.original[-(my.train)[i,],p50[i,]],ncomp=pca.components.p50[i])$posterior[,1]#
    class <- ly.original[-(my.train)[i,]]#
    pca.posterior.p50 <- c(pca.posterior.p50,predict)#
    pca.class.p50 <- c(pca.class.p50,class)#
}#
pca.roc.p50 <- roc(as.vector(pca.posterior.p50),as.factor(pca.class.p50))#
#
pca.posterior.p100 <- NULL#
pca.class.p100 <- NULL#
for(i in 1:nrow(my.train)) {#
    predict <- pca.lda(Xtrain=lx.original[my.train[i,],p100[i,]],Ytrain=ly.original[my.train[i,]],Xtest=lx.original[-(my.train)[i,],p100[i,]],ncomp=pca.components.p100[i])$posterior[,1]#
    class <- ly.original[-(my.train)[i,]]#
    pca.posterior.p100 <- c(pca.posterior.p100,predict)#
    pca.class.p100 <- c(pca.class.p100,class)#
}#
pca.roc.p100 <- roc(as.vector(pca.posterior.p100),as.factor(pca.class.p100))#
#
pca.posterior.p200 <- NULL#
pca.class.p200 <- NULL#
for(i in 1:nrow(my.train)) {#
    predict <- pca.lda(Xtrain=lx.original[my.train[i,],p200[i,]],Ytrain=ly.original[my.train[i,]],Xtest=lx.original[-(my.train)[i,],p200[i,]],ncomp=pca.components.p200[i])$posterior[,1]#
    class <- ly.original[-(my.train)[i,]]#
    pca.posterior.p200 <- c(pca.posterior.p200,predict)#
    pca.class.p200 <- c(pca.class.p200,class)#
}#
pca.roc.p200 <- roc(as.vector(pca.posterior.p200),as.factor(pca.class.p200))
rasch.posterior.p50 <- NULL#
rasch.class.p50 <- NULL#
for(i in 1:nrow(my.train)) {#
    predict <- rasch.lda.predict(Xtrain=lx.original[my.train[i,],p50[i,]],Ytrain=ly.original[my.train[i,]],Xtest=lx.original[-(my.train)[i,],p50[i,]],ncomp=rasch.components.p50[i])$posterior[,1]#
    class <- ly.original[-(my.train)[i,]]#
    rasch.posterior.p50 <- c(rasch.posterior.p50,predict)#
    rasch.class.p50 <- c(rasch.class.p50,class)#
    cat("Complete: ",i,"\n")#
}#
rasch.roc.p50 <- roc(as.vector(rasch.posterior.p50),as.factor(rasch.class.p50))
rasch.posterior.p100 <- NULL#
rasch.class.p100 <- NULL#
for(i in 1:nrow(my.train)) {#
    predict <- rasch.lda.predict(Xtrain=lx.original[my.train[i,],p100[i,]],Ytrain=ly.original[my.train[i,]],Xtest=lx.original[-(my.train)[i,],p100[i,]],ncomp=rasch.components.p100[i])$posterior[,1]#
    class <- ly.original[-(my.train)[i,]]#
    rasch.posterior.p100 <- c(rasch.posterior.p100,predict)#
    rasch.class.p100 <- c(rasch.class.p100,class)#
    cat("Complete: ",i,"\n")#
}#
rasch.roc.p100 <- roc(as.vector(rasch.posterior.p100),as.factor(rasch.class.p100))
rasch.posterior.p200 <- NULL#
rasch.class.p200 <- NULL#
for(i in 1:nrow(my.train)) {#
    predict <- rasch.lda.predict(Xtrain=lx.original[my.train[i,],p200[i,]],Ytrain=ly.original[my.train[i,]],Xtest=lx.original[-(my.train)[i,],p200[i,]],ncomp=rasch.components.p200[i])$posterior[,1]#
    class <- ly.original[-(my.train)[i,]]#
    rasch.posterior.p200 <- c(rasch.posterior.p200,predict)#
    rasch.class.p200 <- c(rasch.class.p200,class)#
    cat("Complete: ",i,"\n")#
}#
rasch.roc.p200 <- roc(as.vector(rasch.posterior.p200),as.factor(rasch.class.p200))
pdf("roc1_bw_temp.pdf",width=20,height=10,paper="special")#
par(cex.axis=1.7,cex.lab=1.7,cex.main=1.7,pty="s",mfrow=c(1,2))#
plot(1-pca.roc.r50[,3],pca.roc.r50[,2],type="l",cex=1.7,main="(a)",xlab=bquote(1-"Specificity"),ylab="Sensitivity",col="gray60",lty=1,lwd=3)#
lines(1-pca.roc.r100[,3],pca.roc.r100[,2],col="gray60",lty=2,lwd=3)#
lines(1-pca.roc.r200[,3],pca.roc.r200[,2],col="gray60",lty=3,lwd=3)#
lines(1-rasch.roc.r50[,3],rasch.roc.r50[,2],col="black",lty=1,lwd=3)#
lines(1-rasch.roc.r100[,3],rasch.roc.r100[,2],col="black",lty=2,lwd=3)#
lines(1-rasch.roc.r200[,3],rasch.roc.r200[,2],col="black",lty=3,lwd=3)#
legend("bottomright",c(expression(RM~(paste(p,"*",sep="")==50)),expression(RM~(paste(p,"*",sep="")==100)),expression(RM~(paste(p,"*",sep="")==200)),expression(PCA~(paste(p,"*",sep="")==50)),expression(PCA~(paste(p,"*",sep="")==100)),expression(PCA~(paste(p,"*",sep="")==200))),lty=c(1,2,3,1,2,3),lwd=c(3,3,3,3,3,3),col=c("black","black","black","gray60","gray60","gray60"),bty="n",cex=1.7)#
plot(1-pca.roc.p50[,3],pca.roc.p50[,2],type="l",cex=1.7,main="(b)",xlab="1-Specificity",ylab="Sensitivity",col="gray60",lty=1,lwd=3)#
lines(1-pca.roc.p100[,3],pca.roc.p100[,2],col="gray60",lty=2,lwd=3)#
lines(1-pca.roc.p200[,3],pca.roc.p200[,2],col="gray60",lty=3,lwd=3)#
lines(1-rasch.roc.p50[,3],rasch.roc.p50[,2],col="black",lty=1,lwd=3)#
lines(1-rasch.roc.p100[,3],rasch.roc.p100[,2],col="black",lty=2,lwd=3)#
lines(1-rasch.roc.p200[,3],rasch.roc.p200[,2],col="black",lty=3,lwd=3)#
legend("bottomright",c(expression(RM~(paste(p,"*",sep="")==50)),expression(RM~(paste(p,"*",sep="")==100)),expression(RM~(paste(p,"*",sep="")==200)),expression(PCA~(paste(p,"*",sep="")==50)),expression(PCA~(paste(p,"*",sep="")==100)),expression(PCA~(paste(p,"*",sep="")==200))),lty=c(1,2,3,1,2,3),lwd=c(3,3,3,3,3,3),col=c("black","black","black","gray60","gray60","gray60"),bty="n",cex=1.7)#
dev.off()#
#
pdf("roc1_cl_temp.pdf",width=20,height=10,paper="special")#
par(cex.axis=1.7,cex.lab=1.7,cex.main=1.7,pty="s",mfrow=c(1,2))#
plot(1-pca.roc.r50[,3],pca.roc.r50[,2],type="l",cex=1.7,main="(a)",xlab=bquote(1-"Specificity"),ylab="Sensitivity",col="blue",lty=1,lwd=3)#
lines(1-pca.roc.r100[,3],pca.roc.r100[,2],col="blue",lty=2,lwd=3)#
lines(1-pca.roc.r200[,3],pca.roc.r200[,2],col="blue",lty=3,lwd=3)#
lines(1-rasch.roc.r50[,3],rasch.roc.r50[,2],col="red",lty=1,lwd=3)#
lines(1-rasch.roc.r100[,3],rasch.roc.r100[,2],col="red",lty=2,lwd=3)#
lines(1-rasch.roc.r200[,3],rasch.roc.r200[,2],col="red",lty=3,lwd=3)#
legend("bottomright",c(expression(RM~(paste(p,"*",sep="")==50)),expression(RM~(paste(p,"*",sep="")==100)),expression(RM~(paste(p,"*",sep="")==200)),expression(PCA~(paste(p,"*",sep="")==50)),expression(PCA~(paste(p,"*",sep="")==100)),expression(PCA~(paste(p,"*",sep="")==200))),lty=c(1,2,3,1,2,3),lwd=c(3,3,3,3,3,3),col=c("red","red","red","blue","blue","blue"),bty="n",cex=1.7)#
plot(1-pca.roc.p50[,3],pca.roc.p50[,2],type="l",cex=1.7,main="(b)",xlab="1-Specificity",ylab="Sensitivity",col="blue",lty=1,lwd=3)#
lines(1-pca.roc.p100[,3],pca.roc.p100[,2],col="blue",lty=2,lwd=3)#
lines(1-pca.roc.p200[,3],pca.roc.p200[,2],col="blue",lty=3,lwd=3)#
lines(1-rasch.roc.p50[,3],rasch.roc.p50[,2],col="red",lty=1,lwd=3)#
lines(1-rasch.roc.p100[,3],rasch.roc.p100[,2],col="red",lty=2,lwd=3)#
lines(1-rasch.roc.p200[,3],rasch.roc.p200[,2],col="red",lty=3,lwd=3)#
legend("bottomright",c(expression(RM~(paste(p,"*",sep="")==50)),expression(RM~(paste(p,"*",sep="")==100)),expression(RM~(paste(p,"*",sep="")==200)),expression(PCA~(paste(p,"*",sep="")==50)),expression(PCA~(paste(p,"*",sep="")==100)),expression(PCA~(paste(p,"*",sep="")==200))),lty=c(1,2,3,1,2,3),lwd=c(3,3,3,3,3,3),col=c("red","red","red","blue","blue","blue"),bty="n",cex=1.7)#
dev.off()
aucRoc(pca.roc.r50)
aucRoc(pca.roc.r100)
aucRoc(pca.roc.r200)
aucRoc(rasch.roc.r50)
aucRoc(rasch.roc.r100)
aucRoc(rasch.roc.r200)
aucRoc(pca.roc.p50)
aucRoc(pca.roc.p100)
aucRoc(pca.roc.p200)
aucRoc(rasch.roc.p50)
aucRoc(rasch.roc.p100)
aucRoc(rasch.roc.p200)
