V2EX = way to explore
V2EX 是一个关于分享和探索的地方
现在注册
已注册用户请  登录
V2EX  ›  imlm  ›  全部回复第 2 页 / 共 2 页
回复总数  29
1  2  
分子++++
个人年前传实验数据,一百多 G,使用种子的方式,你可以参考。
2019-12-29 23:21:50 +08:00
回复了 1oNflow 创建的主题 问与答 个人基因检测产品值得购买吗?
2019-12-29 23:19:07 +08:00
回复了 1oNflow 创建的主题 问与答 个人基因检测产品值得购买吗?
可以考虑自己取样,然后去测序,拿到原始数据,然后分析的。
2019-11-26 10:12:20 +08:00
回复了 imlm 创建的主题 问与答 提取两个文件中的数据,并把完整格式存储到一个新的文件中
@crella 已经可以了,自己修改了一下,嘿嘿,谢啦。
2019-11-25 16:39:44 +08:00
回复了 imlm 创建的主题 问与答 提取两个文件中的数据,并把完整格式存储到一个新的文件中
@crella 修复之后运行,不是很明白为什么一运行直接就结束了,很快的那种,然后也没有什么文件生成或者改变...
2019-11-25 16:23:40 +08:00
回复了 imlm 创建的主题 问与答 提取两个文件中的数据,并把完整格式存储到一个新的文件中
@crella 不清楚为什么,你的两个脚本在同一个地方报错了...
syntax error, unexpected '|', expecting ')'
serials.each( |cur_serial|

syntax error,unexpected tSTRING_DEND, expecting keyword_end
syntax error,unexpected end-of-input, expecting keyword_end
serial_processor()
2019-11-25 07:24:35 +08:00
回复了 imlm 创建的主题 问与答 提取两个文件中的数据,并把完整格式存储到一个新的文件中
话说谢谢你们的帮助。
2019-11-25 07:24:10 +08:00
回复了 imlm 创建的主题 问与答 提取两个文件中的数据,并把完整格式存储到一个新的文件中
语言,这个倒是无所谓,自己用 Perl 弄了一个,不过感觉很差,而且不知道哪里错误了,反正没有能够出来。
#!/usr/bin/perl
use strict;
my %hash=();
open list,"<$ARGV[0]";
while (<list>){
chomp;
$hash{$_}=1; #先把需要的序列 id 放到哈希表中,id 是键所有值都设为 1
}
close list;

open unigene,"<$ARGV[1]";
$/=">";
<unigene>;
while (<unigene>){
chomp;
my $ID=(split /\n/,$_,2)[0]; #分割成两部分,以大于号为分割点,分成 id 和序列.
if (exists $hash{$ID}){ #用 exists 函数判断当前 id 在哈希键是否存在,存在输出 id 与序列,否则下一步
print ">$_\n";
}
else{
next;
}

}
close unigene;
1  2  
关于   ·   帮助文档   ·   博客   ·   API   ·   FAQ   ·   实用小工具   ·   3159 人在线   最高记录 6543   ·     Select Language
创意工作者们的社区
World is powered by solitude
VERSION: 3.9.8.5 · 18ms · UTC 14:08 · PVG 22:08 · LAX 07:08 · JFK 10:08
Developed with CodeLauncher
♥ Do have faith in what you're doing.