opar<-par(no.readonly=TRUE)
par(cex=1.5,lty=5)
# library(sm)
attach(mtcars)
cyl.f<-factor(cyl, levels = c(4, 6, 8), labels = c("4 缸", "6 缸", "8 缸"))
sm.density.compare(mpg, cyl, xlab = "MPG" , lwd=4)
title(main = "不同气缸汽车MPG 分布")
colfill<-c(26,254,552)
legend("topright" , levels(cyl.f), fill = c("blue","green","red") )
detach(mtcars)
par(opar)