lm function on
  transformed data.plm(formula, data, subset, na.action, effect = c("individual","time","twoways"),
    model = c("within","random","ht","between","pooling","fd"),
    random.method = c("swar","walhus","amemiya","nerlove"),
    inst.method = c("bvk","baltagi"), index = NULL, ...)
## S3 method for class 'plm':
summary(object, ...)
## S3 method for class 'summary.plm':
print(x, digits = max(3, getOption("digits") - 2),
    width = getOption("width"), ...)"plm",data.frame,lm,lm,"individual", "time" or "twoways","pooling", "within",
    "between", "random", "fd" and "ht","swar" (the default
    value), "amemiya", "walhus" and "nerlove","bvk" and "baltagi",c("plm","panelmodel").
  A "plm" object  has the following elements :'pFormula' descrbing the model,'pdata.frame' containing the variables used for the
    estimation: the response is in first position and the two indexes in
    the last positions,'ercomp' providing the
    estimation of the components of the errors (for random effect models only),print, summary and print.summary methods.plm is a general function for the estimation of linear
  panel models. It supports the following estimation methods:
  pooled OLS (model="pooling"), fixed effects ("within"),
  random effects ("random"), first--difference ("fd") and
  between ("between"). It supports unbalanced panels and two--ways
  effects (although not with all methods).
  For random effect models, 4 estimators of the transformation
  parameter are available : swar (Swamy and Arora),
  amemiya, walhus (Wallace and Hussain) and nerlove.
  Instrumental variables estimation is obtained using two-parts formula,
  the second part indicating the instrumental variables used. This can
  be a complete list of instrumental variables or an update of the first
  part. If, for example, the model is y~x1+x2+x3, x1,
  x2 are endogenous and z1, z2 are external
  instruments, the model can be estimated with :
  formula=y~x1+x2+x3 | x3+z1+z2,formula=y~x1+x2+x3 | .-x1-x2+z1+z2.inst.method="bvk" or if  inst.method="baltagi".
  
  The Hausman and Taylor estimator is computed if model="ht".data("Produc", package="plm")
zz <- plm(log(gsp)~log(pcap)+log(pc)+log(emp)+unemp, data=Produc, index=c("state","year"))
summary(zz)Run the code above in your browser using DataLab