在贝类选择育种中,我们经常需要对亲本进行选择,如截头选择、判断亲本个体大小分布、截头选择的截点以及选择强度的计算等。因此,我们常需要做出关于某性状频率的直方图以及密度曲线。下面记录用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。
网友评论