library(Ecdat)
data(Produc)
pdata.frame(Produc,"state","year","pprod")
zw <- pvcm(log(gsp)~log(pcap)+log(pc)+log(emp)+unemp,data=pprod,model="within")
zr <- pvcm(log(gsp)~log(pcap)+log(pc)+log(emp)+unemp,data=pprod,model="random")
Run the code above in your browser using DataLab