希腊字母和下标在`pair`图表矩阵的直方图标题中

Question

我得到a great answer a few days ago关于如何自定义pairs函数，但我需要将其与expression合并，我不知道如何在答案的框架内这样做给出。

虚拟例子：

N=1000
x1 = rgamma(N,5,10)
x2 = rnorm(N)+x1
x3 = (x1+x2)/(runif(N)+1)
d = data.frame(x1,x2,x3)
plot(d,col=rgb(1,0,0,.5),pch=19,cex=.5)

我正在寻找低位线和直方图，前面的回答者用以下方法做了：

panel.hist <- function(x, ...)
{
    usr <- par("usr"); on.exit(par(usr))
    par(usr = c(usr[1:2], 0, 1.5) )
    h <- hist(x, plot = FALSE)
    breaks <- h$breaks; nB <- length(breaks)
    y <- h$counts; y <- y/max(y)
    rect(breaks[-nB], 0, breaks[-1], y, ...)
}
panel.loess<-function(x, y, ...) {
    ll <- loess(y~x)
    points(x,y, ...)
    nx<-seq(min(x), max(x), length.out=100)
    lines(nx, predict(ll, nx), col="black")

}


pairs(d,col=rgb(1,0,0,.5),pch=19,cex=.5, 
   diag.panel=panel.hist, 
   lower.panel=panel.loess)

现在我想在标题中加上子/上标和希腊字母。如果我只是致电main = 'whatever'，我会expression(paste(hat('pr'),delta['t-1'],sep=''))。鉴于这不是一个字符串，我不能只添加一个字符向量参数panel.hist

有什么想法吗？