Figures in pdf.

figure1.R

@@ -26,8 +26,10 @@ plotit <- function(m=c(.3, .4, .525, .80), density=dBetam, shape, ylab="", xlab=
 
 }
 
+pdf(file = "pdf/figure1.pdf", width = 6, height = 6)
+
 op <- par(no.readonly = TRUE)
-par(mfcol=c(3,2), mar=c(5, 5.5, 2, 1), cex=.8)
+par(mfcol=c(3,2), mar=c(5, 5.5, 2, 1), cex=.6)
 
 plotit(shape=1, density=dBetam, xlab="Statistical power", ylab="density", main="S=1")
 plotit(shape=1/2, density=dBetam, xlab="Statistical power", ylab="density", main="S=1/2")
@@ -36,5 +38,7 @@ plotit(shape=1, density=dbimodal1090, xlab="Statistical power", ylab="density",
 plotit(shape=2, density=dbimodal1090, xlab="Statistical power", ylab="density", main="S=2 (10/90)")
 plotit(shape=1, density=dbimodal0595, xlab="Statistical power", ylab="density", main="S=1 (05/95)")
 
-dev.copy(device = jpeg, filename = 'jpeg/figure1.jpg', quality=100, width = 700, height = 700)
+#dev.copy(device = jpeg, filename = 'jpeg/figure1.jpg', quality=100, width = 700, height = 700)
+par(op)
 dev.off()
+

figure2.R

@@ -40,6 +40,8 @@ plotit <- function(d, est, ylim=c(0,1), log="x", legpos="topleft", xlab="", ylab
 }
 
 d=data$range
+
+pdf(file = "pdf/figure2.pdf", width = 6, height = 6)
 op <- par(no.readonly = TRUE)
 par(mfrow=c(2,2), mar=c(5, 5.5, 2, 1), cex=.6)
 
@@ -65,5 +67,5 @@ plotit(d, "bias",  legpos="bottomleft", ylim=c(1,250), log="xy",
 
 par(op)
 
-dev.copy(device = jpeg, filename = 'jpeg/figure2.jpg', quality=100, res=300, unit="mm", width = 150, height = 150)
+#dev.copy(device = jpeg, filename = 'jpeg/figure2.jpg', quality=100, res=300, unit="mm", width = 150, height = 150)
 dev.off()

figure3.R

@@ -54,6 +54,7 @@ CIs = c(.95, .75, .50)
 Rs=sort(c(35/97, fv(CIs)))
 d=data$rep
 
+pdf(file = "pdf/figure3.pdf", width = 6, height = 6)
 op <- par(no.readonly = TRUE)
 par(mfrow=c(2,2), mar=c(5, 5.5, 2, 1), cex=.6)
 
@@ -79,5 +80,5 @@ plotit(d, "bias", CIs, Rs, legpos="bottomleft", ylim=c(1,300), log="xy",
 
 par(op)
 
-dev.copy(device = jpeg, filename = 'jpeg/figure3.jpg', quality=100, res=300, unit="mm", width = 150, height = 150)
+#dev.copy(device = jpeg, filename = 'jpeg/figure3.jpg', quality=100, res=300, unit="mm", width = 150, height = 150)
 dev.off()