On 2/7/09, Juliet Hannah <juliet.han...@gmail.com> wrote: > Hi Group, > > Here is some data. > > p <- runif(1000) # sample data > groups <- rep(c(1,2),each=500) #conditioning variable > mydata <- cbind(p,groups) > n <- length(p) > u <- (1:n)/(n + 1) # uniform distribution reference for qqplot > logp <- -log(p,base=10) > logu <- -log(u,base=10) > qqplot(logp,logu) > > How can I make the above qqplot in lattice and/or ggplot2. The sample > is uniform, and I take -log10. In addition, I would like to use > the group variable to condition. I was not able to get the stat_qq > function in ggplot2 to work.
With lattice: library(lattice) qqmath(~ -log(p), data = as.data.frame(mydata), distribution = qexp) qqmath(~ -log(p) | factor(groups), data = as.data.frame(mydata), distribution = qexp) -Deepayan ______________________________________________ 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.
