爱吱声

标题: 字符串匹配 [打印本页]
作者: 喜欢喝冰茶    时间: 2013-10-7 04:13
标题: 字符串匹配
本帖最后由 橡树村 于 2013-10-15 15:48 编辑 , n" z; y+ k" j* w2 g
6 k) ^. @  B( A( v% }+ }
字符串匹配，string match，这个是计算机里面常见的问题，例如：
, F7 L% q' I: U4 o: N/ ~  U
& ?, E* j6 u1 ostring1: TACGGCATGGCTATCGTAGCTAG  H/ u) f8 @5 W' N

/ ~- s2 Z) e! d; M4 h! R0 d# Bstring2: GCTAT8 [; z% Q: u9 c1 }3 u& }/ B7 Y0 o

; q6 L7 j$ p4 s% y9 s0 V要求在string1里找到string2的位置，如果存在多个的话，都要找出来。) b. b/ W; T7 M- {, i! u8 c
3 c0 c, J5 v1 R
可以自己估计一下时间复杂度，真实的例子是，String1长达3billion，或者6个billion。string2长约一、二百，但是数目可以是以billion计的。+ n' d; ~% j4 X4 S3 H( T
# Q3 @6 M2 g& B
先扛着。
作者: 晨枫    时间: 2013-10-7 05:16
怎么感觉这个和基因研究有关？识别基因里的某种特征？
作者: MacArthur    时间: 2013-10-7 06:17
晨枫 发表于 2013-10-6 16:16 ' J  W0 X  L# H$ C
怎么感觉这个和基因研究有关？识别基因里的某种特征？

) J6 |0 c6 C6 F5 f) Y这就是基因BLAST算法的标准定义嘛。。。
) Z0 K9 \. _7 k( H4 f/ m5 l8 B) g
6 e" h( M* q% q$ BGoogle ＂Blast算法＂，一堆开源程序。, F! |' ]7 B, g! `  |
作者: 晨枫    时间: 2013-10-7 07:21
MacArthur 发表于 2013-10-6 16:17
/ `  X/ c  _! w% ]  Q这就是基因BLAST算法的标准定义嘛。。。: P/ b# ^) ~9 s* C/ D' B+ c( F! x
1 w  L. V7 x% `7 j
Google ＂Blast算法＂，一堆开源程序。

) v5 {4 |; B4 @7 w莫非麦帅也是干这一行的？不然怎么会这么熟悉？
作者: 假如十八    时间: 2013-10-7 07:28
晨枫 发表于 2013-10-7 07:21
$ ]+ c+ a  h( {1 u2 t9 V莫非麦帅也是干这一行的？不然怎么会这么熟悉？

% R; S+ E7 T& |' Z, R5 c) G+ W4 X这都不用动手写吧？正则表达式很多语言都支持啊
作者: 晨枫    时间: 2013-10-7 07:29
假如十八 发表于 2013-10-6 17:28 # x! b1 |: h+ |$ t; h2 g3 A
这都不用动手写吧？正则表达式很多语言都支持啊
2 e8 x! F6 D2 s
有可能是普通的bubble sort和quick sort之间的差别？也就是说，普通算法的效率不适合特别大的海量比较计算？
作者: 冰蚁    时间: 2013-10-7 07:33
老兵有个群组，软件人家。
作者: 巴巴爸爸    时间: 2013-10-7 08:38
MacArthur 发表于 2013-10-7 06:17
* ~- [) j" u) a' x这就是基因BLAST算法的标准定义嘛。。。
7 U& K( ?) l& d" e- r3 X/ j9 K" F, e' J. A& @+ L- U
Google ＂Blast算法＂，一堆开源程序。
( V& ?2 u/ x$ n* q% V
blast......恍若隔世，你不说我都忘了曾经有那么一段时间天天跟这个较劲。。。。
作者: 喜欢喝冰茶    时间: 2013-10-7 09:01
最野蛮也是最简单的办法：一个一个比。( c; W8 u, \8 L

9 O( y5 j' i! _8 j4 g
string1: TACGGCATGGCTATCGTAGCTAG
9 W4 j1 ~( h; J6 m6 G/ e9 c) F) x6 m/ R+ G1 i, j
string2: GCTAT
* e  k% e% ]2 z7 u$ W7 ~1 ~5 e: v4 C
要求在string1里找到string2的位置，如果存在多个的话，都要找出来。

* \" i$ W8 w! X/ E  ?% }- b9 v; f% e: N2 M# j
拿string2和string1比，至少需要比string1的长度减去string2的长度再加1次。在实际应用中，如果string1的长度是10^9，而string2只有一二百，那么做一次基本上就是比10^9次。当然如果很幸运，string2在string1开始的地方，那一次就够了。所以平均来说，要比10^9／2次，也就是O(10^9)。
  r" K+ w) g1 f4 ^1 C( L2 T! z) C1 x1 S
但是如果实际情况中，有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年！！！人都死几次还没比完，太郁闷了，所以不可接受。
" E3 V1 e' ^3 z
) N* X: S; Q' W5 {) e5 K& G0 q那如果是这个样子9 y3 G" l; V3 t

2 O, k. C+ U9 {+ e8 I5 J( Y/ ^( k" Rstring1: AAAAAATTTTCCCCCGGGTTTTAAAACCCCCCGG
! e* h% ^, \6 O. _5 `' Ustring2: TTAAA* C  `: p3 p1 E
1 N% D' D1 q" W9 h% R
是不是会快很多？
' j% ]9 ]$ e/ I" d* j$ Z* O! b7 D+ ^! H: }3 c2 O3 O# H4 @; f
继续扛。
作者: qiuwen777    时间: 2013-10-7 09:47
喜欢喝冰茶 发表于 2013-10-7 09:01 8 D& M6 T$ ?/ C6 k
最野蛮也是最简单的办法：一个一个比。

; `8 `" N* @+ @/ @4 W$ g7 ^" ~7 y帮你顶一下。很多问题，看起来很简单，但是在数据量大的情况下，根本不是那么回事。
3 b3 ?" Q) ^  S- n) q" t- a这可以产生很多计算机专家。
1 ~' o2 S$ H. S# Z) i8 ?% @对于算法的问题，一般我都躲着。
作者: 喜欢喝冰茶    时间: 2013-10-7 10:35
本帖最后由 喜欢喝冰茶 于 2013-10-6 20:40 编辑 $ I- |4 {$ ]6 q7 @
如果两个字符串是这个样子  }% X" T3 P: }! i4 ^+ {) @4 Y+ {* M1 S
8 B! r" [& T, |' z
string1: AAAAAATTTTCCCCCGGGTTTTAAAACCCCCCGG& I$ L; D; w7 m" @  t2 s
string2: TTAAA

4 j, a: H* D4 H" z
- O; [! u  x; _! [& E, q当然要省很多时间，因为不需要对string1一个一个比了！！！
( @0 s% P/ a5 R7 S3 m" P- ^& s2 b- I! Z! d7 i/ p1 E
string1可以写成：& R2 G! }' c7 y$ ~
长度 字符 起始位置
. w5 J* T' g% h0 R% y9 Z6      A     1
4 p* o# o$ N% a: i1 L- H  H0 T6 h4      T     7
4 O1 I9 O# B, ^- P5      C     11" u0 u: Z2 [3 [$ y' @/ r: q; O
3      G     16
) }( E# y- [7 N, Y2 h4      T     19
. C# a6 h2 b2 C4 u! T......
/ X* L' z; Z3 G' o; _6 k2 n! J6 t* {, y, t0 Y7 w- c
所以当用string2去比的时候，一开始根本就不用考虑字符为A，G，C的行，因为string2开始是T。因此在这个例子中，不需要去检查string1的每个位置，而是非常有限的几个位置，所以可以省很多时间。, ]& ?8 B7 M% P: \8 C1 ]9 W

. N: q/ n/ q$ |7 @2 D那么如果存在一种这样的转换方法能够将主贴里的字符串转化成这种，势必会省很多时间。有这样一种方法吗？哪里去找？
2 p* l/ ?: [1 T. o0 n! z' b) P; v) {2 I; j) @
如果你是有心人，你觉得这个东西最常用在哪里？
! N7 \( Y4 H5 j2 c
" j# R; L4 ]7 @) [: O2 v) T- Y0 \8 R对了，文件压缩。% ~* c/ ~/ `- i! h6 w9 O7 [9 F

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

2 E2 p7 R# {' E# q那么是否就十全十美，万事大吉了呢？这个需要从实际的具体需求来看。
% y, x* p6 C; g! {, j8 @
6 V. k8 L5 c% F' u* W扛吧，没什么好说的。
5 O! f% E) ^. C3 t9 ^* t  m# ]. ?& \+ B) D5 r* B) ~; _
作者: 喜欢喝冰茶    时间: 2013-10-7 10:38
MacArthur 发表于 2013-10-6 16:17
8 l+ u3 U$ S) u* _2 s8 U这就是基因BLAST算法的标准定义嘛。。。6 w5 g1 k9 G; n* s. ?

  U$ H- O. A% kGoogle ＂Blast算法＂，一堆开源程序。

' F( d' @5 Z' M, {blast并非只是解决基因的问题，蛋白同样适用，只不过相对于DNA／RNA而言，蛋白质要复杂得多。
作者: MacArthur    时间: 2013-10-7 11:35
喜欢喝冰茶 发表于 2013-10-6 21:38 ' U" O1 o0 j; n3 T0 L- F9 K" o6 I( C
blast并非只是解决基因的问题，蛋白同样适用，只不过相对于DNA／RNA而言，蛋白质要复杂得多。 ...
4 E* j( Z3 B3 N- _, ~3 ]
不明觉厉。。。 上Billion字节的操作，想想就头大。。。 这个规模是不是得上并行计算了？$ s" w3 k0 I/ A; x
作者: 喜欢喝冰茶    时间: 2013-10-7 11:45
MacArthur 发表于 2013-10-6 21:35
/ Q5 _! C7 n( h/ g/ S+ s不明觉厉。。。 上Billion字节的操作，想想就头大。。。 这个规模是不是得上并行计算了？& K$ b% I/ l( X$ W
...

0 P- c. U( ?# v* kblast好像支持吧，主要不是分布式的问题，而是blast使用的算法对蛋白的分数的定义是很有讲究的，这些分数的定义大概要涉及到另外一个和进化相关的计算领域。
作者: 老芒    时间: 2013-10-7 11:47
本帖最后由 老芒 于 2013-10-8 00:21 编辑
4 [6 `. ?: }$ z! H4 }2 L. J. Y$ h+ ~6 `% w
这样的短贴适合聊天版。
作者: martian    时间: 2013-10-7 12:27
喜欢喝冰茶 发表于 2013-10-7 10:35
1 u$ n# e, _! |6 b当然要省很多时间，因为不需要对string1一个一个比了！！！
9 [& m' p) O4 p. n+ [) E# t0 [1 d3 ~$ o* y: Y, q9 J6 v
string1可以写成：

" [; g/ z- j# _7 E是要找pattern: Q+ S3 t1 q1 m* \3 @6 f! G" c
& D& ~* ~8 W- U
Ensembl就用了这种数据压缩方法，0 x9 |7 Z1 A7 |9 y* u" ?6 Y
AAATTCCGG
2 g$ G. s9 K0 M. @A3T2C2G29 T- Z% |7 v+ V
) l  }. P+ |8 G+ g4 D7 H9 C
把string1，分成若干片段，多进程查找。
作者: 喜欢喝冰茶    时间: 2013-10-7 21:57
本帖最后由 喜欢喝冰茶 于 2013-10-7 12:05 编辑 ; T. |% z9 |& ]

& l+ ?) [  Z6 R/ T已经有同校指出了blast，一个用于对生物大分子的测序程序。这东西正式诞生于1990年，同年人类基因组计划启动。它不仅可以用于DNA/RNA的测序，还可以对蛋白测序。六年前在河里曾经写过一个搭积木的小玩意儿，其中涉及到利用计算方法来预测实验中很难测量的蛋白质空间结构，最有效的计算方法就是同源模型。在这个模型里面，我们允许寻找的字符串,string2可以不用严格的和string1，也就是模板序列，匹配的。换句话说，就是字符串匹配是容错的。具体到主贴里提到的真实问题的挑战，就是模板字符串是人类的基因组（genome），单链长达3billion结构单位，也就是3billion的字符长度，而DNA是双链的，也就是6B的长度。想一想，每个人都是unique的，所以即使都是“健康人”，每个人的基因组都不会完全一样，因此，这种字符串的匹配一定是允许错误的，否则的话，如果大家都一样，即使算法速度再慢都不是问题，因为只做一次，从时间复杂度上，就是0。
' T3 ^3 W" p* [, V* M/ ]) R! z% g
那么如果考虑容错的匹配，基于bwt方法的算法将面临一个巨大的挑战，因为BWT是严格转换的，可是容错的实际需要，就要考虑转换前和转换后的字符串的错误问题，这会使的问题非常复杂。因此，现在的解决方案就是，当我们需要更多的关注容错的问题时，我们仍然使用传统的，也就是blast所使用的Smith－Waterman算法。具体的细节，感兴趣的可以google／wiki，基本来讲，这是一种被称之为动态编程的计算机算法，可以很好的考虑容错问题，诸如substitution、insertion、deletion或者indel（也就是insertion和deletion的混合体），而这些“错误”，则广泛的存在于病人的基因组中，特别是癌症基因组。当然，现在所使用的smith－waterman基于的方法都是修改了的，主要是计算机算法上的加速，诸如hash的SW算法。
+ E0 {7 ]- e- n$ n9 B* k
5 o7 \6 C) |+ |# q( W3 y2 R( A但是，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序列的检测成为新生儿的常规检查。
作者: 喜欢喝冰茶    时间: 2013-10-9 10:43
本帖最后由 喜欢喝冰茶 于 2013-10-10 15:20 编辑
( T$ Y9 ?* Q$ Z
6 G. t+ I/ a( k6 A6 T这个帖子里的东西，看起来似乎是一个简单的计算问题，但是却很可能是一场改变人类健康革命的基础。正是由于高速有效的匹配算法的发展，才使得新的测序技术得以广泛应用，不仅是基础研究，而更重要的是临床应用，这里简单的提到一些。等有空的时候，写一些有关基于genome和tanscriptome的技术以及实际应用，像分子诊断技术，这将对癌症的诊断和预期提供巨大的帮助。感兴趣的同学，可以看看AML，也就是急性骨髓性白血病的亚型分类标准，特别是WHO的新标准，还有预期，是否能看懂。. |5 V' z# r: b2 R  w2 m, u$ C& }% t
作者: chalet    时间: 2013-10-11 12:06
喜欢喝冰茶 发表于 2013-10-9 10:43
! c9 d, q$ F7 q5 w这个帖子里的东西，看起来似乎是一个简单的计算问题，但是却很可能是一场改变人类健康革命的基础。正是由于 ...

+ R$ A$ C# {; ~欢迎欢迎。
% h/ K9 E, f4 K5 s' A( p( G: }2 m9 G/ T# Y* @! k
如果说硬要分类，这个应该还是属于生物信息学领域吧，在科教学圃应该很对口啊。
6 P) j2 Q8 G2 t/ g9 }( i, p6 R! l6 i5 S( D1 F$ K
使用NCBI网站上的BLAST都是10多年前的事情了，一眨眼人都老了。。。
8 K% O5 @5 B/ i7 V$ Q/ F2 G/ x% q) n0 v0 u# o
关于你下面会介绍的内容，非常期待。很久没有学习这方面的知识了，落伍了，亟待补课
5 v) O7 @  s9 ]* d2 w; O% N8 a! M
作者: 喜欢喝冰茶    时间: 2013-10-11 23:06
chalet 发表于 2013-10-10 22:06 5 \. F; L5 q( k  A
欢迎欢迎。% L- F+ t: N: e; u% \0 p5 r
) d4 _* C# j3 s1 c
如果说硬要分类，这个应该还是属于生物信息学领域吧，在科教学圃应该很对口啊。
0 Z8 o5 P, g' @" v% y/ B
握握手，看起来也是生物计算的啊，现在在做什么？; X5 u6 K7 C) X5 I0 m

( F0 V: l8 H* u8 t2 W" _* W0 e还没想好怎么写，涉及的范围得控制一下，要不太大了，而且需要用些篇幅介绍一下疾病得复杂性，像癌症，这个得加好多分子遗传方面的科普，否则的话会很难理解分子诊断方面所面临的挑战。
. ^* }: J& {0 K% P4 P2 Z4 F$ Q; x" e$ d5 B) j* L- N  P
昨天刚在八卦里贴了一个，http://www.aswetalk.org/bbs/thread-25733-1-1.html，这是一个很好的分子诊断的例子，并且符合现代医学的发展趋势。以后尽可能的在八卦里发一些小东西，最后写的东西应该都能用的上，慢慢来吧。
作者: chalet    时间: 2013-10-12 10:29
喜欢喝冰茶 发表于 2013-10-11 23:06 % C! n  g  b1 Q. k  D
握握手，看起来也是生物计算的啊，现在在做什么？+ s6 T/ \) Q) O2 n

8 e+ E) @  @  j- R  \还没想好怎么写，涉及的范围得控制一下，要不太大了， ...

. V. `' T& }' ?' s3 e: ~4 k2 r我可不是作生物计算的，是学临床医学出身的，而数学正是我的致命伤，哭啊~~~" S' N5 l% o8 z5 |* O3 x
那是上一轮生物技术泡沫的时候，进了国内一家生物技术公司，稍微了解了点BLAST，2003年就离开那个公司，也离开这个行业了。
& V# W9 K6 h4 d+ d2 S对于用基因技术研究人体奥秘的价值，可以说我的认识也经历了起起落落的波折。90年代一直到本世纪初，不管是上课听到的还是看到的资料，都非常乐观，好像用测序把基因序列搞清楚了就能解决一切问题，像找到了终极解决方案一样。可是接触到这个行业后，才理解这个编码其实和天书有的一比，它具体在说什么，还需要非常复杂的分析才可能稍有理解。结合到临床，其实大部分常见病都是多基因病，单靠基因序列数据，人们其实所能做的极其有限。8 H+ M' ~* v1 V( t% w( r6 b
目前基因测序能够帮助到临床的情况，仍然非常局限，还局限在一些基因层面比较简单的情况下。个人之前以为，这是一个系统工程，基因测序技术和数据挖掘方法是其中一块，这块的发展程度可能已经超出其他领域（说实话我不太懂其他领域包括什么，但是起码觉得对序列的临床相关性，目前还所知甚少）。现在看到你说到已经有些实际的临床应用，非常高兴，希望看到这样的进步日积月累，逐渐改变当前有点进退维谷的医学实践。* U/ z8 o; o8 y+ Q6 }
关于AML分类标准，我先去学习下。期待你的佳作!
作者: 喜欢喝冰茶    时间: 2013-10-12 15:10
chalet 发表于 2013-10-11 20:29 & g# `& U  I0 |/ `' K# j+ D0 d. T
我可不是作生物计算的，是学临床医学出身的，而数学正是我的致命伤，哭啊~~~' C3 \) j1 Q2 Y/ N5 c, w
那是上一轮生物技术泡沫的时 ...
5 r0 z  j1 J  ]
兄弟原来是医生，幸会幸会，我的很多合作者都是MD。
2 {# L, D* G$ d1 F0 {# D. i9 X; j7 a; ~1 o% G# Q/ }
呵呵，03年太早了，HGP刚完成，那会儿还没看出个所以然来。到07年，新一代测序技术引入市场，一场“革命”就真得来了。对于癌症，人们了解的越来越多，这方面TCGA功不可没。5 B0 r) G  x' \/ N6 Y
0 F! w! e1 p6 A7 a" [1 u
事实上，大家谈论的是分子诊断，而非基因诊断。虽然大多数生物学家和医生仍然以基因作为基本单元考虑疾病，但在做生物医学信息的人看来，基因是个太大的概念，也是太不准确的概念，例如，临床常用的mRNA的表达，一般是以基因为基本单位的。但是，基因里面有exon，如果一个exon正向表达，另一个负向表达，从基因层次上可能没区别，但是其实已经涉及到了isoform或者splicing的问题，也就是蛋白的编码有可能被改变了，自然功能就会发生改变，那么这种异常就可以成为该种疾病的一种特征，也就是biomarker。8 Q3 A9 d4 E& V# s* N

, [% L! t  X1 q: B幸运的是，越来越多的人开始意识到这个问题。我给医学院代的一门课的听众就是MD student，MD fellow和一些即是医生也是faculty的同学们，大家也越来越意识到新的技术正在颠覆我们的概念。
作者: chalet    时间: 2013-10-12 17:12
喜欢喝冰茶 发表于 2013-10-12 15:10
! r' I0 Q" s0 r& ]兄弟原来是医生，幸会幸会，我的很多合作者都是MD。! M# N( a, Z) |/ o

1 h- f& `" W" _6 y3 m  |呵呵，03年太早了，HGP刚完成，那会儿还没看出个所 ...

, w( b, ]) N9 t1 x' w4 |非常赞同你说的。确实当年我去那家公司的时候，他们拿手的是cDNA表达谱芯片，后面的事实证明，这个层面的发现和病人的实际情况还隔着很远的距离，当时并没有合适的解读手段。. }: f9 Q  f8 [' F8 v' T

: W' D  K, w2 j2年前，和当时认识的、仍旧战斗在该领域的朋友聊天，他谈了一些表型基因组学的进展，也让我感受到一些新的进步。  a: ~  M* U- T8 ~. x! j

8 s; H+ B. K( ^% A0 k! E不过，对于分子诊断在临床的大规模应用，相对而言我没有那么乐观，觉得起码有2大障碍需要克服：: ~! m* f. z9 x
1. 分子层面的异常和临床异常之间的相关性如何确认。对于部分简单的疾病，可能一对一的关系会比较容易发现，但是大多疾病是属于异质性、存在多个环节异常的。如何识别出一群有效的标志物、如何明确一个个体究竟是哪些标志物的异常是真实发挥作用的，似乎都不是很容易（好像是K-RAS吧，很好的一个标志物，可是却没有治疗靶标的价值）。在这方面需要大量的临床信息和检测结果的综合分析。但是对这个阶段，似乎大部分医生的兴趣度还不够高，而光靠academic institute里面的医生似乎还不够收集足够的数据。, j+ K) e5 L  C  A! n
8 B' n+ G$ c8 N7 Q5 \8 {/ y+ b
2. 分子诊断和个体化医学时代，必然需要繁复的基因、蛋白层面的检测，和数量巨大的针对性药物，这与现有的药物开发和营销模式完全不合拍。孤儿药模式不可能成为未来药物治疗的支柱。那么医药企业是否会成为个体化医学的发展阻力呢？对此我是有疑问的。& d( r) ]; _8 ~1 I8 N4 `
% }+ R! g* Q( K  l0 E- f
总之，我相信个体化医学将是最终的方向，但是面对现有研究和诊治体系，个体化医学可能需要现有体系改变目前的游戏规则和流程，我对这一个矛盾的解决途径非常好奇。可以说，有生之年可以目睹这场大戏的上演，也是我们的幸事。
作者: 喜欢喝冰茶    时间: 2013-10-13 15:04
chalet 发表于 2013-10-12 03:12 & A, R5 @" n3 N/ ^; R( w+ S
非常赞同你说的。确实当年我去那家公司的时候，他们拿手的是cDNA表达谱芯片，后面的事实证明，这个层面的 ...

% I5 D2 n, _$ W7 u
; t3 O+ |* m, F" s& z7 a4 F2 Y; J% S( N+ J1 j( d' I- Q% g5 Z/ `
cDNA算是经典的array了，该不会在程氏公司吧~_~。现在mRNA的expression，特别是经典的几十个癌症基因的表达，仍然是对一些癌症进行相对早期诊断的手段之一，只不过受益者比较小众而已。! r: ]  A: f3 I9 Z0 z2 D

! V6 h: V5 [: g1 z8 y分子诊断在临床上的大规模应用现在确实没开始，但是这方面的工作很早就铺开了，而且很多临床机构已经参加进来。这个月cell上刚发了一篇TCGA网络的作者们有关glioblastoma的文章，还没看完，刚看了个开头，不过又看到一个朋友的名字在上面。除了像Washington Univ at St louis(它有和Johnhopkins，harvard齐名的医学院)，broad institute这样的学术机构(这个机构其实是很应用的)，如果对美国的癌症治疗了解的话，会看到很多著名的医疗机构，诸如sloan-kettering, dana-farber,md anderson, mayo clinic, fred hutchinson等等这些名闻遐迩的癌症中心的临床人员也参与进来。大量的临床信息和分子测量信息被综合考虑加以分析，事实上，TCGA的array数据是公开可以下载的。TCGA因为是NCI直接支持的，资源非常充足。无论是这次的GBM还是五月份在New England Journal of Medicine发表的AML，动辄就是几百个病人一下子上五六个平台同时测量分析，使用了当前最先进的设备和分析手段，知道了很多以前不了解的信息，这些东西势必对以后的应用提供很大的帮助。另外一方面，不仅是对第一个问题的补充，同时也算是回答了第二个问题。那就是现在数得着的药厂，都已经或者正在建立分子诊断部门，而他们的实力绝对不可小觑，只不过人家闷声发财罢了。就我所了解的情况，很多医生还是很感兴趣的，不仅是美国的，还有中国的，因为已经有越来越多的病例是依赖新技术而得到治疗的。
4 z! j0 [; E5 G- `; P2 V
- E5 e- @1 y- f, L( ^  i& }然而，诚如你所说，癌症的heterogeneous特性，使得这一领域所面临的挑战远远超过了我们的想象。在AML的研究中，recurrent 的定义是5%，可是看看有多少variants是可以被5%的病人所共享的？这东西不像别的学科，了解的越多，会越来越明朗，而是知道的越多，会越来越困惑。没办法，人自己设计的东西和自然选择出来的东西，论精巧和控制真的不能比。但是，现在的努力仍然对以后是有很大帮助的。例如现在被部分人批评的GWAS研究，不可否认那东西烧了不少钱，而结果相对有限。但是23&me能卖100块钱的疾病风险评估的基础，不就是GWAS鉴定出来的一两千个和疾病关联的SNP吗？
; \) o2 F; V( r4 F1 C3 l) X; u5 x/ p; P& Q% g0 k- R! h4 `5 @/ o
坦白的说，不要说personalized medicine了，就是diagnosis在我们这代人的有生之年都没戏，因为疾病太复杂。不要试图解决所有的问题，但是现在的工作仍然有很大的意义。例如，根据NCI和ACA的数据，美国刚刚过去的财年，确诊了大约一万五千个AML病人，但是也死掉了一万多点儿的病人，当然不都是这一万五里的。你也许知道这种癌症是现在常规手段根本无法诊断和发现的，因为癌细胞在骨髓里。所以不要指望能一下子解决它，但是如果现在的工作能够使一百个病人受益，虽然统计上没有什么意义，毕竟不到1%，但是对这一百个家庭，那又意味着什么？. }0 `- R% L$ \4 A, T

6 P+ O: t8 ^/ o! d. }/ {5 Y' I8 N9 M5 c) e! c' I. Y+ b7 `% V

: t/ w9 j1 X2 M. Z$ q( V
1 M) A  p: t" C9 x$ s6 ~
: [5 K. P  I# d% f9 i
作者: chalet    时间: 2013-10-14 11:20
喜欢喝冰茶 发表于 2013-10-13 15:04
5 O3 U; c2 ?4 ?: j& HcDNA算是经典的array了，该不会在程氏公司吧~_~。现在mRNA的expression，特别是经典的几十个癌症基因的 ...

6 j# T# Q1 F! S7 X* n$ r我对当前这个领域的研究有2个观点：
! P$ Y$ Y2 h: o* }5 o+ U% ]( s( M1. 当前这种研究思路，有点花大本钱做剥丝抽茧的小事的味道。一方面进步是实际存在的，另一方面又是相对片面和局限的，效率不高甚至低下。当然，这也是没有其他选择的，只能从每一个细节处着手，用蚂蚁啃骨头的办法一点点来。之前一个肺癌医生这么看待这个问题：靶向治疗药物每个针对的都是小众，虽然小，但是他们一个接一个在肺癌这个大盘子里切出小片来，积少成多，总有一天，靶向治疗会成为覆盖大多数人的有效方法。
% T1 J2 E/ e. }! |" K' A
4 r" y& \0 d1 R  V2. 个体化医学是方向，但是我们也不能奢望会有一整套的新理论新方法完全取代现在的治疗手段，那肯定是遥遥无期的事情。将现有治疗手段和个体化前沿研究的成果相结合，这个应该是更实际的路径。比如我感觉最近对靶向药物的期待已经不像几年前那么乐观，应该说是更加认识到它们的局限性了。那么现在这个情况下，有必要加强对传统手段的优化研究，比如通过对患者个体化特质的分析对传统化疗进行优化。最近听说有些单位已经在做患者化疗药物的AUC监测，这就是用现成手段进行个体化医学实践的例子，期待他们能脚踏实地，有所收获。
作者: 喜欢喝冰茶    时间: 2013-10-15 14:25
chalet 发表于 2013-10-13 21:20
, p+ n( f0 U2 P* f- I, ?0 f/ {我对当前这个领域的研究有2个观点：* T1 L/ ^" B+ \+ A8 R' G+ P5 I
1. 当前这种研究思路，有点花大本钱做剥丝抽茧的小事的味道。一方面 ...

8 Y/ q3 u: r" e* K! w/ z; A0 F呵呵，那你觉得什么样的思路才是花小钱办大事的呢？HGP刚开始的时候，也有人觉得是浪费钱。没有大量的片面，何来所谓的全面。这样做的动机，其实很简单，就是大量的证据表明，很多疾病，特别是癌症，大部分获得性癌症病人的DNA，RNA序列上有非常明显的异常，有些时候，连形态学上都能看出来。那么大家自然就想既然生出来是“正常的”，发病了就不正常了，所以DNA／RNA序列上的变化一定是时间的“函数”。既然现在没什么别的技术系统地可以用于癌症的早期诊断，那分析DNA/RNA异常大概是最可能也是最现实的方法。; x) G! ?$ w9 C- w3 z+ `( L8 _  C

* R9 H0 F# @$ u从来没有人，再没有确认更好的方法前，去准备取代现在的手段。而且一定要明确得是，基本上，现在得研究给出的是风险，而非确认。毕竟DNA/RNA上观测到的异常只是癌症一个因素，癌症的发展还和“环境”有很大关联的。对现在这种diagnosis方面的工作，相当一部分人存在一种误区。有点像一些实验生物的人看计算生物学一样，开始以为这玩意儿什么都能干，结果发现蛮不是那回事，然后就弃之不用。都不想想，真要是啥都能做，做实验的不早失业了？要真是没用，这门学科也早死了。整天吵吵garbage in garbage out的同学们，你们怎么知道人家input的就是garbage呢？远在这帮做实验的同学们去质疑计算是否正确之前，做biocomputing的早想到了这个问题，并且已经从实验中去寻找间接证据去确认了。把新生的东西当成一种手段就是了，不要排斥，和已有的成熟手段一起使用就行，只要能提高癌症病人的生存率，那每一种方法都有意义。事实上新技术在临床上的成功应用都是和其他手段一起使用的，毕竟这些技术不能治疗，像著名的Nic案例，真正的治疗手段仍然是骨髓移植，但是新技术对最后确定手术起了很大的作用。
& M' H# c1 C: N. o# c3 ?! J/ [) s3 t# C& M
. j$ X, [* |; f1 W" K$ C, B% s1 F0 u
作者: 一叶飞刀    时间: 2014-11-15 21:01
关于字符串匹配，应当已经解决完毕了，大概不会有更高级的算法了。
8 J/ ?. b' J. {: v& o- s1 V' W3 i6 _
. T8 P* B5 z( Y) U1 M% `  X) g6 [0 X从S中找ss简单匹配算法为用ss的第一个字母对齐S，逐一比对，直到不匹配或者完全匹配，然后ss的第一个字母对齐S的第二个字母，如此反复，算法复杂度为S和ss长度的乘积，即无论哪个字符串长度增加一倍，执行时间大体增加一倍。5 R3 _" B2 _0 _- I0 ^

& O5 [" W8 t3 F高级算法是，既然以前已经比较过了，那么后面的比较就可以利用前面失败的匹配所收集的信息，从而下次匹配时直接往前跳，比如一开始从S的第一个字符跟ss的第一个字符比对，比对到第十个字符，失配，那么由于已经进行了九个字符的匹配且能对的上，那么以前已经匹配过的九个字符就不要重复匹配了，第十个字符根据某种规则直接跟ss中的某一个字符匹配。这种算法的复杂度为S的长度。由于要匹配，至少要遍历S一次，所以这个最多在细节上改进，不会有本质的改进。
1 X( s7 B* E* Z' T4 W
7 j& C- p1 T* c, M3 m以上，任何一本数据结构教材，比如清华版的数据结构与算法，均有提及
作者: 喜欢喝冰茶    时间: 2015-1-5 21:18
一叶飞刀 发表于 2014-11-15 07:015 ?& f- X! C- G* q
关于字符串匹配，应当已经解决完毕了，大概不会有更高级的算法了。1 t/ Q2 ]9 X" Z
" R0 W9 q: J2 V, S" o
从S中找ss简单匹配算法为用ss的第一个 ...

3 A. i2 i# I. |4 R嗯，perfect match不是大问题，问题在于容错匹配，甚至有时候只有头上或者尾巴上的部分匹配。这部分在RNA seq或者DNA seq中的translocation中应用很广。



欢迎光临 爱吱声 (http://129.226.69.186/bbs/) Powered by Discuz! X3.2