proc nlmixed data=lab5;
parms mu 0.9 sigma 1 gamma0 0.5 gamma1 -0.5 ;
bounds 0 < sigma;
sqrt2pi = sqrt(2*constant('pi'));
eta = exp((gamma0 + gamma1*y)) / (1 + exp((gamma0 + gamma1 * y)));
f1 = -log(sigma) - log(sqrt2pi*y) - (log(y)-mu)**2 / (2*sigma**2);
f2 = (eta**(d=1) )* ((1 - eta)**(d=2)); /*d is binary variable, 1 or 2*/
ll= log (f1 * f2);
model y~general(ll);
run;
i want to calculate the parameters for joint distribution.
f(y,d) = f(y)f(d|y)
f(y) follows a normal distribution, f(d|y) follow Bernoulli distribution. Then I wrote the syntax above.
however, there was no results and NOTE: Execution error for observation 1.
is there anyone know the problem and how to deal with it?
... View more