Hi, I want to use PROC MCMC to run a random-effects linear negative binomial (NB1) model and estimate ~200,000 random intercepts (as an alternative to NLMIXED). Since the OUTPOST= dataset can only save up to 32,767 variables, what approaches can I take to be able to output the posterior samples for all of the ~200K random-effects parameters? Thanks, Andrea P.S. I have SAS v9.4 (TS1M4) and SAS/STAT 14.2.
... View more
