cmtest: Conditional moments test

View source: R/cmtest.R

cmtestR Documentation

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 class 'tobit'
cmtest(
  x,
  test = c("normality", "reset", "heterosc", "skewness", "kurtosis"),
  powers = 2:3,
  heter_cov = NULL,
  opg = FALSE
)

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

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

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

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

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

Author(s)

Yves Croissant

References

\insertRef

NEWE:85micsr

\insertRef

PAGA:VELL:89micsr

\insertRef

TAUC:85micsr

\insertRef

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)

micsr documentation built on May 29, 2024, 7:32 a.m.