昨天生信技能树发布了学徒作业:学徒作业-在CCLE数据库里面根据指定基因在指定细胞系里面提取表达矩阵 很有意思,任务简单的说就是重复这个图
?...id<- lapply(id1[,9], function(x){
y=strsplit(x,';')[[1]][5]
strsplit(y,' ')[[1]][3]
}) #把gtf的第9列拆一下获得...#CFP=properdin #SERPING1=C1INH #C4BP=C4BPA #VSIG4=CRIg
我手动替换了新的i基因名字,然后我喜欢处理数据框,把它变成了数据框。...),先建立向量,然后变成数据框
xb <- c("BT474","BT549",'MDA-MB-231','HCC1937','MDA-MB-361','MDA-MB-436',"AU565","SK-BR...(x3)))
将细胞的名字全部取出来,变成数据框 因为我喜欢处理数据框
w3<- data.frame(n=colnames(x3),
n2=rep(1,1021)) #建立相匹配的列