* open log file
log using cig_industry.log, replace

* open stata data file
use cig_industry

* get description of data
desc

* generate logs of all variables

gen kl=ln(k)
gen ll=ln(l)
gen ql=ln(q)

* run unrestricted model
reg ql ll kl

* test for constant returns to scale
test kl+ll=1

* generate variables for brute force model
gen ql_kl=ql-kl
gen ll_kl=ll-kl

* run restricted model
reg ql_kl ll_kl

log close