plot.DiSSMod {DiSSMod}

Relative Log Likelihood Plot for Discrete Sample Selection Model Fits

Description

plot method for a class "DiSSMod".

Usage

## S3 method for class 'DiSSMod'
plot(x, ...)

Arguments

x

an object of class "DiSSMod" made by the function DiSSMod.

...

additional control argument is as follows.

• level: an option for controlling the significance level of confidence interval. It has to be given in probability between 0 and 1. Initial level is set to 1 - \alpha = 0.95.

Details

Function plot draws a convex line due to the values of twice relative log likelihoods by using the profile likelihood approach with following the grids of alpha. If confidence interval created from the function confint exists between the maximum and minimum value of the alpha, there will be two points drawn with the color red. Also, the Maximum Likelihood Estimator (MLE) of alpha can be seen easily, if it exists between the maximum and minimum value of the alpha.

See Also

See also DiSSMod and plot.

Examples

# 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)