一始苇、在跑cellranger的過程中出了個(gè)問題踢匣。要琢磨一下怎么回事椿浓。 報(bào)錯(cuò):sequence and quality length mismat...
10X文庫送測序后作烟,從測序公司拿到的測序數(shù)據(jù)是fastq格式的,要經(jīng)過linux上跑cellranger程序,得到表達(dá)矩陣,才能做后面的功能分析...
學(xué)習(xí)及參考教程:https://mp.weixin.qq.com/s/uLt-bYCdVcBH6mqoIG_r6w[https://mp.wei...
接著day51壹置,繼續(xù)作圖: 說明先設(shè)置pvalue和logFC的閾值根據(jù)閾值分別為上調(diào)基因設(shè)置‘up’竞思,下調(diào)基因設(shè)置‘Down’表谊,無差異設(shè)置‘N...
對day49,day50兩次學(xué)習(xí)的總結(jié)盖喷。拿到reads表達(dá)矩陣后爆办,在自己電腦上用DESeq2這個(gè)R包進(jìn)行數(shù)據(jù)分析。 一课梳、先要準(zhǔn)備表達(dá)矩陣和分組表...
接著day49 一距辆、看看整體表達(dá)情況 原始矩陣為變量exprSet,是從txt文件讀入得到的數(shù)據(jù)框歸一化后矩陣為變量exprSet_new暮刃,用D...
學(xué)習(xí)資料:B站視頻:生物技能樹 第八季轉(zhuǎn)錄組測序數(shù)據(jù)分析得到表達(dá)矩陣之后就要考慮如何知道到差異表達(dá)基因(DEG)跨算,這個(gè)要在R中操作,具體是需要a...
轉(zhuǎn)錄組分析要用比對后的bam或sam文件進(jìn)行后續(xù)計(jì)數(shù)椭懊。也就是對基因/轉(zhuǎn)錄本的區(qū)域上比對到的reads進(jìn)行計(jì)數(shù)诸蚕。最后算出來某個(gè)轉(zhuǎn)錄本或者外顯子上比...
學(xué)習(xí)資料:B站視頻:生物技能樹 第八季轉(zhuǎn)錄組測序數(shù)據(jù)分析 一、轉(zhuǎn)錄組比對常用軟件 轉(zhuǎn)錄組比對常用:hisat2、subjunct(subread...