0

ggplot2 を使用して、添付のプロットに表示されるポイントに信頼区間を追加しようとしています。(2 つのデータセット、X 軸 1 つ、Y 軸 2 つ)

ポイントは平均値を表します (データフレームの val 列)。データ フレームには、上限値と下限値も含まれます (L 列と U 列)。

t<-data.frame(rep(c( "AA","AG","GG"),2), c(79.29365,47.67004,22.86688,28.41990,33.29651,36.76852),
           c(71.18651,45.97425,22.04625,27.23392,32.97892,36.57194),
           c(87.40079,49.36582,23.68751,29.60588,33.61411,36.96510),
           c(rep("a",3),c(rep("b",3))))
colnames(t)<-c("x","val","L","U","panel")

d1<-t[1:3,]
d2<-t[4:6,]

p <- ggplot(data = t, mapping = aes(x =t$x, y = t$val))
p <- p + facet_grid(panel ~ ., scale = "free")
p <- p + layer(data = d1, mapping=aes(x = d1$x, y = d1$val), geom = c( "point"), stat = "identity",size=3) 
p <- p + layer(data = d2, mapping =aes(x = d2$x, y = d2$val),geom = "point", stat = "identity",size=3)
p<-p+theme_bw()
p<-p+xlab("Genotypes")+ylab("Quantitative trait")+ggtitle("RS")
p<-p+scale_color_hue()
p

ここに画像の説明を入力

1 つのデータセットの単純な ggplot グラフにエラー バーをプロットすることができましたが、上記で試みたように、エラー バーを両方のプロットに正しく追加できませんでした。

ggplot(d1, aes(x = d1$x, y = d1$val)) +  geom_point(colour="gray", size = 4)+geom_errorbar(aes(ymax = d1$U, ymin = d1$L),width=0.05) +theme_bw()

前もって感謝します。

4

1 に答える 1