python - 将最后一次迭代存储在循环中(python)

Question

所以我有 2 个变量 mean 和 service 由以下值组成:

mean=[5.76,6.5,7,5,4.5,6,9,2.5,3,5.5]
limit=[90,80,95,96,98,85,82,75,90,91]

现在，我对变量 mean 应用了 poisson 函数:

from scipy.stats import poisson
for m in mean:
    r_v=poisson(m)
    r_v1.append(r_v)

接下来我需要做的是获取所有均值(出现次数从0开始)的pmfs之和，并将它们与limit中的相应值进行比较。我需要获取pmfs之和大于或等于limit中对应值的出现次数。例如，对于均值 =5.76，(pmf(0)+pmf(1)...+pmf(9))*100=93.2，在极限变量中超过 90 的次数有 9 次。 我运行了以下似乎工作正常的代码:

for r,s in zip(r_v1, limit):
        l=[]
        z=0
        while True:
            y=r.pmf(z)
            l.append(y)
            z+=1
            if round(sum(l),2)*100>=s:
                break
Store=%z

问题是，当我尝试存储最后一次出现的 z 时 round(sum(l),2)*100>=s，我收到以下错误:

ERROR:root:Line magic function `%z` not found.

有人可以帮我纠正代码中的错误吗？非常感谢。

Answer 1

% 是在 IPython 的交互式 shell 中调用魔法线函数的前缀，并且没有 z 线魔法函数，因此你会得到这个错误。

在 IPython 中键入 %quickref，您将获得可用魔术函数或 readthedocs 的概览。 .

---

你的问题是你在一个步骤中尝试了太多。尝试获得一次输入的正确结果，然后在另一步中累积:

def calc(mean_, limit_):

    def inner(m, s):
        r = poisson(m)
        l = []
        z = 0
        while True:
            y = r.pmf(z)
            l.append(y)
            z += 1
            if round(sum(l), 2) * 100 >= s:
                return z, l

    return [inner(*args) for args in zip(mean_, limit_)]

调用 calc(mean, limit) 给出:

[(10,
  [0.0031511115984444414,
   0.018150402807039979,
   0.052273160084275148,
   0.10036446736180832,
   0.14452483300100394,
   0.16649260761715656,
   0.15983290331247035,
   0.13151964615426115,
   0.094694145231068075,
   0.060604252947883679]),
 (10,
  [0.0015034391929775724,
   0.009772354754354215,
   0.031760152951651209,
   0.068813664728577653,
   0.11182220518393866,
   0.14536886673912017,
   0.1574829389673803,
   0.14623415761256733,
   0.11881525306021086,
   0.08581101609904139]),
 (12,
  [0.00091188196555451624,
   0.0063831737588816145,
   0.022341108156085643,
   0.052129252364199796,
   0.091226191637349643,
   0.1277166682922895,
   0.14900277967433773,
   0.14900277967433773,
   0.1303774322150455,
   0.10140466950059107,
   0.070983268650413558,
   0.045171170959354162]),
 (10,
  [0.006737946999085467,
   0.033689734995427337,
   0.084224337488568321,
   0.1403738958142805,
   0.17546736976785063,
   0.17546736976785068,
   0.1462228081398754,
   0.10444486295705395,
   0.065278039348158651,
   0.036265577415643714]),
 (10,
  [0.011108996538242306,
   0.049990484422090385,
   0.11247858994970336,
   0.168717884924555,
   0.18980762054012446,
   0.17082685848611215,
   0.1281201438645839,
   0.082362949627232548,
   0.046329159165318316,
   0.0231645795826592]),
 (9,
  [0.0024787521766663585,
   0.014872513059998144,
   0.044617539179994441,
   0.089235078359988937,
   0.13385261753998332,
   0.16062314104797995,
   0.16062314104798009,
   0.13767697804112569,
   0.10325773353084421]),
 (13,
  [0.00012340980408667956,
   0.0011106882367801166,
   0.0049980970655105232,
   0.014994291196531574,
   0.033737155192196056,
   0.06072687934595293,
   0.091090319018929264,
   0.1171161244529091,
   0.13175564000952278,
   0.13175564000952278,
   0.11858007600857066,
   0.097020062188830414,
   0.072765046641622894]),
 (4,
  [0.0820849986238988,
   0.20521249655974699,
   0.25651562069968376,
   0.21376301724973648]),
 (6,
  [0.049787068367863944,
   0.14936120510359185,
   0.22404180765538775,
   0.22404180765538775,
   0.16803135574154085,
   0.10081881344492458]),
 (10,
  [0.0040867714384640666,
   0.02247724291155237,
   0.06181241800676901,
   0.1133227663457432,
   0.15581880372539689,
   0.17140068409793663,
   0.15711729375644187,
   0.12344930223720431,
   0.084871395288077939,
   0.051865852676047694])]