我试图使用numpy.random.multivariate_normal生成多个样本,其中每个样本都是从具有不同mean和cov的多元正态分布中提取的。例如,如果我想画两个样本,我试过covs = np.array([np.identity(2) for k in xrange(2)])
rand.m
我试图找到50个样本均值的95%可信区间。样本大小范围从2到600,每个样本中的值界于1到5之间。大小为10或更大的样本具有对数正态分布,其中我使用JAGS对来自John K.Kruschke的对数正态参数进行贝叶斯估计,模型规范如下: model {但是,当3个<=样本< 10时,我得到了上限的极值(例