#绘制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", header
RNA-seq流程学习笔记(17)- PCA图聚类分析
最新推荐文章于 2024-05-18 03:44:14 发布