myscratchpad

x <- rnorm(1000)

## 正規性の検定
hist(x)
qqnorm(x)
qqline(x)
shapiro.test(x) # p > .05なら正規分布

ks.test(x=x, y="pnorm", mean=mean(x), sd=sd(x))

## ks.testは「2群の分布が等しい」という帰無仮説を検定するもの。正規性の検定の場合、yに正規分布を指定している
x2 <- rexp(1000) # 正規分布じゃないやつ
ks.test(x=x, y=x2)


## 多変量正規性
library(mvnormtest)
data(EuStockMarkets)
C <- t(EuStockMarkets[15:29,1:4])
mshapiro.test(C)