首先根據(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ì)控后還是一樣只有一個文件,不知道是哪里出來問題
隨后