全部版块 我的主页
› 论坛 › 数据科学与人工智能 › 数据分析与数据科学 › R语言论坛
1644 5
2017-02-16
怎么统计不同Chromosome类型,Start_Position每100000bp内的Data_Value的平均值
QQ截图20170216212250.png
二维码

扫码加我 拉你入群

请注明:姓名-公司-职位

以便审核进群资格,未注明则拒绝

全部回复
2017-2-17 13:47:16
我理解的是你要做一个分类的统计,比如你的数据输入到R中命名为data,要统计chr1类,
mean(data$Data_Value[which(data$Start_Position<100000 & data$Chromosome=='chr1')])
二维码

扫码加我 拉你入群

请注明:姓名-公司-职位

以便审核进群资格,未注明则拒绝

2017-2-21 09:42:44
统计初学 发表于 2017-2-17 13:47
我理解的是你要做一个分类的统计,比如你的数据输入到R中命名为data,要统计chr1类,
mean(data$Data_Value ...
差不多,但是怎么做到每个chromosome类型中,start—position每100000内的data_value求平均,用R做循环太慢
二维码

扫码加我 拉你入群

请注明:姓名-公司-职位

以便审核进群资格,未注明则拒绝

2017-2-22 09:30:10
囬囬囬 发表于 2017-2-21 09:42
差不多,但是怎么做到每个chromosome类型中,start—position每100000内的data_value求平均,用R做循环太 ...
你的数据有多大
二维码

扫码加我 拉你入群

请注明:姓名-公司-职位

以便审核进群资格,未注明则拒绝

2017-2-22 20:11:06
统计初学 发表于 2017-2-22 09:30
你的数据有多大
60万行
二维码

扫码加我 拉你入群

请注明:姓名-公司-职位

以便审核进群资格,未注明则拒绝

2017-2-23 17:05:34
不需要循环的,有cut函数对位置分组创建新的因子,然后再算,下面代码不知道行不行
library(tidyverse)
step= seq(start最小值,star最大值,by=100000)
data %>%
mutate(start_10000= cut(Start_Position,step) %>%
group_by(chromosome,start_10000) %>%
summarise( mean= mean(data_value))
二维码

扫码加我 拉你入群

请注明:姓名-公司-职位

以便审核进群资格,未注明则拒绝

相关推荐
栏目导航
热门文章
推荐文章

说点什么

分享

扫码加好友,拉您进群
各岗位、行业、专业交流群