生信学习小组Day6笔记-lamyusam_

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包内置函数

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

最后编辑于
©著作权归作者所有,转载或内容合作请联系作者
  • 序言:七十年代末,一起剥皮案震惊了整个滨河市,随后出现的几起案子,更是在滨河造成了极大的恐慌,老刑警刘岩,带你破解...
    沈念sama阅读 204,530评论 6 478
  • 序言:滨河连续发生了三起死亡事件,死亡现场离奇诡异,居然都是意外死亡,警方通过查阅死者的电脑和手机,发现死者居然都...
    沈念sama阅读 86,403评论 2 381
  • 文/潘晓璐 我一进店门,熙熙楼的掌柜王于贵愁眉苦脸地迎上来,“玉大人,你说我怎么就摊上这事。” “怎么了?”我有些...
    开封第一讲书人阅读 151,120评论 0 337
  • 文/不坏的土叔 我叫张陵,是天一观的道长。 经常有香客问我,道长,这世上最难降的妖魔是什么? 我笑而不...
    开封第一讲书人阅读 54,770评论 1 277
  • 正文 为了忘掉前任,我火速办了婚礼,结果婚礼上,老公的妹妹穿的比我还像新娘。我一直安慰自己,他们只是感情好,可当我...
    茶点故事阅读 63,758评论 5 367
  • 文/花漫 我一把揭开白布。 她就那样静静地躺着,像睡着了一般。 火红的嫁衣衬着肌肤如雪。 梳的纹丝不乱的头发上,一...
    开封第一讲书人阅读 48,649评论 1 281
  • 那天,我揣着相机与录音,去河边找鬼。 笑死,一个胖子当着我的面吹牛,可吹牛的内容都是我干的。 我是一名探鬼主播,决...
    沈念sama阅读 38,021评论 3 398
  • 文/苍兰香墨 我猛地睁开眼,长吁一口气:“原来是场噩梦啊……” “哼!你这毒妇竟也来了?” 一声冷哼从身侧响起,我...
    开封第一讲书人阅读 36,675评论 0 258
  • 序言:老挝万荣一对情侣失踪,失踪者是张志新(化名)和其女友刘颖,没想到半个月后,有当地人在树林里发现了一具尸体,经...
    沈念sama阅读 40,931评论 1 299
  • 正文 独居荒郊野岭守林人离奇死亡,尸身上长有42处带血的脓包…… 初始之章·张勋 以下内容为张勋视角 年9月15日...
    茶点故事阅读 35,659评论 2 321
  • 正文 我和宋清朗相恋三年,在试婚纱的时候发现自己被绿了。 大学时的朋友给我发了我未婚夫和他白月光在一起吃饭的照片。...
    茶点故事阅读 37,751评论 1 330
  • 序言:一个原本活蹦乱跳的男人离奇死亡,死状恐怖,灵堂内的尸体忽然破棺而出,到底是诈尸还是另有隐情,我是刑警宁泽,带...
    沈念sama阅读 33,410评论 4 321
  • 正文 年R本政府宣布,位于F岛的核电站,受9级特大地震影响,放射性物质发生泄漏。R本人自食恶果不足惜,却给世界环境...
    茶点故事阅读 39,004评论 3 307
  • 文/蒙蒙 一、第九天 我趴在偏房一处隐蔽的房顶上张望。 院中可真热闹,春花似锦、人声如沸。这庄子的主人今日做“春日...
    开封第一讲书人阅读 29,969评论 0 19
  • 文/苍兰香墨 我抬头看了看天上的太阳。三九已至,却和暖如春,着一层夹袄步出监牢的瞬间,已是汗流浃背。 一阵脚步声响...
    开封第一讲书人阅读 31,203评论 1 260
  • 我被黑心中介骗来泰国打工, 没想到刚下飞机就差点儿被人妖公主榨干…… 1. 我叫王不留,地道东北人。 一个月前我还...
    沈念sama阅读 45,042评论 2 350
  • 正文 我出身青楼,却偏偏与公主长得像,于是被迫代替她去往敌国和亲。 传闻我的和亲对象是个残疾皇子,可洞房花烛夜当晚...
    茶点故事阅读 42,493评论 2 343