将极值箱分组为一个"> x"箱

Question

是否已存在函数/方法来确定大于某个值的数据频率？与Excel频率分布类似,我想将极值组合到最后一个bin中(例如,> 120,如图像中所示).我一直在手动执行此操作,首先使用hist函数,然后将大于给定值的中断计数求和.

Answer 1

这是一个选项:

d <- rlnorm(1000, 3)
d.cut <- cut(d, c(seq(0, 120, 10), Inf))
hist(as.numeric(d.cut), breaks=0:13, xaxt='n', xlab='', 
     col=1, border=0, main='', cex.axis=0.8, las=1)
axis(1, at=0:13, labels=c(seq(0, 120, 10), '>120'), cex.axis=0.8)
box()