4

我想获得多项逻辑回归的预测值（带有置信区间）。我知道这可以通过 predict 来完成，但在我的情况下，我通过以下方式聚集了标准错误：

multinom <- mlogit(Y ~0| X1+ X2 , data)
cl.mlogit   <- function(fm, cluster){
  M <- length(unique(cluster))
  N <- length(cluster)
  K <- length(coefficients(fm))
  dfc <- (M/(M-1))
  uj  <- apply(estfun(fm),2, function(x) tapply(x, cluster, sum));
  vcovCL <- dfc*sandwich(fm, meat.=crossprod(uj)/N)
 coeftest(fm, vcovCL) 
}
cl.mlogit(multinom, data$group)

例如，我如何使用这些结果来获得 X1=1 和 X2=0 的预测概率（带有置信区间），并将其与 X1=2 和 X2=0 的预测概率进行比较。另外，我怎样才能得到这种差异的置信区间？在 stata prvalue 可以做到这一点，但不知道在 R 中是否有一种简单的方法可以做到这一点。）

4

0 回答 0