1. 程式人生 > >R中三種檢驗正態分佈的方式

R中三種檢驗正態分佈的方式

一、畫出密度函式與正態分佈密度圖比較:

library(MASS)
mu<- c(0,0,0)
Sigma<- matrix(c(1,0.5,0.25,0.5,1,0.5, 0.25,0.5,1),3,3)
M<- mvrnorm(1000, mu, Sigma)
d<- density(M[,2])
plot(d$x,dnorm(d$x),col=2,type = "n")
lines(d$x,dnorm(d$x),col=2,lwd=2)
lines(d,col=3,lwd=2)

在這裡:density函式的返回值:x,y給定觀測,用核函式估計出的密度值y。 標準正太和已知資料的圖 二、QQ-PLOT 兩個quantile。quantile是我們的資料和標準正態分佈的很多分位陣列成的陣列畫出來的圖。橫座標是標準的正態分佈的quantile 縱座標是我們資料的圖。如果兩者基本相等,則這個散點圖會非常像y=x

qqnorm(M[,2]);qqline(M[,2],col=2)

qqplot 三、shapiro假設檢驗

shapiro.test(M[,2])

testresult