R软件，分段直方图

setwd("C:\\RBook\\数据研究\\重庆数据研究\\面积\\")

mhist<-function(sData,freq="T",Color="black",lty="solid",brk=NULL)
{
    exf<-hist(sData,plot="F",breaks=brk)
    brk<-exf$breaks
    if(freq)
        cnt<-exf$intensities
    else
        cnt<-exf$counts
        
    for(i in 1:length(cnt))
    {
        n=i+1
        #第一条竖线
        lines(c(brk[i],brk[i]),c(0,cnt[i]),lty=lty,col=Color)
        #第二条竖线
        lines(c(brk[n],brk[n]),c(0,cnt[i]),lty=lty,col=Color)
        #顶线
        lines(c(brk[i],brk[n]),c(cnt[i],cnt[i]),lty=lty,col=Color)
    }
}
##三个数据文件。只有一列 mp
files=c("all.csv","lessT90.csv","90T120.csv","greaterT120.csv")

xMin=NULL
xMax=NULL
yMin=NULL
yMax=NULL
mp_seq=300
brk=seq(0,15000,mp_seq)

for(fil in files)
{
    d<-read.table(fil,header=T)
    exf<-hist(d$mp,plot="F",brk)
    xMin=min(xMin,exf$breaks)
    xMax=max(xMax,exf$breaks)
    yMin=min(yMin,exf$intensities)
    yMax=max(yMax,exf$intensities)
}
plot(x=c(xMin,xMax),xlab="",ylab="",xaxt="n",frame = FALSE,main="2012年上半年不同面积段的均价直方",y=c(yMin,yMax),type="n")

axis(1,brk,pos=0,las=2)
ltys<-c("solid","dashed","dotdash","dotted")
cols<-c("black","green","red","blue")
legends<-c("all","<90","90-120",">=120")

for(i in 1:length(files))
{
    d<-read.table(files[i],header=T)
    mhist(d$mp,freq="T",Color=cols[i],lty=ltys[i],brk=brk)
}
legend(x=12000,y=2e-04,lty=ltys,col=cols,legend=legends)

结果如下：