U
    Kvf                     @   s   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 d dlmZ d dlmZ G dd	 d	ZG d
d deZG dd deZG dd deZG dd deZdS )    N)GLM)families)links)GEE)Independence)assert_allclosec                   @   s   e Zd Zdd Zdd ZdS )CheckGEEGLMc                 C   s(   | j }| j}t|jj|jjddd d S Ngư>g|=)ZrtolZatol)result1result2r   paramsvaluesselfZres1Zres2 r   I/tmp/pip-unpacked-wheel-2v6byqio/statsmodels/genmod/tests/test_gee_glm.py
test_basic   s     zCheckGEEGLM.test_basicc                 C   st   | j }| j}t|j|jddd t|j|jddd t|j|jddd t|j|jddd t|j|jddd d S r	   )r
   r   r   Zresid_responseZresid_pearsonZresid_devianceZresid_anscombeZresid_workingr   r   r   r   
test_resid   s,    
 
 
 
 
 zCheckGEEGLM.test_residN)__name__
__module____qualname__r   r   r   r   r   r   r      s   r   c                   @   s   e Zd Zedd ZdS )TestCompareLogitc                 C   s   t  }t }tjd dtjjdddk  }tjjdd}tjjdd}tjjdd}tjjdddd}t	||||d}t
jd||||d	}	|	 | _tjd||d
}
|
jdd| _d S N    d   sizer      YX1X2X3Y ~ X1 + X2 + X3familyZ
cov_structdatar&   FZdisp)r   r   ZBinomialnprandomseednormalrandintpd	DataFramer   from_formulafitr
   r   r   clsvsr&   r    r!   r"   r#   groupsDZmod1Zmod2r   r   r   setup_class(   s     
 
zTestCompareLogit.setup_classNr   r   r   classmethodr8   r   r   r   r   r   &   s   r   c                   @   s   e Zd Zedd ZdS )TestComparePoissonc                 C   s   t  }t }tjd tdtjjdd }tjjdd}tjjdd}tjjdd}tjjdddd}t	
||||d}tjd||||d	}	|	 | _tjd||d
}
|
jdd| _d S r   )r   r   ZPoissonr*   r+   r,   expr-   r.   r/   r0   r   r1   r2   r
   r   r   r3   r   r   r   r8   ?   s     
 
zTestComparePoisson.setup_classNr9   r   r   r   r   r;   =   s   r;   c                   @   s   e Zd Zedd ZdS )TestCompareGaussianc           
      C   s   t  }t }tjd tjjdd}tjjdd}tjjdd}tjjdd}ttdt	d}t
||||d}tjd||||d}	|	 | _tjd|d	 | _d S )
Nr   r   r         r   r$   r%   )r(   )r   r   ZGaussianr*   r+   r,   r-   ZkronZarangeZonesr/   r0   r   r1   r2   r
   r   r   )
r4   r5   r&   r    r!   r"   r#   r6   r7   Zmdr   r   r   r8   V   s    
 
zTestCompareGaussian.setup_classNr9   r   r   r   r   r=   T   s   r=   c                   @   s   e Zd Zedd ZdS )TestCompareGammac                 C   s   t  }tjt d}tjd tdtjj	dd }tjj	dd}tjj	dd}tjj	dd}tjj
dddd}t||||d}tjd	||||d
}	|	 | _tjd	||d}
|
jdd| _d S )N)linkr   g?r   r   r   r   r   r$   r%   r'   Fr)   )r   r   Gammar   Logr*   r+   r,   r<   r-   r.   r/   r0   r   r1   r2   r
   r   r   r3   r   r   r   r8   m   s     
 
zTestCompareGamma.setup_classNr9   r   r   r   r   r@   k   s   r@   )Znumpyr*   Zpandasr/   Z+statsmodels.genmod.generalized_linear_modelr   Zstatsmodels.genmodr   Zstatsmodels.genmod.familiesr   Z3statsmodels.genmod.generalized_estimating_equationsr   Zstatsmodels.genmod.cov_structr   Znumpy.testingr   r   r   r;   r=   r@   r   r   r   r   <module>   s   