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;
How about a Run; before the ODS GRAPHICS OFF;
How about a Run; before the ODS GRAPHICS OFF;
Thank you. I guess I didn't think of that because it still produced so much output without it. :smileyblush:
Build your skills. Make connections. Enjoy creative freedom. Maybe change the world. Registration is now open through August 30th. Visit the SAS Hackathon homepage.
Register today!ANOVA, or Analysis Of Variance, is used to compare the averages or means of two or more populations to better understand how they differ. Watch this tutorial for more.
Find more tutorials on the SAS Users YouTube channel.