📄 gompit.m
字号:
function [LIK, ssr, yhat, pderiv] = gompit(beta);
global P T;
[rr cc] = size(T);
[rx, cx] = size(P);
y = T; x = P;
beta1 = beta;
nn = P * beta1(1:end-1)' + ones(length(P),1) * beta1(end);
yhat = exp(-exp(-nn));
ssr = (T - yhat)' * (T-yhat);
y = T;
lik = real(y .* log(yhat) + (1-y) .* log(1-yhat));
LIK = -sum(lik);
% lik = yhat .^ y .* (1-yhat) .^ (1-y);
% LIK = prod(lik);
xmean = mean(x);
yhatmean = xmean * beta1(1:end-1)' + beta(end);
for i= 1:cx,
pderiv(i) = mean(exp(yhat) .* exp(-yhatmean) .* beta1(i));
end
⌨️ 快捷键说明
复制代码
Ctrl + C
搜索代码
Ctrl + F
全屏模式
F11
切换主题
Ctrl + Shift + D
显示快捷键
?
增大字号
Ctrl + =
减小字号
Ctrl + -