2014-03-07 35 views
3

我有两种类型,看起来像这样的数据: 1型(http://dpaste.com/1697615/plain/为什么GGPLOT2饼图方面混淆了小标签

Cluster-6 abTcells 1456.74119 
Cluster-6 Macrophages 5656.38478 
Cluster-6 Monocytes 4415.69078 
Cluster-6 StemCells 1752.11026 
Cluster-6 Bcells 1869.37056 
Cluster-6 gdTCells 1511.35291 
Cluster-6 NKCells 1412.61504 
Cluster-6 DendriticCells 3326.87741 
Cluster-6 StromalCells 2008.20603 
Cluster-6 Neutrophils 12867.50224 
Cluster-3 abTcells 471.67118 
Cluster-3 Macrophages 1000.98164 
Cluster-3 Monocytes 712.92273 
Cluster-3 StemCells 557.88648 
Cluster-3 Bcells 599.94109 
Cluster-3 gdTCells 492.61994 
Cluster-3 NKCells 524.42522 
Cluster-3 DendriticCells 647.28811 
Cluster-3 StromalCells 876.27875 
Cluster-3 Neutrophils 1025.24105 

和类型二,(http://dpaste.com/1697602/plain/)。与类型1的数据处理该代码时

abTcells 1456.74119 
Macrophages 5656.38478 
Monocytes 4415.69078 
StemCells 1752.11026 
Bcells 1869.37056 
gdTCells 1511.35291 
NKCells 1412.61504 
DendriticCells 3326.87741 
StromalCells 2008.20603 
Neutrophils 12867.50224 

但是为什么:

library(ggplot2); 
library(RColorBrewer); 
filcol <- brewer.pal(10, "Set3") 
dat <- read.table("http://dpaste.com/1697615/plain/") 
ggplot(dat,aes(x=factor(1),y=dat$V3,fill=dat$V2))+ 
    facet_wrap(~V1)+ 
    xlab("") + 
    ylab("") + 
    geom_bar(width=1,stat="identity",position = "fill") + 
    scale_fill_manual(values = filcol,guide = guide_legend(title = "")) + 
    coord_polar(theta="y")+ 
    theme(strip.text.x = element_text(size = 8, colour = "black", angle = 0)) 

就绪数据:

> dput(dat) 
structure(list(V1 = structure(c(2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 
2L, 2L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L), .Label = c("Cluster-3", 
"Cluster-6"), class = "factor"), V2 = structure(c(1L, 5L, 6L, 
9L, 2L, 4L, 8L, 3L, 10L, 7L, 1L, 5L, 6L, 9L, 2L, 4L, 8L, 3L, 
10L, 7L), .Label = c("abTcells", "Bcells", "DendriticCells", 
"gdTCells", "Macrophages", "Monocytes", "Neutrophils", "NKCells", 
"StemCells", "StromalCells"), class = "factor"), V3 = c(1456.74119, 
5656.38478, 4415.69078, 1752.11026, 1869.37056, 1511.35291, 1412.61504, 
3326.87741, 2008.20603, 12867.50224, 471.67118, 1000.98164, 712.92273, 
557.88648, 599.94109, 492.61994, 524.42522, 647.28811, 876.27875, 
1025.24105)), .Names = c("V1", "V2", "V3"), class = "data.frame", row.names = c(NA, 
-20L)) 

这些值在1型与集群-6相同的上述生成以下数字: enter image description here

请注意,Facet标签放错位置,Cluster-3应为Cluster-6, ,其中嗜中性粒细胞的比例较大。

我该如何解决问题?

处理类型2数据时完全没有问题。

library(ggplot2)     
df <- read.table("http://dpaste.com/1697602/plain/"); 
library(RColorBrewer); 
filcol <- brewer.pal(10, "Set3") 

ggplot(df,aes(x=factor(1),y=V2,fill=V1))+ 
    geom_bar(width=1,stat="identity")+coord_polar(theta="y")+ 
    theme(axis.title = element_blank())+ 
    scale_fill_manual(values = filcol,guide = guide_legend(title = "")) + 
     theme(strip.text.x = element_text(size = 8, colour = "black", angle = 0)) 

就绪数据:

> dput(df) 
structure(list(V1 = structure(c(1L, 5L, 6L, 9L, 2L, 4L, 8L, 3L, 
10L, 7L), .Label = c("abTcells", "Bcells", "DendriticCells", 
"gdTCells", "Macrophages", "Monocytes", "Neutrophils", "NKCells", 
"StemCells", "StromalCells"), class = "factor"), V2 = c(1456.74119, 
5656.38478, 4415.69078, 1752.11026, 1869.37056, 1511.35291, 1412.61504, 
3326.87741, 2008.20603, 12867.50224)), .Names = c("V1", "V2"), class = "data.frame", row.names = c(NA, 
-10L)) 

enter image description here

+0

你能写出你的数据,以便随时复制/粘贴到R? –

+0

我所有的代码都很容易复制粘贴以生成上面的图形。 – pdubois

+0

您的绘图命令很好,但数据不容易放入工作区。请参阅http://stackoverflow.com/questions/5963269/how-to-make-a-great-r-reproducible-example关于如何改进的一些提示。 –

回答

4

这是因为你在aes(...)使用数据帧的名字。这解决了这个问题。

ggplot(dat,aes(x=factor(1),y=V3,fill=V2))+ 
    facet_wrap(~V1)+ 
    xlab("") + 
    ylab("") + 
    geom_bar(width=1,stat="identity",position = "fill") + 
    scale_fill_manual(values = filcol,guide = guide_legend(title = "")) + 
    coord_polar(theta="y")+ 
    theme(strip.text.x = element_text(size = 8, colour = "black", angle = 0)) 

在定义方面,可以通过电平在默认数据集的上下文中引用V1,和ggplot排序字母(因此“群集-3”至上)。在您拨打aes(...)时,您直接参考dat$V3,因此ggplot将从默认数据集的上下文中跳出到原始数据帧。在那里,Cluster-6是第一个。

作为一般评论,应该在aes(...)从未参考数据外侧data=...定义的数据集的上下文中。所以:

ggplot(data=dat, aes(y=V3...))  # good 
ggplot(data=dat, aes(y=dat$V3...)) # bad 

你的问题是为什么第二个选项不好的完美例子。

+0

对不起,为什么我不能用您的代码重现您的结果.. – pdubois

+0

因为我粘贴了错误的代码...抱歉 - 编辑应该修复它。 – jlhoward