基因检测-外显子链接问题
2015-10-27 17:02
274 查看
人类大部分基因组序列都是被垃圾DNA序列分隔成一段段,给定一个已知的目标蛋白质和基因组序列,在该基因组序列中找出一组子字符串(候选外显子),使得其拼接(剪接)与目标蛋白质最匹配(即去掉垃圾DNA序列)。一个强力方法是寻找基因组序列与目标蛋白质序列间的所有局部相似性。若第一个取自基因组序列的子字符串展示了充分相似性于目标蛋白质,那么这个子字符串可被认为是一个推定的外显子。将推定外显子结构化为基因组序列中的赋权区间,它可用三个参数(l、r、w)来描述,l、r分别是推定的外显子的左边、右边的位置,w为其权重。权重w可反该区间是一个外显子的可能性。链是不重叠赋权区间的任一集合,一个链的总权重是该链中所有区间的权重之和。
本博客所有内容是原创,如果转载请注明来源
http://blog.csdn.net/myhaspl/
给定一个推定的外显子集,寻找非重叠的推定的外显子的一个最大集。输入:赋权区间(推定的外显子)集。输出:该集合中区间的最大链。下列算法中,Si代表图G中终止于顶点Vi的最长路线的长度,S2n是外显子链接的解EXONCHAINING(G,n)for i<-1 to2nSi<-0for i<-1 to2nif G中顶点Vi是某个区间I的右端点j<-代表区间I左端点的顶点的序号W<-区间I的权重Si,-max{Sj+w,Si-1}elseSi<-Si-1return S2n相关文章推荐
- 004.self 用来调用方法
- MyBatis学习总结(一)——MyBatis快速入门
- 设置tabbaritem的title的颜色及按钮图片
- SpringMVC与LogBack集成
- JS中的prototype
- UVa 10815 - Andy's First Dictionary【排序+去重】
- jQuery中使用ajax向后台传递数据
- 3.13伪代码
- kafka+storm初探
- 【UFLDL】多层神经网络
- Java:进制转换
- UDP协议的两个主要方法sendto和recvfrom详解
- 基因检测-剪切联配(1)
- Linux下关于热插拔硬盘的指令
- Pawn Storm网络间谍行动再度现身
- (二)MySQL InnoDB锁模型
- python异常处理
- django cookies 和 session
- iOS中SQLite的简单实用
- 在uboot中加入一键nfs启动选项,方便调试工作