R语言mlogit包系数方差-协方差矩阵的获取与计算方法

mlogit包的mlogit()拟合结果本身确实不会直接存储方差-协方差矩阵，但你可以通过以下两种方式获取（本质都是基于Fisher信息矩阵逆矩阵计算，和summary()输出标准误的逻辑完全一致）：

方法1：直接使用vcov()函数 #

vcov()函数可直接作用于mlogit模型对象，返回系数的方差-协方差矩阵：

# 假设你的嵌套Logit模型对象名为nested_model
vcov_matrix <- vcov(nested_model)

方法2：从summary对象中提取 #

如果已经生成了模型的summary对象，可直接从其中提取该矩阵：

sum_model <- summary(nested_model)
vcov_matrix <- sum_model$vcov

验证逻辑 #

summary()输出的标准误，就是该方差-协方差矩阵对角线元素的平方根，你可以通过以下代码验证一致性：

# 提取summary中的标准误
summary_se <- sum_model$coefficients[, "Std. Error"]
# 从方差矩阵手动计算标准误
manual_se <- sqrt(diag(vcov_matrix))
# 对比两者是否一致
all.equal(summary_se, manual_se)

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