Confidence levels for glmfit - glmval or fitglm using the Likelihood Ratio method

5 次查看（过去 30 天）
Craig
Craig 2024-4-27
评论： Craig 2024-5-6
After several days and many hours of researching and googling, I am unable to determine how I could find confidence levels using the Likelihood Ratio for either glmfit - glmval or fitglm. I know it can be done in R, but I don’t see that capability in Matlab. Does it exist? Neither glmval nor coefCI appear to use or have options for using the Likelihood Ratio method. And paramci does not appear to work with a GeneralizedLinearModel object.
Thanks for any help you can provide.

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

回答（1 个）

Akshat
Akshat 2024-5-3
Hi Craig,
Upon looking up the documentation for the relevant functions which can be used to serve your functionality, I unfortunately couldn't find anything that finds the confidence levels using the Likelihood ratio.
On the other hand, I did see that "coefCI" uses the Wald Method to determine the confidence intervals, which as far as I could find, is an approximation of the Likelihood Ratio test.
My reference was : https://stats.oarc.ucla.edu/other/mult-pkg/faq/general/faqhow-are-the-likelihood-ratio-wald-and-lagrange-multiplier-score-tests-different-andor-similar/#:~:text=The%20Wald%20test%20approximates%20the,is%20equal%20to%20some%20value.
Now as a workaround as you said, there is a functionality in R for the same, I can suggest you the following:
https://www.mathworks.com/matlabcentral/fileexchange/5051-matlab-r-link
This is a filexchange developed to create a link between MATLAB and R, using which you can get the variables from R, and also run scripts there. This can help you keep the entire pipeline in MATLAB, with just the confidence intervals being calculated in R.
Hope this helps!
  1 个评论
Craig
Craig 2024-5-6
Hello Akshat, Thank you for looking at this. I had found your reference, and it was helpful, but as you say, the Wald Method is just an approximation to the Likelihood Ratio method.
I was not aware of the MATLAB R-link file exchange package. That might be a work around, it would require me to reverse engineer the R-code I have, but it is a possibility.
I'm going to leave this open for now to see if anyone can provide additional assistance.
Again, thanks for your helpful reply.

请先登录，再进行评论。

类别

在 Help Center 和 File Exchange 中查找有关 Analysis of Variance and Covariance 的更多信息

标签

产品


版本

R2023a

Community Treasure Hunt

Find the treasures in MATLAB Central and discover how the community can help you!

Start Hunting!

Translated by