Create gibbs.R