首页
学习
活动
专区
工具
TVP
发布
精选内容/技术社群/优惠产品,尽在小程序
立即前往
您找到你想要的搜索结果了吗?
是的
没有找到

Tidyverse补充

Tidyverse补充 sunqi 2020/8/13 概述 休息了几天,罪过 tidyverse中长款数据转换函数,类比于之前reshape2包中meltdcast函数 代码 rm(list=...# 尤其是ggplot函数,上述宽数据格式无法满足绘图需要 # 涉及分组绘图 # 对于type1type2 # 因此需要长款转换 # 需要函数 # pivot_longer 转换长 # pivot_wider...转换宽 long <- pivot_longer(df, 2:3, names_to = "type",#用于显示变量名字 values_to = "value"#用于显示值名字 ) long...( cols = -Day, # 这里去掉day,其实2:3是一个意思 names_to = "type", values_to = "value" ) %>% head() # 此时间就可以直接绘图...# 长数据转换宽数据 # 此时又回到了之前数据 long %>% pivot_wider( names_from = "type", values_from = "value" ) %>%

60120

表达矩阵转换为数据框画图

主要介绍使用pivot_longer进行长宽数据转换,这两个函数都是来自于tidyr包 问题背景 现在有一个表达矩阵,要画箱线图 但是,上面表格不满足向ggplot2画箱线图函数传递参数需求,要变换成数据框把所有数字变成一列传递给...首先行列转置 把原来行名变成第一列 把原来列名变成第二列 就变成数据框形式了。也就是把宽数据变成长数据。 代码如何实现?...,用于画图上不同颜色 pdat = dat%>% pivot_longer(cols = starts_with("gene"), # gene开头列 names_to...pdat = dat%>% pivot_longer(cols = 2:4, names_to = "gene", values_to...()函数直接解析 列名中含有多个变量可以用正则表达式拆分成多列 一行有多个观测 列名有重复 详见使用pivot_longerpivot_wider进行长宽数据转换-CSDN博客

7110

R 数据整理(七:使用tidyrdplyr处理数据框 2.0)

pivot_longer/pivot_wider 大部分功能是类似的,这里主要说下pivot_longer 针对下面情况功能: 我们需要 指定切分变量名随访号模式,以解决一行中有多个属性多次观测情形...,在对应 names_to 中用特殊".value" 名字表示切分出来那一部分实际是变量名,这 时不需要 values_to 选项: dwide4 %>% pivot_longer( -id, names_pattern...= "(x|y)([[:digit:]])", names_to = c(".value", "time") ) %>% knitr::kable() 对应长变宽函数有pivot_wider...长宽混合转换 有时候,需要将数据框先转换为宽列表,再转换回长列表,比如: 这个数据问题是 x, y 应该放在两列中却合并成一个了,2018 2019 应该放在一列中却分成了两列。...先合并 2018 2019 这两列,然后再拆分 x y: dlong6 %>% pivot_longer( `2018`:`2019`, names_to = "year", values_to

10.7K30

R中绘制环状聚类热图

「数据代码已经整合上传到2023VIP交流群」,加群观众老爷可自行下载,有需要朋友可关注文末介绍加入VIP交流群。...❞ 关于永久群内容说明 ❝给予长期支持我们忠实读者们一个特别待遇:凡是购买过小编2022年或2023年VIP会员文档朋友们,「将自动获得2024年及以后绘图资料代码更新,无需额外付费。」...目前这两年会员文档已累记卖出1500+,质量方面各位无需担忧**。简要概括就是只要购买任意1年会员内容,2024及后期公众号所更新绘图文档均会在已经加入会员群内分享。...>% group_by(id,Phylum) %>% summarise(across(where(is.numeric), ~ sum(.x, na.rm=TRUE))) %>% pivot_wider...(names_from = "Phylum",values_from = "Abundance") 构建表达信息 exp % pivot_longer(-id) 构建树文件 tree

22920

基础知识 | 踏实做事,不要偷懒,之前偷懒,以后都是要补回来

01 表格之间处理 上一篇文章推出【R语言】基础知识 | 为了偷懒,我不择手段!,想了想,人还是踏实一点比较好,别老想着走捷径,不然有一天会摔很惨,咱还?️...lxl=lxl %>% pivot_wider(names_from = 区域, values_from = 奖励金额, values_fill = 0) names_from: 指定列名来自哪个变量列...思路: 仓山、福清、高新、鼓楼、西湖这几个字符应该作为区域存储在一列 lxl=lxl%>%pivot_longer(-银行,names_to="区域",values_to="奖励金额",values_drop_na...=TRUE) #查看运行结果: 如何将去掉奖励金额中“0”所在行?...将0替换成缺失值NA lxl$奖励金额[lxl$奖励金额==0]<-NA #查看结果 如何将NA所在行删除,合并区域,查看往期文章~

95110

tidyverse数据清洗案例详解

一旦你有了整洁数据一些包提供整洁工具,您将花费很少时间将数据从一种表示转换到另一种,从而可以将更多时间花在分析问题上。 本文将为您提供整理数据实用介绍以及tidyr包中附带工具。...这是一个非常典型现实示例数据集。它包含冗余列,奇数变量代码许多缺失值。我们需要采取多个步骤来对其进行整理。 不是变量列汇集在一起 首先将不是变量列聚集在一起。...所包含列包括: country,iso2iso3是三个指定国家/地区变量。 year是一个变量。...男性(m)女性(f) 其余数字给出了年龄段。...stocks %>% pivot_wider(names_from = year,values_from = return) ? separate() 该函数可将字符进行分割,具体案例如上.

1.5K10

单细胞韧皮部研究代码解析3-comparison_brady.R

,很多研究学者会将以前做显微切割RNA-seq数据与自己研究内容进行比对,则可以得到在单细胞亚细胞水平上数据结果。...,主要是根据自己相关内容进行更改 # read both sheets # Brady提供是probegene对应数据集,是需要将不同表达组织进行对应 # 文章中作者选用数据链接来源:https...sheet = "LONGITUDINAL", skip = 1) %>% pivot_longer(c(-Probe, -Gene), names_to =...gene内容对两个数据集进行合并 # matrix of expression lon_matrix % filter(set == "longitudinal") %>% pivot_wider...数据集进行整合,计算了细胞与组织之间相关性系数,为鉴定细胞亚群也做了相关参考,在细胞层面亚细胞层面上都做了相关分析,也是在以前文章中没有看到内容,同时我自己对自己数据也进行了测试,是可以进行相关分析

18520

UseGalaxy.cn生信云|零代码使用Tiverse优雅地处理数据集

Dplyr Distinct keep unique rows distinct 函数用于去除数据框中重复观测,仅保留唯一观测。它可以基于指定列对数据框进行去重操作,确保每个观测都是唯一。...Dplyr Join two tables join 函数用于根据指定键将两个数据框连接起来,可以根据共同变量将数据框进行合并,支持多种连接操作,如内连接、左连接、右连接外连接等。...Dplyr Rename columns rename 函数用于重命名数据框中变量名,能够快速修改变量名称,使得数据列名更符合用户需求和习惯。...Tidyr Pivot Longer from wide pivot_longer 函数用于将宽格式数据转换为长格式数据,能够根据用户指定列将数据框中多个列整理成一对 “名-值” 对,便于进一步分析处理...Tidyr Pivot Wider from long pivot_wider 函数用于将长格式数据转换为宽格式数据,能够将数据框中一列分成多个列,根据指定列名进行展开,使得数据以更直观宽格式形式呈现

15120

R语言利用转录组基因表达矩阵做基因共表达分析学习资料推荐

参考资料链接 https://github.com/cxli233/SimpleTidy_GeneCoEx/tree/v1.0.1 提供完整示例数据代码,非常好学习材料 做基因共表达比较常用是...WGCNA那个R包,这个链接里提供代码不是用WGCNA这个R包实现,而是利用表达量数据计算不同基因之间相关性,这种方法也挺常用在论文里见过 表达量数据是来源于论文 High-resolution...dim(Metadata) library(tidyverse) Exp_table_long % rename(gene_ID = `...1`) %>% pivot_longer...Exp_table_long) Exp_table_log_wide % select(gene_ID, library, logTPM) %>% pivot_wider...) PCA_by_stage library(patchwork) PCA_by_method+PCA_by_tissue+PCA_by_tissue image.png 以上用到代码示例数据都可以在推文开头提到链接里找到

35910
领券