I'm trying to perform a DEG analysis of two cell types extracted from 8 different patients (data comes from scRNA).
I'm trying to compare the two cells while adjusting for the patients, there are two possible models:
a) ~0+cell.type+patients b) ~cell.type.patients lrt <- glmWeightedF(fit, coef = 2)
If I'm not wrong:
Model a, each coefficient represents the average of the samples for each coefficient.
Model b, the intercept represents the average of coefficient A and the 2º coefficient represents the increase in average of coefficient B over A.
Which model should I use if I want to test for DEG between those two cell types?