多項ロジスティック回帰の予測値 (信頼区間付き) を取得したいと考えています。これは予測で実行できることはわかっていますが、私の場合、次の方法で標準エラーをクラスター化しました。
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 で行う簡単な方法があるかどうかはわかりません。)