||
刚有学生问如何去掉箱线图中的异常值(就是用点表示那些数据),然后求各组的平均值。我简单写个代码,供大家参考!
pixcel=read.csv("pixcel.csv",header=T)
outliers=function(data1)
{n=ncol(data1)
cmean=vector()
for (i in 1:n)
{
a=na.omit(data1[,i])
cmean[i]=mean(a[a> (quantile(a,0.25)-1.5*(quantile(a,0.75)-quantile(a,0.25)))&a< (quantile(a,0.75)+1.5*(quantile(a,0.75)-quantile(a,0.25)))])
}
boxplot(pixcel,outline = F)
return(cmean)
}
outliers(pixcel)
Archiver|手机版|科学网 ( 京ICP备07017567号-12 )
GMT+8, 2024-10-19 23:06
Powered by ScienceNet.cn
Copyright © 2007- 中国科学报社