[R] Adding lines to scatterplot odd result when creating multiple plots

[markm0705]

Dear R helpers

I'm attempting to create a matrix of scatterplots with X-Y mean lines and
regression lines added to each plot in the matrix.  I have managed to create
the first plot I would like using scatterplot but have run into an odd
result when I use par() to set up the page to take multiple plot. 
Specifically, the mean and regression lines appear to plot in the second
plot window, not on top of the scatterplot (when I include the par(mfcol)
command).  

How do I get the scatterplot and added lines to stay together on this plot
and subsequent plots I would like to include in the scatterplot matrix

Data attached

Thanks in advance MarkM

library("car")

par(mfcol=c(1,2))

# Parameters to change

Infile<-"kt3d_Thk1.dat"
X<-"Estimated thickness (mE)"
Y<-"True thickness (mE)"

#load data
#read the data skip then read header

crossval <- read.table(Infile,skip = 9,sep = "")
head<-readLines(Infile,9)

#decode the header lines wanted

head2<-head[3:9]
colnames(crossval)=head2

#Filter out non-estimated

crossval2<-crossval[crossval$Estimate>0,]

# Compute the means

AveEst<-mean(crossval2$Estimate)
AveTru<-mean(crossval2$True)

#Fit the regression line

Fit<-lm(crossval2$True~crossval2$Estimate )

#create plots

scatterplot(True ~ Estimate,
	data=crossval2,
	xlab= X,
	ylab= Y,
	main= "Min 4 Max 8",
	grid=FALSE,
	xlim=c(0,8),
	ylim=c(0,8),
	pch=21,
	cex=1.2,
	smooth=FALSE,
	reg.line=FALSE
)


#Plot mean lines and regression

abline(h=AveTru,col="grey70")
abline(v=AveEst,col="grey70")
abline(0,1,col="red")
abline(lm(crossval2$True~crossval2$Estimate), col="blue", lty=1, lwd=2)

--
View this message in context: http://r.789695.n4.nabble.com/Adding-lines-to-scatterplot-odd-result-when-creating-multiple-plots-tp4014140p4014140.html
Sent from the R help mailing list archive at Nabble.com.