BookmarkSubscribeRSS Feed
ZMX
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.

1 REPLY 1

SAS Innovate 2025: Register Now

Registration is now open for SAS Innovate 2025 , our biggest and most exciting global event of the year! Join us in Orlando, FL, May 6-9.
Sign up by Dec. 31 to get the 2024 rate of just $495.
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
  • 789 views
  • 0 likes
  • 2 in conversation