使用R语言mixcat包npmlt()函数时遇模型矩阵匹配错误求助

我想用R语言的mixcat包中npmlt()函数拟合混合效应多项逻辑回归模型，模型需要包含两个嵌套随机效应：参与者id，以及嵌套在id下的day（变量名为tag）。我几乎没有R语言使用经验，目前遇到了无法解决的错误。

尝试的代码 #

attach(final_merged_data_MLM)

model <- npmlt(formula = NSSI_risky_alone ~ lag_NA_within + centered_NA_global + 
                 lag_URGE_within + centered_URGE_global + geschlecht + stunde_cluster + 
                 bpd_anz + binary_tag,
               random = ~ 1 | id + tag)

收到的错误提示 #

Error in model.matrix.default(random, data = a) : 
  model frame and formula mismatch in model.matrix()

已尝试但无效的解决方法 #

• 参考经验将变量重命名为短名称，无效；
• 使用仅包含相关变量的数据框运行代码（原数据框含90个变量），无效；
• 仅保留单个预测变量运行代码，无效；
• 移除随机效应、更换因变量运行代码，无效；
• 检查变量格式（均正确指定为因子或数值型，因变量为因子，尝试移除缺失值后仍无效）；
• 尝试其他包（如mclogit）但出现其他错误，现专注于mixcat包。

恳请提供可行的解决方法。

内容的提问来源于stack exchange，提问作者Cristina Balemi