ggplot

Question

我正在尝试制作一个ggplot。当我在美学方面塑造时，代码工作得很好。但是，我需要在geom_point（）中设置形状，因为我正在尝试重现一个数字。当我向geom_point（）添加形状时，它给了我以下错误：

美学必须是长度1或与数据（6）相同：形状

我在这里寻找其他答案，但显然，似乎没有什么对我有用。上面我提供了我的数据的图像。有17000个条目。

以下是我的代码：

summarised_data <-ddply(mammals,c('mammals$chr','mammals$Species','mammals$chrMark'), 
function (x) c(median_rpkm = median(x$RPKM), median = median(x$dNdS)))

ggplot(summarised_data,aes(x = summarised_data$median_rpkm, y = summarised_data$median, 
color = summarised_data$`mammals$Species`)) + geom_smooth(se = FALSE, method = "lm") +
geom_point(shape = summarised_data$`mammals$chrMark`) + xlab("median RPKM") + ylab("dNdS")

  "ENSG00000213221", "ENSG00000213341", "ENSG00000213380", "ENSG00000213424", 
   "ENSG00000213533", "ENSG00000213551", "ENSG00000213619", "ENSG00000213626", 
  "ENSG00000213699", "ENSG00000213782", "ENSG00000213949", "ENSG00000214013", 
   "ENSG00000214338", "ENSG00000214357", "ENSG00000214367", "ENSG00000214517", 
    "ENSG00000214814", "ENSG00000215203", "ENSG00000215305",      "ENSG00000215367", 
   "ENSG00000215440", "ENSG00000215897", "ENSG00000221947", "ENSG00000222011", 
   "ENSG00000224051", "ENSG00000225830", "ENSG00000225921", "ENSG00000239305", 
   "ENSG00000239474", "ENSG00000239900", "ENSG00000241058", "ENSG00000242247", 
   "ENSG00000242612", "ENSG00000243646", "ENSG00000244038", "ENSG00000244045"), 
  class = "factor"), Species = structure(c(1L, 1L, 1L, 1L, 1L, 
   1L, 1L, 1L, 1L, 1L), .Label = c("Chimp", "Gori", "Human", "Maca", 
    "Mouse", "Oran"), class = "factor"), labs = structure(c(2L, 2L, 
    2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L), .Label = c("Chimp-A", "Chimp-X", 
    "Gori-A", "Gori-X", "Human-A", "Human-X", "Maca-A", "Maca-X", 
    "Mouse-A", "Mouse-X", "Oran-A", "Oran-X"), class = "factor"), 
     chrMark = structure(c(2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 
      2L), .Label = c("A", "X"), class = "factor"), chr = structure(c(27L, 
      27L, 27L, 27L, 27L, 27L, 27L, 27L, 27L, 27L), .Label = c("1", 
      "10", "11", "12", "13", "14", "15", "16", "17", "18", "19", 
      "2", "20", "21", "22", "2a", "2A", "2b", "2B", "3", "4", 
        "5", "6", "7", "8", "9", "X"), class = "factor"), dN = c(3.00669, 
       3.27182, 7.02044, 1.01784, 3.0363, 2.32786, 4.92959, 3.03753, 
       3.0776, 1.02147), dS = c(3.15631, 5.87147, 3.13716, 2.05438, 
       4.10205, 5.24764, 4.2014, 3.18086, 5.4942, 3.02169), dNdS = c(0.9525965447, 
     0.5572403504, 2.2378329444, 0.4954487485, 0.7401908802, 0.4436013141, 
1.1733207978, 0.954939859, 0.5601543446, 0.3380459279), RPKM = c(31.6, 
13.9, 26.3, 9.02, 11.3, 137, 242, 1.05, 59.4, 10.1), Tau = c(0.7113820598, 
0.8391023102, 0.3185943152, 0.6887167806, 0.9120531859, 0.6254200542, 
0.7165302682, 0.7257435312, 0.2586613298, 0.6493567251), 
GC3 = c(0.615502, 0.622543, 0.393064, 0.490141, 0.461592, 
0.626407, 0.490305, 0.482853, 0.346424, 0.466484)), .Names = c("gene", 
 "Species", "labs", "chrMark", "chr", "dN", "dS", "dNdS", "RPKM", 
  "Tau", "GC3"), row.names = c(NA, 10L), class = "data.frame")

Answer 1

您的代码存在一些问题以及ggplot如何处理非标准评估，我建议您阅读ggplot教程或文档。在名为summarised_data和'mammals$species'的{​​{1}}内调用列会导致很多问题。

如果我们将这些变为更明智的......

'mammals$chrMark'

我们可以使ggplot代码更友好。请注意，当您将其映射到数据时，形状必须在names(summarised_data)[names(summarised_data) == "mammals$species"] <- "mammals_species" names(summarised_data)[names(summarised_data) == "mammals$chrMark"] <- "mammals_chrMark" 范围内。

aes

希望这应该有用，或者至少让你更接近答案。