===正規性検定=== '''qqnormを使う''' ※ データ x は予め平均0分散1に正規化しておく。 qqnorm( x ) abline(0, 1) 直線上に乗ったら正規分布 '''Shapiro-Wilk, Kolmogorov-Smirnov検定''' [http://cse.naro.affrc.go.jp/takezawa/r-tips/r/63.html ここを参照] ===好きな分布に従っていることを調べる。=== '''t分布かどうか''' ※ データ x は予め平均0分散1に正規化しておく。 free <- 2 plot( qt( ppoints(x), df=free ), sort(x) ) abline(0, 1) '''対数正規分布かどうか''' ※ データ x は予め平均0分散1に正規化しておく。 plot( qlnorm( ppoints(x) ), sort(x) ) abline(0, 1)
