首先根據(jù)文章作者提供的metadata文件(圖一),以及作者的count.txt.gz恋腕,隨后讀取成功得到metadata.csv和single cell_raw_scRNAseq.rds抹锄,打開metadata.csv是圖4,發(fā)現(xiàn)sample不是GSM開頭荠藤,所以sample_info.csv文件那里(圖5)數(shù)據(jù)清洗選擇title列為樣本id所在的列名伙单,樣本分組信息所在的列還是tissue:ch1,批次效應(yīng)列選了patient:ch1哈肖,去掉metastatic 組得到sample_info_sub.csv(圖7)吻育,隨后進行質(zhì)控分析,運行完還是只有一個after_qc_sc_list.rds文件淤井,期間試過把原作者的metadata文件的sample列換成GSM名字重新讀取布疼,質(zhì)控后還是一樣只有一個文件,不知道是哪里出來問題
讀取單個單細胞表達矩陣文件count.txt.gz的問題求助2
?著作權(quán)歸作者所有,轉(zhuǎn)載或內(nèi)容合作請聯(lián)系作者
- 文/潘曉璐 我一進店門碍扔,熙熙樓的掌柜王于貴愁眉苦臉地迎上來肴沫,“玉大人,你說我怎么就攤上這事蕴忆〔遥” “怎么了?”我有些...
- 文/不壞的土叔 我叫張陵,是天一觀的道長站蝠。 經(jīng)常有香客問我汰具,道長,這世上最難降的妖魔是什么菱魔? 我笑而不...
- 正文 為了忘掉前任留荔,我火速辦了婚禮,結(jié)果婚禮上澜倦,老公的妹妹穿的比我還像新娘聚蝶。我一直安慰自己,他們只是感情好藻治,可當我...
- 文/花漫 我一把揭開白布碘勉。 她就那樣靜靜地躺著,像睡著了一般桩卵。 火紅的嫁衣襯著肌膚如雪验靡。 梳的紋絲不亂的頭發(fā)上,一...
- 文/蒼蘭香墨 我猛地睜開眼孙技,長吁一口氣:“原來是場噩夢啊……” “哼!你這毒婦竟也來了排作?” 一聲冷哼從身側(cè)響起牵啦,我...
- 正文 年R本政府宣布,位于F島的核電站毡惜,受9級特大地震影響拓轻,放射性物質(zhì)發(fā)生泄漏。R本人自食惡果不足惜经伙,卻給世界環(huán)境...
- 文/蒙蒙 一扶叉、第九天 我趴在偏房一處隱蔽的房頂上張望。 院中可真熱鬧帕膜,春花似錦枣氧、人聲如沸。這莊子的主人今日做“春日...
- 文/蒼蘭香墨 我抬頭看了看天上的太陽。三九已至危纫,卻和暖如春,著一層夾襖步出監(jiān)牢的瞬間乌庶,已是汗流浹背种蝶。 一陣腳步聲響...
推薦閱讀更多精彩內(nèi)容
- 利用文章作者提供的count.txt.gz文件和metadata文件成功運行了第一步1.seurat酗电,現(xiàn)在里面的文...
- 單細胞研究使我們能夠了解細胞發(fā)育的分子過程和病理的發(fā)生機理魄藕。單細胞ATAC-seq(scATAC-seq)通過分析...
- 這次學習的是單細胞測序的第三個比較常用的包:Scater。本文主要參考教程:單細胞轉(zhuǎn)錄組學習筆記-15-利用scR...
- 轉(zhuǎn)自Seurat_Satija 如何使用 Seurat 分析單細胞測序數(shù)據(jù)( Q&A[https://www.ji...
- 1. 關(guān)于 QC Q:過濾線粒體基因是去除低質(zhì)量 cell撵术,核糖體基因用過濾嗎背率? A: Seurat 教程中沒有提...