3 
totalReads=10000
mappedReads=8000
mappedTargetReads=6000
unmappedReads=totalReads-mappedReads
offtargetmappedReads=mappedReads-mappedTargetReads
data<-c(mappedReads,unmappedReads,mappedTargetReads,offtargetmappedReads)
round(data[1]/totalReads*100)->a
round(data[2]/totalReads*100)->b
round(data[3]/mappedReads*100)->c
round(data[4]/mappedReads*100)->d
pct<-c(a,b,c,d)
labels<-c("MappedReads","UnmappedReads","On-exon","Off-exon")
labels<-paste(labels,pct)
labels<-paste(labels,"%", sep="")
pie(c(mappedReads,unmappedReads,mappedTargetReads,offtargetmappedReads),labels=labels)

基本的に、4 つのスライスすべてを個別にプロットします。しかし、ここmappedTargetReadsでは、 はマッピングされた読み取りに関して計算されるため、スライスを細分化してを内で表現offtargetmappedReadsしたいと思います。mappedReadsmappedTargetReadsofftargetmappedReadsmappedReads

または、1 つの Rplot で 2 つの異なる円グラフを並べてプロットすることは可能ですか?

4

1 に答える 1

2

これにより、2つの円グラフを含むプロットが生成されます。

par(mfrow = c(1,2))
pie(c(mappedReads,unmappedReads),labels=labels[1:2], main = "Global")
pie(c(mappedTargetReads,offtargetmappedReads),labels=labels[3:4],
    main = "MappedReads")

ここに画像の説明を入力してください

于 2012-10-28T10:05:57.643 に答える