cl.r <- function(alpha)
{
r2 <- seq(-1, 1, 0.2)
r <- seq(-1, 1, 0.05)
z <- atanh(r)
z0 <- atanh(0)
plot(r, r, type="n", xaxt="n", yaxt="n",
xlab="r", ylab="rho")
abline(h=r, v=r, col="pink")
abline(h=r2, v=r2, col="pink3")
sapply(c(5, 10, 25, 50, 100, 250, 1000),
function(n)
{
temp <- qnorm(0.5+alpha/200)/sqrt(n-3)
lines(r, tanh(z+temp), col="blue")
lines(r, tanh(z-temp), col="blue")
text(0, tanh(z0+temp), paste("n =", n), pos=2)
text(0, tanh(z0-temp), paste("n =", n), pos=4)
}
)
axis(1, r2)
axis(2, r2)
}
# 使用例
cl.r(95)