在GAGE软件包中,您需要使用函数kegg.gsets()生成基因集数据:
库(gage)kg.pfa = kegg.gsets(“ pfa”)kegg.gs = kg.pfa $ kg.sets [kg.pfa $ sigmet.idx]
注意该物种是PFA。PFD不是kegg.gsets()当前支持的6833种。这可能是最近对Kegg的补充。您需要使用类似内容手动将其添加到Korg数据中:
pfd.info = c(“ T01051”,“ 57267”,“ PFD”,“ plasmodium falciparum dd2”,“”,“ 0”,“ pfdg_03934”,pfdg_03934“,na,na,na,na,na)data(korg)data(korg)korg = rbind(korg = rbind(korg),pfd.info)kg.pfd = kegg.gsets(“ pfd”)kegg.gs = kg.pfd $ kg.sets [kg.pfd $ sigmet.idx] head(kegg.gs,3)
请注意,基因ID看起来像:“ PFDG_00033”“ PFDG_00121”等。请确保您在表达式数据中使用了相同类型的基因ID。
您也可以尝试我们的新软件包SBGNVIEW,以了解KEGG之外的更多途径:
github
https://github.com/datapplab/sbgnview
Bioc
//www.andersvercelli.com/packages/sbgnview/
Hths。
你尝试过吗keggrest包裹?似乎有一些所需生物的信息
库(keggrest)> org <-kegglist(“有机体”)> org [grep(org [org [topt''],pattern =“ falciparum”),] t.number生物体[1,]“ T00095”“ pfa”“ pfa”“恶性疟原虫3D7” [2,]“ T01051”“ pfd”“” plasmodium falciparum dd2“ [3,]“ T01052”“ pfh”“ pfh”“ pfh”“ pfh”疟原虫恶性疟原虫hb3“ hb3” hb3“ hb3”系统发育[1,]2,]“真核生物;生物;肺泡; apicomplexans''[3,]”真核生物;生物学家;肺泡; apicomplexans''
StringDB软件包具有可用于的GO-TIMM,KEGG和其他富集分析恶性疟原虫
//www.andersvercelli.com/packages/release/bioc/html/stringdb.html
您也可以尝试基于Web的G:Profiler
恶性疟原虫可作为Go,Kegg,...富集的生物。
本网站的使用构成接受我们的用户协议和隐私政策。
恐怕我是R的初学者
如何上传“恶性疟原虫DD2”数据?
查看Keggrest Vignette,以获取有关如何使用软件包的信息。在这里找到://www.andersvercelli.com/packages/release/bioc/vignettes/keggrest/keggrest/inst/doc/keggrest-vignette.html
我没有使用太多的包裹,所以我鼓励您查看包装插图和文档,但我认为要获得可以做类似的途径
kegglink(“ pathway”,“ pfd”)