ó
q™xTc           @   sZ   d  d l  m Z d  d l m Z d „  Z d „  Z d „  Z d „  Z d „  Z d „  Z	 d	 S(
   iÿÿÿÿ(   t   division(   t   flexc   
      C   s¥   | d k s t  ‚ | d k s$ t  ‚ |  j ƒ  } t j |  ƒ } t | d ƒ | } t | ƒ } | t | ƒ } |  | | d } |  | | } | | | | }	 |	 S(   Ng        g      ð?i   (   t   AssertionErrort   sizeR   t   sort_permutationt   floatt   int(
   t   xt   percentt   nt   ordert   npt   irt   frt   tmp1t   tmp2t   result(    (    sS   /public2/ftp/pub/ccp4/6.5/ccp4-6.5/lib/py2/scitbx/python_utils/robust_statistics.pyt
   percentile   s    c         C   s   t  |  d ƒ } | S(   Ng      à?(   R   (   R   R   (    (    sS   /public2/ftp/pub/ccp4/6.5/ccp4-6.5/lib/py2/scitbx/python_utils/robust_statistics.pyt   median   s    c         C   s(   t  |  d ƒ } t  |  d ƒ } | | f S(   Ng      Ð?g      è?(   R   (   R   t   h1t   h2(    (    sS   /public2/ftp/pub/ccp4/6.5/ccp4-6.5/lib/py2/scitbx/python_utils/robust_statistics.pyt   hinges   s    c         C   s   t  |  ƒ \ } } | | S(   N(   R   (   R   t   lowt   high(    (    sS   /public2/ftp/pub/ccp4/6.5/ccp4-6.5/lib/py2/scitbx/python_utils/robust_statistics.pyt   h_spread   s    c         C   s8   t  |  ƒ \ } } t |  ƒ } | | d | d } | S(   Ng       @g      @(   R   R   (   R   R   R   t   mR   (    (    sS   /public2/ftp/pub/ccp4/6.5/ccp4-6.5/lib/py2/scitbx/python_utils/robust_statistics.pyt   trimean   s    c         C   sm   t  |  | d ƒ } t  |  d | d ƒ } |  j |  | k j ƒ  ƒ } | j | | k  j ƒ  ƒ } t j | ƒ S(   Ng       @g      ð?(   R   t   selectt
   iselectionR   t   mean(   R   R   R   R   t   tmp(    (    sS   /public2/ftp/pub/ccp4/6.5/ccp4-6.5/lib/py2/scitbx/python_utils/robust_statistics.pyt   trimmed_mean%   s
    N(
   t
   __future__R    t   scitbx.array_familyR   R   R   R   R   R   R   (    (    (    sS   /public2/ftp/pub/ccp4/6.5/ccp4-6.5/lib/py2/scitbx/python_utils/robust_statistics.pyt   <module>   s   					