具有行和列预测器的R-线性模型

Question

对创建模型不熟悉。我有一个矩阵（bpt2t），基因作为行（〜2000），样本（〜400）作为列。矩阵可以是0（表示完整的基因）或1（表示缺失的基因）。

我尝试使用glm创建一个nullmodel。假设较长的基因比较小的基因具有更多的断裂，并且某些样本具有比其他基因更多的断裂。我遍历此矩阵中的每个项目并将其命名为“ break”。我还有另一个变量“ genelength”，即每个基因的长度。

因此它使用基因长度来预测断裂：

nullmodel <- glm( break ~ genelength, family = "binomial" )

是否有一种方法还可以包括每个样本的断裂数？我有一个变量，它是每个样本的断裂数，但是它的长度不同于基因长度。