LISTSERV at the University of Georgia
Menubar Imagemap
Home Browse Manage Request Manuals Register
Previous messageNext messagePrevious in topicNext in topicPrevious by same authorNext by same authorPrevious page (April 2008, week 3)Back to main SAS-L pageJoin or leave SAS-L (or change settings)ReplyPost a new messageSearchProportional fontNon-proportional font
Date:   Wed, 16 Apr 2008 12:53:49 -0400
Reply-To:   Sigurd Hermansen <HERMANS1@WESTAT.COM>
Sender:   "SAS(r) Discussion" <SAS-L@LISTSERV.UGA.EDU>
From:   Sigurd Hermansen <HERMANS1@WESTAT.COM>
Subject:   Re: Help Needed in simualting the MCMC
Comments:   To: meesum kazmi <meesumkazmi@gmail.com>
In-Reply-To:   <fa507f530804150244o4a96c306pd134ea560c3255d7@mail.gmail.com>
Content-Type:   text/plain; charset="us-ascii"

Meesum: I wouldn't call what I can tell you an "expert opinion". You'll have to take it for what it's worth.

In the example below I specified in PROC MI the number of MCMC imputations (replicates) with the Nimpute setting. The source dataset, ExperimentU, included some rows with a missing value for the DV. PROC MI imputes missing DV values as well as missing IV values in the five (in this case) MCMC replicates. I then used the imputed dataset, ExperimentUI, and PROC REG to estimate model parameters, and PROC MIANALYZE to adjust estimated variances.

proc mi data=ExperimentU seed=337131 Nimpute=5 round=.01 minmaxiter=20 out=ExperimentUI ; MCMC outest=MCMC /*ACFPlot=WLF*/ ; EM maxiter= 100 /* ; var predM x1 x2 x3; run; quit;

ods trace on; proc reg data=ExperimentUI outest=parmests covout noprint; model predM = x1 x2 x3/ covb ; by _imputation_; output p=ests; ods output covb=covbests; run; quit;

proc mianalyze data=parmests; modeleffects Intercept x1 x2 x3; ; run; ods trace off;

S

-----Original Message----- From: owner-sas-l@listserv.uga.edu [mailto:owner-sas-l@listserv.uga.edu] On Behalf Of meesum kazmi Sent: Tuesday, April 15, 2008 5:45 AM To: SAS-L@listserv.uga.edu Subject: Help Needed in simualting the MCMC

Hi All, Need your help in Simulating the MCMC model. Just to give the background i tried runnning multiple linear regression on different time frame data, as a result of this excercise i got coeffecients for the independent variables. now to get a robust estimates of these variables i want to simulate the MCMC on the these coeffecients. Need your help in validating the methodology, and wants expert opinion on how can we simulate this in SAS. There is experimental proc in SAS 9.2 which simulates MCMC and we can also use proc mi but currently i am not sure on how to use these procedures. Any help would be highly appreciated

Best, Meesum


Back to: Top of message | Previous page | Main SAS-L page