Turn on suggestions

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

Showing results for

- Home
- /
- Analytics
- /
- Stat Procs
- /
- Re: PROC MCMC Multivariate t Prior in Mixed Model

Options

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

- Mark as New
- Bookmark
- Subscribe
- Mute
- RSS Feed
- Permalink
- Report Inappropriate Content

Posted 03-07-2016 02:06 PM
(1588 views)

Hi all,

I was wondering if anyone has tried to specify a multivariate t distribution for the fixed effects in a mixed model. I have some code below that I used to fit a model to some simulated data. In this code, I assumed a multivariate normal prior on the beta's. Does anyone know whether or not I can specify a multivariate t? I looked into the GENERAL statement in PROC MCMC where you can specify custom distributions. Does it look like I'll have to use GENERAL?

```
proc mcmc data = SimulatedData;
array theta[2] b0 b1;
array theta_c[2];
array mu0[2] (0,0);
array sigma_c[2,2];
array sigma_0[2,2] (10000 0 0 10000);
array s[2,2] (0.02 0 0 20);
parms theta_c sigma_c {.1 0 0 .1} var_y;
prior theta_c ~ mvn(mu0, sigma_0);
prior sigma_c ~ iwish(2, s);
prior var_y ~ igamma(0.01, scale = 0.01);
random theta ~ mvn(theta_c, sigma_c) subject = Subject;
mu = b0 + b1*Time;
model Response ~ Normal(mu, var = var_y);
run;
```

Thanks.

3 REPLIES 3

- Mark as New
- Bookmark
- Subscribe
- Mute
- RSS Feed
- Permalink
- Report Inappropriate Content

- Mark as New
- Bookmark
- Subscribe
- Mute
- RSS Feed
- Permalink
- Report Inappropriate Content

Yes, multivariate t is not one of the recognized distributions. I tried specifying a t-distribution on the PARMS statement but that did not lead to the goal of specifying a multivariate t on the fixed effects. I am currently trying to use the GENERAL statement to specify a custom multivariate t. Have you had any experience with this?

- Mark as New
- Bookmark
- Subscribe
- Mute
- RSS Feed
- Permalink
- Report Inappropriate Content

Are you ready for the spotlight? We're accepting content ideas for **SAS Innovate 2025** to be held May 6-9 in Orlando, FL. The call is **open **until September 16. Read more here about **why** you should contribute and **what is in it** for you!

What is ANOVA?

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.