ó
áp7]c           @   s\  d  Z  d d l Z d d l m Z d d l Z d d l m Z d d d „ Z	 e
 d k rXd d l j Z e j j j d e ƒ Z e j e j d	 e ƒe _ e	 e j ƒ Z e	 e j d d … d d
 … f ƒ Z e j j d ƒ rXe j d ƒ Z e	 e j ƒ Z e	 e j e j e j f ƒ ƒ Z e	 e j g  e j  j! D] Z" e e" ^ q6ƒ ƒ Z# qXn  d S(   s,   
Glue for returning descriptive statistics.
iÿÿÿÿN(   t   stats(   t	   sign_testi    c         C   s"  t  j |  ƒ } | d k rY t | t  j ƒ rY | j d k rY | d d … d f } qY n  | j d d k rCd i | d 6d d 6t | ƒ d 6t j	 | ƒ d d d	 6t j	 | ƒ d d d
 6| j
 ƒ  d 6t  j | ƒ d 6d t | j ƒ  ƒ d t | j ƒ  ƒ d d 6| j ƒ  d 6| j ƒ  d 6t j | ƒ d 6t j | ƒ d 6t j | ƒ d 6t  j | d d d ƒd 6t  j | | j
 ƒ  d d d ƒd 6t  j | ƒ d 6} | d t g  d) D] } t j | | ƒ ^ qµƒ 7} t j | d ƒ \ } } t | ƒ \ } }	 t j t  j | ƒ ƒ \ }
 } | d$ | | | |	 |
 | f 7} nÛ | j d d k rd% t j } x¸ t | j d ƒ D]” } | d d … | f } i | d 6t | ƒ d& 6| j
 ƒ  d 6| j ƒ  d 6d t | j ƒ  ƒ d t | j ƒ  ƒ d t j d 6} | d' | 7} qwWn t d( ƒ ‚ | S(*   s¿  
    Prints descriptive statistics for one or multiple variables.

    Parameters
    ----------
    data: numpy array
        `x` is the data

    v: list, optional
        A list of the column number or field names (for a recarray) of variables.
        Default is all columns.

    axis: 1 or 0
        axis order of data.  Default is 0 for column-ordered data.

    Examples
    --------
    >>> descstats(data.exog,v=['x_1','x_2','x_3'])
    i   NsÅ  
    ---------------------------------------------
    Univariate Descriptive Statistics
    ---------------------------------------------

    Var. Name   %(name)12s
    ----------
    Obs.          %(nobs)22i  Range                  %(range)22s
    Sum of Wts.   %(sum)22s  Coeff. of Variation     %(coeffvar)22.4g
    Mode          %(mode)22.4g  Skewness                %(skewness)22.4g
    Repeats       %(nmode)22i  Kurtosis                %(kurtosis)22.4g
    Mean          %(mean)22.4g  Uncorrected SS          %(uss)22.4g
    Median        %(median)22.4g  Corrected SS            %(ss)22.4g
    Variance      %(variance)22.4g  Sum Observations        %(sobs)22.4g
    Std. Dev.     %(stddev)22.4g
    t   names   N/At   sumt   nobsi    t   modet   nmodet   meant   mediant   (s   , t   )t   ranget   variancet   stddevt   coeffvart   skewnesst   kurtosisi   t   axist   usst   sst   sobss  

    Percentiles
    -------------
    1  %%          %12.4g
    5  %%          %12.4g
    10 %%          %12.4g
    25 %%          %12.4g

    50 %%          %12.4g

    75 %%          %12.4g
    90 %%          %12.4g
    95 %%          %12.4g
    99 %%          %12.4g
    i   i
   i   i2   iK   iZ   i_   ic   sc  

    Tests of Location (H0: Mu0=0)
    -----------------------------
    Test                Statistic       Two-tailed probability
    -----------------+-----------------------------------------
    Student's t      |  t %7.5f   Pr > |t|   <%.4f
    Sign             |  M %8.2f   Pr >= |M|  <%.4f
    Signed Rank      |  S %8.2f   Pr >= |S|  <%.4f

    sŽ   
    Var. Name   |     Obs.        Mean    Std. Dev.           Range
    ------------+--------------------------------------------------------t   obss;   %(name)15s %(obs)9i %(mean)12.4g %(stddev)12.4g %(range)20ss   data not understood(	   i   i   i
   i   i2   iK   iZ   i_   ic   (   t   npt   arrayt   Nonet
   isinstancet   recarrayt   ndimt   shapet   lenR    R   R   R   t   strt   mint   maxt   vart   stdt	   variationt   skewR   R   t   tuplet   scoreatpercentilet   ttest_1sampR   t   wilcoxont   squeezet   ost   linesepR   t
   ValueError(   t   datat   colsR   t   xt   desct   pert   tt   p_tt   Mt   p_Mt   St   p_SR!   t   xvt   kwargs(    (    s<   lib/python2.7/site-packages/statsmodels/sandbox/descstats.pyt	   descstats   sH    "!0/ *8	'
 

9t   __main__t	   as_pandast   prependi   s   ./Econ724_PS_I_Data.csv($   t   __doc__t   numpyR   t   scipyR    R*   t"   statsmodels.stats.descriptivestatsR   R   R:   t   __name__t   statsmodels.apit   apit   smt   datasetst   longleyt   loadt   FalseR-   t   add_constantt   exogt   sum1t   sum1at   patht   isfilet
   recfromcsvt   data2t   ahet   sum2t   column_stackt   yrseduct   sum3t   dtypet   namest   _t   sum4(    (    (    s<   lib/python2.7/site-packages/statsmodels/sandbox/descstats.pyt   <module>   s"   	•%!