Multiple Plots eda10.shape

eda10.shape <- function(x,y,t,s)
{
	par(oma = c(10,0,5,0))
	par(fig = c(0.1, 0.4, 0.25, 0.75), mar = c(2,0.5,0.5,0.5), mgp = c(0,0.5,0))
#	par(fig = c(0.1, 0.4, 0.25, 0.75))
	iqd <- summary(x)[5] - summary(x)[2]
	plot(density(x, width = 2 * iqd), xlab = "", type = "l", ylim=c(0, 1.4), pch=15,  
+                        lty = 1, ylab = "", main = "")
	iqd <- summary(y)[5] - summary(y)[2]
	points(density(y, width = 2 * iqd), type = "l", xlab = "", ylab = "", pch=17,  lty = 2)
	iqd <- summary(t)[5] - summary(t)[2]
	points(density(t, width = 2 * iqd), type = "l", pch=18, xlab = "", ylab = "", lty = 3)
	iqd <- summary(s)[5] - summary(s)[2]
	points(density(s, width = 2 * iqd), type = "l", pch=16, xlab = "", ylab = "", lty = 4)
	legend(locator(1), legend=list("Model Studies", "Patients:", "  Examiner 1",
+                       "  Examiner 2", "Inter-examiner"), lty=c(1,-1,2,3,4), pch=c(16,-1,17,18,15), 
+	    cex=0.70)
                      mtext("Density Plot of Angular Differences - Two Trials", side = 1, outer = T,
+                        cex = 0.85, at = 0.2)
       
	mtext("Precision Experiments for 10 Patients and 24 Model Studies",
+                       side = 3, outer = T, cex = 1.5)
	par(fig=c(0.5,1, 0.75, 1),mar = c(2,0.5,0.5,0.5),mgp=c(0,0.5,0), new=TRUE)
 	qqnorm(x, pch = 15, xlab ="", ylab="", cex = 0.75, main="",ylim=c(0, 2.0))
 	qqline(x)
 	invisible()
	par(fig=c(0.5,1, 0.5, 0.75),new=TRUE)
 	qqnorm(y, pch = 17, xlab ="", ylab="", cex = 0.75,main="",ylim=c(0, 2.0))
 	qqline(y)
 	invisible()
	par(fig=c(0.5,1, 0.25, 0.5),new=TRUE)
 	qqnorm(t, pch = 18, xlab ="", ylab="", cex = 0.75, main="",ylim=c(0, 2.0))
 	qqline(t)
	invisible()
	par(fig=c(0.5,1, 0, 0.25),new=TRUE)
 	qqnorm(s, pch = 16, xlab ="", ylab="", cex = 0.75, main="",ylim=c(0, 2.0))
 	qqline(s)
	invisible()
                     mtext("Variation From Normality", side = 1, outer = T, cex = 0.85, at = 0.6)
}

Transcript

Now, one can do multiple plots - trying to summarize this data.  And you can see a long R function here.  Which we just wrote from different things and it goes on and on -