Variables | |
function [dl, d2l] | |
endif | v = g .* (1 - g) ./ p |
v1 = g1 .* (1 - g1) ./ p | |
dlogp = [(diag (v) * z - diag (v1) * z1), (diag (v - v1) * x)] | |
dl = sum (dlogp)' | |
w = v .* (1 - 2 * g) | |
w1 = v1 .* (1 - 2 * g1) | |
d2l |
endif v = g .* (1 - g) ./ p |
w = v .* (1 - 2 * g) |