生信星球 DAY 1 —— 橙子🍊

rm(list=ls())

#读取CGGA_325的两个原始csv
CGGA_325_clinical<-read.csv('CGGA.mRNAseq_325_clinical.20200506.csv')
CGGA_325_genes<-read.csv('CGGA.mRNAseq_325.RSEM-genes.20200506.csv')
class(CGGA_325_genes)
class(CGGA_325_clinical)

#从CGGA_325_genes中选取KLRB1行
a<-CGGA_325_genes[,1]
KLRB1<-CGGA_325_genes[a=='KLRB1',]
KLRB1_new<-KLRB1[,-1]
KLRB1_new<-t(KLRB1_new)
colnames(KLRB1_new)<-'KLRB1'

#将行名赋值为第一列
b<-KLRB1_new
CGGA_ID<-rownames(b)
rownames(b)<-NULL
KLRB1_tr<-cbind(CGGA_ID,b)

#合并2表
CGGA_325_KLRB1<-merge(CGGA_325_clinical,KLRB1_tr,by='CGGA_ID')
View(CGGA_325_KLRB1)
write.csv(CGGA_325_KLRB1,file = 'CGGA_325_KLRB1.csv')