生信星球 DAY 1 —— 橙子

2023-10-21 00:09:36 浏览数 (1)

STEP 1:不会就搜索

1. 好用的搜索引擎

谷歌*、必应、虫部落快搜

snipaste贴图的虫部落快搜snipaste贴图的虫部落快搜

2. 去哪儿找专业教程?

搜狗微信、搜狗知乎、github*

STEP 2:高效学习

1. 高效软件

chrome浏览器|沙拉查词*、scholarscope

snipaste

2. 思维导图

幕布、Xmind(markdown格式)

代码语言:txt复制
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')

生信星球学习小组第159期

0 人点赞