1

我有兴趣将分位数回归模型拟合到具有二元结果响应的数据，因此，我使用GJRMR 中的包。在拟合模型后，我尝试为每个预测变量绘制带有多个分位数的置信区间，但是，我不能吨。在这里，我包含了一个来自包本身的示例。另外，我怎么能找到每个拟合模型的 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)
4

0 回答 0