Have I made a mistake in this code? I don't understand why this isn't printing the plots normally associated with Proc MCMC. (Initially, I had more data, but shortened it for this post. Sorry if it's still too long.) Any suggestions would be appreciated. Thanks. title 'Why Doesnt This Print Plots?'; options validvarname=any; data coded; input X1 X2 X3 @@; datalines; 49.19615 168.9541 339.0285 59.79256 52.22341 239.234 69.32701 428.0111 507.1756 21.82025 26.6339 232.2198 80.69061 130.7855 391.3814 85.07306 334.1087 532.6584 262.8237 78.27474 484.2632 12.34403 26.53437 98.95171 ; ods graphics on; proc mcmc data=coded outpost=bikesout ntu=30000 nmc=500000 thin=100 seed=448 monitor=(_parms_) propcov=nmsimp(itprint) init=pinit; parms Beta1 0 Beta2 0; parms sigma 1; prior Beta1 Beta2 ~ normal(mean = 0, var = 1e6); prior sigma ~ igamma(shape = 3/10, scale = 10/3); mu = Beta1*X1 + Beta2*X2; model X3 ~ n(mu, var = sigma); ods graphics off;
... View more
