如何用ml估計CES生產函數

是 否 +2 论坛币

k人 参与回答

经管之家送您一份

应届毕业生专属福利!

求职就业群

赵安豆老师微信：zhaoandou666

经管之家联合CDA

送您一个全额奖学金名额~ !

立即领取

感谢您参与论坛问题回答

经管之家送您两个论坛币！

+2 论坛币

這是要估計的函數形式，nl可以成功估計beta0,rho,delta等參數(http://www.stata.com/capabilities/nlreg.html)
，stata官方網站有示範以ml的方式估計非線性(http://www.stata.com/support/faqs/stat/nl_ml.html)
試著仿造過程遇到幾個難處：
⑴  delta k^{-rho} 的那個地方有些不解。我看書上類似把k視實際變數，而非代表變數??
(2)  是否可以利用 local logout: word 1 of `varlist' 的語法在輸入 y =k l時
stata 將右手式的第一個式為k(capital)，第二個變數視為l(labor)
可以請各位幫我看看是哪邊的代碼需要修改。

codes:
/*ces dgp*/
clear
setseed 9876
setobs 1000
genu=rnormal()
genk=runiform()
genl=runiform()
scalarrho = 0.2
scalardelta = 0.3
scalar b0 = 0.4

genq = b0 - (1/rho)*log(delta*k^(-1*rho) + (1-delta)*l^(-1*rho)) + u

/*ml ces function */

captureprogram drop mlces
programdefine mlces,rclass

version10

args lnf b0 rho delta

    syntax varlist(min=3 max=3) [aw fw iw] if, at(name)
    local logout: word 1 of `varlist'
    local capital: word 2 of `varlist'
    local labor: word 3 of `varlist'

   // Retrieve parameters out of at
    tempname b0 rho delta
    scalar `b0' = `at'[1,1]
    scalar `rho' = `at'[1,2]
    scalar `delta' = `at'[1,3]

tempvar res
qui gen double `res' = $ML_y1 -(`b0'-(1/`rho’)*(`delta*k^(-`rho')+ (1-`delta’)*l^(-`rho')))
qui replace `lnf' = log(normalden($ML_y1,`res',1))
end

ml model lf mlces (q=l k) /b0 /rho /delta
ml max

扫码加我 拉你入群

请注明：姓名-公司-职位

以便审核进群资格，未注明则拒绝

分享0 收藏2 回帖

关键词：CES 如何用 CAPABILITIES parameters Parameter capital logout local 如何

顶楼主，我也在研究这个问题，不过没有解决