How to model a multivariate multilevel linear model

Questions
6

For the trace it’s just pm.sample().

Also you can clan up the code a little bit if you have something like

sitenames = ['辐照度', '风速', ...]
beta_eff = 0
for beta_idx, site_name in enumerate(sitenames):
  beta_eff += beta_mat[beta_idx][site_idx]*all_site_train[site_name]
mu = pm.Deterministic('mu', alpha[site_idx] + beta_eff)