zoukankan      html  css  js  c++  java
  • clusterProfiler包

    1)enrichGO:(GO富集分析)

    描述:GO Enrichment Analysis of a gene set. Given a vector of genes, this function will return the enrichment GO categories after FDR control

    语法:
    enrichGO(gene, OrgDb, keyType = "ENTREZID", ont = "MF",pvalueCutoff = 0.05, pAdjustMethod = "BH", universe, qvalueCutoff = 0.2,
    minGSSize = 10, maxGSSize = 500, readable = FALSE, pool = FALSE)

    参数介绍:
    gene:a vector of entrez gene id.
    OrgDb:OrgDb
    keyType:keytype of input gene
    ont:One of "MF", "BP", and "CC" subontologies.
    pvalueCutoff:Cutoff value of pvalue.
    pAdjustMethod:one of "holm", "hochberg", "hommel", "bonferroni", "BH", "BY", "fdr", "none"
    universe:background genes
    qvalueCutoff:qvalue cutoff
    minGSSize:minimal size of genes annotated by Ontology term for testing.
    maxGSSize:maximal size of genes annotated for testing
    readable:whether mapping gene ID to gene Name
    pool:If ont=’ALL’, whether pool 3 GO sub-ontologies
    举例:

    data(geneList, package = "DOSE")
    de <- names(geneList)[1:100]
    yy <- enrichGO(de, 'org.Hs.eg.db', ont="BP", pvalueCutoff=0.01)
    head(yy)
    

     

    2)KEGG Enrichment(KEGG富集分析)

    描述:Analysis of a gene set. Given a vector of genes, this function will return the enrichment KEGG categories with FDR control.
    语法:
    enrichKEGG(gene, organism = "hsa", keyType = "kegg", pvalueCutoff = 0.05,
               pAdjustMethod = "BH", universe, minGSSize = 10, maxGSSize = 500,
               qvalueCutoff = 0.2, use_internal_data = FALSE)

    参数:

    gene:a vector of entrez gene id.
    organism:supported organism listed in ’http://www.genome.jp/kegg/catalog/org_list.html’
    keyType:one of "kegg", ’ncbi-geneid’, ’ncib-proteinid’ and ’uniprot’
    pvalueCutoff:Cutoff value of pvalue.
    pAdjustMethod:one of "holm", "hochberg", "hommel", "bonferroni", "BH", "BY", "fdr", "none"
    universe:background genes
    minGSSize:minimal size of genes annotated by Ontology term for testing.
    maxGSSize:maximal size of genes annotated for testing
    qvalueCutoff:qvalue cutoff
    use_internal_data
    logical:use KEGG.db or latest online KEGG data
    举例:

    data(geneList, package='DOSE')
    de <- names(geneList)[1:100]
    yy <- enrichKEGG(de, pvalueCutoff=0.01)
    head(yy)
    

     

     

     

  • 相关阅读:
    海量图片曝光百度新家“搜索框”大厦
    您玩儿转手机通讯录了吗?
    这是给开发者的弥天大谎还是至理名言?
    互联网创业,不要让经验挡住你前进的道路
    永远不要去请示是否应该整理一下你的代码
    LinkedIn开放软件平台 开发者可集成其技术
    马云建新"淘宝" 借传统媒体补课线下消费群
    乔布斯的五大魔法
    全能 Google 服务提醒软件,GoogSysTray
    Twitter用户偏好新闻 Facebook用户更喜欢科技
  • 原文地址:https://www.cnblogs.com/djx571/p/9638880.html
Copyright © 2011-2022 走看看