prs = np.tile(prs, (outcomes.shape[0],1,1,1)).transpose((1,0,2,3))