ó áp7]c @sĂdZddlmZddlZddlmZddlmZddl m Z ddl j Z ddlmZmZddlmZd e fd „ƒYZed d d dddddddƒZd:\ZZZejdƒZdddddgZx‰defdefdefd efd!efgD]T\ZZeeeeƒeeeeƒƒeeejeƒee eeƒƒq+Weejeeed"d#ƒƒeejje eed"d#ƒƒeejeeeƒƒeejje eeƒƒd$„Z d%„Z!d&„Z"d'„Z#dd(d)d*dd+„Z%ee#dd,d)ƒƒee#dd-d)ƒƒej&d.d/dd0d1ƒZ'e%e'd2d;d3dƒZ(ee(dd4ƒee(dd4ƒd5„Z)d6„Z*e)d7ejdd)ƒƒ\Z+Z,Z-ee+ƒee,ƒee-ƒee*ej.d8dƒƒee*ej.d8ej d,ddgƒƒej.j&d0d9ƒZ&e&e&j/ƒZ&ee&j/ƒe&e&d,kj/ƒe&e&dkj/ƒe&e&dkj/ƒƒdS(<se Created on Thu Aug 12 14:59:03 2010 Warning: not tried out or tested yet, Done Author: josef-pktd i˙˙˙˙(tprint_functionN(tstats(tcomb(t rv_continuous(twheretinf(tabstgenpareto2_gencBsGeZd„Zd„Zd„Zd„Zd„Zd„Zd„ZRS(cCsJtj|ƒ}t|dkdt|ƒtƒ|_t|dkddƒS(Nigđ?i(tnptasarrayRtnp_absRtb(tselftc((sJlib/python2.7/site-packages/statsmodels/sandbox/distributions/genpareto.pyt _argchecks%cCs&tjd||dd|ƒ}|S(Nigđżgđ?(Rtpower(R txR tPx((sJlib/python2.7/site-packages/statsmodels/sandbox/distributions/genpareto.pyt_pdfs"cCsdd|tj| |ƒS(Ngđżgđ?(Rtlog1p(R RR ((sJlib/python2.7/site-packages/statsmodels/sandbox/distributions/genpareto.pyt_logpdfscCs dtjd||d|ƒS(Ngđ?i(RR(R RR ((sJlib/python2.7/site-packages/statsmodels/sandbox/distributions/genpareto.pyt_cdfscCs&d|tjd||ƒd}|S(Ngđżi(RR(R tqR tvals((sJlib/python2.7/site-packages/statsmodels/sandbox/distributions/genpareto.pyt_ppf s"cCsntjd|dƒ}d||tjt||ƒd|d||ddƒ}t||dk|tƒS(Niigđ?i˙˙˙˙taxis(RtarangetsumRRR(R tnR tktval((sJlib/python2.7/site-packages/statsmodels/sandbox/distributions/genpareto.pyt_munp#s>cCs5|dkrd|Sd||_tj||ƒSdS(Niigđ?(R Rt_entropy(R R ((sJlib/python2.7/site-packages/statsmodels/sandbox/distributions/genpareto.pyR 's  ( t__name__t __module__RRRRRRR (((sJlib/python2.7/site-packages/statsmodels/sandbox/distributions/genpareto.pyRs      tagtnamet genparetotlongnamesA generalized ParetotshapesR textradocs Generalized Pareto distribution genpareto2.pdf(x,c) = (1+c*x)**(-1-1/c) for c != 0, and for x >= 0 for all c, and x < 1/abs(c) for c < 0. gŕ?iiig{ŽGáz„?gš™™™™™š?gÍĚĚĚĚĚě?gŽGázŽď?tpdftcdftsftppftisftmomentstmvskcCs||||fS(sŻtransform shape scale for peak over threshold y = x-u|x>u ~ GPD(k, sigma-k*u) if x ~ GPD(k, sigma) notation of de Zea Bermudez, Kotz k, sigma is shape, scale ((tthreshtshapetscale((sJlib/python2.7/site-packages/statsmodels/sandbox/distributions/genpareto.pyt paramstopotIscCs||||fS(N((R0R1tscalepot((sJlib/python2.7/site-packages/statsmodels/sandbox/distributions/genpareto.pyt paramsfrompotRscCs|st|dƒndS(Ns does not hold(tprint(tcondtmsg((sJlib/python2.7/site-packages/statsmodels/sandbox/distributions/genpareto.pytwarnifUscCsUt|dkdƒt|dkdƒt|||dkdƒ|||d|S(semean excess function of genpareto assert are inequality conditions in de Zea Bermudez, Kotz i˙˙˙˙s shape > -1is thresh >= 0s(scale - shape*thresh) > 0i(R9(R0R1R2((sJlib/python2.7/site-packages/statsmodels/sandbox/distributions/genpareto.pyt meanexcessYsidi tempc Cs |dkr$|dkrt‚qnÝ|dkrtj|ƒ}|ddd…jƒtjdt|ƒdƒ|ddd…}|ddd…}|rtj|| || ƒ|dk rţ|\}} tj|| | || |d|ƒqţqn||fS(NtestR;i˙˙˙˙igđ?( tNonetNotImplementedErrorRtsorttcumsumRtlentplttplot( tdatatparamstlidxtuidxtmethodRCt datasortedR:R1R2((sJlib/python2.7/site-packages/statsmodels/sandbox/distributions/genpareto.pytmeanexcess_plotds    D  3gŕżiţ˙˙˙gčżR2tsizeičRERCiě˙˙˙cCsźtj|ƒjtƒ}|ddd…jƒtjdt|ƒdƒ|ddd…}|ddd…jƒtjdt|ƒdƒ}|ddd…}|||ddd…fS(Ni˙˙˙˙i(RR?tastypetfloatR@RRA(RDRIR:tmeancont((sJlib/python2.7/site-packages/statsmodels/sandbox/distributions/genpareto.pytmeanexcess_emp…s D3cOsatj|ƒdkr+|jd|dtƒStjg|D]}|jd|dtƒ^q8ƒSdS(Nitlbt conditional(RtndimtexpecttTruetarray(R RPtargstkwdstlbb((sJlib/python2.7/site-packages/statsmodels/sandbox/distributions/genpareto.pytmeanexcess_distŒs gđ?RPi †(gŕ?ii(gčżi(0t__doc__t __future__RtnumpyRtscipyRt scipy.specialRtscipy.stats.distributionsRtmatplotlib.pyplottpyplotRBRRRR Rt genpareto2R1tlocR2RtrvtquantRHRR6tgetattrR%tentropyR3R5R9R:R=RJtrvsRDttmpRORYtdstmetmctnormtmean(((sJlib/python2.7/site-packages/statsmodels/sandbox/distributions/genpareto.pytsd      "*#     %   )IM