1. Seurat對象查看當(dāng)前的Assay 在進(jìn)行了SCTransform操作后,矩陣默認(rèn)會(huì)變成SCT矩陣,如果不加設(shè)置,后續(xù)的PCA等操作都是基于SCT矩陣。修改Defau...
hello,大家好,今天周一,做一個(gè)簡單的匯總,主要針對基因集打分,無論是單細(xì)胞還是空間,基因集打分都很普遍,方法也很多,我們來總結(jié)一下。 研究背景 在單細(xì)胞轉(zhuǎn)錄組測序中,整...
殺殺 (第一步不作贅述,一般會(huì)使用單因素cox或者尋找差異基因等方法,篩選與生存相關(guān)的分子,準(zhǔn)備進(jìn)行下一步的模型構(gòu)建) 本篇舉例是篩選了12個(gè)單因素cox與生存顯著相關(guān)的分子...
1.利用箱線圖比較兩類樣本的某個(gè)細(xì)胞比例差異 比較直觀,但是缺點(diǎn)在于如果單細(xì)胞樣本個(gè)數(shù)過少且異質(zhì)性大,導(dǎo)致很難有統(tǒng)計(jì)學(xué)顯著意義 2.R o/e 比值 好多文章都有用這個(gè),我的...
生存模型的構(gòu)建方法包括:1. Lasso回歸;2. Cox多因素回歸;3. 隨機(jī)森林;4. 支持向量機(jī) 可以把log_rank_test或cox篩選出的基因單獨(dú)做模型構(gòu)建,也...
1. 簡介 隨機(jī)森林是機(jī)器學(xué)習(xí)中的一種分類算法,這個(gè)算法在之前分享的TCGA生存模型的構(gòu)建以及模型預(yù)測和評(píng)估[http://www.itdecent.cn/p/334254...
殺殺報(bào)錯(cuò)信息:Error in download.KEGG.Path(species): 'species' should be one of organisms l 再運(yùn)行...
ggplot2自定義函數(shù)方法 前言:(百年更新,太懶了實(shí)在)最近寫R包和平常使用的時(shí)候,有的時(shí)候想快速畫一個(gè)自己定義好的圖,每一次都要一長串代碼,有點(diǎn)不是很方便,于是想編寫一...
——by不是殺殺Samba是在Linux和UNIX系統(tǒng)上實(shí)現(xiàn)SMB協(xié)議的一個(gè)免費(fèi)軟件,由服務(wù)器及客戶端程序構(gòu)成(來自百度百科),簡單來說它可以實(shí)現(xiàn)服務(wù)器端(Linux)與Wi...
殺殺 發(fā)一些拼圖的小技巧 拼圖之前,先要根據(jù)自己圖片的張數(shù),設(shè)置好畫布大小,這很重要我畫圖的經(jīng)驗(yàn)中,正常一張方形的圖會(huì)設(shè)置為width和height都等于5,如果兩張橫著拼在...
生存分析圖最近我發(fā)現(xiàn)可以用patchwork來拼。 a$plot + b$plot。這樣子就可以了,因?yàn)閜lot是ggplot格式,取出來就可以拼。
R語言拼圖(含生存分析),R語言畫布設(shè)置殺殺 發(fā)一些拼圖的小技巧 拼圖之前,先要根據(jù)自己圖片的張數(shù),設(shè)置好畫布大小,這很重要我畫圖的經(jīng)驗(yàn)中,正常一張方形的圖會(huì)設(shè)置為width和height都等于5,如果兩張橫著拼在...
方法一在MSigDB下載,網(wǎng)址GSEA | MSigDB | MSigDB Collections。下載到 gmt 格式文件使用 GSEABase 包的 getGmt 函數(shù)讀...