micEcon (version 0.4-0)

probit: Probit binary choice model

Description

The functions calculates probit binary choice model

Usage

probit(formula, subset, start = NULL, data = sys.frame(sys.parent()),
       x=FALSE, y = FALSE, model = FALSE, method="ML", ... )

Arguments

formula

a symbolic description of the model to be fit, in the form response ~ explanatory variables. The response variable should be either logical or numeric with values only 0-1 or a factor with two levels.

subset

an optional vector specifying a subset of observations to be used in the fitting process.

start

inital value of parameters.

data

an optional data frame containing the variables in the model. If not found in data, the variables are taken from environment(formula), typically the environment from which probit is called.

x, y, model

logicals. If TRUE the corresponding components of the fit (the model matrix, the response, the model frame) are returned.

method

the method to be used; for fitting, currently only method = "ML" (Maximum Likelihood) is supported; method = "model.frame" returns the model frame (the same as with model = TRUE, see below).

...

further arguments for maxLik.

Value

• An object of class "probit". It is a list with following components:
• LRT
  Likelihood ration test. The full model is tested against H0: the parameters have no effect on the result. This is a list with components LRT{The LRT value} df{Degrees of freedom for LRT (= df of the model - 1)} LRT is distributed by chi2(df) under H0.
• nParam
  Number of parameters of the model including constant
• nObs
  Number of the observations
• N1
  Number of observations with non-zero (true) response
• N0
  Number of observations with zero (false) response
• df
  Number of free parameters
• x
  if requested, the model matrix used.
• y
  if requested, the response used.
• model
  the model frame, only if model = TRUE or method = "model.frame".
• Other components are inherited from maxLik.

Details

The model is estimated using Maximum Likelihood method.

See Also

maxLik for ready-packaged likelihood maximisation routines and methods, glm for generalised linear models, including probit, binomial.

Examples

## A simple MC trial: note probit assumes normal errors
x <- runif(100)
e <- 0.5*rnorm(100)
y <- x + e
summary(probit((y > 0) ~ x))
## female labour force participation probability
data(Mroz87)
Mroz87$kids <- Mroz87$kids5 > 0 | Mroz87$kids618 > 0
Mroz87$age30.39 <- Mroz87$age < 40
Mroz87$age50.60 <- Mroz87$age >= 50
summary(probit(lfp ~ kids + age30.39 + age50.60 + educ + hushrs +
               huseduc + huswage + mtr + motheduc, data=Mroz87))