トップ 差分 一覧 ソース 検索 ヘルプ PDF RSS ログイン

正規性の検定

*disclaimer
49111

R

正規性の検定 というか確認


shapiro.test() シャピロ・ウィルク検定

shapiro.test(x)

(「正規分布している」が帰無仮説。つまり p>0.05なら帰無仮説を捨てない、つまり正規分布している)

  • いくつくらいに使えるか
    • SWは50以上、KSは400以上くらいといううわさ。(要出典)

ks.test() コルモゴロフ・スミルノフ検定

ks.test(x, "pnorm", alternative="two.sided")

(同上)

歪度、尖度を見てみる

(いずれも0に近いほど正規分布に近い)

> install.packages("psych")
> library(psych)
> describe(YNUL4_15$MDD)
   vars   n mean   sd median trimmed  mad min max range skew kurtosis   se
X1    1 978 1.97 0.56   1.86    1.93 0.51   1 3.9   2.9 0.75     0.38 0.02

ヒストグラムを描いてみる

hist(x)

Q-Qプロットして直線になっているか見てみる

qqnorm(x)