0

我正在尝试解决http://rosalind.info/problems/list-view/上的问题，并且有一个名为 IPRB 的问题。简而言之，我们有 3 种不同的生物体 k（纯合子显性）、m（杂合子）和 n （纯合子隐性）。它要求您计算两个随机选择的个体的后代具有显性等位基因的概率。您可以从这里 ( http://rosalind.info/problems/iprb/ ) 检查它。

我的问题是；我用示例输入（2、2、2 表示 k、m、n）编写了一个工作代码，并得到了估计的输出。但是，当输入值彼此不同（即 25、20、18）时，我收到如下错误消息：

回溯（最近一次调用最后）：文件“C:/Users/mNm/PycharmProjects/PySummer/venv/Scripts/IPRB.py”，第 43 行，在 pr3 = float(pr31 + pr32 + pr33) NameError: name 'pr31'没有定义

我在网上研究了该代码的原因。似乎此错误与名称定义有关，例如使用后定义或功能外定义等。但是我检查时没有任何名称。简而言之，我不知道如何解决它，我决定在这里问。

k = 2
m = 2
n = 2

whole = k + m + n
org = [k, m, n]

我有 3 个变量和一个保存它们的列表。该输入工作正常，但...

例如，当它是这样的：

k = 25
m = 25
n = 19

whole = k + m + n
org = [k, m, n]

我得到错误：

回溯（最近一次调用最后）：文件“C:/Users/mNm/PycharmProjects/PySummer/venv/Scripts/IPRB.py”，第 43 行，在 pr3 = float(pr31 + pr32 + pr33) NameError: name 'pr31'没有定义

第 43 行在这里：

    if i == n:
        n_start = float(n / whole)
        for j in org:
            if j == k:
                pr31 = float(n_start * (k / (whole - 1)))
            if j == m:
                pr32 = float(n_start * (m / (whole - 1)) * 0.5)
            if j == n:
                pr33 = float(n_start * ((n - 1) / (whole - 1)) * 0)

43 --> pr3 = 浮动（pr31 + pr32 + pr33）

如果您想查看整个内容并进行检查，我也想粘贴我的完整代码。

k = 2
m = 2
n = 2

whole = k + m + n
org = [k, m, n]

for i in org:
    if i == k:
        k_start = float(k/whole)
        for j in org:
            if j == k:
                pr11 = float(k_start * ((k - 1) / (whole - 1)))
            if j == m:
                pr12 = float(k_start * (m / (whole - 1)))
            if j == n:
                pr13 = float(k_start * (n / (whole - 1)))

    pr1 = float(pr11 + pr12 + pr13)

    if i == m:
        m_start = float(m/whole)
        for j in org:
            if j == k:
                pr21 = float(m_start * (k / (whole - 1)))
            if j == m:
                pr22 = float(m_start * ((m - 1) / (whole - 1)) * 0.75)
            if j == n:
                pr23 = float(m_start * (n / (whole - 1)) * 0.5)

    pr2 = float(pr21 + pr22 + pr23)

    if i == n:
        n_start = float(n / whole)
        for j in org:
            if j == k:
                pr31 = float(n_start * (k / (whole - 1)))
            if j == m:
                pr32 = float(n_start * (m / (whole - 1)) * 0.5)
            if j == n:
                pr33 = float(n_start * ((n - 1) / (whole - 1)) * 0)

    pr3 = float(pr31 + pr32 + pr33)

dom_pr = float(pr1 + pr2 + pr3)

print(dom_pr)

我希望输出在 0-1 之间，因为它是一个概率，但更重要的是了解导致该错误的原因以及我做错了什么。

4

1 回答 1

1

这对我来说看起来不错，但如果您仍然面临错误，请尝试覆盖 try catch 块内的 i 循环。

try:
    k = 25
    m = 20
    n = 19
    whole = k + m + n
    org = [k, m, n]
    for i in org:
        if i == k:
        k_start = float(k/whole)
        for j in org:
            if j == k:
                pr11 = float(k_start * ((k - 1) / (whole - 1)))
            if j == m:
                pr12 = float(k_start * (m / (whole - 1)))
            if j == n:
                pr13 = float(k_start * (n / (whole - 1)))

    pr1 = float(pr11 + pr12 + pr13)

    if i == m:
        m_start = float(m/whole)
        for j in org:
            if j == k:
                pr21 = float(m_start * (k / (whole - 1)))
            if j == m:
                pr22 = float(m_start * ((m - 1) / (whole - 1)) * 0.75)
            if j == n:
                pr23 = float(m_start * (n / (whole - 1)) * 0.5)

    pr2 = float(pr21 + pr22 + pr23)

    if i == n:
        n_start = float(n / whole)
        for j in org:
            if j == k:
                pr31 = float(n_start * (k / (whole - 1)))
            if j == m:
                pr32 = float(n_start * (m / (whole - 1)) * 0.5)
            if j == n:
                pr33 = float(n_start * ((n - 1) / (whole - 1)) * 0)

                pr3 = float(pr31 + pr32 + pr33)

                dom_pr = float(pr1 + pr2 + pr3)
except Exception as e:
    dom_pr=None
    print(e)
print(dom_pr)
于 2019-08-10T16:36:11.370 回答