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
- /
- How to estimate repeated semicontinuous data as in...

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

10-06-2015 11:16 AM

I would like to implement the two-part model of Olsen and Schafer (2001), who estimate coefficients in a model where the dependent variable is semicontinuous (lots of zeroes) and individuals are observed repeatedly over time. I guess my question is, is there a snatdard or easy way to do this, since the random intercepts estimated in the first part affect the covariance matrix of the effects in the second part (at least thats how I read it).

The first part is a logistic model with random effects. The second part numerically approximates a marginal likelihood w/ Gauss-Hermite quadrature, then maximizes it with Newton-Raphson.

Accepted Solutions

Solution

10-12-2015
01:06 PM

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

10-09-2015 10:15 AM

I ended up going with PROC NLMIXED. A very helpful resource was the white paper by McTernan and Blozis, "User-specified likelihood expressions using NLMIXED and the GENERAL statement." Their first example was very close to what i was looking for. For future reference here is a link:

http://wuss.org/Proceedings13/67_Paper.pdf

All Replies

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

10-07-2015 02:34 PM

This sounds like either fitting zero inflated models using PROC GENMOD or a mixture model using PROC FMM. Have you looked at the documentation for these procedures?

Steve Denham

Solution

10-12-2015
01:06 PM

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

10-09-2015 10:15 AM

I ended up going with PROC NLMIXED. A very helpful resource was the white paper by McTernan and Blozis, "User-specified likelihood expressions using NLMIXED and the GENERAL statement." Their first example was very close to what i was looking for. For future reference here is a link:

http://wuss.org/Proceedings13/67_Paper.pdf

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

10-12-2015 08:17 AM

I don't know if you can mark your own response as "Correct" in the new forums, but you certainly should, if able. Good finding.

Steve Denham