看来您要将完整的精度矩阵传递给正态分布：
mu_u = pm.Normal(
'mu_u', mu=0, tau=beta_0 * lambda_u, shape=dim,
testval=np.random.randn(dim) * std)
我假设您只想传递对角线值：
mu_u = pm.Normal(
'mu_u', mu=0, tau=beta_0 * t.diag(lambda_u), shape=dim,
testval=np.random.randn(dim) * std)
这会更改mu_u并mu_v为您修复吗？
602392714
清零编程群