適用背景 之前寫了一篇文章介紹四步完成單細(xì)胞數(shù)據(jù)調(diào)控網(wǎng)絡(luò)流程分析[http://www.reibang.com/p/71ed43163ce1]易迹,但生成的結(jié)果沒有進(jìn)行可視化训貌,...
適用背景 之前寫了一篇文章介紹四步完成單細(xì)胞數(shù)據(jù)調(diào)控網(wǎng)絡(luò)流程分析[http://www.reibang.com/p/71ed43163ce1]易迹,但生成的結(jié)果沒有進(jìn)行可視化训貌,...
轉(zhuǎn)錄因子 (transcription factors, TFs) 是直接作用于基因組斗塘,與特定DNA序列結(jié)合 (TFBS/motif) 胀溺,調(diào)控DNA轉(zhuǎn)錄過程的一類蛋白質(zhì)。轉(zhuǎn)錄...
往期回顧: 在前面的課程中我們已經(jīng)進(jìn)行過“單樣本數(shù)據(jù)分析[https://link.zhihu.com/?target=http%3A//mp.weixin.qq.com/s...
這是ggplot2可視化專題的第二個實(shí)例操作 【ggplot2的基本思路見前文總論】:基于ggplot2的RNA-seq轉(zhuǎn)錄組可視化:總述和分文目錄 【ggplot2繪圖具體...
Q: 得到的差異基因和保守基因蚜厉,用哪個作為 marker 基因好长已? / FindallMarkers 和findconservedmarkers 結(jié)果上有什么區(qū)別嗎? A: ...
前面也說過用CCA昼牛,SCTransform和Harmony等可以消除多個數(shù)據(jù)的批次效應(yīng)术瓮,來實(shí)現(xiàn)多組數(shù)據(jù)的整合。但有時候我們會遇到兩個datasets只有部分重疊贰健,這和之前介紹...
1. 關(guān)于 QC Q:過濾線粒體基因是去除低質(zhì)量 cell胞四,核糖體基因用過濾嗎? A: Seurat 教程中沒有提到用核糖體基因的 UMI 含量來過濾細(xì)胞霎烙,我看到的文獻(xiàn)中基本...
前段時間在生信技能樹單細(xì)胞學(xué)習(xí)小分隊(duì)分享單細(xì)胞文章復(fù)現(xiàn)撬讽,一不小心就踩了一個坑……記錄在這兒和大家一起分享蕊连。 背景 主要是想復(fù)現(xiàn)2020年發(fā)表在 PNAS 上的一篇題為 Sin...
首先華大沒有給我塞錢悬垃,最近有用華大的測序,就學(xué)習(xí)了一下甘苍。 華大的測序技術(shù)起源于 CG(Complete Genomics) 的 cPAL(combinatorial Prob...
1.利用箱線圖比較兩類樣本的某個細(xì)胞比例差異 比較直觀尝蠕,但是缺點(diǎn)在于如果單細(xì)胞樣本個數(shù)過少且異質(zhì)性大,導(dǎo)致很難有統(tǒng)計(jì)學(xué)顯著意義 2.R o/e 比值 好多文章都有用這個载庭,我的...
一文了解單細(xì)胞對象數(shù)據(jù)結(jié)構(gòu)/數(shù)據(jù)格式囚聚,單細(xì)胞數(shù)據(jù)操作不迷茫靖榕。本文內(nèi)容包括 單細(xì)胞seurat對象數(shù)據(jù)結(jié)構(gòu), 內(nèi)容構(gòu)成顽铸,對象的調(diào)用茁计、操作,常見函數(shù)的應(yīng)用等谓松。 1. 單細(xì)胞seu...
前段時間一直在用cellranger count進(jìn)行單細(xì)胞測序數(shù)據(jù)的處理星压,不得不說Cell Ranger作為10X Genomics官方的單細(xì)胞數(shù)據(jù)處理工具,功能還是很強(qiáng)大的...
一鬼譬、ScaleData()簡介 單細(xì)胞基因表達(dá)counts矩陣數(shù)據(jù)經(jīng)過NormalizeData()歸一化處理后娜膘,還需要進(jìn)行scale標(biāo)準(zhǔn)化。ScaleData()[http...
和前面的scCATCH一樣军洼,SCSA(https://github.com/bioinfo-ibms-pumc/SCSA)也是基于已知marker gene進(jìn)行打分的注釋方法...
昨天頭疼了大半天的報錯問題報錯如下: 查遍了資料,無果演怎。 結(jié)果發(fā)現(xiàn)我的路徑里面含有中文歉眷,數(shù)據(jù)就無法載入!可是我下載的網(wǎng)上的GSE數(shù)據(jù)颤枪,路徑無論是中文還是英文都沒問題呀汗捡!而且R...
1扇住、如果你是直接拿到的R1/R2的fastq文件,那么就直接上cellranger count盗胀。氮素艘蹋,如果你是I1/R1/R2的數(shù)據(jù),那就麻煩還要跑個cellranger m...