This function enables the estimation of the variance components of a panel model.
ercomp(object, ...)
# S3 method for class 'plm'
ercomp(object, ...)
# S3 method for class 'pdata.frame'
ercomp(
object,
effect = c("individual", "time", "twoways", "nested"),
method = NULL,
models = NULL,
dfcor = NULL,
index = NULL,
...
)
# S3 method for class 'formula'
ercomp(
object,
data,
effect = c("individual", "time", "twoways", "nested"),
method = NULL,
models = NULL,
dfcor = NULL,
index = NULL,
...
)
# S3 method for class 'ercomp'
print(x, digits = max(3, getOption("digits") - 3), ...)a formula or a plm object,
further arguments.
the effects introduced in the model, see plm() for
details,
method of estimation for the variance components, see
plm() for details,
the models used to estimate the variance components (an alternative to the previous argument),
a numeric vector of length 2 indicating which degree of freedom should be used,
the indexes,
a data.frame,
an ercomp object,
digits,
An object of class "ercomp": a list containing
sigma2 a named numeric with estimates of the variance
components,
theta contains the parameter(s) used for
the transformation of the variables: For a one-way model, a
numeric corresponding to the selected effect (individual or
time); for a two-ways model a list of length 3 with the
parameters. In case of a balanced model, the numeric has length
1 while for an unbalanced model, the numerics' length equal the
number of observations.
Amemiya T (1971). “The Estimation of the Variances in a Variance–Components Model.” International Economic Review, 12, 1–13.
Nerlove M (1971). “Further Evidence on the Estimation of Dynamic Economic Relations from a Time–Series of Cross–Sections.” Econometrica, 39, 359–382.
Swamy PAVB, Arora SS (1972). “The Exact Finite Sample Properties of the Estimators of Coefficients in the Error Components Regression Models.” Econometrica, 40, 261–275.
Wallace TD, Hussain A (1969). “The Use of Error Components Models in Combining Cross Section With Time Series Data.” Econometrica, 37(1), 55–72.
plm() where the estimates of the variance components are
used if a random effects model is estimated
data("Produc", package = "plm")
# an example of the formula method
ercomp(log(gsp) ~ log(pcap) + log(pc) + log(emp) + unemp, data = Produc,
method = "walhus", effect = "time")
#> var std.dev share
#> idiosyncratic 0.0075942 0.0871449 0.985
#> time 0.0001192 0.0109175 0.015
#> theta: 0.2448
# same with the plm method
z <- plm(log(gsp) ~ log(pcap) + log(pc) + log(emp) + unemp,
data = Produc, random.method = "walhus",
effect = "time", model = "random")
ercomp(z)
#> var std.dev share
#> idiosyncratic 0.0075942 0.0871449 0.985
#> time 0.0001192 0.0109175 0.015
#> theta: 0.2448
# a two-ways model
ercomp(log(gsp) ~ log(pcap) + log(pc) + log(emp) + unemp, data = Produc,
method = "amemiya", effect = "twoways")
#> var std.dev share
#> idiosyncratic 0.0011695 0.0341975 0.046
#> individual 0.0238635 0.1544780 0.929
#> time 0.0006534 0.0255613 0.025
#> theta: 0.9464 (id) 0.8104 (time) 0.8084 (total)