赞
踩
Matlab 中的 MK 突变检验代码可以使用现有的工具箱或第三方工具箱进行实现。
一种可行的方法是使用 Bioinformatics Toolbox。Bioinformatics Toolbox 提供了 seqpdist 函数,该函数可以计算序列之间的距离,从而可以用来实现 MK 突变检验。
以下是一段示例代码:
- % 定义序列
- seq1 = 'ATGAGA';
- seq2 = 'ATAGGA';
-
- % 计算距离
- distance = seqpdist(char({seq1; seq2}), 'Method', 'Jukes-Cantor');
-
- % 计算 P 值
- p = 1 -
Copyright © 2003-2013 www.wpsshop.cn 版权所有,并保留所有权利。