Nlmefit covb and PSI
5 次查看(过去 30 天)
显示 更早的评论
Hello
When using nlmefit or nlmefitsa for estimating a mixed effects model, the algorithm gives me 2 covariance matrices as result: PSI and stats.covb. Anyone can tell me what is the difference between these two? The help says this: PSI, an r-by-r estimated covariance matrix for the random effects. covb — The estimated covariance matrix of the parameter estimates
I know that if I want to simulate a population using nlmefitsa's results I should generate random samples using beta and PSI. But, how can I use covb?
Thanks if anyone can help me.
0 个评论
采纳的回答
Star Strider
2012-9-15
编辑:Star Strider
2012-9-15
I haven't done much with nmlefitsa, so I don't have much experience with it. As a general rule however, COVB (the covariance matrix of the parameter estimates), provides one rather important bit of information on the parameters themselves: it allows you to calculate the confidence intervals on the parameters.
With Beta a column vector of the parameter values, and CovB the covariance matrix, the 95% confidence intervals are given by:
CI95 = [Beta-1.96*sqrt(diag(CovB)) Beta+1.96*sqrt(diag(CovB))];
The 1.96 value is the ‘critical value’ (the z-statistic or z-score) corresponding to the 95% confidence interval:
cv = norminv(0.975, 0, 1);
The usual interpretation of the confidence intervals on the parameters is that if the confidence interval for a particular parameter includes zero, that parameter is not needed in the model.
2 个评论
Star Strider
2012-9-16
My pleasure!
Your question stimulated my curiosity as well to see if there was other information in the parameter covariance matrix and other information about the estimated parameters. I discovered a good discussion in ‘A Sensitivity Matrix Based Methodology for Inverse Problem Formulation’. It did not change what I stated in my answer, but it is the most thorough discussion of parameter characteristics and their interpretation that I've seen recently.
I am not certain what you mean by ‘intrinsic variability’, but expressions of the effect of the parameters on the model are variously the Jacobian and the sensitivity matrix, and there may be others. The Jacobian is the partial derivatives of the function with respect to each parameter. This becomes a problem in models that are ‘nonlinear in the parameters’, in which some or all of the partial derivatives of the function with respect to each parameter are functions of at least one other parameter. An interesting discussion of that is ‘Traditional and Generalized Sensitivity Functions and Their Uses in Biomathematical Models’.
更多回答(0 个)
另请参阅
类别
在 Help Center 和 File Exchange 中查找有关 Linear Model Identification 的更多信息
Community Treasure Hunt
Find the treasures in MATLAB Central and discover how the community can help you!
Start Hunting!