Proc MCMC: Why doesn't this print plots?

DavidClark — Tue, 19 Jun 2012 18:50:56 GMT

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;

Re: Proc MCMC: Why doesn't this print plots?

ballardw — Tue, 19 Jun 2012 19:41:20 GMT

How about a Run; before the ODS GRAPHICS OFF;

Re: Proc MCMC: Why doesn't this print plots?

DavidClark — Tue, 19 Jun 2012 19:48:15 GMT

Thank you. I guess I didn't think of that because it still produced so much output without it. :smileyblush:

topic Re: Proc MCMC: Why doesn't this print plots? in Statistical Procedures

Proc MCMC: Why doesn't this print plots?

Re: Proc MCMC: Why doesn't this print plots?

Re: Proc MCMC: Why doesn't this print plots?