BookmarkSubscribeRSS Feed
Fluorite | Level 6 ZMX
Fluorite | Level 6

Hello there,


I set up a mix-logit model with random coefficients in PROC IML. Each random coefficient equation is a function of observed characteristics and unobserved heterogeneity. The unobserved components (error terms) are coming from a multivariate normal distribution with zero means and a variance-covariance matrix and they are added to the observed characteristics. I need to estimate the variance-covariance matrix parameters (variances and correlation coefficients) of the MVN along with the other parameters of the model using simulated maximum likelihood estimator.


I tried to estimate the parameters using Cholesky decomposition by simply adding standard normal errors that are generated outside the optimization procedure, to these coefficient equations. The coefficients of these error terms are then, the lower triangular matrix. However, I’m wondering if there is another way that does not require Cholesky decomposition of the error terms so that I can draw the errors from a MVN and estimate the error variances and correlation parameters directly (along with the other model parameters). I appreciate any suggestions.





The early bird rate has been extended! Register by March 18 for just $695 - $100 off the standard rate.


Check out the agenda and get ready for a jam-packed event featuring workshops, super demos, breakout sessions, roundtables, inspiring keynotes and incredible networking events. 


Register now!

Multiple Linear Regression in SAS

Learn how to run multiple linear regression models with and without interactions, presented by SAS user Alex Chaplin.

Find more tutorials on the SAS Users YouTube channel.

Discussion stats
  • 1 reply
  • 2 in conversation