TCGA_miRNA数据和mRNA数据整合

发表于2018年7月18日由daizao

R语言方法

setwd("D:/TCGA_coexpression")
miRNA <- read.table("diff_miRNAExp.txt",row.names=1,header=T,sep="\t",check.names=F)
mRNA <- read.table("diffmRNAExp.txt",row.names=1,header=T,sep="\t",check.names=F)
colnames(miRNA) <- gsub("(....)?-(..)?-(....)?-(...)?.*","\\1-\\2-\\3-\\4",colnames(miRNA)) #此处的?号用于非贪婪匹配，最后的.*可以让最后的替换成功
colnames(mRNA) <- gsub("(....)?-(..)?-(....)?-(...)?.*","\\1-\\2-\\3-\\4",colnames(mRNA))
sameid <- intersect(colnames(miRNA),colnames(mRNA))
merge <- rbind(id=same,miRNA[,sameid],mRNA[,sameid])
write.table(merge,file="merge.txt",sep="\t",quote=F,col.names=F)

miRNALable <- cbind(rownames(miRNA),"miRNA")
mRNALable <- cbind(rownames(mRNA),"gene")
node <- rbind(c("ID","Classify"),miRNALable,mRNALable)
write.table(node,file="nodelable.txt",sep="\t",quote=F,col.names=F,row.names=F)
此条目发表在TCGA分类目录。将固定链接加入收藏夹。

发表评论

邮箱地址不会被公开。 必填项已用*标注

To create code blocks or other preformatted text, indent by four spaces:

    This will be displayed in a monospaced font. The first four 
    spaces will be stripped off, but all other whitespace
    will be preserved.
    
    Markdown is turned off in code blocks:
     [This is not a link](http://example.com)

To create not a block, but an inline code span, use backticks: 

Here is some inline `code`.

For more help see http://daringfireball.net/projects/markdown/syntax

Protected with IP Blacklist CloudIP Blacklist Cloud