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