参考:https://www.jianshu.com/p/3b9eee8c31cc
优雅R:「Workshop」第三十三期 用iRegulon进行主转录因子的预测
挑圈联靠:转录因子经典数据库,为你的文章增光添彩
挑圈联靠:90%以上标书中有这个生信图!用这几个数据库就够了!能不能成为明年黑马就靠它了!
打开网站https://www.networkanalyst.ca/home.xhtml
点击gene list input,选择物种为人,ID type为entrez ID。
注意这里涉及到基因名和entrez ID的转换,代码如下:
library(org.Hs.eg.db)
hs <- org.Hs.eg.db
my.symbols <- c("gene1", "gene2")
select(hs,
keys = my.symbols,
columns = c("ENTREZID", "SYMBOL"),
keytype = "SYMBOL")
点击上传、下一步,选择TF-gene interaction:
这里选择steiner forest network:
再点击下一步,最后下载sif文件:
将subnetwork1.sif导入cytoscape(network from file),将subnetwork1_Types.NA重新整理成一个两列的文件,第一列为基因名,第二列为属性(基因or转录因子):
然后将这个文件导入cytoscape(table from file),最后调整网络属性即可。
几个小tips:
1.在edge里,把线的宽度和透明度选小一点。
2.将名字调整到图形外,这个可以在node-label position里面调。