
    M/PhE                        d Z ddlZddlZddlmZ ddl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 ej                            ej                            e                    Zd	Zej                            ed
e          Z ej        e          ZdZ ej        ee                                          Zdej         fdej!        fdej"        fdej#        fdej$        fgZ% G d d          Z&dS )zJ
Created on Feb 3, 2022 1:04:22 PM

Author: Josef Perktold
License: BSD-3
    N)assert_allclose)OLS)Probit)TreatmentEffect   )results_teffectszcataneo2.csvresultsz3mbsmoke_ ~ mmarried_ + mage + mage2 + fbaby_ + meduraipwaipwaipw_wlsipw_rac                   r    e Zd Zed             Zd Zej                            de	          d             Z
dS )TestTEffectsc                     d}t          j        |t                    }t          j        t          d                   }t          ||t                    | _        d S )Nz0bweight ~ prenatal1_ + mmarried_ + mage + fbaby_mbsmoke_)results_select)r   from_formuladta_catnpasarrayr   
res_probitteff)clsformula_outcomemodtinds       i/var/www/html/test/jupyter/venv/lib/python3.11/site-packages/statsmodels/treatment/tests/test_teffects.pysetup_classzTestTEffects.setup_class,   sG    L88z'*-.."3ZHHH    c                 R    t                                           }|j        dk    sJ d S )N)i"  )r   predictshape)selfprobs     r   test_auxzTestTEffects.test_aux3   s,    !!##zW$$$$$$r    casec                  	   |\  }}| j         } t          ||          d          }t          |d d         |j        d ddf         d            t          ||          d          }t          ||j        d           |j        }t          |j        |j        d           t          |j        d d         |j        d ddf         d           t          |j        d d         |j        d dd	f         d
           t          |j	        d d         |j        d ddf         d
           t          |j
        d d         |j        d ddf         dd           |                                }t          |d ddf         |j        d ddf         d           t          |d dd	f         |j        d ddf         d           t          |j                  }|dk    rg d}	n*|dk    rg d}	n|dk    rg d}	nt          j        |          }	t          |j        |j        |	df         d           t          |j        |j        |	d	f         d           |                    d          s|j        }
 t          ||          dd	          }t          |d d         |
d ddf         d            t          ||          dd	          }t          ||j        d           |j        }t          |j        |j        d           t          |j        d d         |
d ddf         d           t          |j        d d         |
d dd	f         d
           t          |j	        d d         |
d ddf         d
           t          |j
        d d         |
d ddf         dd           |                                }t          |d ddf         |
d ddf         d           t          |d dd	f         |
d ddf         d            t          ||          dd          } t          ||          dd          }t          ||j        d           t          |j        |j        j        d           d S d S )NF)return_results   r   g-C6?)rtolTgh㈵>r   gMbP?   gV瞯<)r+   atol   gMb@?      )r   r      r*   r,   r.   r/         )r   r   r2   r*   r,   r.   r/      r1   r0   	   
                     r8   )r   r   r2   r*   r,   r.   r/   r4   r1   r0   r5   r6   g?r   )r)   effect_groupgh㈵>g-C6
?g-q=)r   getattrr   tableeffectresults_gmmstart_paramsparamsbsetvaluespvaluesconf_intlenr   arange
startswithtable_t)r$   r'   methres2r   res1res0cik_pidxr?   s              r   test_effectszTestTEffects.test_effects7   s   
dy"wtT""%888RaR$*RaRU"3$???? #wtT""$777dk5555)4;TBBBBBQBBQBE):FFFF!dj!Q&7dCCCCRaR($*RaRU*;$GGGGRaR($*RaRU*;!	/ 	/ 	/ 	/]]__2A2q5	4:bqb!e#44@@@@2A2q5	4:bqb!e#44@@@@ $+!88***CCBYY  CCBYY888CC)C..C 	TZQ%7dCCCC$*S!V"44@@@@ v&& 	(LE&74&&e!LLLDD!HeBQBEl>>>>&74&&dKKKDD$+D9999#DD-t{FFFFDKOU2A2q5\EEEEDHRaRL%A,TBBBBDL!,eBQBElFFFFDL!,eBQBEl!%E3 3 3 3BBrr1uIuRaRU|$????Brr1uIuRaRU|$???? '74&&e!LLLD&74&&dKKKDD$+E::::D-t/?/F!&( ( ( ( ( (3	( 	(r    N)__name__
__module____qualname__classmethodr   r&   pytestmarkparametrizemethodsrS    r    r   r   r   *   sr        I I [I% % % [VW--D( D( .-D( D( D(r    r   )'__doc__osnumpyr   numpy.testingr   pandaspdrX   #statsmodels.regression.linear_modelr   #statsmodels.discrete.discrete_modelr   'statsmodels.treatment.treatment_effectsr   r	   r   res_stpathabspathdirname__file__cur_dir	file_namejoin	file_pathread_csvr   formular   fitr   
results_raresults_ipwresults_aipwresults_aipw_wlsresults_ipwrar[   r   r\   r    r   <module>rw      s    
			     ) ) ) ) ) )      3 3 3 3 3 3 6 6 6 6 6 6      0 / / / / / '//"'//(33
4
4	GLL)Y77	
"+i
 
 
? V '226688
 
6
FV !()v#$R( R( R( R( R( R( R( R( R( R(r    