create virtual observations
wi[:,i] = w + np.random.multivariate_normal(np.zeros(m), R)