mò
GïEc           @   s\   d  k  l Z l Z l Z d „  Z h  Z d „  Z d „  Z d „  Z d d „ Z	 d „  Z
 d S(	   (   s   ceils   sqrts   pic         C   s    t  d „  t d |  d ƒ d ƒ S(   Nc         C   s   |  | S(   N(   t   xt   y(   R    R   (    (    t4   /home/.inka/novalis_dt/novalis.org/cgi/samplesize.pyt   <lambda>   s    i   (   t   reducet   ranget   n(   R   (    (    R   t   fact   s    c         C   s˜   |  d j o d Sn |  t j o t |  Sn d } xS t d |  ƒ D]B } | t t | ƒ t |  d | ƒ ƒ | d d | d 7} qD W| t |  <| S(   Ni    i   i   (   t   kt   inverse_erf_c_tablet   totalR   t   mt   floatt   inverse_erf_c(   R   R   R
   (    (    R   R      s     @
c         C   s   d } d } d } xf t | | ƒ d j oN | } | t | ƒ d | d t t ƒ |  d d | d 7} | d 7} q W| S(   Ni    iÿÿÿÿf1.0000000000000001e-05i   i   (   R   t   lastt   curt   absR   t   sqrtt   piR    (   R    R   R   R   (    (    R   t   inverse_erf   s     :c         C   s   t  d ƒ t d |  d ƒ S(   Ni   i   (   R   R   R    (   R    (    (    R   t   probit_inverse   s    f0.94999999999999996c         C   sš   d | } t d | d ƒ } d |  j  o
 d j  n p t ‚ | | d t |  ƒ d d } | d j o | d | d | :} n t	 t
 | ƒ ƒ S(   Ni   i   i    i2   f0.25f100.0(   t   confidence_levelt   alphaR   t   zt   confidence_intervalt   AssertionErrorR   t   sizet
   populationt   intt   ceil(   R   R   R   R   R   R   (    (    R   t   sample_size"   s    
% c         C   sœ   d } |  d j p t ‚ d | j o
 d j n p t ‚ | o t | |  ƒ | d } n d } t | ƒ d } t | | d |  | ƒ d } | S(   Nf3.8416000000000001i    id   i   i   (
   t   zcR   R   t
   percentageR   R   t   pft   pcR   R   (   R   R   R    R"   R   R!   R   (    (    R   R   -   s    % N(   t   mathR   R   R   R   R	   R   R   R   R   R   (
   R   R	   R   R   R   R   R   R   R   R   (    (    R   t   ?   s   			
	