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

Question

Craig 2024-4-27

0
链接

此问题的直接链接

https://ww2.mathworks.cn/matlabcentral/answers/2112266-confidence-levels-for-glmfit-glmval-or-fitglm-using-the-likelihood-ratio-method

评论： 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.

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

请先登录，再进行评论。

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

Answer 1

Akshat 2024-5-3

0
链接

此回答的直接链接

https://ww2.mathworks.cn/matlabcentral/answers/2112266-confidence-levels-for-glmfit-glmval-or-fitglm-using-the-likelihood-ratio-method#answer_1451901

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 个评论
显示 -1更早的评论隐藏 -1更早的评论

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.

请先登录，再进行评论。