the line: m_test[abs(int(alpha))] is what i want to get.
but this is wrong。
X1 = np.random.randn(size)
basic_model = pm.Model()
m_test=[1,1,1,1,3,3,4,3]
with basic_model:
# Priors for unknown model parameters
alpha = pm.Uniform("alpha", lower=0, upper=5,transform=None)
# Expected value of outcome
mu = m_test[abs(int(alpha))]*X1
# Likelihood (sampling distribution) of observations
Y_obs = pm.Normal("Y_obs", mu=mu, observed=Y)
trace = pm.sample(1000)