PROC MCMC for surveylogistic to resample error and get full CI

saslove — Wed, 06 Mar 2013 16:58:12 GMT

I came across this question at work. I would like to have your input. I have never done monte carlo and i am a newbie. So, please excuse my dumbness!

I want the standard error from model , error estimate and confidence interval on predicted number of response variable. I want to do simulation (i was told multiply predicted probability by case weight and sum it all up and repeat 1000 times). I need 2.5th and 97.5th % confidence interval and errors on percent predictions. PROC SURVEYLOGISTIC was used for analysis. I then used xbeta and stdxbeta to get the linear predictors and standard error of linear predictor. I back transformed the values. I'm trying to use PROC MCMC here. So, I assume standard error is normally distributed and resample error 1000 times. Also, I need a confidence interval on full prediction and sample error on output.

I would like some help here as to use PROC MCMC for getting the desired output. These are just the notes I took from our meeting.

This is the surveylogistic code.

proc surveylogistic data=outboot;

where var1="abc" and var2 in ("abc") and var3 in (11, 13) and var4>15

and var6>1999 and 85>var7>15;

format var varfmt..;

stratum strat;

cluster psu;

weight weight;

class var var2 /param=ref;

model resp (desc)=vara varb varc vard

/ Rsquare ;

output out=outall predprobs=individual xbeta=predval stdxbeta=sepred;

run;

How do I use PROC MCMC to get confidence intervals and standard errors. I need CI on full prediction and need to sample error on output. Please help.

Re: PROC MCMC for surveylogistic to resample error and get full CI

saslove — Wed, 06 Mar 2013 20:47:10 GMT

Anyone done anything similar????

topic Re: PROC MCMC for surveylogistic to resample error and get full CI in Statistical Procedures

PROC MCMC for surveylogistic to resample error and get full CI

Re: PROC MCMC for surveylogistic to resample error and get full CI