在线时间:8:00-16:00
迪恩网络APP
随时随地掌握行业动态
扫描二维码
关注迪恩网络微信公众号
> boxplot(mtcars$mpg, main="box plot", ylab="miles per gallon") > boxplot.stats(mtcars$mpg) $stats [1] 10.40 15.35 19.20 22.80 33.90 $n [1] 32 $conf [1] 17.11916 21.28084 $out numeric(0) 1. 使用并列箱型图进行跨组比较 > boxplot(mpg ~ cyl, data=mtcars, + main="cars mileage data", + xlab="number of cylinders", + ylab="miles per gallon") > > boxplot(mpg ~ cyl, data=mtcars, + notch=TRUE, + varwidth=TRUE, + col="red", + main="Car Mileage Data", + xlab="Number of Cylinders", + ylab="Miles Per Gallon") > mtcars$cyl.f <- factor(mtcars$cyl, + levels=c(4,6,8), + labels=c("4","6","8")) > mtcars$am.f <- factor(mtcars$am, + levels=c(0,1), + labels=c("auto", "standard")) > boxplot(mpg ~ am.f *cyl.f, + data=mtcars, + varwidth=TRUE, + col=c("gold", "darkgreen"), + main="MPG Distribution by Auto Type", + xlab="Auto Type") > 2. 小提琴图 > library(vioplot) > x1 <- mtcars$mpg[mtcars$cyl==4] > x2 <- mtcars$mpg[mtcars$cyl==6] > x3 <- mtcars$mpg[mtcars$cyl==8] > vioplot(x1, x2, x3, + names=c("4 cyl", "6 cyl", "8 cyl"),) Error in vioplot(x1, x2, x3, names = c("4 cyl", "6 cyl", "8 cyl"), ) : argument is missing, with no default > vioplot(x1, x2, x3, + names=c("4 cyl", "6 cyl", "8 cyl"), + col="gold") > title("Violin Plots of Miles Per Gallon") >
|
请发表评论