r - 分位点回帰モデルの信頼区間で 10 番目、50 番目、および 90 番目の分位点をプロットする方法

Question

分位点回帰モデルをバイナリ結果応答のデータに適合させることに興味があるため、GJRMR のパッケージを使用します。モデルを適合させた後、予測変数ごとにいくつかの分位点を使用して信頼区間をプロットしようとしましたが、できませんでした。 t。ここでは、パッケージ自体の例を含めます。また、各適合モデルの AIC をどのように見つけることができますか? 与えられた例の AIC は、両方の適合方程式に対するものです。

set.seed(0) 
n <- 400; Sigma <- matrix(0.5, 2, 2); diag(Sigma) <- 1 
u <- GJRM::rMVN(n, rep(0,2), Sigma) 
x1 <- round(runif(n)); x2 <- runif(n); x3 <- runif(n) 
f1 <- function(x) cos(pi*2*x) + sin(pi*x) 
f2 <- function(x) x + exp(-30*(x - 0.5)^2) 
y1 <- ifelse(-1.55 + 2*x1 + f1(x2) + u[, 1] > 0, 1, 0) 
y2 <- ifelse(-0.25 - 1.25*x1 + f2(x2) + u[, 2] > 0, 1, 0) 
dataSim <- data.frame(y1, y2, x1, x2, x3) ## CLASSIC BIVARIATE PROBIT 
out <- GJRM::gjrm(list(y1 ~ x1 + x2 + x3, y2 ~ x1 + x2 + x3),
                  data=dataSim, margins=c("probit", "probit"), 
                  Model="B")
GJRM::conv.check(out) 
summary(out) 
AIC(out) 
BIC(out)