R语言与直方图与核密度图

 

如果想要研究某一变量的分布情况,则需要使用直方图和密度曲线图。

函数hist(x,breaks=,)绘制直方图。

x:是由一个数据值组成的数值向量

breaks:用于控制组的数量,breaks=seq(220,280,3)从220开始到280,间隔为3

freq=FALSE:控制y轴,使其变为密度​

rug(jitter(mtcar$mpg))​:轴须图

代码:

hist(mtcars$mpg,breaks=12,col="red4",xlab="miles per gallon",main="histogram")

R语言与直方图与核密度图_第1张图片直方图

带有密度曲线的直方图

代码:

​hist(mtcars$mpg,freq=FALSE,breaks=12,col="red4",xlab="miles per gallon",main="histogram")

rug(jitter(mtcars$mpg))

lines(density(mtcars$mpg),col="blue4",lwd=2)

R语言与直方图与核密度图_第2张图片

核密度图:

代码:

(1)d<-density(mtcars$mpg)

plot(d)

R语言与直方图与核密度图_第3张图片

(2)d<-density(mtcars$mpg)

plot(d,main="he mi Du tu ")

​polygon(d,col="red4",border="blue4")

rug(mtcars$mpg,col="brown")

R语言与直方图与核密度图_第4张图片

创建可比较的核密度图

函数:sm.density.compare(x,factor)

x:是一个数值型向量,factor是一个分组变量​

代码:

library​(sm)

attach(mtcars)

cyl.f<-factor(cyl,levels = c(4,6,8),labels = c("4 cylinder","6 cylinder","8 cylinder"))

sm.density.compare(mpg,cyl.f,xlab="miles per gallon")

title(main = "MPG distribution by car cylinders")

legend("topright",pch = 15,legend = c("4 cylinder","6 cylinder","8 cylinder"),col = c("blue","green","red"))

R语言与直方图与核密度图_第5张图片

你可能感兴趣的:(41,学术科研)