24小时热门版块排行榜    

查看: 1917  |  回复: 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 的主题更新
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[考研] 求调剂 +3 熬夜的猫头鹰 2026-03-02 3/150 2026-03-02 11:45 by 刘兵
[考研] 材料工程专硕283求调剂 5+4 ,!? 2026-03-02 6/300 2026-03-02 11:07 by 黑!在干嘛
[考研] 0856材料与化工,270求调剂 +8 YXCT 2026-03-01 9/450 2026-03-02 11:01 by 无际的草原
[考研] 一志愿郑大材料学硕298分,求调剂 +6 wsl111 2026-03-01 6/300 2026-03-02 11:00 by ydudjddnd
[考研] 材料学硕318求调剂 +7 February_Feb 2026-03-01 7/350 2026-03-02 10:57 by 无际的草原
[考研] 264求调剂 +4 巴拉巴拉根556 2026-02-28 4/200 2026-03-02 10:48 by yuchj
[考研] 材料调剂 +6 爱擦汗的可乐冰 2026-02-28 7/350 2026-03-02 10:42 by Jy?
[考研] 275求调剂 +3 L-xin? 2026-03-01 6/300 2026-03-02 10:22 by 热情沙漠
[考研] 材料学调剂 +10 提神豆沙包 2026-02-28 12/600 2026-03-02 09:26 by 李老师!
[考研] 材料工程269求调剂 +3 白刺玫 2026-03-02 3/150 2026-03-02 09:25 by 一休哥FU
[考研] 282求调剂 +3 2103240126 2026-03-02 4/200 2026-03-02 09:25 by 汪!?!
[考研] 306分材料调剂 +4 chuanzhu川烛 2026-03-01 5/250 2026-03-01 19:48 by 无际的草原
[考研] 321求调剂一志愿东北林业大学材料与化工英二数二 +4 虫虫虫虫虫7 2026-03-01 7/350 2026-03-01 16:52 by caszguilin
[考研] 285求调剂 +8 满头大汗的学生 2026-02-28 8/400 2026-03-01 16:47 by caszguilin
[基金申请] 刚录用,没有期刊号,但是在线可看的论文可以放为代表作吗 10+3 arang1 2026-03-01 3/150 2026-03-01 16:43 by babero
[考研] 303求调剂 +4 今夏不夏 2026-03-01 4/200 2026-03-01 14:46 by 嘟嘟小浣熊
[考研] 302材料工程求调剂 +4 Doleres 2026-03-01 5/250 2026-03-01 11:52 by liqiongjy
[考研] 调剂 +3 简木ChuFront 2026-02-28 3/150 2026-03-01 11:46 by 王伟要上岸啊
[考研] 寻找调剂 +4 LYidhsjabdj 2026-02-28 4/200 2026-03-01 10:56 by sunny81
[考研] 304求调剂 +3 52hz~~ 2026-02-28 5/250 2026-03-01 00:00 by 52hz~~
信息提示
请填处理意见