Standar Errors for variance components of non-linear mixed effect models

Question

Alexis Villacis 2018-12-15

0
链接

此问题的直接链接

https://ww2.mathworks.cn/matlabcentral/answers/435944-standar-errors-for-variance-components-of-non-linear-mixed-effect-models

评论： Jeff Miller 2018-12-16

Hello!

I am estimating a non-linear mixed effect model using nlmefit, and although Matlab produces standard errors for the estimated betas, it does not produce standard errors for the variance components (the PSI matrix). How can I obtain/estimate the standard errors for the variance components? Would it also be possible to estimate the s.e. for the covariance components of the PSI matrix?

Any help will be really appreciated from the bottom of my heart.

Thank you in advance,

0 个评论
显示 -2更早的评论隐藏 -2更早的评论

请先登录，再进行评论。

请先登录，再回答此问题。

Answer 1

Jeff Miller 2018-12-16

0
链接

此回答的直接链接

https://ww2.mathworks.cn/matlabcentral/answers/435944-standar-errors-for-variance-components-of-non-linear-mixed-effect-models#answer_352592

One way is to use the bootstrp function. For each bootstrap sample, you would write your own function to call nlmefit and save the output PSI matrix values. After you've done that for many bootstrap samples, you can get s.e.'s using the stored PSI matrix values from the individual bootstrap results (e.g., Wikipedia ).

2 个评论
显示无隐藏无

Alexis Villacis 2018-12-16

Thank you for the tip! Will definitely take it into account.

In R, you can do it by using the confidence interval of these estimates and I was just wondering if there was something similar in Matlab, i.e. Is it possible in Matlab to generate the confidence intervals for the variance components (the PSI matrix)?

I am running 6 different models using LaPlacian approximation and each model takes around 15 minutes in my computer (macbook pro 2.9 GHz intel core i5 8Gb Ram). Generating enough draws using bootstrapping will take a long time.

Any thoughts?

Jeff Miller 2018-12-16

Sorry, I have no helpful thoughts. I don't know enough about these models to know how the R routines get confidence intervals. Some procedure based on Fisher information maybe?

Did you look at nlmefitsa? It might use the same algorithm as the R package.

请先登录，再进行评论。