美文网首页GEO基因表达芯片数据生信分析
(一)GEO分析相关R包下载和加载——DAY2

(一)GEO分析相关R包下载和加载——DAY2

作者: 养猪场小老板 | 来源:发表于2020-01-21 20:55 被阅读0次

    背景知识

    hgu133plus2.db
    1.注释包的作用:基因芯片给出的分析结果都是针对芯片探针组的,每个探针组只有一个公司自己定义的ID,既没有序列信息,也不知道来自于该物种的哪条基因,因此需要在注释包的帮助下对芯片处理结果进行注释。
    2.注释包的组成:构建一个完整的芯片注释信息通常需要三个包,它们有统一的命名方式,例如Affymetrix公司的人类基因组表达芯片hgu133plus2(型号),其对应的三个注释包是hgu133plus2.db、hgu133plus2cdf和hgu133plus2probe。实际使用时,主要用“.db”那个包。

    代码如下

    > options("repos"="https://mirrors.ustc.edu.cn/CRAN/")#设置镜像
    > if(!require("BiocManager")) install.packages("BiocManager",update = F,ask = F)#意思是,如果没有安装BiocManager包,那就安装BiocManager。、.
    > options(BioC_mirror="https://mirrors.ustc.edu.cn/bioc/")
    > cran_packages <- c('tidyr',#包里的主要函数gather和spread数据长宽变,separate和union列的切割和合并
    +                    'tibble',#高效的显示表格数据的结构
    +                    'dplyr',#数据清理
    +                    'stringr',#一个字符串处理工具集
    +                    'ggplot2',#画图
    +                    'ggpubr',#画图
    +                    'factoextra',#主成分分析
    +                    'FactoMineR') #主成分分析
    > Biocductor_packages <- c('GEOquery',
    +                          'hgu133plus2.db',
    +                          "KEGG.db",
    +                          "limma",
    +                          "impute",
    +                          "GSEABase",
    +                          "GSVA",
    +                          "clusterProfiler",
    +                          "org.Hs.eg.db",
    +                          "preprocessCore",
    +                          "hugene10sttranscriptcluster.db",
    +                          "enrichplot")
    > for (pkg in cran_packages){
    +   if (! require(pkg,character.only=T) ) {#它只接受字符串值,它可以接受一个字符串变量;但是第一种不能识别字符串变量。
    +     install.packages(pkg,ask = F,update = F)
    +     require(pkg,character.only=T) 
    +   }
    + }#以上是向量的形式安装包
    ###############################################
    > for (pkg in Biocductor_packages){
    +   if (! require(pkg,character.only=T) ) {
    +     BiocManager::install(pkg,ask = F,update = F)
    +     require(pkg,character.only=T) 
    +   }
    + }#cran和bioconductor安装包的形式不一样
    ##################################################
    > #前面的所有提示和报错都先不要管。主要看这里
    > for (pkg in c(Biocductor_packages,cran_packages)){
    +   require(pkg,character.only=T) 
    + }#对所有的包进行加载
    

    相关文章

      网友评论

        本文标题:(一)GEO分析相关R包下载和加载——DAY2

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