pr0 = np.ones((modelparams.shape[0], expparams.shape[0])) / 2