s*********e
发帖数: 1051 | 1 does anyone here have experience in simplex distribution?
I drafted a piece of nlmixed code based upon the density function but am not
sure it is correct or not.
proc nlmixed data = tmp1 technique = NRRIDG ;
parms b0 = 0 b1 = 0.1 sigma = 5;
bounds sigma > 0;
u = 1 / (1 + exp(-(b0 + b1 * x)));
dd = (y - u) ** 2 / (y * (1 - y) * u ** 2 * (1 - u) ** 2);
lh = (2 * constant('pi') * sigma * (y * (1 - y)) ** 3) ** (-0.5) * exp(-(2
* sigma) ** (-1) * dd);
ll = log(lh);
model y ~ general(ll);
run;
here, the mean parameter is a linear component and the dispersion parameter
is a constant.
any insight is highly appreciated. |
s*********e
发帖数: 1051 | 2 i got the confirmation from matt flynn in travelers that my code is correct.
thx for your attention.
not
(2
【在 s*********e 的大作中提到】
: does anyone here have experience in simplex distribution?
: I drafted a piece of nlmixed code based upon the density function but am not
: sure it is correct or not.
: proc nlmixed data = tmp1 technique = NRRIDG ;
: parms b0 = 0 b1 = 0.1 sigma = 5;
: bounds sigma > 0;
: u = 1 / (1 + exp(-(b0 + b1 * x)));
: dd = (y - u) ** 2 / (y * (1 - y) * u ** 2 * (1 - u) ** 2);
: lh = (2 * constant('pi') * sigma * (y * (1 - y)) ** 3) ** (-0.5) * exp(-(2
: * sigma) ** (-1) * dd);
|
d*******o
发帖数: 493 | 3 statcompute 和 matt是用proc nlmixed的大牛啊 |
s*********e
发帖数: 1051 | 4 he is. iam not.
【在 d*******o 的大作中提到】
: statcompute 和 matt是用proc nlmixed的大牛啊
|
