24小时热门版块排行榜    

查看: 1886  |  回复: 7

lemian

木虫 (正式写手)

[求助] 如何从不同测序深度的样品中随机抽取相同数量的序列

我想请教一个问题,我用高通量测了很多样品,每个样品的reads数差别很大,我想在每个样品中随机抽取相同数量的序列进行后续分析(比如每个样品都随机抽取2万个reads)。请问用什么方法或者软件可以实现随机抽取?还请您给出具体操作步骤,最好容易一点的。
回复此楼
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

lemian

木虫 (正式写手)

顺便说一下,测序公司给我的数据就是一个OTU的丰度表。
2楼2013-10-28 23:34:29
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

lemian

木虫 (正式写手)

没有人做过吗?

[ 发自小木虫客户端 ]
3楼2013-10-29 11:54:00
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

lingonghua

铁杆木虫 (正式写手)

虽然我是菜鸟,弱弱问一句,高通量测序给出的reads不是随机的吗,那么,直接从其中剁一块出来就可以了嘛
4楼2013-11-02 17:14:40
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

lemian

木虫 (正式写手)

引用回帖:
4楼: Originally posted by lingonghua at 2013-11-02 17:14:40
虽然我是菜鸟,弱弱问一句,高通量测序给出的reads不是随机的吗,那么,直接从其中剁一块出来就可以了嘛

问题在于你如何“剁”。用手动的方法太费时间。我现在是用SPSS里面的随机抽样来做,也很费时间。
5楼2013-12-04 23:14:52
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

bullfrog325

木虫 (正式写手)

【答案】应助回帖

为什么会麻烦呢? 在每个样品测序行数的范围内随机生成2万不重复的数字, 然后把这些行的reads提出来。
6楼2013-12-04 23:24:07
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

hongcosb

铜虫 (小有名气)

已经解决了吗?我也遇到同样的问题,急急急!
7楼2014-08-18 10:53:21
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

lemian

木虫 (正式写手)

引用回帖:
7楼: Originally posted by hongcosb at 2014-08-18 10:53:21
已经解决了吗?我也遇到同样的问题,急急急!

用了很笨的方法,SPSS的随机抽取。用muthur qiime应该是可以的,比较快。

[ 发自小木虫客户端 ]
8楼2014-08-18 12:19:14
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
相关版块跳转 我要订阅楼主 lemian 的主题更新
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[考博] 售SCI一区文章,我:8 O5 51O 54,科目齐全 +4 2h7du0nuhk 2026-02-07 4/200 2026-02-07 22:24 by ioveiuleh5
[硕博家园] 售SCI一区文章,我:8 O5 51O 54,科目齐全 +3 2h7du0nuhk 2026-02-07 4/200 2026-02-07 22:11 by ioveiuleh5
[考博] 售SCI一区文章,我:8 O5 51O 54,科目齐全 +3 2h7du0nuhk 2026-02-07 5/250 2026-02-07 22:06 by ioveiuleh5
[教师之家] 售SCI一区文章,我:8 O5 51O 54,科目齐全 +3 2h7du0nuhk 2026-02-07 5/250 2026-02-07 22:04 by ioveiuleh5
[教师之家] 有院领导为了换新车,用横向课题经费买了俩车 +7 瞬息宇宙 2026-02-04 7/350 2026-02-07 21:47 by tfang
[找工作] 售SCI一区文章,我:8 O5 51O 54,科目齐全 +3 2h7du0nuhk 2026-02-07 5/250 2026-02-07 21:46 by ioveiuleh5
[公派出国] 售SCI一区文章,我:8 O5 51O 54,科目齐全 +3 2h7du0nuhk 2026-02-07 6/300 2026-02-07 21:44 by ioveiuleh5
[考博] 售SCI一区文章,我:8 O5 51O 54,科目齐全 +3 2h7du0nuhk 2026-02-07 6/300 2026-02-07 21:31 by ioveiuleh5
[教师之家] 售SCI一区文章,我:8 O5 51O 54,科目齐全 +3 2h7du0nuhk 2026-02-07 6/300 2026-02-07 21:26 by ioveiuleh5
[硕博家园] 售SCI一区文章,我:8 O5 51O 54,科目齐全 +3 2h7du0nuhk 2026-02-07 6/300 2026-02-07 21:11 by ioveiuleh5
[有机交流] 酰胺脱乙酰基 10+5 chibby 2026-02-03 12/600 2026-02-07 19:29 by 江东闲人
[硕博家园] 售SCI一区文章,我:8 O5 51O 54,科目齐全 +3 2h7du0nuhk 2026-02-07 3/150 2026-02-07 18:31 by 6ntj4f1rhx
[基金申请] 同年申请2项不同项目,第1个项目里不写第2个项目的信息,可以吗 +4 hitsdu 2026-02-06 4/200 2026-02-07 13:07 by jurkat.1640
[基金申请] 有时候真觉得大城市人没有县城人甚至个体户幸福 +9 苏东坡二世 2026-02-04 10/500 2026-02-07 12:37 by 小毛球
[考博] 天津大学招2026.09的博士生,欢迎大家推荐交流(博导是本人) +4 a793625982 2026-02-05 5/250 2026-02-07 10:57 by a793625982
[公派出国] CSC & MSCA 博洛尼亚大学能源材料课题组博士/博士后招生|MSCA经费充足、排名优 +4 雨念 2026-02-01 6/300 2026-02-06 23:32 by MelissaPon
[基金申请] 面上项目申报 +3 Tide man 2026-02-01 3/150 2026-02-05 22:56 by god_tian
[硕博家园] 博士延得我,科研能力直往上蹿 +7 偏振片 2026-02-02 7/350 2026-02-04 17:36 by 陈氏帝国
[基金申请] 面上基金申报没有其他的参与者成吗 +5 默默的小虫子 2026-01-31 5/250 2026-02-04 11:26 by 学员fOzRO9
[教师之家] 遇见不省心的家人很难过 +18 otani 2026-02-03 22/1100 2026-02-04 11:06 by tangmnt
信息提示
请填处理意见