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

背景知识

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) 
+ }#对所有的包进行加载
最后编辑于
©著作权归作者所有,转载或内容合作请联系作者
【社区内容提示】社区部分内容疑似由AI辅助生成,浏览时请结合常识与多方信息审慎甄别。
平台声明:文章内容(如有图片或视频亦包括在内)由作者上传并发布,文章内容仅代表作者本人观点,简书系信息发布平台,仅提供信息存储服务。

相关阅读更多精彩内容

  • 今天,完成了第一个目标的一半,其他都没有完成。然后根据分析,我是没有利用好时间,我花费了过多的时间在看视频上面。所...
    hmaccelerate阅读 143评论 0 1
  • michaeljacc阅读 159评论 0 0
  • 今天读完了兄弟下。表示最近两天因为工作量剧增,拖拉没读完的下部,终于在今天终结了。 网上有很多评论说,余华的兄弟,...
    蛋挞米糊阅读 93评论 1 1
  • 今天是西方的情人节,可不知什么时候咱们也过起了这个节,这个节也不知养活了多少花店。朋友圈里,大街上随处可见各...
    卡丽曦阅读 459评论 0 0
  • 今天迅速到某宝给嘉宝买了一双鞋,然后就下线了。没有在某宝和各类购物网站上闲逛。感觉还不错,竟然没有抓心挠肝。 今天...
    会跳舞的长颈鹿阅读 273评论 0 0

友情链接更多精彩内容