plm (version 2.6-4)

pldv: Panel estimators for limited dependent variables

Description

Fixed and random effects estimators for truncated or censored limited dependent variable

Usage

pldv(
  formula,
  data,
  subset,
  weights,
  na.action,
  model = c("fd", "random", "pooling"),
  index = NULL,
  R = 20,
  start = NULL,
  lower = 0,
  upper = +Inf,
  objfun = c("lsq", "lad"),
  sample = c("cens", "trunc"),
  ...
)

Value

For model = "fd", an object of class c("plm", "panelmodel"), for model = "random" and model = "pooling" an object of class c("maxLik", "maxim").

Arguments

formula

a symbolic description for the model to be estimated,

data

a data.frame,

subset

see lm,

weights

see lm,

na.action

see lm,

model

one of "fd", "random", or "pooling",

index

the indexes, see pdata.frame(),

R

the number of points for the gaussian quadrature,

start

a vector of starting values,

lower

the lower bound for the censored/truncated dependent variable,

upper

the upper bound for the censored/truncated dependent variable,

objfun

the objective function for the fixed effect model (model = "fd", irrelevant for other values of the model argument ): one of "lsq" for least squares (minimise sum of squares of the residuals) and "lad" for least absolute deviations (minimise sum of absolute values of the residuals),

sample

"cens" for a censored (tobit-like) sample, "trunc" for a truncated sample,

...

further arguments.

Author

Yves Croissant

Details

pldv computes two kinds of models: a LSQ/LAD estimator for the first-difference model (model = "fd") and a maximum likelihood estimator with an assumed normal distribution for the individual effects (model = "random" or "pooling").

For maximum-likelihood estimations, pldv uses internally function maxLik::maxLik() (from package maxLik).

References

HONO:92plm

Examples

## as these examples take a bit of time, do not run them automatically
if (FALSE) {
data("Donors", package = "pder")
library("plm")
pDonors <- pdata.frame(Donors, index = "id")

# replicate Landry/Lange/List/Price/Rupp (2010), online appendix, table 5a, models A and B
modA <- pldv(donation ~ treatment +  prcontr, data = pDonors,
            model = "random", method = "bfgs")
summary(modA)
modB <- pldv(donation ~ treatment * prcontr - prcontr, data = pDonors,
            model = "random", method = "bfgs")
summary(modB)
}