1

こんにちは、以下のデータで ggplot2 を使用して積み上げバープロットを作成したかったのです

Chr NonSyn_Snps Total_exonic_Snps
A01 9217    13725
A02 6226    9133
A03 14888   21531
A04 5272    7482
A05 4489    6608
A06 8298    12212
A07 6351    9368
A08 3737    5592
A09 12429   18119
A10 7165    10525

基本的に、染色体ごとに NonSyn_Snps と Total_exonic_Snps をスタックしたいのですが、残念ながらできません。

これは私がこれまで運がなかったものです

ggplot(Chr.df_mod, aes(Chr, Total_exonic_Snps, fill = NonSyn_Snps)) + geom_bar(stat = "identity", colour = "white") + xlab("Chromosome") + ylab("Number of SNPs")

プロットを取得していますが、積み上げられていません。

ここに画像の説明を入力

誰かがこれのトラブルシューティングを手伝ってくれませんか。

ありがとうウペンドラ

4

1 に答える 1

6

このggplotイディオムは、幅の広いデータよりも長いデータに最適です。ggplot の多くのオプションを利用するには、ワイド データ フレームをロング フォーマットに変換する必要があります。

# get data
dat <- read.table(text = "Chr NonSyn_Snps Total_exonic_Snps
A01 9217    13725
A02 6226    9133
                  A03 14888   21531
                  A04 5272    7482
                  A05 4489    6608
                  A06 8298    12212
                  A07 6351    9368
                  A08 3737    5592
                  A09 12429   18119
                  A10 7165    10525", header= TRUE)


# load libraries
require(ggplot2)
require(reshape2)

# melt data from wide to long
dat_m <- melt(dat)

# plot
ggplot(dat_m, aes(Chr, value, fill = variable)) + 
  geom_bar(stat = "identity") + 
  xlab("Chromosome") + 
  ylab("Number of SNPs")

ここに画像の説明を入力

于 2013-07-02T02:50:14.440 に答える