前の質問に続いて(これは前の質問とは独立していますが)、次のような素敵な (そうではありませんか?!) 画像が得られました (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"))