我有一组买单和买单的数据,并能写出一般形式的对数似然函数,其中包含5个参数;想通过极大似然呢估计来得到这5个参数的估计值,请教高手我该如何做?
下面尝试写了一个粗浅的程序,但无法正确运行,不知nlmixed 能否达到要求?
data PIN;
input buy sell @@;
datalines;
123 140
78 103
106 100
212 231
187 152
128 192
49 70
0 2
62 147
68 106
62 120
137 204
85 220
11 27
93 113
48 124
36 113
32 138
24 137
35 192
;
proc nlmixed data=PIN;
parms alpha = 0.3 delta = 0.1 eb = 0.9 es =1.1 mu = 0.6;
bounds alpha > 0;
f = (1-alpha) * pdf('POISSON',buy, eb) * pdf('POISSON',sell, es) + alpha * delta * pdf('POISSON',buy, eb) * pdf('POISSON',sell, es+mu) + alpha * (1-delta) * pdf('POISSON',buy, eb+mu) * pdf('POISSON',sell, es);
logf=log(f);
model buy~ general(logf);
run;
恳求热心人士帮忙,十分感谢!