2011-02-17から1日間の記事一覧

なんでn-1なんですか

R

不偏分散を計算する際になぜnではなくn-1で割るのか。理由はともかくnで割ってマズイのは割ってみれば分かります。 やってみよう サイコロ: 1から6までが等しい確率で出る理想的なサイコロ様を用意 分散の計算: サイコロをn回振るごとに分散と不偏分散を計算…

ggplot2のセパレータをドットにしたい

アンダースコアは library(ggplot2) for(o in ls(pos = "package:ggplot2")) assign((gsub("_", ".", o)), eval(as.symbol(o))) p <- qplot(wt, mpg, data = mtcars) p + geom_smooth(aes(group=cyl), method="lm") p + geom.smooth(aes(group=cyl), method=…