turn on suggestions

Auto-suggest helps you quickly narrow down your search results by suggesting possible matches as you type.

Showing results for

Find a Community

- Home
- /
- Analytics
- /
- Stat Procs
- /
- PRIOR statement in mixed procedure

Topic Options

- Subscribe to RSS Feed
- Mark Topic as New
- Mark Topic as Read
- Float this Topic for Current User
- Bookmark
- Subscribe
- Printer Friendly Page

- Mark as New
- Bookmark
- Subscribe
- Subscribe to RSS Feed
- Highlight
- Email to a Friend
- Report Inappropriate Content

01-25-2016 05:01 PM

Dear Everybody,

I need some help about the PRIOR statement in the mixed procedure.

Is thsi option allow bayesian estimation?

The different estimates are exactly the same with and without the option and the posterior distribution mean and 5% and 95% percentiles are very close to the frequentist mean and the 90% confidence interval.

I use the non informative JEFFREYS prior.

SAS doesn't mention this procedure as a procedure able to make bayesian estimation.

Thank in advance for your feedback.

Samir

Accepted Solutions

Solution

01-26-2016
04:56 AM

- Mark as New
- Bookmark
- Subscribe
- Subscribe to RSS Feed
- Highlight
- Email to a Friend
- Report Inappropriate Content

01-25-2016 08:31 PM

What PROC MIXED provide is falling under what is called "Empirical Bayes". Because the parameters or hyperparameters of the prior distributions of the random effects are actually estimated from the data via the likelihood functions.

Not only PROC MIXED offeres that but other procedures also see: Usage Note 23407: Available Bayesian methods

But if we talk about fully Bayesian treatment then you should consider using one of the procedures mentioned as in here: Bayesian Analysis Using SAS/STAT Software

All Replies

Solution

01-26-2016
04:56 AM

- Mark as New
- Bookmark
- Subscribe
- Subscribe to RSS Feed
- Highlight
- Email to a Friend
- Report Inappropriate Content

01-25-2016 08:31 PM

What PROC MIXED provide is falling under what is called "Empirical Bayes". Because the parameters or hyperparameters of the prior distributions of the random effects are actually estimated from the data via the likelihood functions.

Not only PROC MIXED offeres that but other procedures also see: Usage Note 23407: Available Bayesian methods

But if we talk about fully Bayesian treatment then you should consider using one of the procedures mentioned as in here: Bayesian Analysis Using SAS/STAT Software

- Mark as New
- Bookmark
- Subscribe
- Subscribe to RSS Feed
- Highlight
- Email to a Friend
- Report Inappropriate Content

01-26-2016 04:55 AM

Thanks a lot mohamed-zaki for your answer.

It confirms what I suspect.

Have a nice day,

Samir

- Mark as New
- Bookmark
- Subscribe
- Subscribe to RSS Feed
- Highlight
- Email to a Friend
- Report Inappropriate Content

01-26-2016 10:37 AM

I somewhat disagree with the 'solution' post. With he PRIOR statement, one is getting true Bayesian analysis with PROC MIXED. See the very nice presentation in chapter 13 of the book SAS for Mixed Models, 2nd edition (2006). One just does not have much control over the priors. A Jeffrey's prior is used for the variance parameters (these are based on the data, but are considered to be non-informative). A flat (constant) improper prior is used for the fixed effects. There is direct sampling from the posteriors (using a independence chain). With these restrictions, the posterior is conditionally multivariate normal. The statistics paper describing the approach in MIXED is: Wolfinger and Kass (2000). Non-conjugate Bayesian Analysis of Variance Component Models. BIOMETRICS 56:768-774.

For more flexible approaches to Bayesian anlaysis, one has to use the other listed procedures.

- Mark as New
- Bookmark
- Subscribe
- Subscribe to RSS Feed
- Highlight
- Email to a Friend
- Report Inappropriate Content

02-01-2016 12:35 PM - edited 02-01-2016 12:36 PM

I would like to add to Larry's comment. C.R. Henderson showed in the 80s, based on a report from Laird and Ware, that using maximum likelihood estimates of the G and R matrices yields an empirical Bayes estimate of the random effects. So even without sampling, mixed model estimates of the variance components are Bayesian, given the data. This traces back to Harville's work that REML is related to Bayesian marginal inference.{Harville, D. A. (1974). Bayesian inference for variance components using only error contrasts. Biometrika 61, 383–385.}

I think I just proved I'm old...

Steve Denham