美文网首页生信星球培训第六十六期
生信学习小组Day6笔记-lamyusam_

生信学习小组Day6笔记-lamyusam_

作者: lamyusam_ | 来源:发表于2020-06-20 20:41 被阅读0次

R包学习

R语言的一大特点就是拥有许许多多的第三方扩展包,这些包可以从官方网站下载,但由于速度限制,需要在Rstudio上配置镜像下载源,今天主要学习了如何设置镜像以及R包的安装、加载和使用。

1.镜像设置

参考:生信星球:你还在每次配置Rstudio的下载镜像吗?

  • 初级模式-程序设置

ToolsGlobal OptionsPackages

初始配置
  • 升级模式-代码设置

利用options函数-R运行过程中的选项设置

options("repos" = c(CRAN="https://mirrors.tuna.tsinghua.edu.cn/CRAN/")) #设置CRAN的下载镜像为清华源
options(BioC_mirror="https://mirrors.ustc.edu.cn/bioc/") #设置Bioconductor的下载镜像为中科大源
options()$repos #查询CRAN镜像
options()$BioC_mirror#查询Bioconductor镜像

  • 高级模式-配置文件设置

利用file.edit函数修改R的配置文件.Rprofile

file.edit('~/.Rprofile') 

最后保存,重启即可。

2.安装与加载R包

dplyr包是R语言中一个具有强大数据处理能力的Package,今天在花花的指导下,开始初步学习dplyr包的安装、加载,并进行了五个基础函数练习。

  • 使用install.packages()函数安装
install.packages("dplyr")

Ps:一开始我没有注意到安装包的时候,Package名一定要加双引号,结果出现了报错信息...以后需要注意这一细节。

  • 使用libraryrequire函数加载
library (dplyr)
require (dplyr)

3. dplyr函数练习

以R中的内置数据集iris的简化版为示例数据集

#分别取出iris数据集中的
#第1、2、51、52、101、102这6行组成一个新的6行5列数据集test
test <- iris[c(1:2,51:52,101:102),]
如图: 示例数据集test
  • mutate()函数——新增列

## 新增一列,该列数值为Sepal.Length(第1列)和Sepal.Width(第二列)的乘积
mutate(test, new = Sepal.Length * Sepal.Width)
  • select()函数——筛选列

(1)按列号筛选
select(test,1) # 筛选第一列
select(test,c(1,5))# 筛选第1列和第5列
(2)按列名筛选
select(test, Petal.Length, Petal.Width)
  • filtert()函数——筛选行

filter(test, Species == "setosa") #筛选物种名为setosa的行

filter(test, Species == "setosa"&Sepal.Length > 5 )#在上一步的基础上,附加一个第一列数值大于5的条件

filter(test, Species %in% c("setosa","versicolor"))#筛选出包含物种名中serosa和versicolor的行。
  • arrang()函数——按列排序

arrange(test, Sepal.Length)#默认按照Sepal.Length从小到大排序
arrange(test, desc(Sepal.Length))#用desc变成从大到小排序
  • Summarise()函数结合group_by——汇总

# 按照Species分组,计算每组Sepal.Length的平均值和标准差
summarise(group_by(test, Species),mean(Sepal.Length), sd(Sepal.Length))

4. dplyr实用技能

(1)R管道: %>%
  • 快捷键 ctrl+shift+M


    基本含义

参考:
https://blog.csdn.net/luo617/article/details/83451052
简单理解就是将上一个函数的输出作为下一个函数的输入

test %>% 
  group_by(Species) %>% 
  summarise(mean(Sepal.Length), sd(Sepal.Length))

数据集→按物种分组→计算每组的Sepal.Length的平均值和标准差


(2)count()函数统计某列的unique值
count(test,Species)

dplyr处理关系数据

  • 内连inner_join()函数
  • 左连left_join()函数
  • 全连full_join()函数
  • 半连inner_join()函数
  • 反连inner_join()函数

假设其形式均为join(x,y)

  • inner_join 返回所有在y中能查找到的x的行,且包含x和+y的所有列;
  • left_join 返回所有x的行,且包含x和y的所有列,在y中没有查找到的x的行新增的列的值会以NA填充;
  • right_join 同上,只是x和y调换了一下;
  • full_join 返回所有x和y的行和列,未查找的部分同样会被NA填充;
  • anti_join 返回所有未能在y中能查找到的x的行,也只返回x的列
  • semi_join 返回所有在y中能查找到的x的行,也只返回x的列
    参考:https://zhuanlan.zhihu.com/p/27505645

构建两个数据框

#将字符串变量转换为因子,然后构建两个数据框data1,data2
options(stringsAsFactors = F)
data1 <- data.frame(x = c('b','e','f','x'), 
                    z = c("A","B","C",'D'),
                    stringsAsFactors = F)
data2 <- data.frame(x = c('a','b','c','d','e','f'), 
                    y = c(1,2,3,4,5,6),
                    stringsAsFactors = F)
示例数据框data1,data2
inner_join(data1,data2,by = "x") #内连
left_join(data1,data2,by ="x") #左连
full_join(data1,data2,by = "x")#全连
semi_join(x=data1,y=data2,by = "x")#半连接
anti_join(x=data2,data1,by = "x")#反连接

效果:

  • bind_cols( ),bind_rows( )函数——数据合并

bind_rows()函数需要两个表格列数相同,而bind_cols()函数则需要两个数据框有相同的行数

test1 <- data.frame(x = c(1,2,3,4), y = c(10,20,30,40))
test2 <- data.frame(x = c(5,6), y = c(50,60))
test3 <- data.frame(z = c(100,200,300,400))
bind_rows(test1, test2)
bind_cols(test1, test3)

最后,我查了一下这个dplyr包究竟有多少的内置函数

dplyr包内置函数
dplyr包内置函数

好可怕....这比我想象中的还要多.......

相关文章

  • 生信学习小组Day6笔记-lamyusam_

    R包学习 R语言的一大特点就是拥有许许多多的第三方扩展包,这些包可以从官方网站下载,但由于速度限制,需要在Rstu...

  • 生信学习小组Day5笔记-lamyusam_

    R语言学习 今天,了解了R中的数据类型,初步学习R语言的数据结构。 R的数据类型 向量(vector) R语言中的...

  • 生信学习小组Day1笔记-lamyusam_

    前言 本科毕业已经一周有余,每当闭上眼,过往的点点滴滴便开始如电影般清晰地在脑海中放映。回首过往的岁月,有苦也有甜...

  • 生信学习小组Day3笔记-lamyusam_

    是一个管理软件的工具,类似于“应用商店”的概念,目前已成为很好的生物软件管理工具。今天从Miniconda入手,开...

  • 生信学习小组Day4笔记-lamyusam_

    新手上路依然清晰地记得,当初为了在Windows系统中安装R语言学习环境耗费了我多少精力,走了多少弯路最终才得以正...

  • 生信学习小组Day2笔记-lamyusam_

    2020/6/16学习生物信息,Linux是必须掌握的内容正所谓:"万事开头难,后面更难..." 在得到花花给的腾...

  • 生信学习小组Day7笔记-lamyusam_

    在生信星球的学习小组里,每一天的生活可以用两个字形容——“充实”花花给我们制定了很有层次的学习计划,内容有详有略、...

  • 学习小组Day6笔记--kan

    笔记来自生信星球学习小组资料 Day6 学习内容-R包安装使用 1.思维导图镇楼 2.R包基本情况 由社区开发(d...

  • Day-1--小陈

    生信星球-小组学习-day1 ##小陈 1.关于学习方法 学习、索引、小组讨论、正确提问、笔记、总结;心得 2...

  • 生信学习小组 Day6 -- ldj

    R语言学习-R包 R包的下载,在CRAN或 Bioconductor 里,自定义的下载是在CRAN, 配置 Rs...

网友评论

    本文标题:生信学习小组Day6笔记-lamyusam_

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