0

sciplot の bargraph を使用すると、エラー バーを含む棒グラフをプロットできます。また、独立変数 (因子) によるグループ化も可能です。従属変数でグループ化したいのですが、どうすればそれを達成できますか

bargraph.CI(x.factor, response, group=NULL, split=FALSE,
col=NULL, angle=NULL, density=NULL,
lc=TRUE, uc=TRUE, legend=FALSE, ncol=1,
leg.lab=NULL, x.leg=NULL, y.leg=NULL, cex.leg=1,
bty="n", bg="white", space=if(split) c(-1,1),
err.width=if(length(levels(as.factor(x.factor)))>10) 0 else .1,
err.col="black", err.lty=1,
fun = function(x) mean(x, na.rm=TRUE),
ci.fun= function(x) c(fun(x)-se(x), fun(x)+se(x)),
ylim=NULL, xpd=FALSE, data=NULL, subset=NULL, ...)

bargraph.CI の仕様は上記の通りです。通常、応答変数は数値ベクトルです。今回は、同じ独立変数に対して 3 つの応答変数 (A、B、C) をプロットしたいと思います。データ フレーム "mpg" を使用して問題を説明します。次のコードでプロットを正常に取得できます。ここで DV は hwy です

data(mpg)
attach(mpg)

bargraph.CI(
class,  #categorical factor for the x-axis
hwy,    #numerical DV for the y-axis
group=NULL,   #grouping factor
legend=T, 
ylab="Highway MPG",
xlab="Class")

また、DV (hwy から cty に変更) だけを変更して、プロットを正常に取得することもできます。

data(mpg)
attach(mpg)

bargraph.CI(
class,  #categorical factor for the x-axis
cty,    #numerical DV for the y-axis
group=NULL,   #grouping factor
legend=T, 
ylab="Highway MPG",
xlab="Class")

ただし、2 つの DV を同時に使用する場合は、グループごとに 2 つのバー (1 つは cty 用、もう 1 つは hwy 用) を表示したいということです。

data(mpg)
attach(mpg)

bargraph.CI(
class,  #categorical factor for the x-axis
c(cty,hwy),    #numerical DV for the y-axis
group=NULL,   #grouping factor
legend=T, 
ylab="Highway MPG",
xlab="Class")

寸法の不一致により機能しません。どうすればこれを達成できますか？実際には、Boxplot schmoxplot の方法を使用して棒グラフの同様の効果を実現できます: R の係数によって条件付けられた平均値と標準誤差をプロットする方法は? ggplot2で。したがって、ggplot2 でそれを行う方法について何か考えがあれば、それは私にとっても問題ありません。

4

1 に答える 1

0

データを表示するときによくあることですが、最初にデータを操作してから を使用する必要がありますbargraph.CI。あなたの例では、data.frame視覚化したいのは次のとおりです。

df <- data.frame(class=c(mpg$class, mpg$class), 
                 value=c(mpg$cty, mpg$hwy), 
                 grp=rep(c("cty", "hwy"), each=nrow(mpg)))

次にbargraph.CI、この新しいdata.frame.

bargraph.CI(
  class,        #categorical factor for the x-axis
  value,        #numerical DV for the y-axis
  group=grp,    #grouping factor
  data=df, 
  legend=T, 
  ylab="Highway MPG",
  xlab="Class")
于 2013-10-21T15:46:53.953 に答える