生物信息学:了解 ViennaRNA RNAdistance 评分 table
Bioinformatics: Understanding ViennaRNA RNAdistance scoring table
我正在尝试比较 2 种不同的 RNA 结构预测算法的输出
(我的 Nussinov vs RNA-mfold 算法的实现)使用 RNAdistance algorithm that is part of ViennaRNA 包。
我得到 Nussinov 预测和实际 NMR 结构之间的距离得分为 38,m 倍算法和 NMR 结构之间的距离得分为 14,我正在尝试了解它是否是一个大的差异。
为了做到这一点,我必须查看 RNAdistance 使用的评分 table,并在 source code 中找到以下 score table
:
RIVATE CostMatrix UsualCost =
{
/* Null, U, P, H, B, I, M, S, E, R */
{ 0, 1, 2, 2, 2, 2, 2, 1, 1, DIST_INF}, /* Null replaced */
{ 1, 0, 1, DIST_INF, DIST_INF, DIST_INF, DIST_INF, DIST_INF, DIST_INF, DIST_INF}, /* U replaced */
{ 2, 1, 0, DIST_INF, DIST_INF, DIST_INF, DIST_INF, DIST_INF, DIST_INF, DIST_INF}, /* P replaced */
{ 2, DIST_INF, DIST_INF, 0, 2, 2, 2, DIST_INF, DIST_INF, DIST_INF}, /* H replaced */
{ 2, DIST_INF, DIST_INF, 2, 0, 1, 2, DIST_INF, DIST_INF, DIST_INF}, /* B replaced */
{ 2, DIST_INF, DIST_INF, 2, 1, 0, 2, DIST_INF, DIST_INF, DIST_INF}, /* I replaced */
{ 2, DIST_INF, DIST_INF, 2, 2, 2, 0, DIST_INF, DIST_INF, DIST_INF}, /* M replaced */
{ 1, DIST_INF, DIST_INF, DIST_INF, DIST_INF, DIST_INF, DIST_INF, 0, DIST_INF, DIST_INF}, /* S replaced */
{ 1, DIST_INF, DIST_INF, DIST_INF, DIST_INF, DIST_INF, DIST_INF, DIST_INF, 0, DIST_INF}, /* E replaced */
{ DIST_INF, DIST_INF, DIST_INF, DIST_INF, DIST_INF, DIST_INF, DIST_INF, DIST_INF, DIST_INF, 0}, /* R replaced */
};
我也找到了enter link description here下面的解释,但是我不明白table字母U,P,H,B,I, MS,E,R ?
比如分数[H,B] = 2是什么意思?
谁能解释一下table对我的项目会有很大的帮助,谢谢!
好的,答案是正确的here我错过了,张贴它以防其他人将来犯同样的错误。
U = unpaired
P = paired
H = hairpin
B = bulge
I = internal loop
M = multiloop
S = stem (or stack)
E = external loop
我正在尝试比较 2 种不同的 RNA 结构预测算法的输出 (我的 Nussinov vs RNA-mfold 算法的实现)使用 RNAdistance algorithm that is part of ViennaRNA 包。
我得到 Nussinov 预测和实际 NMR 结构之间的距离得分为 38,m 倍算法和 NMR 结构之间的距离得分为 14,我正在尝试了解它是否是一个大的差异。
为了做到这一点,我必须查看 RNAdistance 使用的评分 table,并在 source code 中找到以下 score table
:
RIVATE CostMatrix UsualCost =
{
/* Null, U, P, H, B, I, M, S, E, R */
{ 0, 1, 2, 2, 2, 2, 2, 1, 1, DIST_INF}, /* Null replaced */
{ 1, 0, 1, DIST_INF, DIST_INF, DIST_INF, DIST_INF, DIST_INF, DIST_INF, DIST_INF}, /* U replaced */
{ 2, 1, 0, DIST_INF, DIST_INF, DIST_INF, DIST_INF, DIST_INF, DIST_INF, DIST_INF}, /* P replaced */
{ 2, DIST_INF, DIST_INF, 0, 2, 2, 2, DIST_INF, DIST_INF, DIST_INF}, /* H replaced */
{ 2, DIST_INF, DIST_INF, 2, 0, 1, 2, DIST_INF, DIST_INF, DIST_INF}, /* B replaced */
{ 2, DIST_INF, DIST_INF, 2, 1, 0, 2, DIST_INF, DIST_INF, DIST_INF}, /* I replaced */
{ 2, DIST_INF, DIST_INF, 2, 2, 2, 0, DIST_INF, DIST_INF, DIST_INF}, /* M replaced */
{ 1, DIST_INF, DIST_INF, DIST_INF, DIST_INF, DIST_INF, DIST_INF, 0, DIST_INF, DIST_INF}, /* S replaced */
{ 1, DIST_INF, DIST_INF, DIST_INF, DIST_INF, DIST_INF, DIST_INF, DIST_INF, 0, DIST_INF}, /* E replaced */
{ DIST_INF, DIST_INF, DIST_INF, DIST_INF, DIST_INF, DIST_INF, DIST_INF, DIST_INF, DIST_INF, 0}, /* R replaced */
};
我也找到了enter link description here下面的解释,但是我不明白table字母U,P,H,B,I, MS,E,R ?
比如分数[H,B] = 2是什么意思?
谁能解释一下table对我的项目会有很大的帮助,谢谢!
好的,答案是正确的here我错过了,张贴它以防其他人将来犯同样的错误。
U = unpaired
P = paired
H = hairpin
B = bulge
I = internal loop
M = multiloop
S = stem (or stack)
E = external loop