07_string

作者: 养猪场小老板 | 来源:发表于2020-02-12 22:17 被阅读0次

    制作string的输入数据

    load("step4output.Rdata")
    gene_up= deg[deg$change == 'up','symbol'] 
    gene_down=deg[deg$change == 'down','symbol'] 
    write.table(gene_up,
                file="upgene.txt",
                row.names = F,
                col.names = F,
                quote = F)
    write.table(gene_down,
                file="downgene.txt",
                row.names = F,
                col.names = F,
                quote = F)
    write.table(deg$symbol[1:200],
                file="diffgene.txt",
                row.names = F,
                col.names = F,
                quote = F)
    #--------------
    

    准备cytoscape的输入文件

    tsv = read.table("string_interactionsdiff.tsv",comment.char = "!",header = T)
    tsv2 = tsv[,c(1,2,ncol(tsv))]
    head(tsv2)
    write.table(tsv2,
                file = "cyto.txt",
                sep = "\t",
                quote = F,
                row.names = F)
    
    p = deg[deg$change != "stable",c("symbol","logFC","P.Value")]
    head(p)
    write.table(p,
                file = "deg.txt",
                sep = "\t",
                quote = F,
                row.names = F)
    

    相关文章

      网友评论

        本文标题:07_string

        本文链接:https://www.haomeiwen.com/subject/axdhfhtx.html