问题遇到的现象和发生背景
尝试用pymc3做多元线性回归，这段代码来自《Python贝叶斯分析》这本书第四章，但遇到这个错误You can find the C code in this temporary file: C:\Users\33406\AppData\Local\Temp\theano_compilation_error_y1g7480r library mkl_rt is not found.可能是pymc3相关包的问题，也可能是环境配置的问题，搞几天没整出来
问题相关代码，请勿粘贴截图
代码如下：
%matplotlib inline
import pymc3 as pm
import numpy as np
import pandas as pd
import scipy.stats as stats
import matplotlib.pyplot as plt
import seaborn as sns
palette = 'muted'
sns.set_palette(palette); sns.set_color_codes(palette)
np.set_printoptions(precision=2)
pd.set_option('display.precision', 2)
np.random.seed(314)
N = 100
alpha_real = 2.5
beta_real = [0.9, 1.5]
eps_real = np.random.normal(0, 0.5, size=N)
X = np.array([np.random.normal(i, j, N) for i,j in zip([10, 2], [1, 1.5])])
X_mean = X.mean(axis=1, keepdims=True)
X_centered = X - X_mean
y = alpha_real + np.dot(beta_real, X) + eps_real
def scatter_plot(x, y):
plt.figure(figsize=(10, 10))
for idx, x_i in enumerate(x):
plt.subplot(2, 2, idx+1)
plt.scatter(x_i, y)
plt.xlabel('$x_{}$'.format(idx), fontsize=16)
plt.ylabel('$y$', rotation=0, fontsize=16)
plt.subplot(2, 2, idx+2)
plt.scatter(x[0], x[1])
plt.xlabel('$x_{}$'.format(idx-1), fontsize=16)
plt.ylabel('$x_{}$'.format(idx), rotation=0, fontsize=16)
scatter_plot(X_centered, y)
plt.savefig('B04958_04_25.png', dpi=300, figsize=(5.5, 5.5))
with pm.Model() as model_mlr:
alpha_tmp = pm.Normal('alpha_tmp', mu=0, sd=10)
beta = pm.Normal('beta', mu=0, sd=1, shape=2)
epsilon = pm.HalfCauchy('epsilon', 5)
mu = alpha_tmp + pm.math.dot(beta, X_centered)
alpha = pm.Deterministic('alpha', alpha_tmp - pm.math.dot(beta, X_mean))
y_pred = pm.Normal('y_pred', mu=mu, sd=epsilon, observed=y)
start = pm.find_MAP()
step = pm.NUTS(scaling=start)
trace_mlr = pm.sample(5000, step=step, start=start)
运行结果及报错内容
我的解答思路和尝试过的方法
试过用anaconda3环境安装mkl_rt，但显示没有这个包，试过网上多种方法无效
我想要达到的结果
用pymc3做多元线性回归
