RNA-seq流程学习笔记（17）- PCA图聚类分析

本文链接：https://blog.csdn.net/xiaomotong123/article/details/115195507

本文记录了RNA-seq数据处理中的PCA图制作和聚类分析过程，探讨了如何通过PCA揭示样本间的差异，并利用聚类方法对样本进行分类，为后续生物信息学研究提供关键洞察。

摘要生成于 C知道，由 DeepSeek-R1 满血版支持，前往体验 >

#绘制PCA图
#需要使用各组的FPKM进行绘制
#对FPKM数据进行整理
#清空环境变量
rm(list=ls())

##将StringTie分析得到的含有FPKM数据的TAB文件导入当前工作环境中

#设置工作目录
setwd("G:/kongyu/RNA-seq/2021_02_22/gene_tab/")

DIPG4_1_1.gene.tab <- read.table("G:/kongyu/RNA-seq/2021_02_22/gene_tab/DIPG4_1_1.gene.tab", header = TRUE, sep = "\t" , quote = "\"")
DIPG4_1_2.gene.tab <- read.table("G:/kongyu/RNA-seq/2021_02_22/gene_tab/DIPG4_1_2.gene.tab", header = TRUE, sep = "\t" , quote = "\"")
DIPG4_2_1.gene.tab <- read.table("G:/kongyu/RNA-seq/2021_02_22/gene_tab/DIPG4_2_1.gene.tab", header = TRUE, sep = "\t" , quote = "\"")
DIPG4_2_2.gene.tab <- read.table("G:/kongyu/RNA-seq/2021_02_22/gene_tab/DIPG4_2_2.gene.tab", header = TRUE, sep = "\t" , quote = "\"")
DIPG4_3_1.gene.tab <- read.table("G:/kongyu/RNA-seq/2021_02_22/gene_tab/DIPG4_3_1.gene.tab", header = TRUE, sep = "\t" , quote = "\"")
DIPG4_3_2.gene.tab <- read.table("G:/kongyu/RNA-seq/2021_02_22/gene_tab/DIPG4_3_2.gene.tab", header = TRUE, sep = "\t" , quote = "\"")
DIPG4_4_1.gene.tab <- read.table("G:/kongyu/RNA-seq/2021_02_22/gene_tab/DIPG4_4_1.gene.tab", header = TRUE, sep = "\t" , quote = "\"")
DIPG4_4_2.gene.tab <- read.table("G:/kongyu/RNA-seq/2021_02_22/gene_tab/DIPG4_4_2.gene.tab"