【发布时间】:2021-10-23 16:41:42
【问题描述】:
我想在 python 中使用自定义的似然函数进行贝叶斯分析。下面是我的代码。
我在运行代码时收到
import pandas as pd
import numpy as np
import pymc3 as pm
from gurobipy import quicksum
df = pd.DataFrame(np.random.uniform(1, 125, 100), columns= ['Surface'])
def logp(Surface):
return (-c*quicksum(j**b for i in range(0,len(df)) for j in
range(1,Surface[i].astype(int))))
with pm.Model() as model_g:
c = pm.Uniform('c', lower=0, upper=1)
b = pm.Uniform('b', lower=0, upper=1)
y = pm.DensityDist('y', logp, observed={'Surface':df['Surface']})
trace_g = pm.sample(100, tune=100)
az.plot_trace(trace_g, var_names=['c']);
【问题讨论】:
-
如果您显示堆栈跟踪会有所帮助。错误发生在哪一行?
标签: python tensorflow bayesian pymc3