library(plm)
data(RiceFarms)

#head(RiceFarms)

price <- pdata.frame(RiceFarms, index=c("id"))

fm <- log(goutput) ~ log(seed) + log(urea) + I(phosphate/1000) + log(totlabor) + log(size) + pesticide + varieties + region

olsmod <- plm(fm, price, model="pooling")
remod <- plm(fm, price, model="random")
femod <- plm(fm, price, model="within")

library(texreg)
screenreg(list(OLS=olsmod, RE=remod, FE=femod))

phtest(femod, remod)

plmtest(fm, price)