{"id":452,"date":"2023-04-27T01:32:49","date_gmt":"2023-04-27T01:32:49","guid":{"rendered":"https:\/\/yunshangtulv.com.cn\/?p=452"},"modified":"2023-04-27T01:32:49","modified_gmt":"2023-04-27T01:32:49","slug":"sci%e5%9b%be%e7%89%87%e5%a4%8d%e7%8e%b0%ef%bc%9a%e6%b8%90%e5%8f%98%e7%81%ab%e5%b1%b1%e5%9b%be","status":"publish","type":"post","link":"https:\/\/yunshangtulv.com.cn\/?p=452","title":{"rendered":"SCI\u56fe\u7247\u590d\u73b0\uff1a\u6e10\u53d8\u706b\u5c71\u56fe"},"content":{"rendered":"<pre class=\"corepress-code-pre\"><code># \u52a0\u8f7d\u5305\uff1a\r\nlibrary(ggplot2)\r\n\r\n# \u521b\u5efa\u793a\u4f8b\u6570\u636e\u96c6\uff1a\r\nset.seed(42) # \u4e3a\u4e86\u4fdd\u6301\u968f\u673a\u751f\u6210\u7684\u7ed3\u679c\u4e00\u81f4\u6027\r\nsample_size &lt;- 5000\r\ndata &lt;- data.frame(log2FoldChange = rnorm(sample_size, 0, 1.5),\r\n                   padj = runif(sample_size, 0, 1))\r\n# \u663e\u793a\u6570\u636e\u7684\u524d\u51e0\u884c\r\nhead(data)\r\n# \u4fee\u6539\u53d8\u91cf\u540d\u79f0\uff1a\r\ndata &lt;- data.frame(Log2FoldChange = data$log2FoldChange,\r\n                   PAdj = data$padj)\r\n\r\n# \u65b0\u589e\u4e00\u5217\u7528\u4e8e\u5b58\u50a8label\u4fe1\u606f\uff0c\u5c06\u9700\u8981\u663e\u793a\u7684label\u5217\u51fa\u5373\u53ef\uff1a\r\ndata$Label &lt;- c(1:10, rep(NA, (nrow(data) - 10)))\r\n\r\n# \u521b\u5efa\u6563\u70b9\u56fe\uff1a\r\nplot &lt;- ggplot(data, aes(Log2FoldChange, -log10(PAdj))) +\r\n  # \u6a2a\u5411\u6c34\u5e73\u53c2\u8003\u7ebf\uff1a\r\n  geom_hline(yintercept = -log10(0.05), linetype = \"dashed\", color = \"#999999\") +\r\n  # \u7eb5\u5411\u5782\u76f4\u53c2\u8003\u7ebf\uff1a\r\n  geom_vline(xintercept = c(-1.2, 1.2), linetype = \"dashed\", color = \"#999999\") +\r\n  # \u6563\u70b9\u56fe:\r\n  geom_point(aes(size = -log10(PAdj), color = -log10(PAdj))) +\r\n  # \u6307\u5b9a\u989c\u8272\u6e10\u53d8\u6a21\u5f0f\uff1a\r\n  scale_color_gradientn(values = seq(0, 1, 0.2),\r\n                         colors = c(\"#39489f\", \"#39bbec\", \"#f9ed36\", \"#f38466\", \"#b81f25\")) +\r\n  # \u6307\u5b9a\u6563\u70b9\u5927\u5c0f\u6e10\u53d8\u6a21\u5f0f\uff1a\r\n  scale_size_continuous(range = c(1, 3)) +\r\n  # \u4e3b\u9898\u8c03\u6574\uff1a\r\n  theme_bw() +\r\n  theme(panel.grid = element_blank())\r\n\r\n# \u663e\u793a\u7ed8\u5236\u7684\u56fe\u50cf\r\nprint(plot)\r\n\r\n\r\n<\/code><\/pre>\n<p>&nbsp;<\/p>\n<p>&nbsp;<\/p>\n<p><strong><img decoding=\"async\" src=\"https:\/\/yunshangtulv.com.cn\/wp-content\/uploads\/2023\/04\/1682558689474.webp\" \/><\/strong><\/p>\n<h3>\u4ee3\u7801\u89e3\u91ca\uff1a<\/h3>\n<ol>\n<li>\u52a0\u8f7d<code>ggplot2<\/code>\u5305\u3002<\/li>\n<li>\u521b\u5efa\u793a\u4f8b\u6570\u636e\u96c6\uff0c\u8bbe\u7f6e\u968f\u673a\u6570\u79cd\u5b50\u4ee5\u4fdd\u6301\u7ed3\u679c\u7684\u4e00\u81f4\u6027\u3002<\/li>\n<li>\u4fee\u6539\u6570\u636e\u96c6\u4e2d\u7684\u53d8\u91cf\u540d\u79f0\u3002<\/li>\n<li>\u4e3a\u6570\u636e\u96c6\u6dfb\u52a0\u4e00\u4e2a\u65b0\u5217<code>Label<\/code>\uff0c\u524d10\u884c\u5305\u542b1\u523010\u7684\u503c\uff0c\u5176\u4f59\u884c\u5305\u542bNA\u503c\u3002<\/li>\n<li>\u4f7f\u7528<code>ggplot()<\/code>\u51fd\u6570\u521b\u5efa\u6563\u70b9\u56fe\uff0c\u8bbe\u7f6eX\u8f74\u4e3a<code>Log2FoldChange<\/code>\uff0cY\u8f74\u4e3a<code>-log10(PAdj)<\/code>\u3002<\/li>\n<li>\u6dfb\u52a0\u6a2a\u5411\u6c34\u5e73\u53c2\u8003\u7ebf\uff0c\u4f4d\u7f6e\u4e3a<code>-log10(0.05)<\/code>\u3002<\/li>\n<li>\u6dfb\u52a0\u7eb5\u5411\u5782\u76f4\u53c2\u8003\u7ebf\uff0c\u4f4d\u7f6e\u4e3a<code>-1.2<\/code>\u548c<code>1.2<\/code>\u3002<\/li>\n<li>\u6dfb\u52a0\u6563\u70b9\u56fe\uff0c\u8bbe\u7f6e\u6563\u70b9\u7684\u5927\u5c0f\u548c\u989c\u8272\u4e0e<code>-log10(PAdj)<\/code>\u6210\u6bd4\u4f8b\u3002<\/li>\n<li>\u6307\u5b9a\u989c\u8272\u6e10\u53d8\u6a21\u5f0f\u3002<\/li>\n<li>\u6307\u5b9a\u6563\u70b9\u5927\u5c0f\u6e10\u53d8\u6a21\u5f0f\u3002<\/li>\n<li>\u8c03\u6574\u56fe\u50cf\u4e3b\u9898\u3002<\/li>\n<li>\u6253\u5370\u7ed8\u5236<\/li>\n<\/ol>\n<h3>\u8c03\u6574\u4f18\u5316\uff1a<\/h3>\n<p>\u5728\u706b\u5c71\u56fe\u4e0a\u6807\u6ce8\u6709\u610f\u4e49\u7684\u70b9\uff0c\u6211\u4eec\u53ef\u4ee5\u4f7f\u7528<code>geom_text()<\/code>\u51fd\u6570\u3002\u6211\u4eec\u53ef\u4ee5\u8bbe\u7f6e\u9608\u503c\uff0c\u4f8b\u5982<code>Log2FoldChange<\/code>\u7684\u7edd\u5bf9\u503c\u5927\u4e8e\u67d0\u4e2a\u503c\uff08\u4f8b\u59821.5\uff09\u4e14<code>PAdj<\/code>\u503c\u5c0f\u4e8e0.05\uff0c\u4f5c\u4e3a\u663e\u8457\u7684\u70b9\u3002\u4e0b\u9762\u7684\u4ee3\u7801\u6f14\u793a\u4e86\u5982\u4f55\u5728\u706b\u5c71\u56fe\u4e0a\u6807\u6ce8\u663e\u8457\u7684\u70b9\uff1a<\/p>\n<div class=\"bg-black rounded-md mb-4\">\n<div class=\"flex items-center relative text-gray-200 bg-gray-800 px-4 py-2 text-xs font-sans justify-between rounded-t-md\"><\/div>\n<pre class=\"corepress-code-pre\"><code># \u521b\u5efa\u6563\u70b9\u56fe\uff1a\r\nplot &lt;- ggplot(data, aes(Log2FoldChange, -log10(PAdj))) +\r\n  # \u6a2a\u5411\u6c34\u5e73\u53c2\u8003\u7ebf\uff1a\r\n  geom_hline(yintercept = -log10(0.05), linetype = \"dashed\", color = \"#999999\") +\r\n  # \u7eb5\u5411\u5782\u76f4\u53c2\u8003\u7ebf\uff1a\r\n  geom_vline(xintercept = c(-1.5, 1.5), linetype = \"dashed\", color = \"#999999\") +\r\n  # \u6563\u70b9\u56fe:\r\n  geom_point(aes(size = -log10(PAdj), color = -log10(PAdj))) +\r\n  # \u6307\u5b9a\u989c\u8272\u6e10\u53d8\u6a21\u5f0f\uff1a\r\n  scale_color_gradientn(values = seq(0, 1, 0.2),\r\n                         colors = c(\"#39489f\", \"#39bbec\", \"#f9ed36\", \"#f38466\", \"#b81f25\")) +\r\n  # \u6307\u5b9a\u6563\u70b9\u5927\u5c0f\u6e10\u53d8\u6a21\u5f0f\uff1a\r\n  scale_size_continuous(range = c(1, 3)) +\r\n  # \u4e3b\u9898\u8c03\u6574\uff1a\r\n  theme_bw() +\r\n  theme(panel.grid = element_blank()) +\r\n  # \u5728\u663e\u8457\u7684\u70b9\u4e0a\u6dfb\u52a0\u6807\u7b7e\uff1a\r\n  geom_text(data = subset(data, (abs(Log2FoldChange) &gt; 1.5) &amp; (PAdj &lt; 0.05)),\r\n            aes(label = Label), size = 3, hjust = -0.5, vjust = -0.5)\r\n\r\n# \u663e\u793a\u7ed8\u5236\u7684\u56fe\u50cf\r\nprint(plot)\r\n<\/code><\/pre>\n<\/div>\n<p>\u5728\u8fd9\u4e2a\u4ee3\u7801\u4e2d\uff0c\u6211\u4eec\u6dfb\u52a0\u4e86\u4e00\u4e2a<code>geom_text()<\/code>\u51fd\u6570\uff0c\u5e76\u6839\u636e\u7ed9\u5b9a\u6761\u4ef6\u7b5b\u9009\u6570\u636e\uff08<code>Log2FoldChange<\/code>\u7684\u7edd\u5bf9\u503c\u5927\u4e8e1.5\u4e14<code>PAdj<\/code>\u503c\u5c0f\u4e8e0.05\uff09\u3002\u6211\u4eec\u8fd8\u8bbe\u7f6e\u4e86\u6807\u7b7e\u7684\u5927\u5c0f\u3001\u6c34\u5e73\u548c\u5782\u76f4\u5bf9\u9f50\u53c2\u6570\uff0c\u4ee5\u4fbf\u5728\u663e\u8457\u70b9\u65c1\u8fb9\u8fdb\u884c\u6807\u6ce8\u3002<\/p>\n<p>&nbsp;<\/p>\n","protected":false},"excerpt":{"rendered":"<p># \u52a0\u8f7d\u5305\uff1a library(ggplot2) # \u521b\u5efa\u793a\u4f8b\u6570\u636e\u96c6\uff1a set.seed(42) # \u4e3a\u4e86\u4fdd\u6301\u968f [&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-452","post","type-post","status-publish","format-standard","hentry","category-r"],"_links":{"self":[{"href":"https:\/\/yunshangtulv.com.cn\/index.php?rest_route=\/wp\/v2\/posts\/452","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=452"}],"version-history":[{"count":0,"href":"https:\/\/yunshangtulv.com.cn\/index.php?rest_route=\/wp\/v2\/posts\/452\/revisions"}],"wp:attachment":[{"href":"https:\/\/yunshangtulv.com.cn\/index.php?rest_route=%2Fwp%2Fv2%2Fmedia&parent=452"}],"wp:term":[{"taxonomy":"category","embeddable":true,"href":"https:\/\/yunshangtulv.com.cn\/index.php?rest_route=%2Fwp%2Fv2%2Fcategories&post=452"},{"taxonomy":"post_tag","embeddable":true,"href":"https:\/\/yunshangtulv.com.cn\/index.php?rest_route=%2Fwp%2Fv2%2Ftags&post=452"}],"curies":[{"name":"wp","href":"https:\/\/api.w.org\/{rel}","templated":true}]}}