R语言画经验分布函数，直方图，QQ图，实现正态检验（非参数统计，习题1.11）

习题1.11：（1）画出beenswax数据的经验分布函数，直方图和QQ图

data1.11<-read.table("beenswax.txt",header=TRUE)   #读取数据
Hy<-data1.11$Hydrocarbon

par(mfrow=c(2,2))    #画两行两列四个图

hist(Hy,breaks=6,freq=TRUE,main="频数直方图")

hist(Hy,breaks=6,freq=FALSE,main="频率密度直方图")

plot(ecdf(Hy),verticals=TRUE,do.points=FALSE,main="经验分布函数图",col="red")

qqnorm(Hy,main="QQ图")

qqline(Hy)

（2）找出0.9，0.75，0.5，0.25和0.10分位数.

q1<-quantile(Hy,0.90)

q2<-quantile(Hy,0.75)

.....

(3)这个分布是正态分布吗？

使用Shapiro—Wilk检验：原假设 :  该随机变量服从正态分布；

                                         备择假设：该随机变量不服从正态分布.

 shapiro.test(Hy)

输出结果：

Shapiro-Wilk normality test

data:  Hy
W = 0.9766, p-value = 0.3123

P值>0.05, 不能拒绝原假设，该分布是正态分布.

数据：

 Hy
14.27 14.80 12.28 17.09 15.10 12.92 15.56 15.38 15.15 13.98 14.90 15.91 14.52 15.63 13.83
13.66 13.98 14.47 14.65 14.73 15.18 14.49 14.56 15.03 15.40 14.68 13.33 14.41 14.19 15.21
14.75 14.41 14.04 13.68 15.31 14.32 13.64 14.77 14.30 14.62 14.10 15.47 13.73 13.65 15.02
14.01 14.92 15.47 13.75 14.87 15.28 14.43 13.96 14.57 15.49 15.13 14.23 14.44 14.57