ó áp7]c@sèdZddlmZddlZddd„ZeZedkräej j ddd ƒZ ee ƒZ e e ƒgeee ƒƒD]Ze ejƒ^q„Ze e dd…d fed ƒee d d gƒZe eƒndS(s5 Created on Wed May 04 06:09:18 2011 @author: josef iÿÿÿÿ(tprint_functionNgš™™™™™©?cCsd}tj|ƒ}|j|}tj|d|ƒ}|dkrL|}n |tj||ƒjtƒ}tj||ddƒ}g}x|t t |ƒdƒD]d} || } || j ƒ} || j ƒ} | tj || ƒ} |j| || | | fƒq¡Wtj|ƒ}d}|dd…dd…f||dd…d d…ftjd dggƒ}tj||fƒS( sÒemprirical mean residual life or expected shortfall Parameters ---------- x : 1-dimensional array-like frac : list[float], optional All entries must be between 0 and 1 alpha : float, default 0.05 FIXME: not actually used. TODO: check formula for std of mean doesn't include case for all observations last observations std is zero vectorize loop using cumsum frac doesn't work yet itaxistsidetrightig\Âõ(\ÿ?Niiÿÿÿÿ(tnptasarraytshapetsorttNonetfloortastypetintt searchsortedtrangetlentmeantstdtsqrttappendtarrayt column_stack(txtfractalphaRtnobstxsortedt xthresholdt xlargerindextresulttitk_indtrmeantrstdtrmstdtrestcrittconfint((sHlib/python2.7/site-packages/statsmodels/sandbox/distributions/try_pot.pytmean_residual_life s&     !Pt__main__itsizei iiRgà?(t__doc__t __future__RtnumpyRRR%texpected_shortfallt__name__trandomt standard_ttrvsR"tprintR RRRRtres_frac(((sHlib/python2.7/site-packages/statsmodels/sandbox/distributions/try_pot.pyts -   /"