[R] R package to automatically produce combination plot?

[Jason Rupert]

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.murrell at auckland.ac.nz> wrote:

> From: Paul Murrell <p.murrell at auckland.ac.nz>
> Subject: Re: [R] R package to automatically produce combination plot?
> To: jasonkrupert at yahoo.com
> Cc: R-help at 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
> <jasonkrupert at yahoo.com> wrote:
> > 
> >> From: Jason Rupert <jasonkrupert at yahoo.com>
> Subject: [R] R package
> >> to automatically produce combination plot? To:
> R-help at 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 at 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 at 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
> paul at stat.auckland.ac.nz
> http://www.stat.auckland.ac.nz/~paul/