24小时热门版块排行榜    

查看: 1919  |  回复: 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 的主题更新
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[考研] 江苏省农科院招调剂1名 +4 Qwertyuop 2026-03-01 4/200 2026-03-02 14:27 by 升格阿达
[考研] 282求调剂 +4 2103240126 2026-03-02 5/250 2026-03-02 13:45 by littlehu66
[考研] 材料学硕318求调剂 +9 February_Feb 2026-03-01 9/450 2026-03-02 13:31 by njzyff
[考研] 材料与化工328求调剂 +3 。,。,。,。i 2026-03-02 3/150 2026-03-02 13:09 by houyaoxu
[考研] 275求调剂 +4 明远求学 2026-03-01 4/200 2026-03-02 12:45 by 无际的草原
[基金申请] 成果系统访问量大,请15分钟后再尝试。由此给您造成的不便,敬请谅解。 +5 xhuama 2026-03-02 5/250 2026-03-02 12:34 by stidwellNK
[考研] 考研复试调剂,过国家线的同学都可报名 +3 黑!在干嘛 2026-02-28 4/200 2026-03-02 12:20 by 花鱼不是鱼
[考研] 281求调剂 +5 2026计算机_诚心 2026-03-01 8/400 2026-03-02 11:05 by 汪!?!
[考研] 求调剂 +8 yunziaaaaa 2026-03-01 9/450 2026-03-02 11:03 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
[考研] 264求调剂 +4 巴拉巴拉根556 2026-02-28 4/200 2026-03-02 10:48 by yuchj
[考研] 272求调剂 +6 田智友 2026-02-28 6/300 2026-03-01 21:40 by 公瑾逍遥
[考研] 299求调剂 +3 Y墨明棋妙Y 2026-02-28 5/250 2026-03-01 21:01 by tangxiaotian
[考研] 298求调剂 +6 axyz3 2026-02-28 6/300 2026-03-01 19:00 by 18137688336
[考研] 328求调剂 +3 aaadim 2026-03-01 5/250 2026-03-01 17:29 by njzyff
[考研] 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
[考研] 材料284求调剂,一志愿郑州大学英一数二专硕 +10 想上岸的土拨鼠 2026-02-28 10/500 2026-03-01 14:12 by yc258
[考研] 调剂 +3 简木ChuFront 2026-02-28 3/150 2026-03-01 11:46 by 王伟要上岸啊
信息提示
请填处理意见