>hSrSSKrSSKJr SSKJr SSKJrJ r J r J r J r J r "SS5r"SS \5r"S S \5rg) zM Created on Wed Feb 17 15:35:23 2021 Author: Josef Perktold License: BSD-3 N)stats)cache_readonly)_Grid cdf2prob_grid prob2cdf_grid_eval_bernstein_dd_eval_bernstein_2d_eval_bernstein_1dcV\rSrSrSrSr\S5r\S5r Sr Sr Sr S r S rg ) BernsteinDistributionaZDistribution based on Bernstein Polynomials on unit hypercube. Parameters ---------- cdf_grid : array_like cdf values on a equal spaced grid of the unit hypercube [0, 1]^d. The dimension of the arrays define how many random variables are included in the multivariate distribution. Attributes ---------- cdf_grid : grid of cdf values prob_grid : grid of cell or bin probabilities k_dim : (int) number of components, dimension of random variable k_grid : (tuple) shape of cdf_grid k_grid_product : (int) total number of bins in grid _grid : Grid instance with helper methods and attributes c0[R"U5=UlnURUlUR Ul[R"UR Vs/sHo"S- PM sn5Ul[UR 5Ul gs snf)N) npasarraycdf_gridndimk_dimshapek_gridprodk_grid_productr_grid)selfris vC:\Users\julio\OneDrive\Documentos\Trabajo\Ideas Frescas\venv\Lib\site-packages\statsmodels/distributions/bernstein.py__init__BernsteinDistribution.__init__&si#%::h#77 ]] nn  ggDKK&@KqsK&@A4;;' 'AsBc |[R"U5n[R"US:5(d[R"US:5(a [S5eURS:Xa USS2S4nUR Sn[R "U5S:XaU/U-nUVs/sH!n[R"SU- SUS-5PM# nn[R"XSS9upg[UVs/sH oSS:HPM sn5(deU[U5-n[U5n U"U 5$s snfs snf) aCreate distribution instance from data using histogram binning. Classmethod to construct a distribution instance. Parameters ---------- data : array_like Data with observation in rows and random variables in columns. Data can be 1-dimensional in the univariate case. k_bins : int or list Number or edges of bins to be used in numpy histogramdd. If k_bins is a scalar int, then the number of bins of each component will be equal to it. Returns ------- Instance of a Bernstein distribution rrzdata needs to be in [0, 1]NF)binsdensity) rrany ValueErrorrrsizelinspace histogramddalllenr) clsdatak_binsrnir"ceeirs r from_dataBernsteinDistribution.from_data-s(zz$ 66$(  rvvdQh//9: : 99>4=D 1  776?a X%F:@A&B BGQQ/&A~~du=+2qEQJ+,,,, SY #8}B,s *(D44D9c*[URSS9$)N)prepend)rr)rs r prob_gridBernsteinDistribution.prob_gridUsT]]D99c[R"U5nURS:XaURS:Xa USS2S4n[ XR 5nU$)ancdf values evaluated at x. Parameters ---------- x : array_like Points of multivariate random variable at which cdf is evaluated. This can be a single point with length equal to the dimension of the random variable, or two dimensional with points (observations) in rows and random variables in columns. In the univariate case, a 1-dimensional x will be interpreted as different points for evaluation. Returns ------- pdf values Notes ----- Warning: 2-dim x with many points can be memory intensive because currently the bernstein polynomials will be evaluated in a fully vectorized computation. rN)rrrrrrrxcdf_s rcdfBernsteinDistribution.cdfYsF. JJqM 66Q;4::?!T' A!!]]3 r8c[R"U5nURS:XaURS:Xa USS2S4nUR[ XR 5-nU$)anpdf values evaluated at x. Parameters ---------- x : array_like Points of multivariate random variable at which pdf is evaluated. This can be a single point with length equal to the dimension of the random variable, or two dimensional with points (observations) in rows and random variables in columns. In the univariate case, a 1-dimensional x will be interpreted as different points for evaluation. Returns ------- cdf values Notes ----- Warning: 2-dim x with many points can be memory intensive because currently the bernstein polynomials will be evaluated in a fully vectorized computation. rN)rrrrrrr6rr;pdf_s rpdfBernsteinDistribution.pdfvsR. JJqM 66Q;4::?!T' A""%7>>%JJ r8cURS:XaU$S/UR-n[R"U5S:XaU/nUHn[SSS5X#'M UR[ U5n[ U5nU$)zGet marginal BernsteinDistribution. Parameters ---------- idx : int or list of int Index or indices of the component for which the marginal distribution is returned. Returns ------- BernsteinDistribution instance for the marginal distribution. rr N)rrrslicertupler )ridxsliicdf_m bpd_marginals r get_marginal"BernsteinDistribution.get_marginalsw ::?KTDJJ  88C=B %CB4t,BF eBi(,U3 r8c [RRXRR 55nUR n/n[ [U55HnX%S:wdM [R"XPRR5n/n[ U5HnnURUn URRUXhn UR[RR!X-S-U SU - -S-X%S95 Mp UR""U55 M [R$"U5n U $)z|Generate random numbers from distribution. Parameters ---------- nobs : int Number of random observations to generate. rr)r&)rrandom multinomialr6flattenrranger* unravel_indexrrr x_marginalappendrbetarvs column_stack concatenate) rnobsrvs_mnlk_comprvs_mrrHrvsijnxgirvsms rrXBernsteinDistribution.rvss))''nn.D.D.FGs7|$AzQ&&q..*>*>?vA AA**//236:CKK qw{A3K!O4;J!/!@A ' R__T23%~~e$ r8)rrrrrN)__name__ __module__ __qualname____firstlineno____doc__r classmethodr2rr6r=rBrMrX__static_attributes__rEr8rr r sI&(%%N:::<6r8r c \rSrSrSrSrSrg)BernsteinDistributionBVc0[XR5nU$N)r rr:s rr=BernsteinDistributionBV.cdfs!!]]3 r8cJUR[XR5-nU$rp)rr r6r@s rrBBernsteinDistributionBV.pdfs!""%7>>%JJ r8rENrerfrgrhr=rBrkrEr8rrmrms r8rmc(\rSrSrSSjrSSjrSrg)BernsteinDistributionUVc.[XRUS9nU$N)method)r r)rr;rzr<s rr=BernsteinDistributionUV.cdfs!!]]6B r8cHUR[XRUS9-nU$ry)rr r6)rr;rzrAs rrBBernsteinDistributionUV.pdfs*""%7>>?E&GG r8rEN)binomrtrEr8rrvrvs  r8rv)rinumpyrscipyrstatsmodels.tools.decoratorsrstatsmodels.distributions.toolsrrrrr r r rmrvrEr8rrsH7DD vvr 3  3 r8