作者s0935018 (ko)
看板Statistics
標題[問題] stata計算最適化問題
時間Wed Jul 25 15:39:29 2018
大家好，我的資料型態舉例如下:
id year sales cogs sga ppe rd gw oita
1101 2002 1000 2000 3000 4000 5000 6000 7000
1101 2003 1500 2500 3000 4000 5500 6500 7500
....
1102 2002 3200 4300 4800 5400 3400 5600 7300
1102 2003 2900 2300 2300 2300 6300 4300 3300
....
我想利用文獻的公式，求出每公司每年的營運效率
公式: max=sales/(v1*cogs+v2*sga+v3*ppe+v4*rd+v5*gw+v6*oita)
v1~v6為非負之權重會使整個式子最大化，
我一開始想說求整個式子最大化，是不是可以等於求分母最小化，
嘗試用optimize，先求第一家第一年的營運效率:
mata:
void plleval(todo,b, val, g, H)
{
real vector zb,Z
real matrix X
X=st_data(.,"cogs sga ppe rd gw oita")
Z=X[1,(1::6)]
zb=Z*b'
val=zb
}
S=optimize_init()
optimize_init_evaluator(S,&plleval())
optimize_init_params(S,J(1,6, 0.001))
optimize_init_which(S,"min")
b=optimize(S)
optimize_result_params(S)
b
end
結果stata回報"3499 S not found"
網路上的例子比較多是解幾元幾次方程式，
不太清楚權重在optimize設定是否一樣，還請大家指正，
先謝謝大家。
--
※ 發信站: 批踢踢實業坊(ptt.cc), 來自: 1.160.56.20
※ 文章網址: https://webptt.com/m.aspx?n=bbs/Statistics/M.1532504371.A.2E6.html