ó
áp7]c           @  s[  d  Z  d d l m Z d d l Z d „  Z d „  Z e j d ƒ Z e j d ƒ Z	 e j d ƒ Z
 e j d	 ƒ Z e j d
 ƒ Z e j d ƒ Z e j d ƒ Z e j e e e e	 ƒ ƒ j e ƒ Z e e j e ƒ ƒ e e j e j e j e d d ƒ ƒ e ƒ ƒ e e j e d d e j d e e d ƒ e ƒ ƒ e j e j e d d ƒ ƒ e j e j e d ƒ ƒ d e j e e j ƒ Z e d d e j d e e d e e
 ƒ d e j e
 ƒ Z e d d e j d e d e ƒ d Z e e Z e e j e ƒ ƒ e e j e ƒ ƒ e e j e ƒ Z e e ƒ e e j e ƒ ƒ e d ƒ x1 e j e ƒ d D] \ Z Z e e d e ƒ qvWe e j e ƒ d d ƒ e d ƒ e e Z e e j e ƒ ƒ e e j e ƒ ƒ e d ƒ e e j e  d
 d d d d d d d ƒ ƒ ƒ e e j e  d
 d d d d d d d ƒ ƒ j! ƒ  ƒ d S(   s9   
Created on Sat Mar 13 07:56:22 2010

Author: josef-pktd
iÿÿÿÿ(   t   print_functionNc         C  s>   d | t  j d t  j ƒ t  j |  | d d | d ƒ S(   s=   Return the probability density function as an expression in xi   i   (   t   syt   sqrtt   pit   exp(   t   xt   mut   sigma(    (    sH   lib/python2.7/site-packages/statsmodels/sandbox/regression/sympy_diff.pyt   pdf   s    c         C  s*   d t  j |  | | t  j d ƒ ƒ d S(   s<   Return the cumulative density function as an expression in xi   i   (   R   t   erfR   (   R   R   R   (    (    sH   lib/python2.7/site-packages/statsmodels/sandbox/regression/sympy_diff.pyt   cdf   s    R   R   t   sigma2R   t   yt   dft   si   i   g       @g      à?g      ð?s0   
 derivative of loglike of t distribution wrt dfi    t   =s(   
standard t distribution, dll_df, dll_dys   
g      ø?gŽuq $@("   t   __doc__t
   __future__R    t   sympyR   R   R
   t   SymbolR   R   R   R   R   R   R   t   logt   difft   dldxnormt   printt   simplifyt   gammaR   t   tllf1t   tllf2t   tllf2stdt   tllft   dlddft   cset   kt   vt   tllfstdt   subst   dictt   evalf(    (    (    sH   lib/python2.7/site-packages/statsmodels/sandbox/regression/sympy_diff.pyt   <module>   sF   		$07R@+


  


.