確率分布を図示する簡単な関数の作成。
> normd <- function (z1,z2)
+ {
+ x <- seq(-3,3,0.01)
+ plot(x,dnorm(x,mean=0,sd=1.0),type="n")
+ curve(dnorm(x,mean=0,sd=1.0),type="l",add=T)
+ xaxis <- seq(z1, z2, length=100)
+ yaxis <- c(dnorm(xaxis, mean=0, sd=1.0), 0, 0)
+ xaxis <- c(xaxis, z2, z1)
+ polygon(xaxis, yaxis, density=25)
+ ypos <- dnorm(z1, mean=0, sd=1.0)
+ text(z1, ypos,round(pnorm(z2)-pnorm(z1),digits=4), adj=0)
+ }
>
> normd(0,1)
> normd <- function (z1,z2)
+ {
+ x <- seq(-3,3,0.01)
+ plot(x,dnorm(x,mean=0,sd=1.0),type="n")
+ curve(dnorm(x,mean=0,sd=1.0),type="l",add=T)
+ xaxis <- seq(z1, z2, length=100)
+ yaxis <- c(dnorm(xaxis, mean=0, sd=1.0), 0, 0)
+ xaxis <- c(xaxis, z2, z1)
+ polygon(xaxis, yaxis, density=25)
+ ypos <- dnorm(z1, mean=0, sd=1.0)
+ text(z1, ypos,round(pnorm(z2)-pnorm(z1),digits=4), adj=0)
+ }
>
> normd(0,1)