B ZB @sdZddlmZddlmZmZddlZddZddd Z d d Z d d Z e dkrddl mZddZeeddgddgeddeeddgddgeddeddddgddddgddddgddddgddddggZe edddedddeZeeededgdgdgdggZe edddedddeZeeededdddggZe edddedddeZeeeddS)z?Quantizing a continuous distribution in 2d Author: josef-pktd )print_function)rangelmapNcCsD||}||d|d}||d|d}||}||||S)ahelper function for probability of a rectangle in a bivariate distribution Parameters ---------- lower : array_like tuple of lower integration bounds upper : array_like tuple of upper integration bounds cdf : callable cdf(x,y), cumulative distribution function of bivariate distribution how does this generalize to more than 2 variates ? r)loweruppercdfZprobuuZprobulZprobluZprobllrrIlib/python3.7/site-packages/statsmodels/sandbox/distributions/quantize.pyprob_bv_rectangle s r c Cst|tjszttj|}t|}g}tttj|t|krxHt |D],}dg|}t d||<| |||qHWn|j d}|}t t|||}|}xt |D]}tj||d}qW|S)zhelper function for probability of a rectangle grid in a multivariate distribution how does this generalize to more than 2 variates ? bins : tuple tuple of bin edges, currently it is assumed that they broadcast correctly Nr)axis) isinstancenpZndarrayrasarraylenallndimonesrsliceappendshapeprintcopyZdiff) Zbinsr r Zn_dimZbins_dZsl cdf_valuesprobsrrr prob_mv_grids"      rc st|}t|}t|d}t|d}tjt||f}||dddf|fdd}xbtd|dD]P}xJtd|dD]8}||f} |d|df} t| | |||d|df<qWqrWt|rt |S)zquantize a continuous distribution given by a cdf Parameters ---------- binsx : array_like, 1d binedges rNcs ||fS)Nr)xy)rrr Nsz#prob_quantize_cdf..) rrrnanrrr isnananyAssertionError) binsxbinsyr nxnyrZcdf_funcxindyindrrr)rr prob_quantize_cdf?s     $r+c Cst|}t|}t|d}t|d}tjt||f}xrtd|dD]`}xZtd|dD]H}||||f}||d||df} t| ||||d|df<qdWqPWt|rt |S)zquantize a continuous distribution given by a cdf old version without precomputing cdf values Parameters ---------- binsx : array_like, 1d binedges r) rrrr!rrr r"r#r$) r%r&r r'r(rr)r*rrrrr prob_quantize_cdf_oldYs    $r,__main__)assert_almost_equalcCs||S)Nr)rrrrr r xsr rg?g?g?)r )__doc__Z __future__rZstatsmodels.compat.pythonrrZnumpyrr rr+r,__name__Z numpy.testingr.Zunif_2dZarrayZarr1bZlinspaceZarr1aZarr2bZarr2aZarr3bZarr3arrrr s8  !