添付の画像に似たグラフを再現する必要があります。この図を使用して、比率の差の信頼区間を比較したいと考えています。Rを使用して添付のグラフィックを作成するにはどうすればよいですか? 正しい方向への指針をいただければ幸いです。
質問する
1210 次
2 に答える
8
コンテキストと再現可能な例がなければ、適切な答えを出すのは困難です。でもストーリーは面白いと思います。
ここで私のggplot2の試み。alpha layer にはまだ問題がありますが、プロットの主なアイデアはここにあります。
一部のデータ
structure(list(scen = 1:6,
name = c("I", "II", "III", "IV", "V","VI"),
ymin = c(0.06, -0.102, 0.487, 0.116, -0.436, 0.021),
ymax = c(-0.231,0.135, 0.117, 0.338, -0.347, -0.025)),
.Names = c("scen", "name", "ymin", "ymax"),
row.names = c(NA, 6L),
class = "data.frame")
データはこんな感じ
データ
scen name ymin ymax y
1 1 I 0.060 -0.231 I
2 2 II -0.102 0.135 II
3 3 III 0.487 0.117 III
4 4 IV 0.116 0.338 IV
5 5 V -0.436 -0.347 V
6 6 VI 0.021 -0.025 VI
これが結果です
theme_new <- theme_set(theme_bw())
p <- ggplot(data=dat) +
geom_segment(aes(x=ymin,y=scen,xend=ymax,yend=scen),
arrow=arrow(length=unit(0.3,"cm"),
ends='both'),size=1)
p <- p+ geom_rect(xmin=max(dat$ymin)/2,
xmax=min(dat$ymax)/2,
ymin=0,
ymax=max(dat$scen)+1,
alpha=0.2,fill='grey')
p <- p + geom_text(aes(x=(ymin+ymax)/2,
y=scen+0.2,label =name),size=6)
p<- p + coord_cartesian(ylim=c(0,max(dat$scen)+3))+
xlab(expression(P[1]-P[0]))+
theme(
axis.ticks = element_blank(),
axis.text.y = element_blank(),
axis.text.x = element_blank(),
axis.title.x = element_text(face="bold", size=20))
p <- p + geom_vline(linetype ='dashed',
xintercept = mid.dash)
p <- p + geom_text(aes(x= mid.dash,
y = max(dat$scen)+2,
label="Zone of Indifference",
color="NA*"),rotate=180)
p <- p + theme(legend.position = "none")
p
于 2012-12-08T06:08:18.060 に答える