概率分布模型

Question

我试图在这个数据科学站点上运行一个用作示例的python脚本：https://www.datascience.com/blog/introduction-to-bayesian-inference-learn-data-science-tutorials

import numpy as np
from scipy.misc import factorial
import matplotlib.pyplot as plt
%matplotlib inline

plt.rcParams['figure.figsize'] = (16,7)

def likelihood(theta, n, x):
    """
    likelihood function for a binomial distribution

    n: [int] the number of experiments
    x: [int] the number of successes
    theta: [float] the proposed probability of success
    """
    return (factorial(n) / (factorial(x) * factorial(n - x))) \
            * (theta ** x) * ((1 - theta) ** (n - x))

# The number of impressions for our facebook-yellow-dress campaign
n_impressions = 10

# The number of clicks for our facebook-yellow-dress campaign
n_clicks = 7

# Observed click through rate
ctr = n_clicks / n_impressions

# 0 to 1, all possible click through rates
possible_theta_values = map(lambda x: x/100., range(100))

# Evaluate the likelihood function for possible click through rates
likelihoods = map(lambda theta: likelihood(theta, n, x)\
                                , possible_theta_values)

# Pick the best theta
mle = possible_theta_values[np.argmax(likelihoods)]

# Plot
f, ax = plt.subplots(1)
ax.plot(possible_theta_values, likelihoods)
ax.axvline(mle, linestyle = "--")
ax.set_xlabel("Theta")
ax.set_ylabel("Likelihood")
ax.grid()
ax.set_title("Likelihood of Theta for New Campaign")
plt.show()

我在用Spyder。我所犯的错误包括：

    %matplotlib inline
    ^
SyntaxError: invalid syntax

以及

mle = possible_theta_values[np.argmax(likelihoods)]
TypeError: 'map' object is not subscriptable

和

likelihoods = list(map(lambda theta: likelihood(theta, ctr, x)\
NameError: name 'x' is not defined

我尝试了一个我在这里找到的解决方案：Python map object is not subscriptable

# Evaluate the likelihood function for possible click through rates
likelihoods = list(map(lambda theta: likelihood(theta, n_clicks, ctr)\
                                , possible_theta_values))

但它给了我们错误

TypeError: object of type 'map' has no len()

请问如何修复在Spyder中运行的脚本？

Answer 1

对代码进行一些小的修改以使其正常工作。出于您提到的原因，我不得不将map可迭代对象转换为列表，以使代码工作：TypeError: 'map' object is not subscriptable。我还必须显式地定义x数组，以便在计算likelihoods时能够使用它。

# 0 to 1, all possible click through rates
possible_theta_values = list(map(lambda x: x/100., range(100)))

n = n_impressions
x = n_clicks

# Evaluate the likelihood function for possible click through rates
likelihoods = list(map(lambda theta: likelihood(theta, n, x)\
                                , possible_theta_values))

​

​