felogit_lik.m

来自「时间序列分析中常用到的matlab代码」· M 代码 · 共 18 行

function [P,lnL] = felogit_lik(y,x,b,c,d)
% PURPOSE: Compute probabilities and value of log-likelihood
%   for fixed effects logit model. Called from felogit()

% written by:
% Simon D. Woodcock
% CISER / Economics
% Cornell University
% 201 Caldwell Hall
% Ithaca, NY 14850
% [e-mail]

index = x*b + d*c;
e_index = exp(index);
P = e_index ./ ( 1 + e_index);
lnL = sum( y.*index - log(1+e_index) );