micsr (version 0.1-1)

cmtest: Conditional moments test

Description

Conditional moments tests for maximum likelihood estimators, particularly convenient for the probit and the tobit model to test relevance of functional form, omitted variables, heteroscedasticity and normality.

Usage

cmtest(
  x,
  test = c("normality", "reset", "heterosc", "skewness", "kurtosis"),
  powers = 2:3,
  heter_cov = NULL,
  opg = FALSE
)
# S3 method for tobit
cmtest(
  x,
  test = c("normality", "reset", "heterosc", "skewness", "kurtosis"),
  powers = 2:3,
  heter_cov = NULL,
  opg = FALSE
)
# S3 method for micsr
cmtest(
  x,
  test = c("normality", "reset", "heterosc", "skewness", "kurtosis"),
  powers = 2:3,
  heter_cov = NULL,
  opg = FALSE
)
# S3 method for censReg
cmtest(
  x,
  test = c("normality", "reset", "heterosc", "skewness", "kurtosis"),
  powers = 2:3,
  heter_cov = NULL,
  opg = FALSE
)
# S3 method for glm
cmtest(
  x,
  test = c("normality", "reset", "heterosc", "skewness", "kurtosis"),
  powers = 2:3,
  heter_cov = NULL,
  opg = FALSE
)

Value

an object of class "htest" containing the following components:

• data.mane: a character string describing the fitted model

• statistic: the value of the test statistic

• parameter: degrees of freedom

• p.value: the p.value of the test

• method: a character indicating what type of test is performed

Arguments

x

a fitted model, currently a tobit model either fitted by AER::tobit, censReg::censReg or micsr::tobit1 or a probit model fitted by glm with family = binomial(link = "probit") or by micsr::binomreg with link = "probit"

test

the kind of test to be performed, either a normality test (or separately a test that the skewness or kurtosis are 0 and 3), a heteroscedasticity test or a reset test,

powers

the powers of the fitted values that should be used in the reset test,

heter_cov

a one side formula that indicates the covariates that should be used for the heteroscedasticity test (by default all the covariates used in the regression are used),

opg

a boolean, if FALSE (the default), the analytic derivatives are used, otherwise the outer product of the gradient formula is used

Author

Yves Croissant

References

NEWE:85micsr

PAGA:VELL:89micsr

TAUC:85micsr

WELL:03micsr

Examples

charitable$logdon <- with(charitable, log(donation) - log(25))
ml <- tobit1(logdon ~ log(donparents) + log(income) + education +
             religion + married + south, data = charitable)
cmtest(ml, test = "heterosc")
cmtest(ml, test = "normality", opg = TRUE)