R语言如何估计带参数约束的方程？

如附件，我要估计附件中的3.1节的这几个方程。数据已经准备好了，我不会联合估计两个方程，且附带β的约束条件（方程3）

在R论坛上找到一个translogCostEst函数专门估计Translog 的成本函数，可惜它的{homPrice = TRUE }功能尚未开发出来。
estResult2 <- translogCostEst( cName = "COST", yName = "GDP",
                              pNames = c( "PL", "PE","PK" ),
                              data = dat, fNames = NULL,
                              shifterNames = NULL, dataLogged = FALSE, homPrice = TRUE )
总之，希望大家能够指点一下，怎么能估计出βke，βkl，βlk 等参数，它们处于不同的方程，还要满足约束条件。

万分感谢！

Elasticities of substitution between energy an.pdf (302.87 KB)

2015-12-15 11:25:52 上传

Final Data -2015-12-15-09-35-13.xlsx (12.94 KB)

2015-12-15 11:25:27 上传