Programming the statistical procedures from SAS

Proc MCMC: Why doesn't this print plots?

Accepted Solution Solved
Reply
New Contributor
Posts: 2
Accepted Solution

Proc MCMC: Why doesn't this print plots?

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;


Accepted Solutions
Solution
‎06-19-2012 03:41 PM
Super User
Posts: 10,871

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

How about a Run;  before the ODS GRAPHICS OFF;

View solution in original post


All Replies
Solution
‎06-19-2012 03:41 PM
Super User
Posts: 10,871

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

How about a Run;  before the ODS GRAPHICS OFF;

New Contributor
Posts: 2

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

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


🔒 This topic is solved and locked.

Need further help from the community? Please ask a new question.

Discussion stats
  • 2 replies
  • 180 views
  • 0 likes
  • 2 in conversation