'YHhmdZddlZddlmZddlmZddlmZdZ dZ Gdd eZ Gd d eZ Gd d eZ y)a Created on Mon Jul 26 08:34:59 2010 Author: josef-pktd changes: added offset and zero-inflated version of Poisson - kind of ok, need better test cases, - a nan in ZIP bse, need to check hessian calculations - found error in ZIP loglike - all tests pass with Issues ------ * If true model is not zero-inflated then numerical Hessian for ZIP has zeros for the inflation probability and is not invertible. -> hessian inverts and bse look ok if row and column are dropped, pinv also works * GenericMLE: still get somewhere (where?) "CacheWriteWarning: The attribute 'bse' cannot be overwritten" * bfgs is too fragile, does not come back * `nm` is slow but seems to work * need good start_params and their use in genericmle needs to be checked for consistency, set as attribute or method (called as attribute) * numerical hessian needs better scaling * check taking parts out of the loop, e.g. factorial(endog) could be precalculated N)stats) factorial)GenericLikelihoodModelcXtjtj||z S)Nnpmaxabsarr1arr2s U/var/www/html/planif/env/lib/python3.12/site-packages/statsmodels/miscmodels/count.pymaxabsr$s 66"&&% &&c^tjtj||z dz S)Nrr s r maxabsrelr's" 66"&&q) **rceZdZdZdZdZy) PoissonGMLEafMaximum Likelihood Estimation of Poisson Model This is an example for generic MLE which has the same statistical model as discretemod.Poisson. Except for defining the negative log-likelihood method, all methods and results are generic. Gradients and Hessian and all resulting statistics are based on numerical differentiation. ctj|j|}|j}tj|||zz tj t |zS| Loglikelihood of Poisson model Parameters ---------- params : array_like The parameters of the model. Returns ------- The log likelihood of the model evaluated at `params` Notes ----- .. math:: \ln L=\sum_{i=1}^{n}\left[-\lambda_{i}+y_{i}x_{i}^{\prime}\beta-\ln y_{i}!\right] )rdotexogendogexplogr)selfparamsXBrs r nloglikeobszPoissonGMLE.nloglikeobs:sK"VVDIIv & vvbzU2X%y/?(@@@rct|dst|j}|j}t j t j ||}tj|dS)zOreturn frozen scipy.stats distribution with mu at estimated prediction resultr)loc) hasattr ValueErrorr#rrrrrpoisson)rrr#rmus rpredict_distributionz PoissonGMLE.predict_distributionOsRtX& [[F]]FtV,-B==+ +rN)__name__ __module__ __qualname____doc__r!r)rrrr,s A* ,rrc*eZdZdZdfd ZdZxZS)PoissonOffsetGMLEauMaximum Likelihood Estimation of Poisson Model This is an example for generic MLE which has the same statistical model as discretemod.Poisson but adds offset Except for defining the negative log-likelihood method, all methods and results are generic. Gradients and Hessian and all resulting statistics are based on numerical differentiation. c |.|jdk(r |dddf}|j|_nd|_t|||fd|i|y)Nrmissing)ndimraveloffsetsuper__init__rrrr6r3kwds __class__s rr8zPoissonOffsetGMLE.__init__lsV  {{a$ ,,.DKDK  g  rc|jtj|j|z}|j}tj |||zz tj t|z}|Sr)r6rrrrrrr)rrr rnlogliks rr!zPoissonOffsetGMLE.nloglikeobs|sX$[[266$))V4 4 &&*b(266)E2B+CCrNNnoner*r+r,r-r8r! __classcell__r;s@rr0r0_s   rr0c*eZdZdZdfd ZdZxZS) PoissonZiGMLEaMaximum Likelihood Estimation of Poisson Model This is an example for generic MLE which has the same statistical model as discretemod.Poisson but adds offset and zero-inflation. Except for defining the negative log-likelihood method, all methods and results are generic. Gradients and Hessian and all resulting statistics are based on numerical differentiation. There are numerical problems if there is no zero-inflation. c d|_t|||f|dgd||.|jdk(r |dddf}|j |_nd|_|&t j|jdf|_ |jjd|_ t jt j|jdf|_ |xjdz c_ dg|_y)Nrzi)r3extra_params_namesr2r start_params)k_extrar7r8r4r5r6ronesnobsrshapenparamshstackrH cloneattrr9s rr8zPoissonZiGMLE.__init__s   3g$(6 3-1 3  {{a$ ,,.DKDK <1 .DIyyq) IIrwwt||'/BFF2J.%(:RVVIeDT=UU ffTBFFGE1H4E3E,F%FGGqrr>r@rBs@rrDrDs *6rrD)r-numpyrscipyr scipy.specialrstatsmodels.base.modelrrrrr0rDr.rrrXsM:#9'+ /,(/,f2.2hD*Dr