paramerarization.R

minaka — Oct 30, 2013, 4:15 PM

# =====【標準正規分布からのサンプリング】=====

# 標準正規乱数500個のヒストグラムを描く
   data <- rnorm(500, mean=0, sd=1); 
     hist(data, freq=F);

# 位置パラメーター（平均）を図示する
     cat("平均 =", mean(data), "\n")

平均 = -0.005797 

     abline(v=mean(data), col="red", lwd=5)

#分散パラメーター（標準偏差）を図示する
     cat("標準偏差 =", sd(data), "\n")

標準偏差 = 0.959 

     xaxis <- c(mean(data)-sd(data), mean(data)+sd(data));
     yaxis <- c(0.2, 0.2);
     lines(xaxis, yaxis, col="blue", lwd=5);

# 標準正規分布と乱数500個のヒストグラムを描く（まとめて）
  data <- rnorm(500, mean=0, sd=1); 
    hist(data, freq=F);
    curve(dnorm(x, mean=0, sd=1), add=T);

  cat("平均 =", mean(data), "\n")

平均 = 0.05223 

  abline(v=mean(data), col="red", lwd=5)

  cat("標準偏差 =", sd(data), "\n")

標準偏差 = 1.043 

  xaxis <- c(mean(data)-sd(data), mean(data)+sd(data));
  yaxis <- c(0.2, 0.2);
  lines(xaxis, yaxis, col="blue", lwd=5);