系列文章目录

生信分析第一步:R语言基础应用以及数据前处理

文章目录

  • R包下载
  • 使用GEOquery包下载原始数据
  • 芯片数据读取
  • GEOquery 下载并读取数据
  • 提取GEO表达矩阵
  • 提取GEO注释信息
  • 提取GEO芯片样本及临床信息

前言

用最精简的代码展示生信分析的全部流程


一、如何一步安装生信所有R包

生信分析过程中会用到各种各样的R包。当我们在全身心投入到数据分析的工作中,突然发现有个R包尚未安装,而安装过程又费时费力还往往出现安装错误,这会很影响我们的工作效率。本节代码可以一步安装我们生信分析中遇到的95%的R包,做好前期准备工作,会让我们以后的分析畅通无阻。

首先安装BiocManager,若已安装则跳过if (!requireNamespace("BiocManager", quietly = TRUE)) { install.packages("BiocManager",ask = F, update = F)}安装stringr,若已安装则跳过if (!requireNamespace("stringr", quietly = TRUE)) { install.packages("stringr",ask = F, update = F)}去重,识别有无github格式的安装包pkgs <- unique(pkgs)pkgs2 <- pkgslogi <- stringr::str_detect(pkgs2, "/")pkgs2[logi] <- stringr::str_match(pkgs2[logi], ".*/(.*)$")[,2]安装pkgs中尚未安装的包new <- !(sapply(pkgs2, requireNamespace, quietly = T))显示需安装的包if (sum(new) > 0) {cat("pkgs to install: ", pkgs[new], "\n")} else {cat("All pkgs already installed \n")}install pkgsif(any(new)) BiocManager::install(pkgs[new], ask = F, update = F)
}需要安装的pkgs ----pkgs <- c("tidyverse", "limma", "affy", "oligo", "lumi","beadarray", "GEOquery", "simpleaffy", "gcrma", "readxl","impute", "genefilter", "pd.hugene.1.0.st.v1", "pd.hg.u133.plus.2","tkWidgets", "illuminaHumanv4.db", "AnnotationDbi", "org.Hs.eg.db","hgug4112a.db", "AgiMicroRna", "sva", "DESeq2", "edgeR","lumiHumanIDMapping", "remotes", "pheatmap", "shiny", "aggregation","tidyverse/dplyr", "limma", "hwriter", "devtools")安装pkgs中的R包 ----pkgs_in(pkgs)  

# 最后运行pkgs_in()时应无提示,表明安装成功;

# library(),不报error,提示安装成功。

# 若仍有安装不成功的包,可参考下节的方法。


生信分析之R语言常用R包一步下载相关推荐

  1. r语言python print包_python : 批量下载R语言库包

    soupR.py 代码如下# -*- coding: cp936 -*- import urllib import urllib2 import os, re from BeautifulSoup i ...

  2. r语言python0基础_万丈高楼平地起,从零开始学R语(生信分析向,R语言基础)...

    一.R语言的安装 我所做的R语言学习主要是为应对工作上的生信分析,如果你是往这个方向走,那么这篇文章很大可能会帮助你.如果你也是有python基础的话,那我写的可能会更対你口味. 不过现在看到这篇文章 ...

  3. 生信分析用python还是r_生信分析中常用的R语句

    一.前言 好难受,今天的R崩了.昨天才安装的包怎么报错了.为啥我装了这个包没反应,什么鬼怎么又乱码了.这玩泥巴(╯‵□′)╯︵┻━┻ 不行要佛系生活,好久没写简书了,还是来整理下代码吧┬-┬ノ('-' ...

  4. 找不到r低版本_R 语言与数据挖掘直播班开始招生,生信分析帮你发高分文章

    数据单薄很难支撑文章内容?数据有了不知道怎么处理作出高级的图片?这个时候需要的是生信分析--深度的数据挖掘和分析处理,可以帮助临床医生通过数据处理得到自己想要的信息,更快速地发文章. 学习哪种生信分析 ...

  5. R语言---生信分析---count转换成TPM、FPKM

    R语言---生信分析---count转换成TPM.FPKM 背景介绍 代码 0. 设置工作目录,加载需要的包 1. 读取 reads count 的数据 2. 下载基因长度的数据,并读取 3. cou ...

  6. 生信分析R语言助力作图----单基因批量相关性分析

    单基因批量相关性分析 嘻嘻嘻~~~,晚上秒变生信分析小白,一个游走在生物学和计算机变成之间的小白,享受着里面的快乐和痛苦.不停的挣扎,不停的成长,多学习,多尝试,一定会有意想不到的收获.加油!!! 首 ...

  7. 生信分析常用编程语言汇总

    0525 Cloudy 说到编程语言,其实大致可以分成两种,一种是计算机编译器语言,比方说C,C++,优点在于基于计算机硬件特点编写,可以使代码有效高速的运行,特别适合大规模复杂的计算.日常用到的程序 ...

  8. perl语言入门第七版 电子版_百迈客带您走近生信分析【入门篇】

    年末促销倒计时:59天 百迈客推出年末活动促销,发文有礼.推广有礼以及多种产品钜惠来袭,百迈客为您倾情打造科研福利,您还在等什么?快快行动起来领取您的超级奖励吧!(详情请见"决战2020!品 ...

  9. 生信分析常用软件记录

    20190727,在学习二代分析的过程中,只是根据别人已经建好的轮子照抄照搬,并不能真正理解每一步为什么要用这个软件,以及软件之间的区别.因此今天记录一些生信分析过程(主要是二代测序)中常用的软件,若 ...

最新文章

  1. Lucid Sight宣布即将推出“VR云霄飞车之星际之旅”
  2. Oracle-使用awrrpt.sql生成AWR报告
  3. 发布订阅之fanout
  4. ps怎么清屏_黑洞PS大赛刷屏!最后一张扎心了……
  5. 立足于飞秋和即时通讯软件市场
  6. 系统架构师学习笔记-系统性能评价
  7. 七人のオンラインゲーマーズ 全年龄正式汉化补丁
  8. 第3章 一切基于pom
  9. Segment Advisor
  10. java 网页 录音_Java+FlashWavRecorder实现网页录音并上传【新】
  11. java 观察者模式类图_设计模式——观察者模式
  12. android 带刻度的滑动条_Android实现滚动刻度尺效果
  13. xcode7的那些坑-“Your binary is not optimized for iPhone 5” (ITMS-90096) when submitting
  14. 第一章 操作系统引论之多道批处理系统
  15. 百度地图api根据详细地址反查坐标
  16. GJB6600-装备交互式电子技术手册
  17. 网络安全技术与黑客攻击威胁
  18. 华为的强大再次得到证明,高通为它定制功耗更低的芯片
  19. python提取word表格中数据
  20. 移动电源当中的自动检测负载方式讲解

热门文章

  1. C++——Hook教程[1]:虚函数表(VMT)Hook
  2. 计算机视觉基础——3D空间坐标点的重建(三角测量)
  3. scrapy抓斗鱼主播的图片
  4. 【小学】小学汉语拼音知识复习汇总
  5. pytorch版本下的yolov3训练实现火焰检测
  6. 基于搜索的贝叶斯网络结构学习算法-K2
  7. 关于人机智能的几点思考
  8. android 机器人gif,看完这15张GIF图 就再也不用担心机器人叛乱了
  9. 【收藏夹2020-2021】
  10. HTC Desire获取ROOT方法--HTC Desire rooting guide (now with new improved / easier root method!)