作者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/cn.aspx?n=bbs/Statistics/M.1532504371.A.2E6.html