关于如何使用PyMC将两个Normal分布拟合为数据,有a question on CrossValidated。的Cam.Davidson.Pilon答案是使用伯努利分布数据分配到两条法线之一: size = 10
p = Uniform("p", 0 , 1) #this is the fraction that come from mean1 vs mean2
ber = Bernoul
我写了一个PyMC模型拟合3个法线使用数据(类似于一个在this question)。 import numpy as np
import pymc as mc
import matplotlib.pyplot as plt
n = 3
ndata = 500
# simulated data
v = np.random.randint(0, n, ndata)
data = (