理学療法士がまとめたノート

統計学備忘録 since2016

Rを使って統計学を勉強するブログです

正規分布の色塗り

plot(dnorm, -4, 4, xaxt="n")
xvals <- seq(-4, -1.96, length=10)        # -4以上-1.33以下 領域をx軸方向に10個の多角形(台形)に等分割
dvals <- dnorm(xvals)    # 対応するグラフの高さ
polygon(c(xvals,rev(xvals)),
        c(rep(0,10),rev(dvals)),col=5)       # 塗りつぶす

xvals2 <- seq(1.96, 4, length=10)        
dvals2 <- dnorm(xvals2)   
polygon(c(xvals2,rev(xvals2)),
        c(rep(0,10),rev(dvals2)),col=5) 

name<-c("-1.96","0","1.96")
axis(side=1,at=c(-1.96,0,1.96),labels=name)  
f:id:yoshida931:20180719112607p:plain:w500

N(45,102)におけるz<40の確率を示す図の描き方

pn <- pnorm(40, mean=45, sd=10)

par(ps=15)
plot(dnorm, -4, 4, xaxt="n", xlab="N(45,10^2)")
xvals <- seq(-4, qnorm(pn), length=100) # -4以上-1.33以下 領域をx軸方向に10個の多角形(台形)に等分割
dvals <- dnorm(xvals)    # 対応するグラフの高さ

polygon(c(xvals,rev(xvals)),
        c(rep(0,100),rev(dvals)),col=7)       # 塗りつぶす

name <- c("40")
axis(side=1,at=c(qnorm(pn)),labels=name)  
f:id:yoshida931:20180904134826p:plain