EiSSKrSSKJr SSKJr /SQr/SQr\"\SSS25r\"\SSS25rSr /SQr /S Qr \"\ SSS25r \"\ SSS25r S rg) N)poly1d)beta)gSbQ @go|-a?g_3 L/g|A"?gCUG)?g<*x @gy@g`B{dA?g-~?ct[R"U5nURn[R"U5R [R 5nSnUS:nUS:Hn[R "U5U:nX4-U-)nXnXn[R"U5n [RX'[RX'URS:a[U5[U5- X'URS:a)SUS-- SSSU--- [US-US-5- -X'U RU5n U $)aVariance of the Tukey Lambda distribution. Parameters ---------- lam : array_like The lambda values at which to compute the variance. Returns ------- v : ndarray The variance. For lam < -0.5, the variance is not defined, so np.nan is returned. For lam = 0.5, np.inf is returned. Notes ----- In an interval around lambda=0, this function uses the [4,4] Pade approximation to compute the variance. Otherwise it uses the standard formula (https://en.wikipedia.org/wiki/Tukey_lambda_distribution). The Pade approximation is used because the standard formula has a removable discontinuity at lambda = 0, and does not produce accurate numerical results near lambda = 0. g333333?grg@r)npasarrayshape atleast_1dastypefloat64abs empty_likenaninfsize_tukeylambda_var_p_tukeylambda_var_qrreshape) lamshp thresholdlow_mask neghalf_mask small_maskreg_masksmallregvs [/home/mande/repo/quber/.venv/lib/python3.13/site-packages/scipy/stats/_tukeylambda_stats.pytukeylambda_variancer#+s$. **S/C ))C --  # #BJJ /CI TzH$;Ly(J(:56H OE -C cA&&AKffAO zzA~*514Fu4MM  xx!|S!V|sQW}(=(,S1WcAg(>)?@  #A H)g333333?g6|igeSH6gѐ환^?g˝)kPd@)rg?ݻA@gID@)@gPr?g`2fQc[R"U5nURn[R"U5R [R 5nSnUS:nUS:Hn[R "U5U:nX4-U-)nXnXn[R"U5n [RX'[RX'URS:a[U5[U5- X'URS:amSSU-S-- S[SU-S-US-5-- S[SU-S-SU-S-5--n SSSU-S-- [US-US-5- S--n X- S- X'U RU5n U $) a*Kurtosis of the Tukey Lambda distribution. Parameters ---------- lam : array_like The lambda values at which to compute the variance. Returns ------- v : ndarray The variance. For lam < -0.25, the variance is not defined, so np.nan is returned. For lam = 0.25, np.inf is returned. g)\(?gпrrr r)r r r r rrrrrrr_tukeylambda_kurt_p_tukeylambda_kurt_qrr) rrrr negqrtr_maskrrrr knumerdenoms r"tukeylambda_kurtosisr.s| **S/C ))C --  # #BJJ /CIU{H%r6s8-? /"56/"569 z?A 1$B$781$B$784 r$