Preliminary factor_gibbs model