我正在使用gnmR 中的函数。我试图拟合一个具有积分部分的模型，我通过使用六个矩阵对其进行离散化，这些矩阵的列是每个时间间隔的协变量的二进制指标。对于矩阵 1，这些列中的每一列都标有a1, a2等，b1, b2对于矩阵 2，等等。对于这些矩阵中的每一个，我都有大约 250 列。
我将这些矩阵作为字符串传递给 gnm 命令：
mod<-gnm(as.formula(expr),family=binomial)
其中expr是一个字符串，其中包含所有矩阵列变量名称的列表（即a1, a2, ..., b1, b2, ....）。
当我使用 5 个矩阵时，我的代码会产生合理的结果，但是当我尝试使用 6 个矩阵时，当我尝试运行gnm命令时会出现以下错误。
variable names are limited to 10000 bytes
有没有办法绕过这个？