即使我不包含任何观察到的值,pymc3仍将运行并为我提供结果。这只是从先前的样本中抽取而来吗?
'''导入pymc3作为pm模型= pm.Model()
带型号:
# Define the prior of the parameter lambda.
lam = pm.Gamma('lambda', alpha=3.5, beta=2)
带型号:trace = pm.sample(draws = 20,chains = 3)
pm.traceplot(trace)'''
是的,您只是从先前的样本中取样。如果需要,可以通过将样本绘制为直方图或核密度估计值并将其与从scipy.stats.gamma获得的pdf进行比较来进行检查。