Ë
    £�Dj©)  ã            	       óF  — d Z ddlZddlmZ ddlmZ ddlmZ ddl	m
Z ddlmZmZ ddlmZ  G d„ d	e«      Z ed
ddd¬«      Zd\  ZZZ ej*                  d«      Zg d¢ZdefdefdefdefdeffD ]@  \  ZZ e  eee«      eeee«      «        e  eej8                  e«      ee ee«      «       ŒB  eej	                  eeed¬«      «        eej8                  j	                  e eed¬«      «        eej;                  eee«      «        eej8                  j;                  e ee«      «       d„ Zd„ Zd„ Z d„ Z!d/d„Z" e e!ddd«      «        e e!dd d«      «       ejG                  d!dd"¬#«      Z$ e"e$d$d%¬&«      Z% ee%d%   d'd «        ee%d   d'd «       d(„ Z&d)„ Z' e&d* ej*                  d%d«      z  «      \  Z(Z)Z* ee(«        ee)«        ee*«        e e'ejV                  d+¬,«      «        e e'ejV                  ej                   ddd+g¬,«      «       ejV                  jG                  d-¬.«      Z#e#e#jY                  «       z
  Z# ee#jY                  «       e#e#dkD     jY                  «       e#e#dkD     jY                  «       e#e#d+kD     jY                  «       «       y)0ze
Created on Thu Aug 12 14:59:03 2010

Warning: not tried out or tested yet, Done

Author: josef-pktd
é    N)Ústats)Úcomb)Úrv_continuous)ÚwhereÚinf)Úabsc                   ó6   — e Zd Zd„ Zd„ Zd„ Zd„ Zd„ Zd„ Zd„ Z	y)	Úgenpareto2_genc                 ó–   — t        j                  |«      }t        |dkD  dt        |«      z  t        «      | _        t        |dk(  dd«      S )Nr   ç      ð?é   )ÚnpÚasarrayr   Únp_absr   Úb©ÚselfÚcs     úoC:\Crop_Prediction\Backend\crop-ai-system\venv\Lib\site-packages\statsmodels/sandbox/distributions/genpareto.pyÚ	_argcheckzgenpareto2_gen._argcheck   s>   € Ü�J‰J�q‹MˆÜ�q˜1‘u˜c¤F¨1£I™o¬sÓ3ˆŒÜ�Q˜!‘V˜Q Ó"Ð"ó    c                 óJ   — t        j                  d||z  z
  dd|z  z   «      }|S )Nr   ç      ð¿r   ©r   Úpower)r   Úxr   ÚPxs       r   Ú_pdfzgenpareto2_gen._pdf   s'   € Ü�X‰X�a˜!˜a™%‘i ¨¨a©¡Ó0ˆØˆ	r   c                 óF   — dd|z  z   t        j                  | |z  «      z  S )Nr   r   )r   Úlog1p©r   r   r   s      r   Ú_logpdfzgenpareto2_gen._logpdf   s$   € Ø�s˜Q‘w‘¤"§(¡(¨A¨2°©6Ó"2Ñ2Ð2r   c                 óF   — dt        j                  d||z  z
  d|z  «      z
  S )Nr   r   r   r!   s      r   Ú_cdfzgenpareto2_gen._cdf    s$   € Ø”R—X‘X˜a ! a¡%™i¨¨q©Ó1Ñ1Ð1r   c                 óJ   — d|z  t        j                  d|z
  |«      dz
  z  }|S )Nr   r   r   )r   Úqr   Úvalss       r   Ú_ppfzgenpareto2_gen._ppf#   s)   € Ø�a‰xœ2Ÿ8™8 A¨¡E¨1Ó-°Ñ1Ñ2ˆØˆr   c                 óÔ   — t        j                  d|dz   «      }d|z  |z  t        j                  t        ||«      d|z  z  d||z  z   z  d¬«      z  }t	        ||z  dkD  |t
        «      S )Nr   r   r   éÿÿÿÿ)Úaxis)r   ÚarangeÚsumr   r   r   )r   Únr   ÚkÚvals        r   Ú_munpzgenpareto2_gen._munp'   sg   € Ü�I‰I�a˜˜Q™ÓˆØ�Q‰w˜1‰nœrŸv™v¤d¨1¨a£j°B¸1±9Ñ&<ÀÀaÈ!ÁeÁÑ&LÐSTÔUÑUˆÜ�Q˜‘U˜R‘Z ¤cÓ*Ð*r   c                 óV   — |dk  rd|z
  S d|z  | _         t        j                  | |«      S )Nr   r   r   )r   r   Ú_entropyr   s     r   r3   zgenpareto2_gen._entropy,   s0   € Ø�ŠEØ�q‘5ˆLà˜1‘WˆDŒFÜ ×)Ñ)¨$°Ó2Ð2r   N)
Ú__name__Ú
__module__Ú__qualname__r   r   r"   r$   r(   r1   r3   © r   r   r
   r
      s%   „ ò#ò
ò3ò2òò+ó
3r   r
   g        Ú	genparetozA generalized Paretor   )ÚaÚnameÚlongnameÚshapes)ç      à?r   r   é   )g{®Gáz„?gš™™™™™¹?r=   gÍÌÌÌÌÌì?g®Gáz®ï?ÚpdfÚcdfÚsfÚppfÚisfÚmvsk)Úmomentsc                 ó   — |||| z  z
  fS )z¯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
    r7   ©ÚthreshÚshapeÚscales      r   ÚparamstopotrK   Q   s   € ð �%˜% &™.Ñ(Ð(Ð(r   c                 ó   — |||| z  z   fS )Nr7   )rH   rI   Úscalepots      r   ÚparamsfrompotrN   [   s   € Ø�(˜U V™^Ñ+Ð+Ð+r   c                 ó"   — | st        |d«       y y )Nzdoes not hold)Úprint)ÚcondÚmsgs     r   ÚwarnifrS   _   s   € ÙÜˆc�?Õ#ð r   c                 ó„   — t        |dkD  d«       t        | dk\  d«       t        ||| z  z
  dkD  d«       ||| z  z
  d|z   z  S )zemean excess function of genpareto

    assert are inequality conditions in de Zea Bermudez, Kotz
    r*   z
shape > -1r   zthresh >= 0z(scale - shape*thresh) > 0r   )rS   rG   s      r   Ú
meanexcessrU   d   sQ   € ô
 ˆ5�2‰:�|Ô$Ü
ˆ6�Q‰;˜Ô&Ü
ˆE�E˜F‘NÑ" aÑ'Ð)EÔFØ�E˜F‘NÑ" q¨5¡yÑ1Ð1r   é
   c                 ó˜  — |dk(  r|€t         ‚	 fS |dk(  r°t        j                  | «      }|d d d…   j                  «       t        j                  dt        | «      dz   «      z  |d d d…   z
  }|d d d…   }|rOt        j                  |d |  |d |  «       |�/|\  }}	t        j                  |d |  |	|d |  |z  z
  d|z   z  «       fS )NÚestÚempr*   r   r   )ÚNotImplementedErrorr   ÚsortÚcumsumr,   ÚlenÚpltÚplot)
ÚdataÚparamsÚlidxÚuidxÚmethodr_   Ú
datasortedrU   rI   rJ   s
             r   Úmeanexcess_plotrf   o   sþ   € Ø�‚ð ˆ>Ü%Ð%àð �zÐ!Ð!ð 
�5Šä—W‘W˜T“]ˆ
Ø ¡ 2 Ñ&×-Ñ-Ó/´2·9±9¸QÄÀDÃ	ÈAÁÓ3NÑNÐQ[Ñ\`Ð^`Ð\`ÑQaÑaˆ
Ø¡ " Ñ%ˆ
ÙÜ�H‰H�Z  $ Ð'¨°F°d°UÐ);Ô<ØÐ!Ø%‘��uÜ—‘˜ F d UÐ+¨e°jÀÀ4À%Ð6HÈ5Ñ6PÑ.PÐUWÐZ_ÑU_Ñ-`ÔaØ�zÐ!Ð!r   g      à¿éþÿÿÿç      è¿iè  )rJ   Úsize)rh   r>   r   )ra   r_   iìÿÿÿc                 ón  — t        j                  | «      j                  t        «      }|d d d…   j	                  «       t        j
                  dt        | «      dz   «      z  |d d d…   z
  }|d d d…   j	                  «       t        j
                  dt        | «      dz   «      z  }|d d d…   }|||d d d…   fS )Nr*   r   )r   r[   ÚastypeÚfloatr\   r,   r]   )r`   re   rU   Úmeanconts       r   Úmeanexcess_emprn   ’   s¬   € Ü—‘˜“×%Ñ%¤eÓ,€JØ™T˜r˜TÑ"×)Ñ)Ó+¬r¯y©y¸¼CÀ»IÈ¹MÓ/JÑJÈZÑX\ÐZ\ÐX\ÑM]Ñ]€JØ™4˜R˜4Ñ ×'Ñ'Ó)¬R¯Y©Y°q¼#¸d»)Àa¹-Ó-HÑH€HØ™D˜b˜DÑ!€JØ�z 8©D¨b¨D¡>Ð1Ð1r   c           	      óÊ   — t        j                  |«      dk(  r| j                  |d¬«      S t        j                  |D �cg c]  }| j                  |d¬«      ‘Œ c}«      S c c}w )Nr   T)ÚlbÚconditional)r   ÚndimÚexpectÚarray)r   rp   ÚargsÚkwdsÚlbbs        r   Úmeanexcess_distrx   š   s`   € ô 
‡w�wˆrƒ{�aÒØ�{‰{˜b¨dˆ{Ó3Ð3ä�x‰xØ "ö$Øð Ÿ™¨¸˜Õ>ò $ó %ð 	%ùò $s   ¿A r   r=   )rp   i † )ri   )Néd   rV   rY   r   )-Ú__doc__Únumpyr   Úscipyr   Úscipy.specialr   Úscipy.stats.distributionsr   Úmatplotlib.pyplotÚpyplotr^   r   r   r   r   r
   Ú
genpareto2rI   ÚlocrJ   r,   ÚrvÚquantrd   r   rP   Úgetattrr8   ÚentropyrK   rN   rS   rU   rf   Úrvsr`   Útmprn   rx   ÚdsÚmeÚmcÚnormÚmeanr7   r   r   ú<module>rŽ      sÄ  ðñó Ý Ý Ý 3Ý ç Ý ô3�]ô 3ñB ˜c¨Ø%;Ø#&ô
€
ð Ñ €€sˆEØ€R‡Y�Yˆqƒ\€Ú#€Ø˜"�+Ø˜"�+Ø˜�*Ø˜%�.Ø˜%�.ð	"ò C�I€FˆAñ
 
Ð
%‰'�*˜fÓ
% a¨°°UÓ
;Ô<Ù	Ð
*‰'�%—/‘/ 6Ó
*¨1¨u¨f°c¸5Ó
AÕBðCñ €j×Ñ�u˜c 5°&ÐÓ9Ô :Ù €e‡o�o×Ñ˜U˜F C¨¸ÐÓ?Ô @Ù €j×Ñ˜  UÓ+Ô ,Ù €e‡o�o×Ñ˜u˜f c¨5Ó1Ô 2ò)ò,ò$ò
2ó"ñ* �j��D˜"ÓÔ Ù �j��B˜ÓÔ à‡~�~�e 1¨4€~Ó0€ñ �d :°AÔ6€Ù €cˆ!�fˆSˆT€lÔ Ù €cˆ!�fˆSˆT€lÔ ò
2ò%ñ ˜B  §¡¨1¨bÓ!1Ñ1Ó2�
€€BˆÙ €b„	Ù €b„	Ù €b„	á �o�e—j‘j SÔ)Ô *Ù �o�e—j‘j r§v¡v g¨t°Q¸Ð%<Ô=Ô >Ø‡j�j‡n�n˜&€nÓ!€Ø	ˆC�H‰H‹JÑ€Ù €c‡h�hƒj�#�c˜D‘j‘/×&Ñ&Ó(¨#¨c°A©g©,×*;Ñ*;Ó*=¸sÀ3ÈÁ9¹~×?RÑ?RÓ?TÔ UðIr   