Rに次の表があります(クランヘルプデータシートに触発されました):
> dfx <- data.frame(
+ group = c(rep('A', 108), rep('B', 115), rep('C', 106)),
+ sex = sample(c("M", "F","U"), size = 329, replace = TRUE),
+ age = runif(n = 329, min = 18, max = 54)
+ )
> head(dfx)
group sex age
1 A U 47.00788
2 A M 32.40236
3 A M 21.95732
4 A F 19.82798
5 A F 30.70890
6 A M 30.00830
エラーバーを含む棒グラフを使用して、各グループの男性 (M)、女性 (F)、および「不明」(U) のパーセンテージをプロットすることに興味があります。このグラフを作成するには、panel.ci/prepanel.ci コマンドを使用する予定です。
prop.table コマンドを使用して、各グループの比率テーブルを簡単に作成できます。
> with(dfx, prop.table(table(group,sex), margin=1)*100)
sex
group F M U
A 29.62963 28.70370 41.66667
B 35.65217 35.65217 28.69565
C 37.73585 33.01887 29.24528
しかし今、エラーバーを含む同様のテーブルを作成し、これら 2 つのテーブルを使用して棒グラフを作成したいと考えています。可能であれば、ddply
同様の目的で使用するコマンドを使用したいと思います(ただし、パーセンテージでも手段でもありませんでした)。