RGLab / MAST

Tools and methods for analysis of single cell assay data in R
224 stars 57 forks source link

Deviance #155

Open sofiapuvogelvittini opened 3 years ago

sofiapuvogelvittini commented 3 years ago

Hello MAST team, I am using MAST for my single cell data. I would like to make models comparisons based on daviance or AIC measurements, because I have different patients ralated variables that I should test if they should be included as fixed effect in the model.

My formula is the following

zlmCond<-MAST::zlm(~condition + cenegeneson + test_var1 + test_var2 + test_varx + (1|donor_id), data, method=´glmer´,ebayes=FALSE)

How could I obtain deviance or AIC measurements? I try to use deviance_residuals_hook(zlmCond) but it is deprecated in R 4.0

Could you help me please? Kind regards, Sofía

amcdavid commented 3 years ago

There's not a method for it at the moment, but @.` and @.Residual` has what you need. Not the interpretation is a bit more complicated than the classical case because of the penalization but I expect they would still serve to compare between zlm models.

On Thu, Apr 15, 2021, 11:44 AM sofiapuvogelvittini @.***> wrote:

Hello MAST team, I am using MAST for my single cell data. I would like to make models comparisons based on daviance or AIC measurements, because I have different patients ralated variables that I should test if they should be included as fixed effect in the model.

My formula is the following

zlmCond<-MAST::zlm(~condition + cenegeneson + test_var1 + test_var2 + test_varx + (1|donor_id), data, method=´glmer´,ebayes=FALSE)

How could I obtain deviance or AIC measurements? I try to use deviance_residuals_hook(zlmCond) but it is deprecated in R 4.0

Could you help me please? Kind regards, Sofía

— You are receiving this because you are subscribed to this thread. Reply to this email directly, view it on GitHub https://github.com/RGLab/MAST/issues/155, or unsubscribe https://github.com/notifications/unsubscribe-auth/AALLAHWAFCOOUL4OYXMWLR3TI4CVVANCNFSM427YJKKA .

sofiapuvogelvittini commented 3 years ago

Thanks!