I would like to have a figure with two graphs. This is easily accomplished
using mfcol:
oldpar <- par(mfcol=c(1,2))
plot(x,y)
plot(z,x)
par(oldpar)
I run into trouble if I try to use identify with the two plots. If, after
identifying points on my first graph I hit the ESC key, or hitting stop menu
bar of my R session, the system stops the identification process, but fails to
give me my second graph. Is there a way to allow for the identification of
points when one is plotting to graphs in a single graph window? My code follows.
plotter <- function(first,second) {
# Allow for two plots in on graph window.
oldpar<-par(mfcol=c(1,2))
#Bland-Altman plot.
plot((second+first)/2,second-first)
abline(0,0)
# Allow for indentification of extreme values.
BAzap<-identify((second+first)/2,second-first,labels = seq_along(data$Line))
print(BAzap)
# Plot second as a function of first value.
plot(first,second,main="Limin vs. Limin",xlab="First (cm^2)",ylab="Second
(cm^3)")
# Add identity line.
abline(0,1,lty=2,col="red")
# Allow for identification of extreme values.
zap<-identify(first,second,labels = seq_along(data$Line))
print(zap)
# Add regression line.
fit1<-lm(first~second)
print(summary(fit1))
abline(fit1)
print(summary(fit1)$sigma)
# reset par to default values.
par(oldpar)
}
plotter(first,second)
Thanks,
John
John David Sorkin M.D., Ph.D.
Chief, Biostatistics and Informatics
University of Maryland School of Medicine Division of Gerontology
Baltimore VA Medical Center
10 North Greene Street
GRECC (BT/18/GR)
Baltimore, MD 21201-1524
(Phone) 410-605-7119
(Fax) 410-605-7913 (Please call phone number above prior to faxing)
Confidentiality Statement:
This email message, including any attachments, is for th...{{dropped:6}}
______________________________________________
R-help@r-project.org mailing list
https://stat.ethz.ch/mailman/listinfo/r-help
PLEASE do read the posting guide http://www.R-project.org/posting-guide.html
and provide commented, minimal, self-contained, reproducible code.
