1、安装bioconductor及go分析涉及的相关包

source("http://bioconductor.org/biocLite.R")
options(BioC_mirror="http://mirrors.ustc.edu.cn/bioc/")
biocLite("DO.db", type = "source")
biocLite("BiocUpgrade")
biocLite('clusterProfiler')
biocLite("fastmatch", type = "source")
biocLite("gridExtra")
biocLite("ggraph")
biocLite("UpSetR")
biocLite("cowplot")
biocLite("ggridges")
biocLite("org.Hs.eg.db")
biocLite("topGO")
biocLite("Rgraphviz")  

在这一步经常出现以下报错:

报错类型1:Warning message: In file(filename, "r", encoding = encoding) : InternetOpenUrl failed: '操作超时'。

这类报错的解决方法就是修改镜像。R修改镜像步骤:

打开RStudio,依次点击:Tools → Packages → CRAN mirror

换镜像的原则是,哪个镜像方便安装哪个包就用哪个镜像

还有第二种修改镜像的方法就是查看R安装在哪儿,修改Rprofile.site.

R.home()

  根据上述命令指出来的路径修改镜像

查找X:\XXX\R-X.X.X\etc\ 路径下的Rprofile.site文件,用记事本的方式打开,写入如下命令:

# set a CRAN mirror
local({r <- getOption("repos")
r["CRAN"] <- "https://mirrors.ustc.edu.cn/CRAN/"
options(repos=r)})

 

报错类型2:package ‘XXX’ is not available (for R version X.X.X)

解决这类报错常见思路:1、查看安装包名字有没有写错;2、版本是太新还是太旧,一般而言,都是版本太旧导致的;3、这个安装包是不是一个依赖包,比如"cowplot"是一个依赖于bioconductor的包,就不能直接用“install.packages("cowplot")”这种方式安装,要先安装好bioconductor,再安装biocLite("cowplot");

2、go基因富集分析

如果你的原始数据为ENTREZID格式,则不需要通过bitr转化,如果不是的话,需要通过bitr转化为ENTREZID格式。如果你不知道啥是ENTREZID和SYMBOL,见下面的例子。

ENTREZID      SYMBOL                Gene name                                                                             location

1501               CTNND2                  catenin delta 2                                                                        5p15.2
154664           ABCA13                   ATP binding cassette subfamily A member 13                       7p12.3
1585              CYP11B2                  cytochrome P450 family 11 subfamily B member 2                8q24.3

go基因富集分析正式开始:

install.packages(xlsx)
install.packages(readxl)
library(xlsx)
library(readxl)
model<-read.xlsx("F:/download/38gene.xlsx",1,header=F)
model<-data.frame(model)
library(clusterProfiler)
GENEID_model = bitr(model, fromType="SYMBOL", toType="ENTREZID", OrgDb="org.Hs.eg.db")
#GENID_LIST=as.vector(t(GENEID_model))
go_model <- enrichGO(GENEID_model, OrgDb = org.Hs.eg.db, ont='ALL',pAdjustMethod = 'BH',pvalueCutoff = 0.05,
qvalueCutoff = 0.2,keyType = 'ENTREZID')
write.csv((as.data.frame(go_model)),"F:/download/GENEID_model.csv",row.names =F)
###go富集结果barplot图
barplot(go_model,showCategory=20,drop=T)
####go富集结果点图
dotplot(go_model,showCategory=50)
###绘制GO的网络关系图
go.BP <- enrichGO(go_model, OrgDb = org.Hs.eg.db, ont='CC',pAdjustMethod = 'BH', pvalueCutoff = 0.05, qvalueCutoff = 0.2, keyType = 'ENTREZID')
plotGOgraph(go.BP)
###ont='CC'也可以改为ont='BP'或ont='MF'

  

R: 修改镜像、bioconductor安装及go基因富集分析的更多相关文章

  1. 使用docker安装centos6.10镜像并安装新版gcc

    使用docker安装centos6.10镜像并安装新版gcc 环境:Linux Ubuntu 16.04.7 LTS 目录 使用docker安装centos6.10镜像并安装新版gcc 使用docke ...

  2. win10 下安装docker,创建镜像,push镜像到私有仓库,创建私有仓库,修改镜像仓库地址

    通过连接下载window docker安装文件,https://download.docker.com/win/stable/Docker%20for%20Windows%20Installer.ex ...

  3. R(二): http与R脚本通讯环境安装

    结合实际的工作环境,在开始R研究的时候,首先着手收集的就是能以Web方式发布R运行结果的基础框架,无耐的是,R一直以来常使用于个人电脑的客户端程序上,大家习惯性的下载R安装包,在自己的电脑上安装 -- ...

  4. 定制LFS镜像及安装过程

    定制LFS镜像及安装过程 http://blog.csdn.net/decload/article/details/7407698 一.定制LFS镜像        定制LFS镜像的思想是在已构建完成 ...

  5. yarn 国内加速,修改镜像源

    为什么慢 由于默认情况下执行 yarn 各种命令是去国外的 yarn 官方镜像源获取需要安装的具体软件信息,所以在不使用代理.不翻墙的情况下,从国内访问国外服务器的速度相对比较慢 可以通过以下命令快速 ...

  6. JavaScript 使用yrm修改镜像源

    安装yrm npm install -g yrm 列出当前可用的镜像源 # yrm ls npm ----- https://registry.npmjs.org/ cnpm ---- http:// ...

  7. 修改镜像文件EI.CFG

    一.EI.cfg说明 Windows 7 安装光盘中存在着 SOURCES\EI.CFG 这样一个配置文件.EI.cfg 是特定于 Windows 安装程序的配置文件,用于确定在安装过程中应该使用哪种 ...

  8. docker 部署 flask(三)高级编写及生成镜像,安装requirements.txt

    简介: 上一篇,我写了如何使用别人的docker基础镜像,生成我们的docker镜像. 也就最基本的flask,没有别的库(包)支持.连数据库支持都没有. 也就让大家了解一下怎么生成镜像而已. 本篇介 ...

  9. centos lamp/lnmp阶段复习 以后搬迁discuz论坛不需要重新安装,只需修改配置文件即可 安装wordpress 安装phpmyadmin 定时备份mysql两种方法 第二十五节课

    centos  lamp/lnmp阶段复习 以后搬迁discuz论坛不需要重新安装,只需修改配置文件即可 安装wordpress  安装phpmyadmin  定时备份mysql两种方法  第二十五节 ...

随机推荐

  1. 莫烦keras学习自修第四天【分类问题】

    1.代码实战 #!/usr/bin/env python #! _*_ coding:UTF-8 _*_ # 导入numpy import numpy as np np.random.seed(133 ...

  2. maven(win10)配置完环境变量后无法识别mvn -v命令

    第一步:http://maven.apache.org/download.cgi官网下载 第二步:把压缩包解压缩到不含中文和空格的目录下 第三步:新建MAVEN_HOME环境变量,值为maven解压缩 ...

  3. hdu-1176(dp)

    解题思路:用dp做的,dp[i][j]表示在i时刻,j点的最大馅饼.a[i][j]表示在i这个时刻j点同时掉落的馅饼: 每个点除了0和10之外,都有三种状态: 1.没有移动,这样值就为dp[i][j] ...

  4. Redis宕机的问题

    在主从模式下宕机要分为区分来看: slave从redis宕机 ​ 在Redis中从库重新启动后会自动加入到主从架构中,自动完成同步数据: ​ 如果从数据库实现了持久化,只要重新假如到主从架构中会实现增 ...

  5. 未能写入输出文件“c:\Windows\Microsoft.NET\Framework64\v4.0.30319\Temporary ASP.NET Files\root\2da42acc\ab2935

    未能写入输出文件“c:\Windows\Microsoft.NET\Framework64\v4.0.30319\Temporary ASP.NET Files\root\2da42acc\ab293 ...

  6. 洛谷P2055假期的宿舍

    题目 此题主要是考察二分图匹配,而二分图匹配最主要的就是建图,而图一般都是要分成两个部分来分,比如该题就需要先将在学校住的人和床连在一起,因为在学校住就会与一个床.然后每两个人之间假如他们相互认识就可 ...

  7. lemon special judge模板

    /* argv[1]:输入文件 argv[2]:选手输出文件 argv[3]:标准输出文件 argv[4]:单个测试点分值 argv[5]:输出最终得分的文件 argv[6]:输出错误报告的文件 */ ...

  8. 【AGC002E】Candy Piles 博弈论

    题目大意 有\(n\)堆糖果,第\(i\)堆有\(a_i\)个. 两个人轮流决策,决策分为两种: 1.选择糖果数最多的一堆糖果,并把这堆糖全吃了. 2.在每堆非空的糖果堆里拿一颗糖吃掉. 吃掉最后一颗 ...

  9. MySQL中 如何查询表名中包含某字段的表

    查询tablename 数据库中 以"_copy" 结尾的表 select table_name from information_schema.tables where tabl ...

  10. fullcalendar 日历插件3.9.0 -- 基本插件使用

    以下主要结构,直接执行即可以使用 ,仅用参考: html: <!DOCTYPE html> <html> <head> <title>test</ ...