用R画直方图

用R画直方图和核密度图

条形图、饼图：分类型名义型数据展示 直方图：描述连续型数据
一.首先展示一下直方图：

接下来将展现如何画直方图的代码：

>par(mfrow=c(2,2))
>hist(mtcars$mpg)
>hist(mtcars$mpg,breaks=12,clo="red",xlab="miles per gallon",main="Colored histogram with 12 bins")

以上语句画了两个直方图，第一个直方图没有任何注释，第二个直方图指定了组数和颜色，并且添加了标题。以下是添加了密度曲线和轴须图：为数据的分布提供了一种平滑的描述。

>hist(mtcars$mpg,freq=FALSE,breaks=12,col="red",xlab="Miles per Gallon",main="Histogram,rug plot,density curve")
>rug(jitter(mtcars$mpg))
>lines(density(mtcars$mpg),col="blue",lwd=2)

>x<-mtcars$mpg
>h<-hist(x,breaks=12,col="red",xlab="Miles per Gallon",main="histogram with normal curve and box")
>xfit<-seq(min(x),max(x),length=40)
>yfit<-dnorm(xfit,mean=mean(x),sd=sd(x))
>lines(xfit,yfit,col="blue",lwd=2)
>box()

二.核密度估计是用于估计随机变量概率密度函数的一种非参数方法。核密度图不失为一种用来观察连续型变量分布的有效方法。

> par(mfrow=c(2,1))
> d<-density(mtcars$mpg)
> plot(d)
Error in plot.new() : figure margins too large
> d<-density(mtcars$mpg)---使用默认设置创建最简图形
> plot(d,main="abd")--添加标题
> polygon(d,col="red",border="blue")---将曲线改为蓝色，使用实心红色填充下方区域
> rug(mtcars$mpg,col="brown")--添加棕色的轴须图

如果遇到上面报的Error in plot.new() : figure margins too large错误，将plot窗口调大即可。核密度图如下：

结束：感谢R语言实战。