姓名 模式 大小
.github文件 040000
R(右) 040000
数据 040000
安装 040000
男人 040000
测验 040000
小插曲 040000
.R建筑物 100644 0千字节
.gitignore(.git忽略) 100644 0千字节
.travis.yml文件 100644 1千字节
描述 100644 1千字节
执照 100644 0千字节
NAMESPACE公司 100644 1千字节
自述。房间 100644 2千字节
自述.md 100644 4千字节
自述.md
[![R-CMD-检查](https://github.com/ctlab/gatom/actions/workflows/R-CMD-check.yaml/bode.svg)](https://github.com/ctlab/gatom/actions/workflows/R-CMD-check.yaml)#关贸总协定寻找原子跃迁中活性代谢模块的R包网络。可以找到完整的渐晕图[此处](https://rpubs.com/asergushichev/gatom-tutorial网站).###安装```第页库(devtools)安装github(“ctlab/gatom”)```###快速入门```第页图书馆(gatom)库(data.table)库(igraph)图书馆(mwcsr)```首先,让我们用原子映射(`network`对象)、酶加载数据小鼠注释(`org.Mm.eg.gatom`)和代谢物注释(`met.kegg.db.rda`):```对数据(“networkEx”)数据(“org.Mm.eg.gatom.annoEx”)数据(“met.kegg.dbEx”)```加载输入数据:```第页数据(“met.de.rawEx”)数据(“gene.de.rawEx”)```获取原子图:```第页g<-makeMetabolicGraph(网络=网络Ex,拓扑=“atoms”,org.gatom.anno=org.Mm.eg.gatom.annoEx,gene.de=基因.de.rawEx,met.db=met.kegg.dbEx,met.de=met.de.rawEx)```##找到具有RefSeq ID的基因的DE表##找到HMDB ID代谢物的DE表```第页打印(g)```##IGRAPH b23a8ad联合国——176 190——##+属性:名称(v/c)、代谢物(v/c##|(v/c),pval(v/n),原点##|(v/n)、logPval(v/n)、信号(v/c)、标签(e/c),pval(e/n),来源##|(e/n),RefSeq(e/c),基因,##|reaction_equation(e/c),url(e/c),reaction(e/c),log2FC(e/n),##|baseMean(e/n),logPval##+来自b23a8ad的边(顶点名称):##[1]C00025_-0.3248_2.8125-C00026_-0.3248 _2.8125##[2]C00025_-1.6238_3.5625-C00026_-1.628_3.5625##[3]C00025_-2.9228_2.8125-C00026_-2.9225_2.8125## + ... 省略了几个边得分图,获得SGMWCS(信号广义最大权重子图)问题实例:```第页gs<-评分图(g,k.gene=25,k.met=25)```初始化SMGWCS解算器(启发式松弛与剪切解算器`rnc_solver `用于简化,请查看`mwcsr`包更多选项的文档):```第页解算器<-rnc_solver()```查找模块:```第页res<-solve_mwcsp(解算器,gs)m<-res$图形``````第页打印(m)```##IGRAPH e74bc24联合国——37 36——##+属性:信号(g/n)、名称(v/c)、代谢物,##|标签(v/c)、url(v.c)、pval(v/n##|(v/n),baseMean(v/n.),logPval(v/n:),signal(v/c),score(v/n-),##|标签(e/c),pval(e/n,##|酶(e/c),反应_名称(e/c),反作用_等式(e&c),url##|(e/c),反应(e/c),log2FC(e/n),基本平均值(e/n),logPval(e/m),##|signal(e/c),signalRank(e/n),score(e/n)##+来自e74bc24的边(顶点名称):##[1]C00025_-2.9228_2.8125-C00026_-2.9225_2.8125##[2]C00024_15.0644_27.8518-C00033_-1.6238_0.5625## + ... 省略了几个边```第页头部(E(m)$标签)```##[1]“Psat1”“Acss2”“Gpt2”“Got2”“Pkm”“Tpi1”```第页头(V(m)$标签)```##[1]“丙酮酸”“乙酰-CoA”“L-谷氨酸”“2-氧化戊二酸”##[5]“醋酸盐”“草酸酯”我们可以将模块保存为不同的格式(dot、xgmml、svg、pdf):```第页saveModuleToPdf(m,file=“M0.vs.M1.pdf”,name=“M0.vs.M1”,n_iter=100,force=1e-5)```![模块](https://rawgit.com/ctlab/gatom/master/inst/M0.vs.M1.png)