[R] Same scale on different q-q plots

Dr. Enrique Hernández Lemus ehernandez at inmegen.gob.mx
Mon Jun 4 20:28:36 CEST 2007


Hi there, I'm new on the list. I am wondering as how to make different q-q
plots and display them on the same scale in order to attain an easy visual
comparison.

I have tried to use xlim and ylim as arguments to qq.plot in a similar
fashion to what is done on plot, also I have tried to use the plot.window
utility to implement xlim and ylim. I have even modified the qq.plot.R
source and re-compiled but none of these efforts has given any result.

I am doing a sequence of steps like the following:

library(car)
 p<-read.csv("inputP",header=T)
 p <- apply(p,c(1,2),as.numeric)
 q<-qchisq(1-p,2)

 png(file="outputP.png", bg="transparent")
plot.window(xlim=c(0,50), ylim=c(0,50))
qq.plot(q, dist="chisq", df=2)
 y<-rchisq(length(p),2)
 lines(y,y)
 dev.off()


Is there anything that I can do to produce several q-q plots from
different data and put them in the same scale, say x-values in (0-50) and
y-values in (0-50)?

Thank you to all in advance

-- 
Dr. Enrique Hernández Lemus
Medical Sciences Researcher "B"
Computational Genomics
Phone: 5350-1970
Fax: 5350-1999
Email: ehernandez at inmegen.gob.mx



More information about the R-help mailing list