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
- /
- modeling autoregressive covariance structure for u...

- 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

06-07-2016 01:41 PM

Hi,

I would like to analyse a longitudinal data set with relatively many missing values using proc mcmc. Before moving to the "missingness issue", however, I am not sure how to implement an autoregessive covariance structure in proc mcmc. I have seen a few examples on unbalanced longitudinal data with proc mcmc, but they do not model the covariance structure explicitly.

My data is in long format, of course, as the number of observations per subject varies. Now, to implement a specific covariance structure, I would need to define a multivariate distribution spanning multiple ROWS, would I not?

None of the examples I could find tackles this problem. Either the data is balanced (then it can be brought to wide format and it is clear how to specify the covariance structure) or the covariance structure was not modeled explicitly.

It would be great to get some hints as to whether this is possible with proc mcmc at all, and if so, which resources could help me understand how to implement this.

Best,

Kevin

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

06-07-2016 05:30 PM

This will take some programming. This article will help you get started.

http://support.sas.com/resources/papers/proceedings16/SAS5601-2016.pdf

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

06-08-2016 03:20 AM

Thanks for that excellent article! However, they again only describe the implementation of a variance structure for balanced data, but mine is unbalanced both in number as well as time difference between observations x)

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

06-08-2016 11:23 AM

can anybody comment on the use of the by statement in proc mcmc? It is supported, but I am unsure about how to utilize it for constructing the likelihood groupwise with sharerd e.g. variance parameters.

Havent seen it in any of the mcmc examples so far and the documentation is "sparse"...

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

06-08-2016 11:34 AM

The BY statement is the same for all SAS procedures. If there are k levels to the BY variables, then you get k independent analyses. It is equivalent to running the procedure k times using WHERE=(ByVar=k_th_level), except it is more efficient.

From your description of this problem, I don't think you want to use a BY variable.

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

06-08-2016 12:24 PM

Rick is right, the by statement is not what you want. As you know, MCMC is essentially a programming language for Bayesian analysis. The programs can be quite complex, and it is hard to give specific advice ahead of time. THe following article might help you with the "missing value" part of your problem.

http://support.sas.com/resources/papers/proceedings13/436-2013.pdf