{"id":28899,"date":"2024-02-19T18:03:28","date_gmt":"2024-02-19T10:03:28","guid":{"rendered":"http:\/\/www.biocloudservice.com\/wordpress\/?p=28899"},"modified":"2024-02-19T18:03:30","modified_gmt":"2024-02-19T10:03:30","slug":"%e5%88%a9%e7%94%a8ggplot2%e7%bb%98%e5%88%b6%e5%a4%9a%e8%a7%92%e5%9b%be%e7%bb%98%e5%88%b6","status":"publish","type":"post","link":"http:\/\/www.biocloudservice.com\/wordpress\/?p=28899","title":{"rendered":"\u5229\u7528ggplot2\u7ed8\u5236\u591a\u89d2\u56fe\u7ed8\u5236"},"content":{"rendered":"

\u4eca\u5929\u5229\u7528ggplot2\u6765\u8fdb\u884c\u591a\u89d2\u56fe\u7684\u7ed8\u5236\uff0c\u4ee3\u7801\u5982\u4e0b\uff1a<\/p>\n

    \n
  1. \u5b89\u88c5\u9700\u8981\u7684R\u5305<\/li>\n<\/ol>\n

    install.packages(\u201cggplot2\u201d)<\/p>\n

      \n
    1. \u5bfc\u5165\u9700\u8981\u7684R\u5305<\/li>\n<\/ol>\n

      library(ggplot2)<\/p>\n

        \n
      1. \u793a\u4f8b\u4ee3\u7801<\/li>\n<\/ol>\n

        raw.data <- read.table(\u201cgroup.txt\u201d, header = T, stringsAsFactors = T)<\/p>\n

        \"Dingtalk_20230213122803\"<\/p>\n

        # \u8bb0\u5f55\u672b\u5217\u7684\u540d\u79f0<\/p>\n

        txt <- colnames(raw.data)[ncol(raw.data)]<\/p>\n

        ## \u8ba1\u7b97\u5750\u6807\u8f74\u89d2\u5ea6<\/p>\n

        sides <- ncol(raw.data) – 1 # \u8fb9\u7684\u4e2a\u6570\uff08\u5217\u6570\u51cf1\uff0c\u56e0\u4e3a\u6709\u672b\u5217\uff09<\/p>\n

        start_angle <- 90 # \u5f00\u59cb\u89d2\u5ea6<\/p>\n

        basis <- t(sapply((1:sides)-1, function(i){<\/p>\n

        c(cos(2 * i * pi \/ sides + start_angle * pi \/ 180),<\/p>\n

        sin(2 * i * pi \/ sides + start_angle * pi \/ 180))<\/p>\n

        }))<\/p>\n

        ## \u5c0612\u7ef4\u7684\u6570\u636e\u5750\u6807\u6295\u5c04\u5230\u4e8c\u7ef4<\/p>\n

        data <- as.matrix(raw.data[, -match(txt, colnames(raw.data))]) %*% basis<\/p>\n

        # \u51c6\u5907\u7ed8\u56fe\u6570\u636e<\/p>\n

        ### \u70b9\u7684\u6570\u636e<\/p>\n

        plot.data <- data.frame(<\/p>\n

        “x” = data[, 1],<\/p>\n

        “y” = data[, 2],<\/p>\n

        “Group” = raw.data$celltype<\/p>\n

        )<\/p>\n

        ### \u591a\u8fb9\u5f62\u7684\u6570\u636e<\/p>\n

        basis <- data.frame(<\/p>\n

        “x” = basis[, 1],<\/p>\n

        “y” = basis[, 2]<\/p>\n

        )<\/p>\n

        ### \u7c7b\u56fe\u4f8b\u7684\u6570\u636e<\/p>\n

        label <- data.frame(<\/p>\n

        “x” = 1.1 * basis[, 1],<\/p>\n

        “y” = 1.1 * basis[, 2],<\/p>\n

        “label” = colnames(raw.data)[-ncol(raw.data)]<\/p>\n

        )<\/p>\n

        ### \u8bbe\u5b9a\u989c\u8272\uff08\u4f7f\u7528\u7684\u662fRGB\u989c\u8272\uff09<\/p>\n

        color.palette <- list(<\/p>\n

        ‘hEndo’= c(190, 10, 10),’hPeric’= c(225, 160, 30),’hMgl’= c(217, 245, 7),<\/p>\n

        ‘hDA1’= c(170, 180, 170),’hDA2’= c(130, 140, 140),’hNbM’= c(180, 140, 130),<\/p>\n

        ‘hNbML1’= c(100, 100, 240),’hProgM’= c( 80, 235, 255),’hProgFPM’=c(190, 235, 255),<\/p>\n

        ‘hProgFPL’=c(210, 255, 215),’hProgBP’=c(230, 140, 120),’hNProg’= c(255, 195, 28),<\/p>\n

        ‘hNbML5’= c(139, 101, 100),’hRgl1’= c(252, 183, 26),’hRgl3’= c(214, 194, 39),<\/p>\n

        ‘hRgl2c’= c(255, 120, 155),’hRgl2b’= c(250, 145, 45),’hRgl2a’= c(250, 125, 25),<\/p>\n

        ‘hDA0’= c(190, 200, 190),’hOPC’= c(255, 35, 155),’hRN’= c(199, 121, 41),<\/p>\n

        ‘hNbGaba’= c( 40, 55, 130),’hGaba’= c( 7, 121, 61),’hOMTN’= c( 95, 186, 70),<\/p>\n

        ‘hSert’= c( 50, 180, 180),’eSCa’= c(245, 205, 170),’eSCb’= c(205, 245, 170),<\/p>\n

        ‘eSCc’= c(205, 205, 220)<\/p>\n

        )<\/p>\n

        color.palette <- lapply(color.palette, function(x){<\/p>\n

        rgb(red = x[1]\/255, green = x[2]\/255, blue = x[3]\/255)<\/p>\n

        })<\/p>\n

        color.palette <- unlist(color.palette) # \u6700\u7ec8\u7684\u989c\u8272\u683c\u5f0f\u4e3a\u4e00\u4e2a\u5e26\u540d\u79f0\u768416\u8fdb\u5236\u5b57\u7b26\u4e32<\/p>\n

        ###\u5f00\u59cb\u7ed8\u56fe<\/p>\n

        ggplot(plot.data, aes(x = x, y = y)) +<\/p>\n

        geom_point(aes(col = Group, size = .5)) +<\/p>\n

        geom_polygon(mapping = aes(x = x, y = y), data = basis,<\/p>\n

        fill = “transparent”, color = “black”, size = 1) +<\/p>\n

        geom_text(mapping = aes(x = x, y = y, label = label), data = label) +<\/p>\n

        scale_color_manual(values = color.palette) +<\/p>\n

        theme_classic() +<\/p>\n

        theme(axis.line = element_blank(),<\/p>\n

        axis.title = element_blank(),<\/p>\n

        axis.ticks = element_blank(),<\/p>\n

        axis.text = element_blank())<\/p>\n

        ggsave(“wheelplot.pdf”, width = 10, height = 8) # \u4fdd\u5b58\u56fe\u5f62<\/p>\n

        \"Dingtalk_20230213122417\"<\/p>\n","protected":false},"excerpt":{"rendered":"

        \u4eca\u5929\u5229\u7528ggplot2\u6765\u8fdb\u884c\u591a\u89d2\u56fe\u7684\u7ed8\u5236\uff0c\u4ee3\u7801\u5982\u4e0b\uff1a \u5b89\u88c5\u9700\u8981\u7684R\u5305 install.packages(\u201cgg […]<\/p>\n","protected":false},"author":1,"featured_media":0,"comment_status":"closed","ping_status":"closed","sticky":false,"template":"","format":"standard","meta":{"jetpack_post_was_ever_published":false,"_jetpack_newsletter_access":""},"categories":[1],"tags":[],"jetpack_featured_media_url":"","_links":{"self":[{"href":"http:\/\/www.biocloudservice.com\/wordpress\/index.php?rest_route=\/wp\/v2\/posts\/28899"}],"collection":[{"href":"http:\/\/www.biocloudservice.com\/wordpress\/index.php?rest_route=\/wp\/v2\/posts"}],"about":[{"href":"http:\/\/www.biocloudservice.com\/wordpress\/index.php?rest_route=\/wp\/v2\/types\/post"}],"author":[{"embeddable":true,"href":"http:\/\/www.biocloudservice.com\/wordpress\/index.php?rest_route=\/wp\/v2\/users\/1"}],"replies":[{"embeddable":true,"href":"http:\/\/www.biocloudservice.com\/wordpress\/index.php?rest_route=%2Fwp%2Fv2%2Fcomments&post=28899"}],"version-history":[{"count":1,"href":"http:\/\/www.biocloudservice.com\/wordpress\/index.php?rest_route=\/wp\/v2\/posts\/28899\/revisions"}],"predecessor-version":[{"id":28902,"href":"http:\/\/www.biocloudservice.com\/wordpress\/index.php?rest_route=\/wp\/v2\/posts\/28899\/revisions\/28902"}],"wp:attachment":[{"href":"http:\/\/www.biocloudservice.com\/wordpress\/index.php?rest_route=%2Fwp%2Fv2%2Fmedia&parent=28899"}],"wp:term":[{"taxonomy":"category","embeddable":true,"href":"http:\/\/www.biocloudservice.com\/wordpress\/index.php?rest_route=%2Fwp%2Fv2%2Fcategories&post=28899"},{"taxonomy":"post_tag","embeddable":true,"href":"http:\/\/www.biocloudservice.com\/wordpress\/index.php?rest_route=%2Fwp%2Fv2%2Ftags&post=28899"}],"curies":[{"name":"wp","href":"https:\/\/api.w.org\/{rel}","templated":true}]}}