在上一个问题之后(不过,这个问题独立于上一个问题),我最终得到了以下漂亮的(不是吗?!)图片(R 代码如下所示)。
但是,当我以 pdf 格式导出图片时,质量有点令人失望;右上角出现虚假的白线:(
你有解决方案吗?
library(gplots)
f <- function(S1, S2, theta)
{
num <- (S1^(-theta) + S2^(-theta) - 1)^(-1 / theta)
denom <- S1 * S2
return(num / denom)
}
S1 <- S2 <- seq(0.05, 0.95, 0.01)
z <- outer(S1, S2, FUN=f, theta=2/3)
levels=seq(1.1, 8, length.out=25)
windows(5,5)
par(mgp=c(2, 0.5, 0))
plot(NA,
xlim=rev(c(0, 1)), xlab=expression(S[1](t)), xaxt="n",
ylim=rev(c(0, 1)), ylab="", yaxt="n",
frame=FALSE)
mtext(text=bquote(S[2](t)), side=2, line=1.8, las=1)
axis(1, c(1, 0), label=TRUE, tcl=-0.5)
axis(1, seq(1, 0, -0.25), label=FALSE, tcl=-0.25)
axis(2, c(1, 0), label=TRUE, tcl=-0.5)
axis(2, seq(1, 0, -0.25), label=FALSE, tcl=-0.25)
par(mgp=c(3, 1, 0))
.filled.contour(x=S1, y=S2, z=z,
levels=levels,
col=colorpanel(length(levels), "gray95", "gray20", "black"))