概率分布模型

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))