字符串匹配

喜欢喝冰茶

字符串匹配，string match，这个是计算机里面常见的问题，例如：
# g; k: `* ?4 e% W
string1: TACGGCATGGCTATCGTAGCTAG

" W; {" d0 t+ X) X3 H4 g# R6 fstring2: GCTAT

要求在string1里找到string2的位置，如果存在多个的话，都要找出来。
: v% U, g3 D: F' P
+ r9 ^8 L. h. ]+ E可以自己估计一下时间复杂度，真实的例子是，String1长达3billion，或者6个billion。string2长约一、二百，但是数目可以是以billion计的。

% }& G  F; i2 V+ N$ j5 w先扛着。

晨枫

怎么感觉这个和基因研究有关？识别基因里的某种特征？

MacArthur

晨枫 发表于 2013-10-6 16:16
怎么感觉这个和基因研究有关？识别基因里的某种特征？

' Q# Z0 J' h) ?- ~' @! S) L/ Y) s这就是基因BLAST算法的标准定义嘛。。。

Google ＂Blast算法＂，一堆开源程序。

晨枫

MacArthur 发表于 2013-10-6 16:17
& Z/ V# C, t2 a8 s) s这就是基因BLAST算法的标准定义嘛。。。

Google ＂Blast算法＂，一堆开源程序。

' }' j5 J1 t, d  Q& J莫非麦帅也是干这一行的？不然怎么会这么熟悉？

假如十八

晨枫 发表于 2013-10-7 07:21
* ^# F, I6 j8 x1 M( c5 N$ }莫非麦帅也是干这一行的？不然怎么会这么熟悉？

这都不用动手写吧？正则表达式很多语言都支持啊

晨枫

假如十八 发表于 2013-10-6 17:28
# s+ K+ B4 J$ ?5 g( P* H0 [这都不用动手写吧？正则表达式很多语言都支持啊

3 g/ K) r& z  H4 r! Q有可能是普通的bubble sort和quick sort之间的差别？也就是说，普通算法的效率不适合特别大的海量比较计算？

冰蚁

老兵有个群组，软件人家。

巴巴爸爸

MacArthur 发表于 2013-10-7 06:17
这就是基因BLAST算法的标准定义嘛。。。
3 s9 W- x7 {/ x) l
Google ＂Blast算法＂，一堆开源程序。

% Z; N+ L- u+ J+ ]- ~& L/ K- }blast......恍若隔世，你不说我都忘了曾经有那么一段时间天天跟这个较劲。。。。

喜欢喝冰茶

最野蛮也是最简单的办法：一个一个比。
+ |4 `  [( w0 t; R  l
. T, f. E* @1 W; S$ _8 o" v8 F7 F

string1: TACGGCATGGCTATCGTAGCTAG
. ]! P5 c" F2 k# m
' S2 P& V4 S# H2 Jstring2: GCTAT
& e1 o; X; R* C% |
' \1 q9 Q5 \2 o+ E1 g要求在string1里找到string2的位置，如果存在多个的话，都要找出来。

$ ^) s& @, _! t
0 `7 n. n2 e. n# u+ I) D拿string2和string1比，至少需要比string1的长度减去string2的长度再加1次。在实际应用中，如果string1的长度是10^9，而string2只有一二百，那么做一次基本上就是比10^9次。当然如果很幸运，string2在string1开始的地方，那一次就够了。所以平均来说，要比10^9／2次，也就是O(10^9)。

! l7 q, B& b* ]# b但是如果实际情况中，有10^6到10^9个string2s，那总共要比多少次？10^15到10^18次。这什么概念？不考虑所有的overhead，比一次只需一个时钟，那3G的CPU，意味着一秒可以比10^9次，要完成这样一个工作，需要10^6到10^9秒，1年＝365天 x 24小时 x 3600秒＝31Millon秒。也就是说，最短大约需要12天，最长需要30年。如果这样的操作做十次，一台CPU要算至少120天到300年！！！人都死几次还没比完，太郁闷了，所以不可接受。
$ }* ]& [9 r1 \* R
) L0 k! q4 s( y2 U3 w那如果是这个样子

string1: AAAAAATTTTCCCCCGGGTTTTAAAACCCCCCGG
# H4 c0 v& @  E3 x+ j0 f% dstring2: TTAAA
6 t7 ]/ D. f3 \6 t" s% i
是不是会快很多？

2 N! x. U" r; }. ^, s( m继续扛。

qiuwen777

喜欢喝冰茶 发表于 2013-10-7 09:01
最野蛮也是最简单的办法：一个一个比。

帮你顶一下。很多问题，看起来很简单，但是在数据量大的情况下，根本不是那么回事。
这可以产生很多计算机专家。
8 @4 i8 R: A& O9 L- ?; Q% O对于算法的问题，一般我都躲着。

喜欢喝冰茶

如果两个字符串是这个样子
4 M6 ^8 [( @( g
* N* T$ J* q' d2 A6 x/ Qstring1: AAAAAATTTTCCCCCGGGTTTTAAAACCCCCCGG
5 [) ]7 e1 O' Z6 Z# astring2: TTAAA

当然要省很多时间，因为不需要对string1一个一个比了！！！
9 u' M! M! \1 O- E  b
( a$ S6 i0 Y# m: P6 Sstring1可以写成：
长度 字符 起始位置
6      A     1
4      T     7
7 c) f; a5 ^! s  N6 a2 d/ C6 Z5      C     11
3      G     16
# q- g$ ~: R# }3 j2 s. Q+ u4      T     19
  s: }; b' Y4 l: Q  \$ G......
' y2 L9 s- u9 c# f+ M" {. W/ v
) t: x7 B5 @: f8 K* R; _. g+ w, V所以当用string2去比的时候，一开始根本就不用考虑字符为A，G，C的行，因为string2开始是T。因此在这个例子中，不需要去检查string1的每个位置，而是非常有限的几个位置，所以可以省很多时间。
/ x7 G: H% I' f5 ~
& Y: M. v( x9 a0 g3 L那么如果存在一种这样的转换方法能够将主贴里的字符串转化成这种，势必会省很多时间。有这样一种方法吗？哪里去找？

如果你是有心人，你觉得这个东西最常用在哪里？

对了，文件压缩。

事实上，真正的解决方法就是借鉴了最早用于文件压缩的一种算法，称为Burrows－Wheeler Transform，又称block-sorting compression。这是当年在DEC工作的Michael Burrows和David Wheeler发明的，所以以他们的名字命名，bzip2的压缩文件就是基于该算法的。它的转换其实很简单，如果感兴趣大家可以google／wiki（wikipedia上很详细的操作细节）上去看细节，但简单的来说，就是把一个字符串头围相接，不停的移动一位，然后排序，最后取出最后一列就行了。Burrows－Wheeler Transform的特性就是转换后的字符串相对于原始字符串含有大量的重复字符片段，所以就可以使得我们的问题变的相对快捷。

那么是否就十全十美，万事大吉了呢？这个需要从实际的具体需求来看。
& L6 G7 p4 E" ~
. o' J! ?* w" q; V扛吧，没什么好说的。

喜欢喝冰茶

MacArthur 发表于 2013-10-6 16:17
这就是基因BLAST算法的标准定义嘛。。。
( c6 J% k7 C7 ~
Google ＂Blast算法＂，一堆开源程序。

blast并非只是解决基因的问题，蛋白同样适用，只不过相对于DNA／RNA而言，蛋白质要复杂得多。

MacArthur

喜欢喝冰茶 发表于 2013-10-6 21:38
blast并非只是解决基因的问题，蛋白同样适用，只不过相对于DNA／RNA而言，蛋白质要复杂得多。 ...

不明觉厉。。。 上Billion字节的操作，想想就头大。。。 这个规模是不是得上并行计算了？

喜欢喝冰茶

MacArthur 发表于 2013-10-6 21:35
不明觉厉。。。 上Billion字节的操作，想想就头大。。。 这个规模是不是得上并行计算了？
0 v, e7 @* O1 [' V, ?; G- w ...

blast好像支持吧，主要不是分布式的问题，而是blast使用的算法对蛋白的分数的定义是很有讲究的，这些分数的定义大概要涉及到另外一个和进化相关的计算领域。

老芒

+ F2 M# L2 S1 v# C# q
& ^& x( N0 ?* k, m& h# ?+ ~" C( v这样的短贴适合聊天版。

martian

喜欢喝冰茶 发表于 2013-10-7 10:35
& K/ f, k( y- j' C1 d当然要省很多时间，因为不需要对string1一个一个比了！！！
; j2 L/ p0 l5 }- x
string1可以写成：

是要找pattern

Ensembl就用了这种数据压缩方法，
AAATTCCGG
A3T2C2G2

; n( T/ K# t& f6 Q把string1，分成若干片段，多进程查找。

喜欢喝冰茶

已经有同校指出了blast，一个用于对生物大分子的测序程序。这东西正式诞生于1990年，同年人类基因组计划启动。它不仅可以用于DNA/RNA的测序，还可以对蛋白测序。六年前在河里曾经写过一个搭积木的小玩意儿，其中涉及到利用计算方法来预测实验中很难测量的蛋白质空间结构，最有效的计算方法就是同源模型。在这个模型里面，我们允许寻找的字符串,string2可以不用严格的和string1，也就是模板序列，匹配的。换句话说，就是字符串匹配是容错的。具体到主贴里提到的真实问题的挑战，就是模板字符串是人类的基因组（genome），单链长达3billion结构单位，也就是3billion的字符长度，而DNA是双链的，也就是6B的长度。想一想，每个人都是unique的，所以即使都是“健康人”，每个人的基因组都不会完全一样，因此，这种字符串的匹配一定是允许错误的，否则的话，如果大家都一样，即使算法速度再慢都不是问题，因为只做一次，从时间复杂度上，就是0。
+ B4 r1 K3 W! S( c7 W: y" s
" {. f- u& D- N* t5 s' M8 J: k( _那么如果考虑容错的匹配，基于bwt方法的算法将面临一个巨大的挑战，因为BWT是严格转换的，可是容错的实际需要，就要考虑转换前和转换后的字符串的错误问题，这会使的问题非常复杂。因此，现在的解决方案就是，当我们需要更多的关注容错的问题时，我们仍然使用传统的，也就是blast所使用的Smith－Waterman算法。具体的细节，感兴趣的可以google／wiki，基本来讲，这是一种被称之为动态编程的计算机算法，可以很好的考虑容错问题，诸如substitution、insertion、deletion或者indel（也就是insertion和deletion的混合体），而这些“错误”，则广泛的存在于病人的基因组中，特别是癌症基因组。当然，现在所使用的smith－waterman基于的方法都是修改了的，主要是计算机算法上的加速，诸如hash的SW算法。
- g: n8 F* e3 W+ H
0 ~* G$ w1 i. P# q9 ?但是，bwt基于的方法，运算速度非常快。曾经有朋友几年前做过测试，拿上一代mac book pro，对于5万个长达100字符的输入string2s，string1是人类第一号染色体的DNA序列。blast要跑a couple of minutes，但是第一代bwt based的应用程序，在用手捋了一下头发之后，结果就出来了，当时朋友还以为自己错了。所以相对于blast所使用的算法而言，以bwt为基础的应用，诸如bowtie／bowtie2，bwa（short read version），soap（主要运行于HPC上，由位于深圳的华大基因开发，他们大概是曙光的最大用户），这些闻名遐迩的应用程序（从文章引用次数上，bowtie是最早的，发表于2009年，引用数已经接近3千，bwa大约2千五，soap约为700。引用它们的文章，不仅包括Nature，Science，Cell还有医学领域里的一些著名杂志，像新英格兰医学杂志），使得新型测量技术得以广泛的应用。如果不考虑监管、法律和伦理方面的问题，在可以预见的将来，你的mobile device里存储自己的DNA序列将不再是梦想而是现实，而人们同时希望，DNA序列的检测成为新生儿的常规检查。

喜欢喝冰茶

- p2 b/ D  [) V5 r. C这个帖子里的东西，看起来似乎是一个简单的计算问题，但是却很可能是一场改变人类健康革命的基础。正是由于高速有效的匹配算法的发展，才使得新的测序技术得以广泛应用，不仅是基础研究，而更重要的是临床应用，这里简单的提到一些。等有空的时候，写一些有关基于genome和tanscriptome的技术以及实际应用，像分子诊断技术，这将对癌症的诊断和预期提供巨大的帮助。感兴趣的同学，可以看看AML，也就是急性骨髓性白血病的亚型分类标准，特别是WHO的新标准，还有预期，是否能看懂。

chalet

喜欢喝冰茶 发表于 2013-10-9 10:43
& }9 K8 _1 o5 W4 E5 y7 M- A这个帖子里的东西，看起来似乎是一个简单的计算问题，但是却很可能是一场改变人类健康革命的基础。正是由于 ...

' v" ]' {: ]1 x4 e. W8 l' f欢迎欢迎。
* A3 j) `: L- j7 C
% r  y( X+ J( M$ l$ \) E如果说硬要分类，这个应该还是属于生物信息学领域吧，在科教学圃应该很对口啊。
* p$ p& O$ M7 v$ u4 y- y8 M
, J# R( W; J* J) V7 @使用NCBI网站上的BLAST都是10多年前的事情了，一眨眼人都老了。。。

2 k+ d. F8 ^+ p1 A  T关于你下面会介绍的内容，非常期待。很久没有学习这方面的知识了，落伍了，亟待补课

! b  D% y& m. S: Y+ e! [9 Y0 W$ [' u

喜欢喝冰茶

chalet 发表于 2013-10-10 22:06
: A1 V6 A! G2 T欢迎欢迎。

如果说硬要分类，这个应该还是属于生物信息学领域吧，在科教学圃应该很对口啊。

9 n6 Y6 f* v; k0 U. H; k, F握握手，看起来也是生物计算的啊，现在在做什么？

2 ^/ a5 s5 `0 m+ ]* L9 C还没想好怎么写，涉及的范围得控制一下，要不太大了，而且需要用些篇幅介绍一下疾病得复杂性，像癌症，这个得加好多分子遗传方面的科普，否则的话会很难理解分子诊断方面所面临的挑战。

昨天刚在八卦里贴了一个，http://www.aswetalk.org/bbs/thread-25733-1-1.html，这是一个很好的分子诊断的例子，并且符合现代医学的发展趋势。以后尽可能的在八卦里发一些小东西，最后写的东西应该都能用的上，慢慢来吧。