直方图+密度曲线

       在贝类选择育种中，我们经常需要对亲本进行选择，如截头选择、判断亲本个体大小分布、截头选择的截点以及选择强度的计算等。因此，我们常需要做出关于某性状频率的直方图以及密度曲线。下面记录用R做直方图及密度曲线的简单方法。

hist（）函数画直方图

hist()函数的常用格式为：hist(x,break=,freq=)

其中x是一个数值型向量，breaks用于控制组的数量，freq用于选择概率密度（freq=F）或频数（freq=T），只有当选择概率密度（freq=F）时，才能在直方图上加密度曲线。

例：数据为data，做其中一列关于壳长的直方图并画上密度曲线。

代码：

>data=read.table("clipboard",header=T)········复制后读取数据

> hist(data$壳长, freq=FALSE, breaks=15, col=c("white", "white", "white", "white", "white", "white", "white", "white", "white", "blue","blue","blue"))·······画data中关于壳长的直方图，并将截点设置为蓝色

>lines(density(data$壳长),lwd=3,col="red")·······直方图上添加密度曲线，并设置线条粗细为3，颜色为红色

参考文献

1、林元震，R与ASReml-R统计学。中国林业出版社。

2、王庆志，长牡蛎品种选育与生长性状的遗传参数估计。中国海洋大学研究生学位论文。2011。