ERGM新手求助：ergm.AME()报non-conformable arguments及结构零/一警告

1. 矩阵不兼容报错（non-conformable arguments） #

这个错误的核心是协变量矩阵和系数向量维度不匹配，结合你使用的复杂模型，可能的原因包括：

• edgecov网络与主网络不匹配：edgecov()要求传入的辅助网络（network1-network4）和主模型的network节点完全一致——节点数量、ID、顺序都不能错。如果辅助网络的节点结构和主网络有差异，会导致计算边际效应时生成的协变量矩阵维度异常，触发矩阵乘法错误。
• 模型项冗余或共线性：你加入了大量节点属性相关项（如nodeifactor/nodeofactor/nodematch的组合，以及nodefactor+nodematch），这类组合可能导致多重共线性，或者在ERGM内部计算AME时，变量矩阵的维度出现冲突。
• MPLE估计的局限性：MPLE（最大伪似然估计）适合简单模型的快速估计，但复杂模型下系数估计可能不稳定，进而导致后续AME计算的系数向量与协变量矩阵维度不匹配。

2. 结构零/一警告 #

这个警告说明网络中存在结构上必然存在或不存在边的二元组，复杂模型加入更多节点属性后更容易触发：

• 主网络本身的结构约束：比如创建主网络时设置了二分图（bipartite）、禁止自环（默认loops=FALSE），或者手动指定了禁止/强制的二元组，这些都会被识别为结构零/一。
• 辅助edgecov网络的结构影响：如果network1-network4本身存在结构零/一，会被主模型的ERGM过程识别为约束条件，进而在计算预测概率时触发警告。

排查与解决建议 #

• 检查edgecov网络一致性：用以下代码对比主网络和辅助网络的节点信息：

# 对比节点数量
network.size(network)
network.size(network1)

# 对比节点名称/ID
all(network.vertex.names(network) == network.vertex.names(network1))

如果不一致，重新对齐辅助网络的节点后再传入edgecov()。

• 逐步简化模型排查问题：先去掉所有edgecov项，运行ergm.AME看是否正常；再逐个加入edgecov，定位是哪个辅助网络导致的问题。同理，逐步去掉部分节点属性项，找到触发矩阵错误的冗余项。
• 改用MLE估计：将模型的estimate = "MPLE"替换为estimate = "MLE"（注意MLE计算时间更长，需要确保模型收敛），复杂模型下MLE的系数估计更稳定，可能解决矩阵维度问题。
• 处理结构零/一：用ergm.getmodel(model5)$constraints查看模型的约束项，或者用summary(network)查看网络的结构属性。如果是合法的结构约束（比如自环），可以忽略警告；如果是异常约束，需要在创建网络或模型时调整约束条件，或者手动标记这些二元组的概率为0/1后再计算AME。

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