Deviance

[sofiapuvogelvittini]

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]

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 <#155>, or unsubscribe <https://github.com/notifications/unsubscribe-auth/AALLAHWAFCOOUL4OYXMWLR3TI4CVVANCNFSM427YJKKA> .

[sofiapuvogelvittini]

Thanks!