Paul, This is great! Like you said it is really close. I made a few changes and for some reason the y-axis label magically came back. I tried to remove it but it wouldn't stay away.
Also, for some reason the title exceeds the margins of the layout. I am going ot mess around with this a bit more, but any suggestions for fixing this are also greatly appreciated. Thank you again. y<-rnorm(1000) yrange <- range(y) #histdata <- hist(y, breaks=breaks, plot=FALSE) histdata <- hist(y, plot=FALSE) densitydata <- density(y) par(oma=rep(3, 4), mar=rep(0, 4)) layout(matrix(1:3, ncol=3), widths=c(.2, .1, 1), respect=TRUE) barplot(histdata$density, horiz=TRUE, space=0, xlim=c(0, max(histdata$density)*1.1), axes=FALSE) par(new=TRUE) plot(densitydata$y, densitydata$x, xlim=c(0, max(histdata$density)*1.1), ylim=yrange, type="l", ann=FALSE, axes=FALSE) axis(2) box() boxplot(y, ylim=yrange, axes=FALSE) box() # qqnorm(y, axes=FALSE, ylim=yrange, xlim=yrange, main="") qqnorm(y, , main="") axis(2, labels = FALSE) #abline(0, 1) qqline(y, col="red", lwd=3) box() kurtosis_val <-(sum((y - mean(y))^4 ))/(var(y)*var(y))/length(y)-3 skewness<-function(x) { m_skew=mean(x) me_skew=median(x) s_skew=sqrt(var(x)) sk_skew=(m_skew-me_skew)/s_skew return(sk_skew) } title_text<-paste("Title Text", "\n Mean = ", format(mean(y), digits=4, scientific=F), " Standard Deviation = ", format(sd(y), digits=4, scientific=F), "\n Skewness = ", format(skewness(y), digits=4, scientific=F), "Kurtosis =", format(kurtosis_val, digits=4, scientific=F), sep="") mtext(title_text, NORTH<-3, line=0, adj=0.5, cex=1.2, col="red", outer=TRUE) --- On Tue, 3/17/09, Paul Murrell <p.murr...@auckland.ac.nz> wrote: > From: Paul Murrell <p.murr...@auckland.ac.nz> > Subject: Re: [R] R package to automatically produce combination plot? > To: jasonkrup...@yahoo.com > Cc: R-help@r-project.org > Date: Tuesday, March 17, 2009, 2:51 PM > Hi > > > Jason Rupert wrote: > > I guess no reply means there is not an existing > package to produce > > the plot? > > > > I will post the results of my script to hopefully help > others who are > > trying to formulate the same plot. > > > > Thanks again. > > > I don't know of an existing function that does that > particular > arrangement of plots, but here's a start ... > > > breaks <- seq(-3.5, 3.5, by=.5) > yrange <- range(breaks) > > histdata <- hist(y, breaks=breaks, plot=FALSE) > densitydata <- density(y) > > par(oma=rep(3, 4), mar=rep(0, 4)) > > layout(matrix(1:3, ncol=3), widths=c(.2, .1, 1), > respect=TRUE) > > barplot(histdata$density, horiz=TRUE, space=0, > xlim=c(0, max(histdata$density)*1.1), axes=FALSE) > par(new=TRUE) > plot(densitydata$y, densitydata$x, > xlim=c(0, max(histdata$density)*1.1), ylim=yrange, > type="l", ann=FALSE, axes=FALSE) > axis(2) > box() > boxplot(y, ylim=yrange, axes=FALSE) > box() > qqnorm(y, axes=FALSE, ylim=yrange, xlim=yrange, > main="") > abline(0, 1) > axis(1) > box() > > > ... the important bits are the call to layout() to get the > plots at > different sizes and the numerous efforts made to make sure > that the > coordinate systems of the various plots are coherent. > > Paul > > > > --- On Mon, 3/16/09, Jason Rupert > <jasonkrup...@yahoo.com> wrote: > > > >> From: Jason Rupert <jasonkrup...@yahoo.com> > Subject: [R] R package > >> to automatically produce combination plot? To: > R-help@r-project.org > >> Date: Monday, March 16, 2009, 8:14 PM By any > chance is there an R > >> package that automatically produces the plot shown > at the following > >> link: > http://n2.nabble.com/Can-R-produce-this-plot--td2489288.html > >> > >> That is an R package to produce on plot that has > the following: (a) > >> a vertically oriented histogram, (b) associated > barplot, and (c) > >> quantile-quantile plot (Q-Q Plot). > >> > >> This is based on a class lecture from University > of Pennsylvania: > >> > stat.wharton.upenn.edu/~mcjon/stat-431/lecture-02.pdf > >> > >> I am pretty confident I can put one together, but > just wanted to > >> check that there does not already exist an R > package to output such > >> a plot. > >> > >> Thanks again. > >> > >> ______________________________________________ > 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. > > > > ______________________________________________ > 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. > > -- > Dr Paul Murrell > Department of Statistics > The University of Auckland > Private Bag 92019 > Auckland > New Zealand > 64 9 3737599 x85392 > p...@stat.auckland.ac.nz > http://www.stat.auckland.ac.nz/~paul/ ______________________________________________ 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.