# NOT RUN {
# example continued from DiSSMod
set.seed(45)
data(DoctorRWM, package = "DiSSMod")
n0 <- 600
set.n0 <- sample(1:nrow(DoctorRWM), n0)
reduce_DoctorRWM <- DoctorRWM[set.n0,]
result0 <- DiSSMod(response = as.numeric(DOCVIS > 0) ~ AGE + INCOME_SCALE + HHKIDS + EDUC + MARRIED,
selection = PUBLIC ~ AGE + EDUC + FEMALE,
data = reduce_DoctorRWM, resp.dist="bernoulli", select.dist = "normal",
alpha = seq(-5.5, -0.5, length.out = 21), standard = TRUE)
plot(result0, level = 0.90)
data(CreditMDR, package = "DiSSMod")
n1 <- 600
set.n1 <- sample(1:nrow(CreditMDR), n1)
reduce_CreditMDR <- CreditMDR[set.n1,]
result1 <- DiSSMod(response = MAJORDRG ~ AGE + INCOME + EXP_INC,
selection = CARDHLDR ~ AGE + INCOME + OWNRENT + ADEPCNT + SELFEMPL,
data = reduce_CreditMDR, resp.dist="poi", select.dist = "logis",
alpha = seq(-0.3, 0.3,length.out = 21), standard = FALSE, verbose = 1)
plot(result1)
# }
Run the code above in your browser using DataLab