使用R计算多项预测周围误差的方法

Question

首先，抱歉重新发布，但我添加了一个示例代码来解释我的请求。 希望这更清楚。

使用＆＃34; multinom＆＃34;将多项模型拟合到我的数据之后function（package nnet），我想显示控制其他变量值的选定变量的效果。我知道＆＃34;效果＆＃34;包主要做我想要的，但我希望能够自己计算预测误差（置信区间）。 有人可以告诉我方法和R代码吗？ 我认为我们应该使用delta方法，但我不确定如何在这种情况下应用它。

这是一个小示例代码（基于效果包中可用的数据）

library(nnet)
library(effects)
mod <- multinom(vote ~ age + gender, data=BEPS)
summary(mod)

# Call:
# multinom(formula = vote ~ age + gender, data = BEPS)

# Coefficients:
#                  (Intercept)         age gendermale
# Labour             1.2241862 -0.01562320  0.1682676
# Liberal Democrat   0.4979706 -0.01551381  0.1240998

# Std. Errors:
#                  (Intercept)         age gendermale
# Labour             0.2277826 0.003830006  0.1204621
# Liberal Democrat   0.2694373 0.004578836  0.1436882

# Residual Deviance: 3186.266 
# AIC: 3198.266 

plot(allEffects(mod))

我唯一需要的是能够计算此图中显示的错误值！

提前谢谢你，