[R] Need a little help setting the upper median using "layout"...

[Jason Rupert]

The code I'm using is shown below.  

I would like to have a larger median at the top of the plot so that I can show the entirity of "title_text". 

Several times I tried messing with "par(mar", but that seemed to make matters worse.  

By any chance can anyone provide any insight as to the best way to increase the top/upper/northern margin, so the entirity of the title is visible?  

Thank you again for any help that can be provided.  


x <- rnorm(10000)
# op <- par(mar=c(2,2,2,2), oma=c(0,0,0,0)+.8)

op <- par(mar=c(2, 2, 2, 2), oma=c(0,0,0,0)+.8)
layout(matrix(c(2, 2, 2, 3, 3, 3, 
                1, 1, 1, 3, 3, 3, 
                1, 1, 1, 3, 3, 3, 
                1, 1, 1, 3, 3, 3), 4, 6, byrow=TRUE))
#  2  2   2  3  3  3
#  1  1   1  3  3  3
#  1  1   1  3  3  3
#  1  1   1  3  3  3

# 1st Plot
hist(x, probability=TRUE, breaks=20, col="light blue", main="", xlab="Density" )
rug(jitter(x, 5))
points(density(x), type='l', lwd=3, col='red')
box()

# 2nd Plot
#par(mar=c(0,0,0,0), oma=c(0,0,0,0)+.2)
#boxplot( x, horizontal=T, col="white", xaxt='n', yaxt='n')
boxplot(x, horizontal = TRUE, col = "pink", lwd=1, xaxt='n', yaxt='n')  
box(col="white")

#par(op)
# 3rd Plot
#par(mar=c(1,1,1,1), oma=c(0,0,0,0)+.8)
qqnorm(x, , main="")
qqline(x, col="red", lwd=3)

kurtosis_val <-(sum((x-mean(x))^4 ))/(var(x)*var(x))/length(x)-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<-c("Title Text", "\nMean = ", format(mean(x), digits=4, scientific=F), " Standard Deviation = ", format(sd(x), digits=4, scientific=F),"\n Skewness = ", format(skewness(x), digits=4, scientific=F), "Kurtosis =", format(kurtosis_val, digits=4, scientific=F))

mtext(title_text, NORTH<-3, line=0, adj=0.5, cex=1.2, col="red", outer=TRUE)