{"id":61543,"date":"2024-10-30T15:35:27","date_gmt":"2024-10-30T07:35:27","guid":{"rendered":"http:\/\/www.biocloudservice.com\/wordpress\/?p=61543"},"modified":"2024-10-30T15:35:27","modified_gmt":"2024-10-30T07:35:27","slug":"%e5%8d%95%e7%bb%86%e8%83%9e%e7%b1%bb%e5%9e%8b%e9%89%b4%e5%ae%9a%e5%92%8c%e8%a7%a3%e9%87%8a%e7%9a%84%e5%87%86%e7%a1%ae%e6%80%a7%e6%ac%a0%e4%bd%b3%ef%bc%9fcascc%e8%87%aa%e9%80%82%e5%ba%94%e5%90%b8","status":"publish","type":"post","link":"http:\/\/www.biocloudservice.com\/wordpress\/?p=61543","title":{"rendered":"\u5355\u7ec6\u80de\u7c7b\u578b\u9274\u5b9a\u548c\u89e3\u91ca\u7684\u51c6\u786e\u6027\u6b20\u4f73\uff1fCASCC\u81ea\u9002\u5e94\u5438\u5f15\u5b50\u7b97\u6cd5\u7ed9\u4f60\u63d0\u4f9b\u989d\u5916\u89c6\u89d2\uff01"},"content":{"rendered":"\n<p><\/p>\n\n\n\n<p><strong>\u524d\u8a00<\/strong><strong><\/strong><\/p>\n\n\n\n<p>\u5c0f\u679c\u89c1\u9762\u4f1a\u53c8\u8981\u5f00\u59cb\u4e86\uff0c\u867d\u7136\u5c0f\u679c\u5b9e\u9a8c\u5f88\u5fd9\uff0c\u4f46\u662f\u770b\u5230\u597d\u7528\u7684\u5355\u7ec6\u80de\u805a\u7c7bR\u5305\u8fd8\u662f\u5fcd\u4e0d\u4f4f\u7b2c\u4e00\u65f6\u95f4\u5411\u5927\u5bb6\u5206\u4eab\uff01\u60f3\u5fc5\u5927\u591a\u6570\u5c0f\u4f19\u4f34\u4eec\u7684\u7ec6\u80de\u7c7b\u578b\u6807\u8bb0\u7269\u7531\u805a\u7c7b\u540e\u8ba1\u7b97\u7684\u5dee\u5f02\u57fa\u56e0\u6765\u9274\u5b9a\u7684\u5427\uff01\u662f\u4e0d\u662f\u6709\u65f6\u5019\u5355\u7ec6\u80de\u7c7b\u578b\u9274\u5b9a\u548c\u89e3\u91ca\u7684\u51c6\u786e\u6027\u6b20\u4f73\uff1f\u4e0d\u59a8\u8ddf\u7740\u5c0f\u679c\u6765\u770b\u770bCASCC\u8fd9\u4e2a\u5171\u8868\u8fbe\u8f85\u52a9\u5355\u7ec6\u80de\u805a\u7c7b\u65b9\u6cd5\u5427\u3002CASCC\u7684\u4e00\u4e2a\u5173\u952e\u7279\u70b9\u662f\u4f7f\u7528\u4e86\u4e00\u79cd\u65b0\u7684\u81ea\u9002\u5e94\u5438\u5f15\u5b50\u7b97\u6cd5\uff0c\u8be5\u7b97\u6cd5\u6709\u52a9\u4e8e\u6b63\u786e\u7684\u7279\u5f81\u9009\u62e9\u3001\u805a\u7c7b\u548c\u7ec6\u80de\u7c7b\u578b\u9274\u5b9a\uff0c\u4ece\u800c\u4e0e\u73b0\u6709\u65b9\u6cd5\u76f8\u6bd4\u63d0\u9ad8\u4e86\u751f\u7269\u5b66\u51c6\u786e\u6027\u3002CASCC\u4f18\u4e8e\u5176\u4ed6\u65b9\u6cd5\uff0c\u8fd9\u5df2\u7ecf\u591a\u4e2a\u8bc4\u4f30\u6307\u6807\u5f97\u5230\u8bc1\u660e\u5566\uff01\u6240\u4ee5\u5c0f\u4f19\u4f34\u4eec\u653e\u5fc3\u4f7f\u7528\u54e6\uff01\u6ca1\u6709\u65f6\u95f4\u5b66\u4e60\u7684\u5c0f\u4f19\u4f34\u4eec\u4e5f\u4e0d\u8981\u7740\u6025\u54e6\uff01\u6709\u9700\u8981\u751f\u4fe1\u5206\u6790\u7684\u5c0f\u4f19\u4f34\u4eec\u4e5f\u53ef\u4ee5\u627e\u5c0f\u679c\u54e6\uff01\u7ec3\u4e86\u5341\u5e74\u751f\u4fe1\u5206\u6790\u7684\u5c0f\u679c\u5bf9\u4e8e\u751f\u4fe1\u5206\u6790\u77e5\u8bc6\u5df2\u7ecf\u5982\u9c7c\u5f97\u6c34\u4ece\u5206\u6790\u5230\u53ef\u89c6\u5316\u76f4\u5230\u4f60\u6ee1\u610f\u4e3a\u6b62\uff01<\/p>\n\n\n\n<p>\u5355\u7ec6\u80de\u6570\u636e\u5904\u7406\u8d77\u6765\u5360\u7528\u5185\u5b58\u5b9e\u5728\u592a\u5927\u4e86\uff0c\u653e\u8fc7\u81ea\u5df1\u7684\u7535\u8111\u5427\uff01\u5c0f\u679c\u5728\u8fd9\u91cc\u7ed9\u5927\u5bb6\u9001\u4e0a\u798f\u5229\u4e86\uff0c\u6709\u9700\u8981\u670d\u52a1\u5668\u7684\u5c0f\u4f19\u4f34\u4eec\uff0c\u6b22\u8fce\u5927\u5bb6\u8054\u7cfb\u5c0f\u679c\uff0c\u4fdd\u8bc1\u670d\u52a1\u5668\u7684\u6027\u4ef7\u6bd4\u6700\u9ad8\u54e6\uff01<\/p>\n\n\n\n<p><strong>\u5b9e\u6218\u6559\u7a0b<\/strong><strong><\/strong><\/p>\n\n\n\n<p><strong>1.<\/strong><strong>\u5b89\u88c5CASCC<\/strong><strong><\/strong><\/p>\n\n\n\n<p>1.1\u4f7f\u7528anaconda\u5b89\u88c5R\u548cCASCC<\/p>\n\n\n\n<p>conda create -n CASCC_env<\/p>\n\n\n\n<p>source activate CASCC_env<\/p>\n\n\n\n<p>conda install conda-forge::r-fs #&nbsp;devtools\u662f\u5fc5\u5907\u7684<\/p>\n\n\n\n<p>conda install conda-forge::r-devtools<\/p>\n\n\n\n<p>R<\/p>\n\n\n\n<p>.libPaths(&#8220;.\/anaconda3\/envs\/CASCC_env\/lib\/R\/library&#8221;) # \u8bbe\u7f6econda R\u5e93\u8def\u5f84<\/p>\n\n\n\n<p>install.packages(&#8220;rlang&#8221;,version=&#8221;1.1.0&#8243;) # \u4f9d\u8d56\u5305:\u66f4\u65b0rlang\u5305<\/p>\n\n\n\n<p>remotes::install_version(&#8220;Seurat&#8221;, &#8220;4.3.0&#8221;) # \u4f9d\u8d56\u5305<\/p>\n\n\n\n<p>devtools::install_github(&#8220;weiyi-bitw\/cafr&#8221;) # \u4f9d\u8d56\u5305<\/p>\n\n\n\n<p>devtools::install_github(&#8220;LingyiC\/CASCC&#8221;)<\/p>\n\n\n\n<p>CASCC::checkDependencies()<\/p>\n\n\n\n<p>1.2\u9009\u98792:\u5728R\u4e2d\u5b89\u88c5<\/p>\n\n\n\n<p>require(devtools)<\/p>\n\n\n\n<p>remotes::install_version(&#8220;Seurat&#8221;, &#8220;4.3.0&#8221;) # dependency<\/p>\n\n\n\n<p>install_github(&#8220;weiyi-bitw\/cafr&#8221;) # dependency<\/p>\n\n\n\n<p>install_github(&#8220;LingyiC\/CASCC&#8221;)<\/p>\n\n\n\n<p>CASCC::checkDependencies()<\/p>\n\n\n\n<p>2.&nbsp;CASCC\u4f7f\u7528\u6559\u7a0b<\/p>\n\n\n\n<p>\u5728\u672c\u6559\u7a0b\u4e2d\uff0c\u6211\u4eec\u5c06\u5206\u6790\u6765\u81ea(Peng et al. 2019)\u7814\u7a76\u7684\u80f0\u817a\u5bfc\u7ba1\u817a\u764c(PDAC)\u60a3\u8005\u7684\u793a\u4f8b\u6570\u636e\u96c6\u3002\u8fd9\u4e2a\u4f8b\u5b50\u5305\u62ec8\u79cd\u4e0d\u540c\u7684\u7ec6\u80de\u7c7b\u578b\u3002<\/p>\n\n\n\n<p>2.1\u5feb\u901f\u5165<\/p>\n\n\n\n<p>2.1.1\u7528\u4e00\u884c\u4ee3\u7801\u8fd0\u884cCASCC<\/p>\n\n\n\n<p>\u8fd9\u91cc\u6211\u4eec\u4f7f\u7528main\u51fd\u6570run.CASCC\u6765\u8fd0\u884cCASCC\u3002\u9010\u6b65\u8fd0\u884cCASCC\u8bf7\u53c2\u89c1\u201c\u9010\u6b65\u8fd0\u884cCASCC\u201d\u4e00\u8282\u3002<\/p>\n\n\n\n<p>rm(list = ls()); gc()<\/p>\n\n\n\n<p>library(CASCC)<\/p>\n\n\n\n<p>library(Seurat)<\/p>\n\n\n\n<p>set.seed(1)<\/p>\n\n\n\n<p>## \u52a0\u8f7d\u793a\u4f8b\u6570\u636e<\/p>\n\n\n\n<p>data(&#8220;Data_PDAC_peng_2k&#8221;)<\/p>\n\n\n\n<p>dim(Data_PDAC_peng_2k) # 400 cells<\/p>\n\n\n\n<p>## \u8fd0\u884cCASCC<\/p>\n\n\n\n<p>CASCC::checkDependencies()<\/p>\n\n\n\n<p>##inputDataType &#8216; data &#8216;\u88ab\u89c4\u8303\u5316\u548c\u5bf9\u6570\u8f6c\u6362\u3002<\/p>\n\n\n\n<p>res &lt;- CASCC::run.CASCC(Data_PDAC_peng_2k)<\/p>\n\n\n\n<p>labels &lt;- res$mainType.output$clusteringResults<\/p>\n\n\n\n<p>2.1.2UMAP\u53ef\u89c6\u5316<\/p>\n\n\n\n<p>## \u7ed8\u56fe<\/p>\n\n\n\n<p>adata &lt;- res$res.predictK$adata # UMAP<\/p>\n\n\n\n<p>labels &lt;- res$mainType.output$clusteringResults<\/p>\n\n\n\n<p>&nbsp;# \u805a\u7c7b\u8f93\u51fa:\u5177\u6709\u5438\u5f15\u5b50\u4e2d\u5fc3\u7684k\u5747\u503c<\/p>\n\n\n\n<p>Idents(adata) &lt;- labels<\/p>\n\n\n\n<p>DimPlot(adata)<\/p>\n\n\n\n<figure class=\"wp-block-image size-full\"><img decoding=\"async\" loading=\"lazy\" width=\"494\" height=\"395\" src=\"https:\/\/i0.wp.com\/www.biocloudservice.com\/wordpress\/wp-content\/uploads\/2024\/10\/1730268771902_4B55C4FE-3FDF-4d0e-A49C-D34E0CF09012.png?resize=494%2C395\" alt=\"\" class=\"wp-image-61544\" srcset=\"https:\/\/i0.wp.com\/www.biocloudservice.com\/wordpress\/wp-content\/uploads\/2024\/10\/1730268771902_4B55C4FE-3FDF-4d0e-A49C-D34E0CF09012.png?w=494 494w, https:\/\/i0.wp.com\/www.biocloudservice.com\/wordpress\/wp-content\/uploads\/2024\/10\/1730268771902_4B55C4FE-3FDF-4d0e-A49C-D34E0CF09012.png?resize=300%2C240 300w\" sizes=\"(max-width: 494px) 100vw, 494px\" data-recalc-dims=\"1\" \/><\/figure>\n\n\n\n<p>2.1.3\u5438\u5f15\u5b50\u5171\u8868\u8fbe\u5f0f\u8f93\u51fa<\/p>\n\n\n\n<p>\u5206\u6790\u786e\u5b9a\u4e86\u516b\u4e2a\u4e0d\u540c\u96c6\u7fa4\u7684\u5b58\u5728\uff0c\u6bcf\u4e2a\u96c6\u7fa4\u5bf9\u5e94\u4e00\u4e2a\u7279\u5b9a\u7684\u5438\u5f15\u5b50\uff0c\u5982\u4e0b\u8868\u6240\u793a\u3002\u6211\u4eec\u5c55\u793a\u4e86\u6bcf\u4e2a\u5438\u5f15\u5b50\u7684\u524d10\u4e2a\u57fa\u56e0\uff0c\u5b83\u4eec\u4e0e\u5df2\u77e5\u7684\u7ec6\u80de\u7c7b\u578b\u6807\u8bb0\u5bc6\u5207\u76f8\u5173\u3002\u4ece\u5de6\u5230\u53f3\u4f9d\u6b21\u4e3a\u5de8\u566c\u7ec6\u80de\u30011\u578b\u5bfc\u7ba1\u7ec6\u80de\u3001B\u7ec6\u80de\u3001\u5185\u76ae\u7ec6\u80de\u30012\u578b\u5bfc\u7ba1\u7ec6\u80de\u3001T\u7ec6\u80de\u3001\u661f\u72b6\u7ec6\u80de\u548c\u6210\u7ea4\u7ef4\u7ec6\u80de\u3002<\/p>\n\n\n\n<p>## \u7ed3\u679c\u5217\u8868<\/p>\n\n\n\n<p>attr.raw &lt;- res$fs.res$attr.raw<\/p>\n\n\n\n<p># \u539f\u59cb\u7684\u5438\u5f15\u5668\u626b\u63cf\u7ed3\u679c\uff0c\u5982\u679c\u91cd\u65b0\u8fd0\u884cCASCC\u53ef\u4ee5\u4f7f\u7528\u3002<\/p>\n\n\n\n<p># res$fs.res$attrs # \u5220\u9664\u76f8\u540c\u7684\u5438\u5f15\u5b50<\/p>\n\n\n\n<p>finalAttrs &lt;- res$fs.res$finalAttrs # \u8fc7\u6ee4\u5438\u5f15\u5b50&nbsp;<\/p>\n\n\n\n<p>## \u5438\u5f15\u5b50\u7528\u4f5c\u805a\u7c7b\u4e2d\u5fc3<\/p>\n\n\n\n<p>attrs &lt;- res$mainType.output$attrs.center<\/p>\n\n\n\n<p>df &lt;- as.data.frame(do.call(cbind, lapply(attrs, function(x) {data.frame(&#8220;gene&#8221; = names(x))})))<\/p>\n\n\n\n<p>titles &lt;- c()<\/p>\n\n\n\n<p>for (attr.ID in paste0(&#8220;attr_&#8221;, unlist(lapply(attrs, function(x){names(head(x)[1])})))) {<\/p>\n\n\n\n<p>&nbsp;&nbsp;titles &lt;- c(titles, c(attr.ID))<\/p>\n\n\n\n<p>}<\/p>\n\n\n\n<p>colnames(df) &lt;- titles<\/p>\n\n\n\n<p>df[1:10, ] # \u6bcf\u4e2a\u5438\u5f15\u5b50\u7684\u524d10\u4e2a\u57fa\u56e0(\u884c)(\u5217)<\/p>\n\n\n\n<figure class=\"wp-block-image size-full\"><img decoding=\"async\" loading=\"lazy\" width=\"554\" height=\"268\" src=\"https:\/\/i0.wp.com\/www.biocloudservice.com\/wordpress\/wp-content\/uploads\/2024\/10\/1730268781998_7C2C8D5A-30A7-41da-B5AA-7739351F73E6.png?resize=554%2C268\" alt=\"\" class=\"wp-image-61545\" srcset=\"https:\/\/i0.wp.com\/www.biocloudservice.com\/wordpress\/wp-content\/uploads\/2024\/10\/1730268781998_7C2C8D5A-30A7-41da-B5AA-7739351F73E6.png?w=554 554w, https:\/\/i0.wp.com\/www.biocloudservice.com\/wordpress\/wp-content\/uploads\/2024\/10\/1730268781998_7C2C8D5A-30A7-41da-B5AA-7739351F73E6.png?resize=300%2C145 300w\" sizes=\"(max-width: 554px) 100vw, 554px\" data-recalc-dims=\"1\" \/><\/figure>\n\n\n\n<p>2.1.4 UMAP\uff0c\u7531\u5438\u5f15\u5b50\u6807\u8bb0<\/p>\n\n\n\n<p>\u5c06\u786e\u5b9a\u76848\u4e2a\u5438\u5f15\u5b50\u7528\u4f5c\u6bcf\u4e2a\u4e9a\u7c07\u7684\u4e2d\u5fc3\u3002<\/p>\n\n\n\n<p>Idents(adata) &lt;- res$mainType.output$clusteringResults<\/p>\n\n\n\n<p>nClust &lt;- length(unique(Idents(adata)))<\/p>\n\n\n\n<p>mycolor &lt;- grDevices::colorRampPalette(RColorBrewer::brewer.pal(n = 8, name = &#8220;Dark2&#8221;))(nClust)<\/p>\n\n\n\n<p>attractorDimplot(adata, attr.list = res$mainType.output$attrs.center, col = mycolor)<\/p>\n\n\n\n<figure class=\"wp-block-image size-full\"><img decoding=\"async\" loading=\"lazy\" width=\"552\" height=\"368\" src=\"https:\/\/i0.wp.com\/www.biocloudservice.com\/wordpress\/wp-content\/uploads\/2024\/10\/1730268792102_DC3F4DFF-9DDB-4eca-927C-109A65C46835.png?resize=552%2C368\" alt=\"\" class=\"wp-image-61546\" srcset=\"https:\/\/i0.wp.com\/www.biocloudservice.com\/wordpress\/wp-content\/uploads\/2024\/10\/1730268792102_DC3F4DFF-9DDB-4eca-927C-109A65C46835.png?w=552 552w, https:\/\/i0.wp.com\/www.biocloudservice.com\/wordpress\/wp-content\/uploads\/2024\/10\/1730268792102_DC3F4DFF-9DDB-4eca-927C-109A65C46835.png?resize=300%2C200 300w\" sizes=\"(max-width: 552px) 100vw, 552px\" data-recalc-dims=\"1\" \/><\/figure>\n\n\n\n<p>2.1.5 Attractor plot<\/p>\n\n\n\n<p>\u6bcf\u4e2a\u5438\u5f15\u5b50\u4e2d\u524d5\u4e2a\u57fa\u56e0\u7684\u5e73\u5747\u8868\u8fbe\u91cf\u3002<\/p>\n\n\n\n<p>attractorPlot(adata, attrs[1:4], n.col = 2)<\/p>\n\n\n\n<figure class=\"wp-block-image size-full\"><img decoding=\"async\" loading=\"lazy\" width=\"522\" height=\"392\" src=\"https:\/\/i0.wp.com\/www.biocloudservice.com\/wordpress\/wp-content\/uploads\/2024\/10\/1730268801390_31A2306E-4B6A-4f6e-81A1-4531BB875732.png?resize=522%2C392\" alt=\"\" class=\"wp-image-61547\" srcset=\"https:\/\/i0.wp.com\/www.biocloudservice.com\/wordpress\/wp-content\/uploads\/2024\/10\/1730268801390_31A2306E-4B6A-4f6e-81A1-4531BB875732.png?w=522 522w, https:\/\/i0.wp.com\/www.biocloudservice.com\/wordpress\/wp-content\/uploads\/2024\/10\/1730268801390_31A2306E-4B6A-4f6e-81A1-4531BB875732.png?resize=300%2C225 300w\" sizes=\"(max-width: 522px) 100vw, 522px\" data-recalc-dims=\"1\" \/><\/figure>\n\n\n\n<p>2.1.6\u5bfb\u627e\u5dee\u5f02\u8868\u8fbe\u57fa\u56e0<\/p>\n\n\n\n<p>\u6211\u4eec\u53d1\u73b0\u4e9a\u7fa4\u4e4b\u95f4\u5b58\u5728\u5dee\u5f02\u8868\u8fbe\u57fa\u56e0(DEGs)\u3002<\/p>\n\n\n\n<p>DEGs &lt;- Seurat::FindAllMarkers(adata, only.pos = T)<\/p>\n\n\n\n<p>head(DEGs<\/p>\n\n\n\n<figure class=\"wp-block-image size-full\"><img decoding=\"async\" loading=\"lazy\" width=\"554\" height=\"182\" src=\"https:\/\/i0.wp.com\/www.biocloudservice.com\/wordpress\/wp-content\/uploads\/2024\/10\/1730268813228_DDC46894-915B-42ed-9CEB-F19DFD67854F.png?resize=554%2C182\" alt=\"\" class=\"wp-image-61548\" srcset=\"https:\/\/i0.wp.com\/www.biocloudservice.com\/wordpress\/wp-content\/uploads\/2024\/10\/1730268813228_DDC46894-915B-42ed-9CEB-F19DFD67854F.png?w=554 554w, https:\/\/i0.wp.com\/www.biocloudservice.com\/wordpress\/wp-content\/uploads\/2024\/10\/1730268813228_DDC46894-915B-42ed-9CEB-F19DFD67854F.png?resize=300%2C99 300w\" sizes=\"(max-width: 554px) 100vw, 554px\" data-recalc-dims=\"1\" \/><\/figure>\n\n\n\n<p>2.2\u5feb\u67e5\u627e\u5171\u8868\u8fbe\u7279\u5f81\u57fa\u56e0<\/p>\n\n\n\n<p>\u6211\u4eec\u63d0\u4f9b\u4e86\u4e00\u79cd\u81ea\u9002\u5e94\u7248\u672c\u7684\u5438\u5f15\u5b50\u7b97\u6cd5\uff0c\u57fa\u4e8e\u6240\u63d0\u4f9b\u7684\u79cd\u5b50\u57fa\u56e0\u6765\u5bfb\u627e\u6536\u655b\u7684\u5438\u5f15\u5b50\u3002\u5982\u679c\u7528\u6237\u53ea\u5bf9\u5171\u540c\u8868\u8fbe\u5f0f\u5438\u5f15\u5668\u611f\u5174\u8da3\uff0c\u5219\u53ef\u4ee5\u4f7f\u7528\u4ee5\u4e0b\u4ee3\u7801\u67e5\u627e\u5171\u540c\u8868\u8fbe\u5f0f\uff0c\u800c\u65e0\u9700\u8fd0\u884c\u6574\u4e2aCASCC\u7ba1\u9053\u3002<\/p>\n\n\n\n<p>library(&#8220;CASCC&#8221;)<\/p>\n\n\n\n<p>library(&#8220;cafr&#8221;)<\/p>\n\n\n\n<p>data(&#8220;Data_PDAC_peng_2k&#8221;)<\/p>\n\n\n\n<p># \u627e\u5230AIF1\u79cd\u5b50\u5438\u5f15\u5b50\u3002AIF1\u57fa\u56e0\u662f\u5de8\u566c\u7ec6\u80de\u7684\u6807\u5fd7\u7269\u3002<\/p>\n\n\n\n<p>attr.res &lt;- findAttractor.adaptive(Data_PDAC_peng_2k, &#8220;AIF1&#8221;)<\/p>\n\n\n\n<p># \u5171\u8868\u8fbe\u7684\u5438\u5f15\u5b50<\/p>\n\n\n\n<p>head(names(attr.res$attractor.final)) &nbsp;&nbsp;<\/p>\n\n\n\n<p>#&gt; [1] &#8220;FCER1G&#8221; &#8220;AIF1&#8221; &nbsp;&nbsp;&#8220;TYROBP&#8221; &#8220;IGSF6&#8221; &#8220;CD163&#8221; &#8220;C1QA&#8221;<\/p>\n\n\n\n<p>2.3\u9010\u6b65\u8fd0\u884cCASCC<\/p>\n\n\n\n<p>\u5728\u672c\u8282\u4e2d\uff0c\u6211\u4eec\u5c06\u9010\u6b65\u5c55\u793a\u5982\u4f55\u8fd0\u884cCASCC\u3002<\/p>\n\n\n\n<p>rm(list = ls()); gc()<\/p>\n\n\n\n<p>library(&#8220;CASCC&#8221;)<\/p>\n\n\n\n<p>data(&#8220;Data_PDAC_peng_2k&#8221;)<\/p>\n\n\n\n<p># 1-3\u6b65<\/p>\n\n\n\n<p>fs.res &lt;- CASCC.featureSelection(Data_PDAC_peng_2k, inputDataType = &#8220;well-normalized&#8221;, attr.raw = NULL, exponent.max = 10, exponent.min = 2,<\/p>\n\n\n\n<p>&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;generalSeedList = NULL, mc.cores = 1, topDEGs = 10, topAttr = 50,<\/p>\n\n\n\n<p>&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;removeMT.RP.ERCC = TRUE, removeNonProtein = FALSE,<\/p>\n\n\n\n<p>&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;topN.DEG.as.seed = 1,<\/p>\n\n\n\n<p>&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;overlapN = 10)<\/p>\n\n\n\n<p>adata &nbsp;&lt;- fs.res$adata;<\/p>\n\n\n\n<p>SeuratVersionCheck = as.numeric(strsplit(as.character(packageVersion(&#8220;Seurat&#8221;)), &#8220;\\\\.&#8221;)[[1]][1])<\/p>\n\n\n\n<p>if (SeuratVersionCheck == 4) {<\/p>\n\n\n\n<p>&nbsp;&nbsp;data &lt;- adata[[&#8220;RNA&#8221;]]@data<\/p>\n\n\n\n<p>}else if (SeuratVersionCheck == 5) {<\/p>\n\n\n\n<p>&nbsp;&nbsp;data &lt;- adata[[&#8220;RNA&#8221;]]$data<\/p>\n\n\n\n<p>}<\/p>\n\n\n\n<p>data &lt;- as.matrix(data)<\/p>\n\n\n\n<p>#\u7b2c\u56db\u6b65<\/p>\n\n\n\n<p>res.predictK &lt;- predictClusterK(data, fs.res, method = &#8220;ward.D2&#8221;, index = &#8220;silhouette&#8221;, min.nc_fix = FALSE)<\/p>\n\n\n\n<p># \u7ed3\u679c\u5217\u8868<\/p>\n\n\n\n<p>res &lt;- list()<\/p>\n\n\n\n<p>res$fs.res &lt;- fs.res<\/p>\n\n\n\n<p>res$res.predictK &lt;- res.predictK<\/p>\n\n\n\n<p># \u7b2c\u4e94\u6b65<\/p>\n\n\n\n<p>res$mainType.output &lt;- kmeansCenterModify(res)<\/p>\n\n\n\n<p>3.\u5176\u4ed6\u811a\u672c<\/p>\n\n\n\n<p>3.1\u751f\u6210\u793a\u4f8b\u6570\u636e\u96c6<\/p>\n\n\n\n<p>set.seed(1)<\/p>\n\n\n\n<p>rm(list = ls()); gc()<\/p>\n\n\n\n<p>## \u539f\u59cb\u6570\u636e\u6765\u81ea\u4e8e&nbsp;PRJCA001063<\/p>\n\n\n\n<p>load(&#8220;.\/adata.T15.RData&#8221;)<\/p>\n\n\n\n<p>load(&#8220;.\/cell.type.RData&#8221;)<\/p>\n\n\n\n<p>cell.type &lt;- cell.type[which(cell.type$cell.name %in% colnames(adata)), ]<\/p>\n\n\n\n<p>table(cell.type$cluster)<\/p>\n\n\n\n<p>library(dplyr)<\/p>\n\n\n\n<p>sampled_cells &lt;- cell.type %&gt;%<\/p>\n\n\n\n<p>&nbsp;&nbsp;group_by(cluster) %&gt;%<\/p>\n\n\n\n<p>&nbsp;&nbsp;## if a cell type does not have more than 50 cells, we removed this cell type<\/p>\n\n\n\n<p>&nbsp;&nbsp;sample_n(size = if_else(n() &gt;= 50, 50, 0)) %&gt;%<\/p>\n\n\n\n<p>&nbsp;&nbsp;ungroup()<\/p>\n\n\n\n<p>print(sampled_cells)<\/p>\n\n\n\n<p>genes &lt;- Seurat::VariableFeatures(adata)<\/p>\n\n\n\n<p>cells &lt;- sampled_cells$cell.name<\/p>\n\n\n\n<p>SeuratVersionCheck = as.numeric(strsplit(as.character(packageVersion(&#8220;Seurat&#8221;)), &#8220;\\\\.&#8221;)[[1]][1])<\/p>\n\n\n\n<p>if (SeuratVersionCheck == 4) {<\/p>\n\n\n\n<p>&nbsp;&nbsp;data &lt;- as.matrix(adata[[&#8220;RNA&#8221;]]@data[genes, cells])<\/p>\n\n\n\n<p>}else if (SeuratVersionCheck == 5) {<\/p>\n\n\n\n<p>&nbsp;&nbsp;data &lt;- as.matrix(adata[[&#8220;RNA&#8221;]]$data[genes, cells])<\/p>\n\n\n\n<p>}<\/p>\n\n\n\n<p>save(data, file = &#8220;.\/CASCC\/data\/Data_PDAC_peng_2k.RData&#8221;)<\/p>\n\n\n\n<p>start_time &lt;- Sys.time()<\/p>\n\n\n\n<p>res &lt;- CASCC::run.CASCC(data)<\/p>\n\n\n\n<p>end_time &lt;- Sys.time()<\/p>\n\n\n\n<p>end_time &#8211; start_time<\/p>\n\n\n\n<p>save(res, file = &#8220;.\/CASCC\/data\/Res_PDAC_peng_2k.RData&#8221;)<\/p>\n\n\n\n<p>3.2\u5728Seurat\u5de5\u4f5c\u6d41\u4e2d\u8fd0\u884cCASCC<\/p>\n\n\n\n<p>\u8fd9\u91cc\u6211\u4eec\u5c06\u5c55\u793a\u5982\u4f55\u5728Seurat\u5de5\u4f5c\u6d41\u4e2d\u8fd0\u884cCASCC\u3002<\/p>\n\n\n\n<p>SeuratVersionCheck = as.numeric(strsplit(as.character(packageVersion(&#8220;Seurat&#8221;)), &#8220;\\\\.&#8221;)[[1]][1])<\/p>\n\n\n\n<p>if (SeuratVersionCheck == 4) {<\/p>\n\n\n\n<p>&nbsp;&nbsp;data &lt;- SeuratObject[[&#8220;RNA&#8221;]]@data<\/p>\n\n\n\n<p>}else if (SeuratVersionCheck == 5) {<\/p>\n\n\n\n<p>&nbsp;&nbsp;data &lt;- SeuratObject[[&#8220;RNA&#8221;]]$data<\/p>\n\n\n\n<p>}<\/p>\n\n\n\n<p>res &lt;- CASCC::run.CASCC(data)<\/p>\n\n\n\n<p>Idents(SeuratObject) &lt;- res$mainType.output$clusteringResults<\/p>\n\n\n\n<p>SeuratObject<\/p>\n\n\n\n<p><strong>\u5c0f\u7ed3<\/strong><strong><\/strong><\/p>\n\n\n\n<p>\u4eca\u5929\u5c0f\u679c\u5e26\u7740\u5927\u5bb6\u5b66\u4e60\u4e86\u4e00\u79cd\u7528\u4e8escRNA-seq\u5206\u6790\u7684\u5171\u8868\u8fbe\u8f85\u52a9\u5355\u7ec6\u80de\u805a\u7c7b\uff08CASCC\uff09\u65b9\u6cd5\u3002CASCC\u7684\u4e00\u4e2a\u5173\u952e\u7279\u70b9\u662f\u4f7f\u7528\u4e86\u4e00\u79cd\u65b0\u7684\u81ea\u9002\u5e94\u5438\u5f15\u5b50\u7b97\u6cd5\uff0c\u8be5\u7b97\u6cd5\u6709\u52a9\u4e8e\u6b63\u786e\u7684\u7279\u5f81\u9009\u62e9\u3001\u805a\u7c7b\u548c\u7ec6\u80de\u7c7b\u578b\u9274\u5b9a\uff0c\u4ece\u800c\u4e0e\u73b0\u6709\u65b9\u6cd5\u76f8\u6bd4\u63d0\u9ad8\u4e86\u751f\u7269\u5b66\u51c6\u786e\u6027\u3002\u6700\u540e\u5c0f\u679c\u7ed9\u5927\u5bb6\u4ecb\u7ecd\u4e00\u4e2a\u4e91\u5de5\u5177\uff01\u540c\u5b66\u4eec\u5982\u679c\u89c9\u5f97\u81ea\u5df1\u7684\u4ee3\u7801\u6c34\u5e73\u4e00\u822c\uff0c\u5bf9\u4e8e\u5f88\u591a\u7684\u53c2\u6570\u4e0d\u77e5\u9053\u600e\u4e48\u6539\uff0c\u53ef\u4ee5\u4f53\u9a8c\u4e00\u4e0b\u6211\u4eec\u7684\u4e91\u751f\u4fe1\u5c0f\u5de5\u5177\uff0c\u53ea\u9700\u8f93\u5165\u6570\u636e\uff0c\u5373\u53ef\u8f7b\u677e\u751f\u6210\u6240\u9700\u56fe\u8868\uff0c\u5b57\u4f53\u5927\u5c0f\u3001\u6807\u9898\u7b49\u4e5f\u53ef\u4e00\u952e\u66f4\u6539\u3002\u611f\u5174\u8da3\u7684\u5c0f\u4f19\u4f34\u53bb\u4e91\u751f\u4fe1\uff08http:\/\/www.biocloudservice.com\/home.html\uff09\u4f53\u9a8c\u4e00\u4e0b\u5427\uff01<\/p>\n","protected":false},"excerpt":{"rendered":"<p>\u524d\u8a00 \u5c0f\u679c\u89c1\u9762\u4f1a\u53c8\u8981\u5f00\u59cb\u4e86\uff0c\u867d\u7136\u5c0f\u679c\u5b9e\u9a8c\u5f88\u5fd9\uff0c\u4f46\u662f\u770b\u5230\u597d\u7528\u7684\u5355\u7ec6\u80de\u805a\u7c7bR\u5305\u8fd8\u662f\u5fcd\u4e0d\u4f4f\u7b2c\u4e00\u65f6\u95f4\u5411\u5927\u5bb6\u5206\u4eab\uff01\u60f3\u5fc5\u5927 [&hellip;]<\/p>\n","protected":false},"author":1,"featured_media":0,"comment_status":"closed","ping_status":"","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\/61543"}],"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=61543"}],"version-history":[{"count":1,"href":"http:\/\/www.biocloudservice.com\/wordpress\/index.php?rest_route=\/wp\/v2\/posts\/61543\/revisions"}],"predecessor-version":[{"id":61549,"href":"http:\/\/www.biocloudservice.com\/wordpress\/index.php?rest_route=\/wp\/v2\/posts\/61543\/revisions\/61549"}],"wp:attachment":[{"href":"http:\/\/www.biocloudservice.com\/wordpress\/index.php?rest_route=%2Fwp%2Fv2%2Fmedia&parent=61543"}],"wp:term":[{"taxonomy":"category","embeddable":true,"href":"http:\/\/www.biocloudservice.com\/wordpress\/index.php?rest_route=%2Fwp%2Fv2%2Fcategories&post=61543"},{"taxonomy":"post_tag","embeddable":true,"href":"http:\/\/www.biocloudservice.com\/wordpress\/index.php?rest_route=%2Fwp%2Fv2%2Ftags&post=61543"}],"curies":[{"name":"wp","href":"https:\/\/api.w.org\/{rel}","templated":true}]}}