data("Wages",package="plm")
ht <- plm(lwage~wks+south+smsa+married+exp+I(exp^2)+
bluecol+ind+union+sex+black+ed |
sex+black+bluecol+south+smsa+ind,
data=Wages,model="ht",index=595)
summary(ht)
Run the code above in your browser using DataLab