{"id":493,"date":"2023-06-30T02:04:49","date_gmt":"2023-06-30T02:04:49","guid":{"rendered":"https:\/\/yunshangtulv.com.cn\/?p=493"},"modified":"2023-06-30T02:04:49","modified_gmt":"2023-06-30T02:04:49","slug":"nc%e5%9b%be%e7%89%87%e5%a4%8d%e7%8e%b0%ef%bc%9a%e5%a4%9a%e7%bb%84%e5%ad%a6%e6%95%a3%e7%82%b9%e5%9b%be%ef%bc%88%e8%b1%a1%e9%99%90%ef%bc%89","status":"publish","type":"post","link":"https:\/\/yunshangtulv.com.cn\/?p=493","title":{"rendered":"NC\u56fe\u7247\u590d\u73b0\uff1a\u591a\u7ec4\u5b66\u6563\u70b9\u56fe\uff08\u8c61\u9650\uff09"},"content":{"rendered":"<p><span>\u91c7\u7528\u7684\u662fMASS\u5305\u4e2d\u7684mvrnorm\u51fd\u6570\u751f\u6210\u4e8c\u5143\u6b63\u6001\u5206\u5e03\u7684\u6570\u636e\uff0c\u6570\u636e\u7684\u5217\u540d\u5206\u522b\u4e3a\"mRNA_FC\"\u548c\"RPF_FC\"\u3002<\/span><\/p>\n<p>\u5b8c\u6574\u4ee3\u7801:<\/p>\n<pre class=\"corepress-code-pre\"><code># \u52a0\u8f7d\u5fc5\u8981\u7684\u5e93\r\nlibrary(MASS)\r\nlibrary(ggplot2)\r\nlibrary(cowplot)\r\n\r\n# \u751f\u6210\u6570\u636e\r\nset.seed(123) # \u8bbe\u7f6e\u968f\u673a\u79cd\u5b50\u4ee5\u83b7\u5f97\u53ef\u590d\u5236\u7684\u7ed3\u679c\r\ncovariance &lt;- matrix(c(1,0.8,0.8,1),nrow=2,byrow=TRUE)\r\ndata &lt;- mvrnorm(n=1000, mu=c(0,0), covariance)\r\ndata &lt;- as.data.frame(data)\r\ncolnames(data) &lt;- c(\"mRNA_FC\", \"RPF_FC\")\r\n\r\n# \u5206\u7ec4\u5904\u7406\r\ndata$group &lt;- ifelse(abs(data$mRNA_FC) &gt; log10(2) &amp; abs(data$RPF_FC) &gt; log10(2),\r\n                     \"mRNA+RPF_both\", \r\n                     ifelse(abs(data$mRNA_FC) &gt; log10(2),\r\n                            \"mRNA_only\", \r\n                            ifelse(abs(data$RPF_FC) &gt; log10(2),\r\n                                   \"RPF_only\", NA)))\r\n\r\n# \u521b\u5efa\u4e00\u4e2a\u51fd\u6570\u6765\u7ed8\u5236\u81ea\u5b9a\u4e49\u7684\u5750\u6807\u8f74\r\ndraw_axis_line &lt;- function(length_x, length_y, \r\n                           tick_step = NULL, lab_step = NULL) {\r\n  axis_x_begin &lt;- -1 * length_x\r\n  axis_x_end &lt;- length_x\r\n  axis_y_begin &lt;- -1 * length_y\r\n  axis_y_end &lt;- length_y\r\n  if (missing(tick_step))\r\n    tick_step &lt;- 1\r\n  if (missing(lab_step))\r\n    lab_step &lt;- 2\r\n\r\n  # \u521b\u5efa\u5750\u6807\u8f74\u523b\u5ea6\u548c\u6807\u7b7e\u7684\u6570\u636e\u6846\r\n  tick_x_frame &lt;- data.frame(ticks = seq(axis_x_begin, axis_x_end, by = tick_step))\r\n  tick_y_frame &lt;- data.frame(ticks = seq(axis_y_begin, axis_y_end, by = tick_step))\r\n  lab_x_frame &lt;- subset(data.frame(lab = seq(axis_x_begin, axis_x_end, by = lab_step), zero = 0), lab != 0)\r\n  lab_y_frame &lt;- subset(data.frame(lab = seq(axis_y_begin, axis_y_end, by = lab_step), zero = 0), lab != 0)\r\n\r\n  tick_x_length = 0.05\r\n  tick_y_length = 0.05\r\n\r\n  # \u5f00\u59cb\u7ed8\u5236\u56fe\u5f62\r\n  p &lt;- ggplot() +\r\n    # \u7ed8\u5236\u5750\u6807\u8f74\r\n    geom_segment(aes(x = axis_x_begin, xend = axis_x_end, y = 0, yend = 0), size = 0.5) +\r\n    geom_segment(aes(x = 0, xend = 0, y = axis_y_begin, yend = axis_y_end), size = 0.5) +\r\n    # \u6dfb\u52a0\u523b\u5ea6\r\n    geom_segment(data = tick_x_frame, aes(x = ticks, xend = ticks, y = 0, yend = -tick_x_length)) +\r\n    geom_segment(data = tick_y_frame, aes(x = 0, xend = -tick_y_length, y = ticks, yend = ticks)) +\r\n    # \u6dfb\u52a0\u523b\u5ea6\u6807\u7b7e\r\n    geom_text(data = lab_x_frame, aes(x = lab, y = zero, label = lab), vjust = 1.5) +\r\n    geom_text(data = lab_y_frame, aes(x = zero, y = lab, label = lab), hjust = 1.5) +\r\n    theme_minimal() +\r\n    theme(panel.grid = element_blank(),\r\n          axis.text = element_blank())\r\n  \r\n  return(p)\r\n}\r\n\r\n# \u751f\u6210\u4e5d\u5bab\u683c\u56fe\r\np &lt;- draw_axis_line(4, 4)\r\np1 &lt;- p + \r\n  geom_point(data = data, aes(mRNA_FC, RPF_FC, color = group))+\r\n  scale_color_manual(values = c(\"mRNA+RPF_both\" = \"#dd8653\", \"mRNA_only\" = \"#59a5d7\", \"RPF_only\" = \"#aa65a4\", \"#878787\")) +\r\n  labs(x = \"mRNA:FC(P42\/E15.5)\", y = \"RPF:FC(P42\/E15.5)\") +\r\n  theme(legend.position = \"bottom\") +\r\n  annotate(\"text\", label = \"bolditalic(Brain)\", parse = TRUE, x = -2, y = 2, size = 4, colour = \"black\") +\r\n  guides(color = guide_legend(title = \"\", ncol = 1, byrow = TRUE))\r\n\r\n# \u4fdd\u5b58\u56fe\u5f62\u4e3aPDF\r\nggsave(\"plot.pdf\", plot = p1, height = 6, width = 5)\r\n\r\n# \u62fc\u63a5\u591a\u4e2a\u56fe\u5f62\r\np_list &lt;- list(p1=p1, p2=p1, p3=p1)\r\np_all &lt;- plot_grid(plotlist = p_list, ncol = 3)\r\n\r\n# \u4fdd\u5b58\u62fc\u63a5\u540e\u7684\u56fe\u5f62\u4e3aPDF\r\nggsave(\"plot2.pdf\", plot = p_all, height = 4, width = 9)\r\n\r\n# \u4fdd\u5b58\u6570\u636e\u4e3aCSV\u683c\u5f0f\r\nwrite.csv(data, file = \"data.csv\", row.names = FALSE)\r\n<\/code><\/pre>\n<p><span>\u5bf9\u4e8e\u4f7f\u7528\u81ea\u5df1\u7684\u6570\u636e\uff0c\u6570\u636e\u5e94\u8be5\u662f\u4e00\u4e2acsv\u6587\u4ef6\uff0c\u5305\u542b\u4e24\u5217\uff0c\u5217\u540d\u5206\u522b\u4e3a\"mRNA_FC\"\u548c\"RPF_FC\"\u3002\u5176\u4e2d\u6bcf\u4e00\u884c\u4ee3\u8868\u4e00\u4e2a\u89c2\u6d4b\uff0c\"mRNA_FC\"\u548c\"RPF_FC\"\u662f\u8fd9\u4e2a\u89c2\u6d4b\u7684\u4e24\u4e2a\u7279\u5f81\u503c\u3002<\/span><\/p>\n<p><img loading=\"lazy\" decoding=\"async\" src=\"https:\/\/yunshangtulv.com.cn\/wp-content\/uploads\/2023\/06\/\u622a\u56fe_20230630095924-847x1024.png\" alt=\"\" width=\"847\" height=\"1024\" class=\"alignnone size-large wp-image-494\" srcset=\"https:\/\/yunshangtulv.com.cn\/wp-content\/uploads\/2023\/06\/\u622a\u56fe_20230630095924-847x1024.png 847w, https:\/\/yunshangtulv.com.cn\/wp-content\/uploads\/2023\/06\/\u622a\u56fe_20230630095924-600x725.png 600w, https:\/\/yunshangtulv.com.cn\/wp-content\/uploads\/2023\/06\/\u622a\u56fe_20230630095924-248x300.png 248w, https:\/\/yunshangtulv.com.cn\/wp-content\/uploads\/2023\/06\/\u622a\u56fe_20230630095924-768x928.png 768w, https:\/\/yunshangtulv.com.cn\/wp-content\/uploads\/2023\/06\/\u622a\u56fe_20230630095924.png 1074w\" sizes=\"auto, (max-width: 847px) 100vw, 847px\" \/><\/p>\n<p><img loading=\"lazy\" decoding=\"async\" src=\"https:\/\/yunshangtulv.com.cn\/wp-content\/uploads\/2023\/06\/\u622a\u56fe_20230630095944-1024x453.png\" alt=\"\" width=\"1024\" height=\"453\" class=\"alignnone size-large wp-image-495\" srcset=\"https:\/\/yunshangtulv.com.cn\/wp-content\/uploads\/2023\/06\/\u622a\u56fe_20230630095944-1024x453.png 1024w, https:\/\/yunshangtulv.com.cn\/wp-content\/uploads\/2023\/06\/\u622a\u56fe_20230630095944-600x266.png 600w, https:\/\/yunshangtulv.com.cn\/wp-content\/uploads\/2023\/06\/\u622a\u56fe_20230630095944-300x133.png 300w, https:\/\/yunshangtulv.com.cn\/wp-content\/uploads\/2023\/06\/\u622a\u56fe_20230630095944-768x340.png 768w, https:\/\/yunshangtulv.com.cn\/wp-content\/uploads\/2023\/06\/\u622a\u56fe_20230630095944-1536x680.png 1536w, https:\/\/yunshangtulv.com.cn\/wp-content\/uploads\/2023\/06\/\u622a\u56fe_20230630095944.png 1548w\" sizes=\"auto, (max-width: 1024px) 100vw, 1024px\" \/><\/p>\n","protected":false},"excerpt":{"rendered":"<p>\u91c7\u7528\u7684\u662fMASS\u5305\u4e2d\u7684mvrnorm\u51fd\u6570\u751f\u6210\u4e8c\u5143\u6b63\u6001\u5206\u5e03\u7684\u6570\u636e\uff0c\u6570\u636e\u7684\u5217\u540d\u5206\u522b\u4e3a\"mRNA_FC\"\u548c\"RPF_ [&hellip;]<\/p>\n","protected":false},"author":111,"featured_media":0,"comment_status":"open","ping_status":"open","sticky":false,"template":"","format":"standard","meta":{"footnotes":""},"categories":[9],"tags":[],"class_list":["post-493","post","type-post","status-publish","format-standard","hentry","category-r"],"_links":{"self":[{"href":"https:\/\/yunshangtulv.com.cn\/index.php?rest_route=\/wp\/v2\/posts\/493","targetHints":{"allow":["GET"]}}],"collection":[{"href":"https:\/\/yunshangtulv.com.cn\/index.php?rest_route=\/wp\/v2\/posts"}],"about":[{"href":"https:\/\/yunshangtulv.com.cn\/index.php?rest_route=\/wp\/v2\/types\/post"}],"author":[{"embeddable":true,"href":"https:\/\/yunshangtulv.com.cn\/index.php?rest_route=\/wp\/v2\/users\/111"}],"replies":[{"embeddable":true,"href":"https:\/\/yunshangtulv.com.cn\/index.php?rest_route=%2Fwp%2Fv2%2Fcomments&post=493"}],"version-history":[{"count":0,"href":"https:\/\/yunshangtulv.com.cn\/index.php?rest_route=\/wp\/v2\/posts\/493\/revisions"}],"wp:attachment":[{"href":"https:\/\/yunshangtulv.com.cn\/index.php?rest_route=%2Fwp%2Fv2%2Fmedia&parent=493"}],"wp:term":[{"taxonomy":"category","embeddable":true,"href":"https:\/\/yunshangtulv.com.cn\/index.php?rest_route=%2Fwp%2Fv2%2Fcategories&post=493"},{"taxonomy":"post_tag","embeddable":true,"href":"https:\/\/yunshangtulv.com.cn\/index.php?rest_route=%2Fwp%2Fv2%2Ftags&post=493"}],"curies":[{"name":"wp","href":"https:\/\/api.w.org\/{rel}","templated":true}]}}