ビンエッジにラベルを合わせてRで可変幅のヒストグラムを作成するにはどうすればよいですか？

Question

私はggplot2を使用しています。デフォルトでは、ビン幅が固定されたヒストグラムが作成され、ビンのラベルは各ビンの中央にプロットされます。私が代わりに欲しい

は、そのビンのラベルこのプロットのように、各ビンのエンドポイントの代表である可変幅のヒストグラムである。この例のプロットを作成するには

、私は手動で入力しましたビンのパラメータを変更し、ビンをエンドポイントに揃えます。

income=data.frame(lx=c(0,10,25,50,100),rx=c(10,25,50,100,150),y=c(20,28,27,18,7)) 
income$width = income$rx-income$lx 


ggplot(income, aes(lx+width/2,y/width)) + geom_bar(aes(width=rx-lx), color='black', stat='identity') + 
    scale_x_continuous(breaks=unique(c(income$lx,income$rx))) + labs(x='Income (thousands of $)', y='% per thousand $') 

元のデータから自動的にこれを実行します。 （元のデータが、次のコードを使用して近似することができる）：

Answer 1

incomes=unlist(sapply(1:nrow(income), function(i) sample(income$lx[i]:(income$rx[i]-1),income$y[i],replace=TRUE))) 
widths=unlist(sapply(1:nrow(income), function(i) rep(income$rx[i]-income$lx[i],income$y[i]))) 
incomes=data.frame(incomes, widths) 

あなたgeom_histogramで所望breaksを指定することにより、可変幅のヒストグラムを生成することができます。 y=..density..（カウントに基づくデフォルトではなく）を使用します。これにより、バーは、バーの総面積の割合に正規化されます。

breaks = c(0,10,25,50,100,150) 

ggplot(incomes, aes(incomes)) + 
    geom_histogram(aes(y=..density..), 
       color="black", fill="grey40", breaks=breaks) + 
    scale_x_continuous(breaks=breaks)