0Ph#dZddlZddlZddlmZddlmZmZddlZddl m Z ddl m Z ddlmZddlmZmZdd lmZdd lmZmZmZdd lmZd Zeeed ddgeed ddgeed ddgeedddgeedddgeed ddgeed ddgddgeedd dgeedddgdgeedddddgeedddddgdgdgdd d]dddddddddddddd dZeeed ddgeed ddgeed ddgeedddgeed ddgdgdgedd hdgdgdgd! d d]d"dd#dd$dd$dd%d&Zeeed ddgdgd'dd^dd)d*Zeeed ddgeed ddgeedddgeed ddgeedddgeed dddgeedd dgeedddgdgdgdgd+ d d_d,d ddd-ddd$dd. d/Zeeeddde gdgeeddddgdgeedd dgd0d d`dddd1d2d3Z!eeed dde gdgeeddddgdgd4dd`dddd5d6Z"eeed dddgeed ddgeed ddddgeeddddge gdgdgdgd7d daddd8ddd$d9d:Z#eeed ddgeed"ddgeedddgdgd;ddbddd<d=Z$eeed ddgeedddgdgd>dd`ddd<d?Z%eeed ddgeedddgdgd>dd`ddd<d@Z&eeed ddgeed ddgeed ddgeedd dgdgdAd d_d,d-ddBdCZ'eeed ddgeed ddgeed ddgeed ddgdgdDddd)dEZ(eeed ddgeed ddgdgdFddbdd)dGZ)eeed ddgdgdHddd)dIZ*eeed ddgeedd dgdgeedd dgeedd dgehdJdgdgdKd dcdLd$dMdNdddOdPZ+eeed ddgeedddgdgdgdQdd`ddd$dRdSZ,eeed ddgeedddgdgd>dd`ddd<dTZ-eddgeedddgeed ddgeed ddgeed ddgdgdgdUddddddVdddUdWZ.dddXZ/ee geed ddgeedddgeedddgeedddgdgdgdYddd,ddddZd[Z0ee geed dddgeedddgeedddgeedddgdgdgdYddd,ddddZd\Z1dS)ez* Generate samples of synthetic data sets. N)Iterable)IntegralReal)linalg)MultiLabelBinarizer) check_arraycheck_random_state)shuffle)Interval StrOptionsvalidate_params)sample_without_replacementcr|dkr?tj|d||dz ft|d|gSt d|z||dd}tj|dd d d | d f}|S) z5Returns distinct binary samples of length dimensions.rsize random_statez>u4F)dtypecopyz>u1) N) nphstackrandint_generate_hypercuberastype unpackbitsviewreshape)samples dimensionsrngouts c/var/www/html/test/jupyter/venv/lib/python3.11/site-packages/sklearn/datasets/_samples_generator.pyrrsBy AWj2o$> ??#GR55    %Q ]G# N N N U U% V  C - ( ( 0 0 : :111zkll? KC Jleft)closedz array-likebothneitherbooleanr) n_samples n_features n_informative n_redundant n_repeated n_classesn_clusters_per_classweightsflip_y class_sep hypercubeshiftscaler rT)prefer_skip_nested_validationdg{Gz??) r0r1r2r3r4r5r6r7r8r9r:r rc  t|}||z|z|krtd|tjzkr/d}|dz }t||d|zt dz fvrtdt dz kr^t trdtz gzn`. Parameters ---------- n_samples : int, default=100 The number of samples. n_features : int, default=20 The total number of features. These comprise ``n_informative`` informative features, ``n_redundant`` redundant features, ``n_repeated`` duplicated features and ``n_features-n_informative-n_redundant-n_repeated`` useless features drawn at random. n_informative : int, default=2 The number of informative features. Each class is composed of a number of gaussian clusters each located around the vertices of a hypercube in a subspace of dimension ``n_informative``. For each cluster, informative features are drawn independently from N(0, 1) and then randomly linearly combined within each cluster in order to add covariance. The clusters are then placed on the vertices of the hypercube. n_redundant : int, default=2 The number of redundant features. These features are generated as random linear combinations of the informative features. n_repeated : int, default=0 The number of duplicated features, drawn randomly from the informative and the redundant features. n_classes : int, default=2 The number of classes (or labels) of the classification problem. n_clusters_per_class : int, default=2 The number of clusters per class. weights : array-like of shape (n_classes,) or (n_classes - 1,), default=None The proportions of samples assigned to each class. If None, then classes are balanced. Note that if ``len(weights) == n_classes - 1``, then the last class weight is automatically inferred. More than ``n_samples`` samples may be returned if the sum of ``weights`` exceeds 1. Note that the actual class proportions will not exactly match ``weights`` when ``flip_y`` isn't 0. flip_y : float, default=0.01 The fraction of samples whose class is assigned randomly. Larger values introduce noise in the labels and make the classification task harder. Note that the default setting flip_y > 0 might lead to less than ``n_classes`` in y in some cases. class_sep : float, default=1.0 The factor multiplying the hypercube size. Larger values spread out the clusters/classes and make the classification task easier. hypercube : bool, default=True If True, the clusters are put on the vertices of a hypercube. If False, the clusters are put on the vertices of a random polytope. shift : float, ndarray of shape (n_features,) or None, default=0.0 Shift features by the specified value. If None, then features are shifted by a random value drawn in [-class_sep, class_sep]. scale : float, ndarray of shape (n_features,) or None, default=1.0 Multiply features by the specified value. If None, then features are scaled by a random value drawn in [1, 100]. Note that scaling happens after shifting. shuffle : bool, default=True Shuffle the samples and the features. random_state : int, RandomState instance or None, default=None Determines random number generation for dataset creation. Pass an int for reproducible output across multiple function calls. See :term:`Glossary `. Returns ------- X : ndarray of shape (n_samples, n_features) The generated samples. y : ndarray of shape (n_samples,) The integer labels for class membership of each sample. See Also -------- make_blobs : Simplified variant. make_multilabel_classification : Unrelated generator for multilabel tasks. Notes ----- The algorithm is adapted from Guyon [1] and was designed to generate the "Madelon" dataset. References ---------- .. [1] I. Guyon, "Design of experiments for the NIPS 2003 variable selection benchmark", 2003. Examples -------- >>> from sklearn.datasets import make_classification >>> X, y = make_classification(random_state=42) >>> X.shape (100, 20) >>> y.shape (100,) >>> list(y[:5]) [np.int64(0), np.int64(0), np.int64(1), np.int64(1), np.int64(0)] ziNumber of informative, redundant and repeated features must sum to less than the number of total featuresz0n_classes({}) * n_clusters_per_class({}) must bez) smaller or equal 2**n_informative({})={}rNr(z:Weights specified but incompatible with number of classes.r>rcLg|] }t|zzz !S)int).0kr3r4r.r5s r& z'make_classification..sE  II . .1E EFFr'rFrrr.?r?r<r)r ValueErrorrlog2formatlen isinstancelistsumresizerangezerosrCrrfloatuniformstandard_normal enumeratedotintpr util_shufflearanger )"r.r/r0r1r2r3r4r5r6r7r8r9r:r r generatormsg n_useless n_clustersn_samples_per_clusteriXy centroidsstoprEcentroidstartX_kABnindices flip_masks"` ``` r&make_classificationrn(sJ#<00I{"Z/*<<    rwy+??@@@@@ :: JJy"6 q-GW X X    w<< 9q=9 9 9L  w<<9q= ( ('4(( 6!S3w<<%7$88)GY77!C $5$55 ?#i/]*[8:EI11Jz"" 9s#8999 : :33a*n---2---- )Z())A #&&&A$J yIIPP EQIYI I @Y&&ZO&<<< Y&&Q ,>&??? %449m:T4UUAaaa- D ++ 8D#8#;;tI %* d N]N*+ !! }'E!FF F J6#q>>C xQ !! {'C!DD Dq H<>F aaa- != = !!!]][8 8 89 A~~ K 'EY..J.???#EMMbgVV#$QQQZ=!!!QZ   1}}%55Iy;Q5RR!!!iZ[[.}}%%9%55>  (((II)  }Y&&J&777!;yHJA }C)+++<<<<JA Aqy9991)J'''"""AAAwJ-!!!QQQ$ a4Kr'densesparse) r.r/r3n_labelslengthallow_unlabeledrpreturn_indicatorreturn_distributionsr2F)r3rqrrrsrprtrurct| |} | | z} tj| |ftjdzfd} t jd} t jddg} g}t |D][}| \}}| || t| ||\tj t| tj }tj || | f|f}||s|}|dvrJt!|d k }|t |g|}|r||| fS||fS) a` Generate a random multilabel classification problem. For each sample, the generative process is: - pick the number of labels: n ~ Poisson(n_labels) - n times, choose a class c: c ~ Multinomial(theta) - pick the document length: k ~ Poisson(length) - k times, choose a word: w ~ Multinomial(theta_c) In the above process, rejection sampling is used to make sure that n is never zero or more than `n_classes`, and that the document length is never zero. Likewise, we reject classes which have already been chosen. For an example of usage, see :ref:`sphx_glr_auto_examples_datasets_plot_random_multilabel_dataset.py`. Read more in the :ref:`User Guide `. Parameters ---------- n_samples : int, default=100 The number of samples. n_features : int, default=20 The total number of features. n_classes : int, default=5 The number of classes of the classification problem. n_labels : int, default=2 The average number of labels per instance. More precisely, the number of labels per sample is drawn from a Poisson distribution with ``n_labels`` as its expected value, but samples are bounded (using rejection sampling) by ``n_classes``, and must be nonzero if ``allow_unlabeled`` is False. length : int, default=50 The sum of the features (number of words if documents) is drawn from a Poisson distribution with this expected value. allow_unlabeled : bool, default=True If ``True``, some instances might not belong to any class. sparse : bool, default=False If ``True``, return a sparse feature matrix. .. versionadded:: 0.17 parameter to allow *sparse* output. return_indicator : {'dense', 'sparse'} or False, default='dense' If ``'dense'`` return ``Y`` in the dense binary indicator format. If ``'sparse'`` return ``Y`` in the sparse binary indicator format. ``False`` returns a list of lists of labels. return_distributions : bool, default=False If ``True``, return the prior class probability and conditional probabilities of features given classes, from which the data was drawn. random_state : int, RandomState instance or None, default=None Determines random number generation for dataset creation. Pass an int for reproducible output across multiple function calls. See :term:`Glossary `. Returns ------- X : ndarray of shape (n_samples, n_features) The generated samples. Y : {ndarray, sparse matrix} of shape (n_samples, n_classes) The label sets. Sparse matrix should be of CSR format. p_c : ndarray of shape (n_classes,) The probability of each class being drawn. Only returned if ``return_distributions=True``. p_w_c : ndarray of shape (n_features, n_classes) The probability of each feature being drawn given each class. Only returned if ``return_distributions=True``. Examples -------- >>> from sklearn.datasets import make_multilabel_classification >>> X, y = make_multilabel_classification(n_labels=3, random_state=42) >>> X.shape (100, 20) >>> y.shape (100, 5) >>> list(y[:3]) [array([1, 1, 0, 1, 0]), array([0, 1, 1, 1, 0]), array([0, 1, 0, 0, 0])] rraxisc 2j\}}|dz}s|dks||kr#  }s|dk||k#t}t||krat j |t|z }||t||kat|}d}|dkr  }|dkt|dkr |}||fS |d d }||dz}t j| |}||fS)Nr(rrryr) shapepoissonsetrMr searchsortedrUupdaterOrtakerPcumsum)_r3y_sizerccn_wordswordscumulative_p_w_samplerscumulative_p_cr\rrr/rqp_w_cs r&sample_examplez6make_multilabel_classification..sample_examples{ 9Q" 1v{{v 7I7I&&x00F# 1v{{v 7I7I EE!ff 0A0AvPSTUPVPV0A0W0WXXA HHQKKK!ff GGll''//Gll q66Q;;%%jw%??E!8O!& 11 5 5 9 9q 9 A A H H J J!6r!:: 5y7H7Hg7H7V7VWWaxr'rarG)r|)Trprorp) sparse_output)r rUrPrrarrayrRextendappendrMonesfloat64sp csr_matrixsum_duplicatestoarrayrfit transform)r.r/r3rqrrrsrprtrurp_cr X_indicesX_indptrYrarrcX_datarblbrr\rs ` ``` @@@r&make_multilabel_classificationr@sn#<00I     + +C37799CYs^^N   J #:  ; ;E RVE " " ""EB C  I{3$$H A 9  !>##qI'''  WS^^2: 6 6 6F vy(3Iz;RSSSA  IIKK444 0@H0L N N N FFE)$$% & & 0 0 3 3 !S% a4Kr')r.r.rc t|}|df}|||}|dzddkt jd}d ||d k<||fS) aiGenerate data for binary classification used in Hastie et al. 2009, Example 10.2. The ten features are standard independent Gaussian and the target ``y`` is defined by:: y[i] = 1 if np.sum(X[i] ** 2) > 9.34 else -1 Read more in the :ref:`User Guide `. Parameters ---------- n_samples : int, default=12000 The number of samples. random_state : int, RandomState instance or None, default=None Determines random number generation for dataset creation. Pass an int for reproducible output across multiple function calls. See :term:`Glossary `. Returns ------- X : ndarray of shape (n_samples, 10) The input samples. y : ndarray of shape (n_samples,) The output values. See Also -------- make_gaussian_quantiles : A generalization of this dataset approach. References ---------- .. [1] T. Hastie, R. Tibshirani and J. Friedman, "Elements of Statistical Learning Ed. 2", Springer, 2009. Examples -------- >>> from sklearn.datasets import make_hastie_10_2 >>> X, y = make_hastie_10_2(n_samples=24000, random_state=42) >>> X.shape (24000, 10) >>> y.shape (24000,) >>> list(y[:5]) [np.float64(-1.0), np.float64(1.0), np.float64(-1.0), np.float64(1.0), np.float64(-1.0)] r@r(rygGz"@FrHr?)r normalr!rPrrr)r.rrsr|rbrcs r&make_hastie_10_2rsp L ) )B OE u %%e,,A S&1   $,,RZe,DDAAa3hK a4Kr') r.r/r0 n_targetsbiaseffective_rank tail_strengthnoiser coefrrrI) r0rrrrrr rrc t||}t| } || ||f} nt||||| } t j||f} d| ||fz| d|ddf<t j| | |z}|dkr|| ||j z }|r[t| || \} }t j |}| || dd|f| ddddf<| |} t j |}| r| |t j | fS| |fS)aS Generate a random regression problem. The input set can either be well conditioned (by default) or have a low rank-fat tail singular profile. See :func:`make_low_rank_matrix` for more details. The output is generated by applying a (potentially biased) random linear regression model with `n_informative` nonzero regressors to the previously generated input and some gaussian centered noise with some adjustable scale. Read more in the :ref:`User Guide `. Parameters ---------- n_samples : int, default=100 The number of samples. n_features : int, default=100 The number of features. n_informative : int, default=10 The number of informative features, i.e., the number of features used to build the linear model used to generate the output. n_targets : int, default=1 The number of regression targets, i.e., the dimension of the y output vector associated with a sample. By default, the output is a scalar. bias : float, default=0.0 The bias term in the underlying linear model. effective_rank : int, default=None If not None: The approximate number of singular vectors required to explain most of the input data by linear combinations. Using this kind of singular spectrum in the input allows the generator to reproduce the correlations often observed in practice. If None: The input set is well conditioned, centered and gaussian with unit variance. tail_strength : float, default=0.5 The relative importance of the fat noisy tail of the singular values profile if `effective_rank` is not None. When a float, it should be between 0 and 1. noise : float, default=0.0 The standard deviation of the gaussian noise applied to the output. shuffle : bool, default=True Shuffle the samples and the features. coef : bool, default=False If True, the coefficients of the underlying linear model are returned. random_state : int, RandomState instance or None, default=None Determines random number generation for dataset creation. Pass an int for reproducible output across multiple function calls. See :term:`Glossary `. Returns ------- X : ndarray of shape (n_samples, n_features) The input samples. y : ndarray of shape (n_samples,) or (n_samples, n_targets) The output values. coef : ndarray of shape (n_features,) or (n_features, n_targets) The coefficient of the underlying linear model. It is returned only if coef is True. Examples -------- >>> from sklearn.datasets import make_regression >>> X, y = make_regression(n_samples=5, n_features=2, noise=1, random_state=42) >>> X array([[ 0.4967..., -0.1382... ], [ 0.6476..., 1.523...], [-0.2341..., -0.2341...], [-0.4694..., 0.5425...], [ 1.579..., 0.7674...]]) >>> y array([ 6.737..., 37.79..., -10.27..., 0.4017..., 42.22...]) Nrr.r/rrrr<r?r:rr)minr rVmake_low_rank_matrixrrSrUrXrr|rZr[r squeeze)r.r/r0rrrrrr rrr\rb ground_truthrcrls r&make_regressionr8sh M22M"<00I  % %Iz+B % C C !!)'"    8Z344L&)I,=,=Y '->--'L-"# q,$&A s{{ Y  E  8 88-Aqy9991)J'''"""AAAwJ-!!!QQQ$#G,  1 A !RZ ----!t r')r.r rrfactorg?)r rrrct|tjr |dz}||z }n't|dkrt d|\}}t |}t jddt jz|d}t jddt jz|d} t j |} t j |} t j | |z} t j | |z} t j t j | | t j | | gj }t jt j|t jt j|t jg}|rt%|||\}}|||||j z }||fS) ahMake a large circle containing a smaller circle in 2d. A simple toy dataset to visualize clustering and classification algorithms. Read more in the :ref:`User Guide `. Parameters ---------- n_samples : int or tuple of shape (2,), dtype=int, default=100 If int, it is the total number of points generated. For odd numbers, the inner circle will have one point more than the outer circle. If two-element tuple, number of points in outer circle and inner circle. .. versionchanged:: 0.23 Added two-element tuple. shuffle : bool, default=True Whether to shuffle the samples. noise : float, default=None Standard deviation of Gaussian noise added to the data. random_state : int, RandomState instance or None, default=None Determines random number generation for dataset shuffling and noise. Pass an int for reproducible output across multiple function calls. See :term:`Glossary `. factor : float, default=.8 Scale factor between inner and outer circle in the range `[0, 1)`. Returns ------- X : ndarray of shape (n_samples, 2) The generated samples. y : ndarray of shape (n_samples,) The integer labels (0 or 1) for class membership of each sample. Examples -------- >>> from sklearn.datasets import make_circles >>> X, y = make_circles(random_state=42) >>> X.shape (100, 2) >>> y.shape (100,) >>> list(y[:5]) [np.int64(1), np.int64(1), np.int64(1), np.int64(0), np.int64(0)] rz7When a tuple, n_samples must have exactly two elements.rF)endpointrGrNr)rNnumbersrrMrJr rlinspacepicossinvstackrTrrSrYrrZrr|)r.r rrr n_samples_out n_samples_inr\ linspace_out linspace_in outer_circ_x outer_circ_y inner_circ_x inner_circ_yrbrcs r& make_circlesrsB)W-..0!Q  =0 y>>Q  VWW W&/# |"<00I;q!be)]UKKKL+aRUL5IIIK6,''L6,''L6+&&/L6+&&/L < . . , 0U0UV   -rw / / /RW1U1U1UV  A:Aqy9991  Y  E  8 88 a4Kr')r.r rr)r rrct|tjr |dz}||z }n) |\}}n"#t$r}td|d}~wwxYwt |}t jt jdt j|}t j t jdt j|} dt jt jdt j|z } dt j t jdt j|z dz } t j t j || t j | | gj } t j t j|t jt j|t jg} |rt#| | |\} } || ||| j z } | | fS) a%Make two interleaving half circles. A simple toy dataset to visualize clustering and classification algorithms. Read more in the :ref:`User Guide `. Parameters ---------- n_samples : int or tuple of shape (2,), dtype=int, default=100 If int, the total number of points generated. If two-element tuple, number of points in each of two moons. .. versionchanged:: 0.23 Added two-element tuple. shuffle : bool, default=True Whether to shuffle the samples. noise : float, default=None Standard deviation of Gaussian noise added to the data. random_state : int, RandomState instance or None, default=None Determines random number generation for dataset shuffling and noise. Pass an int for reproducible output across multiple function calls. See :term:`Glossary `. Returns ------- X : ndarray of shape (n_samples, 2) The generated samples. y : ndarray of shape (n_samples,) The integer labels (0 or 1) for class membership of each sample. Examples -------- >>> from sklearn.datasets import make_moons >>> X, y = make_moons(n_samples=200, noise=0.2, random_state=42) >>> X.shape (200, 2) >>> y.shape (200,) rz8`n_samples` can be either an int or a two-element tuple.Nrr(rIrGrr)rNrrrJr rrrrrrrrrrSrYrrZrr|)r.r rrrrer\rrrrrbrcs r& make_moonsr>sj)W-.. !Q  =0  *3 'M<<   J   #<00I6"+a >>??L6"+a >>??Lrvbk!RULAABBBLrvbk!RULAABBBSHL < . . , 0U0UV   -rw / / /RW1U1U1UV  A:Aqy9991  Y  E  8 88 a4Ks- A AA )r.r/centers cluster_std center_boxr rreturn_centers)g$g$@)rrrr rrc:t|}t|tjrq|d}t|tjr)|} ||d|d| |f}nt |}|jd}|jd} nt|} |&||d|d| |f}t|ts"td |t|| krtd|d|t |}|jd}t|d r6t|| kr#td ||t|tj r"tjt||}t|tr|} n;t|| zg| z} t!|| zD]} | | xxdz cc<tj| } tjt'| |ftj } tjt'| ft }t+t-| |D]O\} \}}| dkr | | dz nd}| | }||| |||f | ||<| |||<P|rt1| || \} }|r| ||fS| |fS)a Generate isotropic Gaussian blobs for clustering. Read more in the :ref:`User Guide `. Parameters ---------- n_samples : int or array-like, default=100 If int, it is the total number of points equally divided among clusters. If array-like, each element of the sequence indicates the number of samples per cluster. .. versionchanged:: v0.20 one can now pass an array-like to the ``n_samples`` parameter n_features : int, default=2 The number of features for each sample. centers : int or array-like of shape (n_centers, n_features), default=None The number of centers to generate, or the fixed center locations. If n_samples is an int and centers is None, 3 centers are generated. If n_samples is array-like, centers must be either None or an array of length equal to the length of n_samples. cluster_std : float or array-like of float, default=1.0 The standard deviation of the clusters. center_box : tuple of float (min, max), default=(-10.0, 10.0) The bounding box for each cluster center when centers are generated at random. shuffle : bool, default=True Shuffle the samples. random_state : int, RandomState instance or None, default=None Determines random number generation for dataset creation. Pass an int for reproducible output across multiple function calls. See :term:`Glossary `. return_centers : bool, default=False If True, then return the centers of each cluster. .. versionadded:: 0.23 Returns ------- X : ndarray of shape (n_samples, n_features) The generated samples. y : ndarray of shape (n_samples,) The integer labels for cluster membership of each sample. centers : ndarray of shape (n_centers, n_features) The centers of each cluster. Only returned if ``return_centers=True``. See Also -------- make_classification : A more intricate variant. Examples -------- >>> from sklearn.datasets import make_blobs >>> X, y = make_blobs(n_samples=10, centers=3, n_features=2, ... random_state=0) >>> print(X.shape) (10, 2) >>> y array([0, 0, 1, 0, 2, 2, 2, 1, 1, 0]) >>> X, y = make_blobs(n_samples=[3, 3, 4], centers=None, n_features=2, ... random_state=0) >>> print(X.shape) (10, 2) >>> y array([0, 1, 2, 0, 2, 2, 2, 1, 1, 0]) Nrr(rz8Parameter `centers` must be array-like. Got {!r} insteadzMLength of `n_samples` not consistent with number of centers. Got n_samples = z and centers = __len__zeLength of `clusters_std` not consistent with number of centers. Got centers = {} and cluster_std = {}r|rlocr:rr)r rNrrrUr r|rMrrJrLhasattrrrfullrCrRremptyrPrrWziprrZ)r.r/rrrr rrr\ n_centersn_samples_per_centerracum_sum_n_samplesrbrcrkstd start_idxend_idxs r& make_blobsrsH#<00I)W-..#& ?G gw/ 0 0 )I''1 z!}Iz3J(GG "'**G q)J a(II NN ?''1 z!}Iz3J(G'8,, JQQ  w<<9 $ $P,5PPFMPP g&&]1% {I&& 3{+;+;y+H+H ##)6';#?#?   +w|,,9gc'llK88 )X&&)( #I$: ; ;bjQQQA 01133???A %9;!G!GHH!! 8As01A%a!e,,1 #A&(// #Q O 0  )G  !)G :Aqy9991!W}!t r')r.r/rr)rrcnt|}|||f}dtjtj|dddfz|dddfzzd|dddfdz dzzzd|ddd fzzd |ddd fzz|||zz}||fS) aCGenerate the "Friedman #1" regression problem. This dataset is described in Friedman [1] and Breiman [2]. Inputs `X` are independent features uniformly distributed on the interval [0, 1]. The output `y` is created according to the formula:: y(X) = 10 * sin(pi * X[:, 0] * X[:, 1]) + 20 * (X[:, 2] - 0.5) ** 2 + 10 * X[:, 3] + 5 * X[:, 4] + noise * N(0, 1). Out of the `n_features` features, only 5 are actually used to compute `y`. The remaining features are independent of `y`. The number of features has to be >= 5. Read more in the :ref:`User Guide `. Parameters ---------- n_samples : int, default=100 The number of samples. n_features : int, default=10 The number of features. Should be at least 5. noise : float, default=0.0 The standard deviation of the gaussian noise applied to the output. random_state : int, RandomState instance or None, default=None Determines random number generation for dataset noise. Pass an int for reproducible output across multiple function calls. See :term:`Glossary `. Returns ------- X : ndarray of shape (n_samples, n_features) The input samples. y : ndarray of shape (n_samples,) The output values. References ---------- .. [1] J. Friedman, "Multivariate adaptive regression splines", The Annals of Statistics 19 (1), pages 1-67, 1991. .. [2] L. Breiman, "Bagging predictors", Machine Learning 24, pages 123-140, 1996. Examples -------- >>> from sklearn.datasets import make_friedman1 >>> X, y = make_friedman1(random_state=42) >>> X.shape (100, 10) >>> y.shape (100,) >>> list(y[:3]) [np.float64(16.8...), np.float64(5.8...), np.float64(9.4...)] rrNrr(r=rrIrrv)r rUrrrrV)r.r/rrr\rbrcs r&make_friedman1rIsL#<00I :677A RVBEAaaadGOa1g- . .. !!!Q$# !# # $ qAw,  a1g+  )++)+== =  > a4Kr')r.rrc t|}||df}|dddfxxdzcc<|dddfxxdtjzzcc<|dddfxxdtjzz cc<|ddd fxxd zcc<|ddd fxxdz cc<|dddfd z|dddf|ddd fzd|dddf|ddd fzz z d zzd z|||zz}||fS) aGenerate the "Friedman #2" regression problem. This dataset is described in Friedman [1] and Breiman [2]. Inputs `X` are 4 independent features uniformly distributed on the intervals:: 0 <= X[:, 0] <= 100, 40 * pi <= X[:, 1] <= 560 * pi, 0 <= X[:, 2] <= 1, 1 <= X[:, 3] <= 11. The output `y` is created according to the formula:: y(X) = (X[:, 0] ** 2 + (X[:, 1] * X[:, 2] - 1 / (X[:, 1] * X[:, 3])) ** 2) ** 0.5 + noise * N(0, 1). Read more in the :ref:`User Guide `. Parameters ---------- n_samples : int, default=100 The number of samples. noise : float, default=0.0 The standard deviation of the gaussian noise applied to the output. random_state : int, RandomState instance or None, default=None Determines random number generation for dataset noise. Pass an int for reproducible output across multiple function calls. See :term:`Glossary `. Returns ------- X : ndarray of shape (n_samples, 4) The input samples. y : ndarray of shape (n_samples,) The output values. References ---------- .. [1] J. Friedman, "Multivariate adaptive regression splines", The Annals of Statistics 19 (1), pages 1-67, 1991. .. [2] L. Breiman, "Bagging predictors", Machine Learning 24, pages 123-140, 1996. Examples -------- >>> from sklearn.datasets import make_friedman2 >>> X, y = make_friedman2(random_state=42) >>> X.shape (100, 4) >>> y.shape (100,) >>> list(y[:3]) [np.float64(1229.4...), np.float64(27.0...), np.float64(65.6...)] rrNrr<r((rrrrI)r rUrrrVr.rrr\rbrcs r&make_friedman2rsbH#<00I 1~..AaaadGGGsNGGGaaadGGGsRU{GGGaaadGGGrBEzGGGaaadGGGrMGGGaaadGGGqLGGG !!!Q$1 !!!Q$!AAAqD')A111a41QQQT71B,CCII  y00y0BBB CA a4Kr'c2t|}||df}|dddfxxdzcc<|dddfxxdtjzzcc<|dddfxxdtjzz cc<|ddd fxxd zcc<|ddd fxxdz cc<tj|dddf|ddd fzd|dddf|ddd fzz z |dddfz |||zz}||fS) aGenerate the "Friedman #3" regression problem. This dataset is described in Friedman [1] and Breiman [2]. Inputs `X` are 4 independent features uniformly distributed on the intervals:: 0 <= X[:, 0] <= 100, 40 * pi <= X[:, 1] <= 560 * pi, 0 <= X[:, 2] <= 1, 1 <= X[:, 3] <= 11. The output `y` is created according to the formula:: y(X) = arctan((X[:, 1] * X[:, 2] - 1 / (X[:, 1] * X[:, 3])) / X[:, 0]) + noise * N(0, 1). Read more in the :ref:`User Guide `. Parameters ---------- n_samples : int, default=100 The number of samples. noise : float, default=0.0 The standard deviation of the gaussian noise applied to the output. random_state : int, RandomState instance or None, default=None Determines random number generation for dataset noise. Pass an int for reproducible output across multiple function calls. See :term:`Glossary `. Returns ------- X : ndarray of shape (n_samples, 4) The input samples. y : ndarray of shape (n_samples,) The output values. References ---------- .. [1] J. Friedman, "Multivariate adaptive regression splines", The Annals of Statistics 19 (1), pages 1-67, 1991. .. [2] L. Breiman, "Bagging predictors", Machine Learning 24, pages 123-140, 1996. Examples -------- >>> from sklearn.datasets import make_friedman3 >>> X, y = make_friedman3(random_state=42) >>> X.shape (100, 4) >>> y.shape (100,) >>> list(y[:3]) [np.float64(1.5...), np.float64(0.9...), np.float64(0.4...)] rrNrr<r(rrrrr)r rUrrarctanrVrs r&make_friedman3rs]H#<00I 1~..AaaadGGGsNGGGaaadGGGsRU{GGGaaadGGGrBEzGGGaaadGGGrMGGGaaadGGGqLGGG 111a41QQQT7 Q!AAAqD'AaaadG"34 4!!!Q$?   )) );;; `. Parameters ---------- n_samples : int, default=100 The number of samples. n_features : int, default=100 The number of features. effective_rank : int, default=10 The approximate number of singular vectors required to explain most of the data by linear combinations. tail_strength : float, default=0.5 The relative importance of the fat noisy tail of the singular values profile. The value should be between 0 and 1. random_state : int, RandomState instance or None, default=None Determines random number generation for dataset creation. Pass an int for reproducible output across multiple function calls. See :term:`Glossary `. Returns ------- X : ndarray of shape (n_samples, n_features) The matrix. Examples -------- >>> from numpy.linalg import svd >>> from sklearn.datasets import make_low_rank_matrix >>> X = make_low_rank_matrix( ... n_samples=50, ... n_features=25, ... effective_rank=5, ... tail_strength=0.01, ... random_state=0, ... ) >>> X.shape (50, 25) reconomicF)mode check_finiterGr(rrg) r rrqrrVrr[rexpidentityrXr)r.r/rrrr\rkurv singular_indlow_ranktailss r&rrEs*`#<00I Iz""A 9!! 1~!66    DAq 9!! A!77    DAq9Qbj111LM!RVDL>4QVW3W,W%X%XXH 26$"5"FGG GD A(T/*A 6"&A,, $ $$r')r. n_componentsr/n_nonzero_coefsrc&t|}|||f}|tjtj|dzdz}tj||f}t |D]P}tj|} || | d|} |||| |f<Qtj ||} | j |j |j }}} ttj | ||fS)aGenerate a signal as a sparse combination of dictionary elements. Returns matrices `Y`, `D` and `X` such that `Y = XD` where `X` is of shape `(n_samples, n_components)`, `D` is of shape `(n_components, n_features)`, and each row of `X` has exactly `n_nonzero_coefs` non-zero elements. Read more in the :ref:`User Guide `. Parameters ---------- n_samples : int Number of samples to generate. n_components : int Number of components in the dictionary. n_features : int Number of features of the dataset to generate. n_nonzero_coefs : int Number of active (non-zero) coefficients in each sample. random_state : int, RandomState instance or None, default=None Determines random number generation for dataset creation. Pass an int for reproducible output across multiple function calls. See :term:`Glossary `. Returns ------- data : ndarray of shape (n_samples, n_features) The encoded signal (Y). dictionary : ndarray of shape (n_components, n_features) The dictionary with normalized components (D). code : ndarray of shape (n_samples, n_components) The sparse code such that each column of this matrix has exactly n_nonzero_coefs non-zero items (X). Examples -------- >>> from sklearn.datasets import make_sparse_coded_signal >>> data, dictionary, code = make_sparse_coded_signal( ... n_samples=50, ... n_components=100, ... n_features=10, ... n_nonzero_coefs=4, ... random_state=0 ... ) >>> data.shape (50, 10) >>> dictionary.shape (100, 10) >>> code.shape (50, 100) rrrryN) r rVrsqrtrPrSrRr[r rXrmapr) r.rr/rrr\Drbraidxrs r&make_sparse_coded_signalrsT#<00I !! L'A!BBAAQ''' ( ((A , *++A 9  DDi %%#"?"#--?-CC#q&  q! Ac13!qA rzAq!9 % %%r')r.r/rc"t|}|dd||f}||dddfd|dddfzzd|dddfzz d|dddfzz tj|}||fS) aiGenerate a random regression problem with sparse uncorrelated design. This dataset is described in Celeux et al [1]. as:: X ~ N(0, 1) y(X) = X[:, 0] + 2 * X[:, 1] - 2 * X[:, 2] - 1.5 * X[:, 3] Only the first 4 features are informative. The remaining features are useless. Read more in the :ref:`User Guide `. Parameters ---------- n_samples : int, default=100 The number of samples. n_features : int, default=10 The number of features. random_state : int, RandomState instance or None, default=None Determines random number generation for dataset creation. Pass an int for reproducible output across multiple function calls. See :term:`Glossary `. Returns ------- X : ndarray of shape (n_samples, n_features) The input samples. y : ndarray of shape (n_samples,) The output values. References ---------- .. [1] G. Celeux, M. El Anbari, J.-M. Marin, C. P. Robert, "Regularization in regression: comparing Bayesian and frequentist methods in a poorly informative situation", 2009. Examples -------- >>> from sklearn.datasets import make_sparse_uncorrelated >>> X, y = make_sparse_uncorrelated(random_state=0) >>> X.shape (100, 10) >>> y.shape (100,) rr(rNr?r)rr:)r rrr)r.r/rr\rbrcs r&make_sparse_uncorrelatedrsr#<00IQay*.EFFA qqq!tWq1QQQT7{ "Q111a4[ 03111a4= @gi     A a4Kr')n_dimrc Zt|}|||f}tjt j|j|d\}}}t jt j|dt j||z|}|S)akGenerate a random symmetric, positive-definite matrix. Read more in the :ref:`User Guide `. Parameters ---------- n_dim : int The matrix dimension. random_state : int, RandomState instance or None, default=None Determines random number generation for dataset creation. Pass an int for reproducible output across multiple function calls. See :term:`Glossary `. Returns ------- X : ndarray of shape (n_dim, n_dim) The random symmetric, positive-definite matrix. See Also -------- make_sparse_spd_matrix: Generate a sparse symmetric definite positive matrix. Examples -------- >>> from sklearn.datasets import make_spd_matrix >>> make_spd_matrix(n_dim=2, random_state=42) array([[2.09..., 0.34...], [0.34..., 0.21...]]) rF)rr>)r rUrsvdrrXrdiag)rrr\riUrVtrbs r&make_spd_matrixr TsL#<00Iu~..Az"&a..u===HAq" rvarwy'8'8e'8'D'DEEEFFKKA Hr'>bsrcoocsccsrdiadoklil)ralpha norm_diag smallest_coef largest_coef sparse_formatrgffffff?g?g?)rrrrrrcttj| }tj||d|z fd}tj|dd}|} || j| }||z }|j|z} |rCtjdtj | z } | | z| z} || S| |S)a4Generate a sparse symmetric definite positive matrix. Read more in the :ref:`User Guide `. Parameters ---------- n_dim : int, default=1 The size of the random matrix to generate. .. versionchanged:: 1.4 Renamed from ``dim`` to ``n_dim``. alpha : float, default=0.95 The probability that a coefficient is zero (see notes). Larger values enforce more sparsity. The value should be in the range 0 and 1. norm_diag : bool, default=False Whether to normalize the output matrix to make the leading diagonal elements all 1. smallest_coef : float, default=0.1 The value of the smallest coefficient between 0 and 1. largest_coef : float, default=0.9 The value of the largest coefficient between 0 and 1. sparse_format : str, default=None String representing the output sparse format, such as 'csc', 'csr', etc. If ``None``, return a dense numpy ndarray. .. versionadded:: 1.4 random_state : int, RandomState instance or None, default=None Determines random number generation for dataset creation. Pass an int for reproducible output across multiple function calls. See :term:`Glossary `. Returns ------- prec : ndarray or sparse matrix of shape (dim, dim) The generated matrix. If ``sparse_format=None``, this would be an ndarray. Otherwise, this will be a sparse matrix of the specified format. See Also -------- make_spd_matrix : Generate a random symmetric, positive-definite matrix. Notes ----- The sparsity is actually imposed on the cholesky factor of the matrix. Thus alpha does not translate directly into the filling fraction of the matrix itself. Examples -------- >>> from sklearn.datasets import make_sparse_spd_matrix >>> make_sparse_spd_matrix(n_dim=4, norm_diag=False, random_state=42) array([[1., 0., 0., 0.], [0., 1., 0., 0.], [0., 0., 1., 0.], [0., 0., 0., 1.]]) r(c4|S)N)lowhighr)rU)xrrrs r&z(make_sparse_spd_matrix..s#<//Lq0  r')mrkdensitydata_rvsrrr )rErLr>) r reyerandomtril permutationrdiagsrrdiagonalrasformat) rrrrrrrcholauxr"precds `` ` r&make_sparse_spd_matrixr*sn&l33L F5MM>D )  E       "   C '#E * * *C**511K k  [ )CCKD 6D=D HS274==??333 4 44x!|||~~}}]+++r')r.rrhole)rrr+ct|}|sFdtjzdd||zzz}d||z}ntjdt dD}tj|dd }|d |}|d|ftjtjgd ggz} ||j| z\}}|tj |z} |tj |z} tj | || f} | || d|fzz } | j} tj |}| |fS) aNGenerate a swiss roll dataset. Read more in the :ref:`User Guide `. Parameters ---------- n_samples : int, default=100 The number of sample points on the Swiss Roll. noise : float, default=0.0 The standard deviation of the gaussian noise. random_state : int, RandomState instance or None, default=None Determines random number generation for dataset creation. Pass an int for reproducible output across multiple function calls. See :term:`Glossary `. hole : bool, default=False If True generates the swiss roll with hole dataset. Returns ------- X : ndarray of shape (n_samples, 3) The points. t : ndarray of shape (n_samples,) The univariate position of the sample according to the main dimension of the points in the manifold. Notes ----- The algorithm is from Marsland [1]. References ---------- .. [1] S. Marsland, "Machine Learning: An Algorithmic Perspective", 2nd edition, Chapter 6, 2014. https://homepages.ecs.vuw.ac.nz/~marslast/Code/Ch6/lle.py Examples -------- >>> from sklearn.datasets import make_swiss_roll >>> X, t = make_swiss_roll(noise=0.05, random_state=0) >>> X.shape (100, 3) >>> t.shape (100,) rr(rrc`g|]+}tdD]}tjd|zz|dzg,S)rr)rRrr)rDrajs r&rFz#make_swiss_roll..<s? L L LA588 L LabesQwQ ' L L L Lr'rrrryr/)r rrrUrrRdeletechoicerrrrrVr) r.rrr+r\trccorners corner_index parametersrzrbs r&make_swiss_rollr9szt#<00I  4 "%K1q9#4#4)#4#D#DDD E "" "22 2( L Lq L L L  )GQQ/// ''955 &&Q N&;;bhRSQT~>V>VV |$&31 BF1II A BF1II A 1a)A**I*?? ??A A 1 A a4Kr'c t|}dtjz|d|fdz z}tj|dftj}tj||dddf<d||z|dddf<tj|tj|dz z|ddd f<||| d|fj zz }tj |}||fS) aGenerate an S curve dataset. Read more in the :ref:`User Guide `. Parameters ---------- n_samples : int, default=100 The number of sample points on the S curve. noise : float, default=0.0 The standard deviation of the gaussian noise. random_state : int, RandomState instance or None, default=None Determines random number generation for dataset creation. Pass an int for reproducible output across multiple function calls. See :term:`Glossary `. Returns ------- X : ndarray of shape (n_samples, 3) The points. t : ndarray of shape (n_samples,) The univariate position of the sample according to the main dimension of the points in the manifold. Examples -------- >>> from sklearn.datasets import make_s_curve >>> X, t = make_s_curve(noise=0.05, random_state=0) >>> X.shape (100, 3) >>> t.shape (100,) rr(rrIrNrrr) r rrrUrrrsignrrVrr)r.rrr\r4rbs r& make_s_curver<NsX#<00I BE Y&&Q N&;;cABA  1~RZ888AfQiiAaaadGI%%9%555AaaadGgajjBF1IIM*AaaadG**I*??A AAA 1 A a4Kr')meancovr.r/r3r rrc ||krtdt|}|tj|}ntj|}|||tj|z|f}tjtj||tj ddfz dzd} || ddf}||z} tj tj tj || tj |dz || |zz g} |rt|| |\}} || fS)aGenerate isotropic Gaussian and label samples by quantile. This classification dataset is constructed by taking a multi-dimensional standard normal distribution and defining classes separated by nested concentric multi-dimensional spheres such that roughly equal numbers of samples are in each class (quantiles of the :math:`\chi^2` distribution). Read more in the :ref:`User Guide `. Parameters ---------- mean : array-like of shape (n_features,), default=None The mean of the multi-dimensional normal distribution. If None then use the origin (0, 0, ...). cov : float, default=1.0 The covariance matrix will be this value times the unit matrix. This dataset only produces symmetric normal distributions. n_samples : int, default=100 The total number of points equally divided among classes. n_features : int, default=2 The number of features for each sample. n_classes : int, default=3 The number of classes. shuffle : bool, default=True Shuffle the samples. random_state : int, RandomState instance or None, default=None Determines random number generation for dataset creation. Pass an int for reproducible output across multiple function calls. See :term:`Glossary `. Returns ------- X : ndarray of shape (n_samples, n_features) The generated samples. y : ndarray of shape (n_samples,) The integer labels for quantile membership of each sample. Notes ----- The dataset is from Zhu et al [1]. References ---------- .. [1] J. Zhu, H. Zou, S. Rosset, T. Hastie, "Multi-class AdaBoost", 2009. Examples -------- >>> from sklearn.datasets import make_gaussian_quantiles >>> X, y = make_gaussian_quantiles(random_state=42) >>> X.shape (100, 2) >>> y.shape (100,) >>> list(y[:5]) [np.int64(2), np.int64(0), np.int64(1), np.int64(0), np.int64(2)] z$n_samples must be at least n_classesNrr(ryr)rJr rrSrmultivariate_normalrargsortrPnewaxisrrepeatr[rZ) r=r>r.r/r3r rr\rbrsteprcs r&make_gaussian_quantilesrEsFj9?@@@"<00I |x ##x~~ %%dC"+j2I2I,II<XXA *RVQbj!!!m!44:CCC D DC #qqq& A  !D Ibi **D 1 1 Ii!mY 1A%A B B   A:Aqy9991 a4Kr'ct|}|j\}}||}||}||dd|f}|||fSN)r r|r")datarr\n_rowsn_colsrow_idxcol_idxresults r&_shufflerNsc"<00IZNFF##F++G##F++G ']111g: &F 7G ##r')r|r_rminvalmaxvalr r)rrOrPr rct|}|\}} ||||} ||tjd|z |} || tjd|z |} tjdt t|| Dtjdt t|| Dtj|tj } t|D]5}tj |k|k}| |xx| |z cc<6|dkr| | || j z } |r$t| |\} }}||tjfdt|D}tjfdt|D}| ||fS) aSGenerate a constant block diagonal structure array for biclustering. Read more in the :ref:`User Guide `. Parameters ---------- shape : tuple of shape (n_rows, n_cols) The shape of the result. n_clusters : int The number of biclusters. noise : float, default=0.0 The standard deviation of the gaussian noise. minval : float, default=10 Minimum value of a bicluster. maxval : float, default=100 Maximum value of a bicluster. shuffle : bool, default=True Shuffle the samples. random_state : int, RandomState instance or None, default=None Determines random number generation for dataset creation. Pass an int for reproducible output across multiple function calls. See :term:`Glossary `. Returns ------- X : ndarray of shape `shape` The generated array. rows : ndarray of shape (n_clusters, X.shape[0]) The indicators for cluster membership of each row. cols : ndarray of shape (n_clusters, X.shape[1]) The indicators for cluster membership of each column. See Also -------- make_checkerboard: Generate an array with block checkerboard structure for biclustering. References ---------- .. [1] Dhillon, I. S. (2001, August). Co-clustering documents and words using bipartite spectral graph partitioning. In Proceedings of the seventh ACM SIGKDD international conference on Knowledge discovery and data mining (pp. 269-274). ACM. Examples -------- >>> from sklearn.datasets import make_biclusters >>> data, rows, cols = make_biclusters( ... shape=(10, 20), n_clusters=2, random_state=42 ... ) >>> data.shape (10, 20) >>> rows.shape (2, 10) >>> cols.shape (2, 20) r>c>g|]\}}tj||SrBrrCrDvalreps r&rFz#make_biclusters..g(OOOc3  OOOr'c>g|]\}}tj||SrBrSrTs r&rFz#make_biclusters..jrWr'rGrrcg|]}|k SrBrB)rDr row_labelss r&rFz#make_biclusters..zAAA!jAoAAAr'cg|]}|k SrBrB)rDr col_labelss r&rFz#make_biclusters..{r[r')r rU multinomialrrCrrrRrSrouterrr|rNr)r|r_rrOrPr rr\rIrJconsts row_sizes col_sizesrMraselectorrKrLrowscolsr]rZs @@r&make_biclustersrfsp#<00INFF   vvz : :F%%fbij8H*.U.UVVI%%fbij8H*.U.UVVIOOSz1B1BI-N-NOOOJOOSz1B1BI-N-NOOOJXe2: . . .F :  &&8J!OZ1_==xF1I% qyy)""V\"BBB)#+FL#A#A ( ( 9AAAAuZ/@/@AAA B BD 9AAAAuZ/@/@AAA B BD 4 r'c Zt|}t|dr|\}n|x}|\} } || tjd|z |} || tjdz } tjdt t|| Dtjdt t| Dtj|tj } t|D]U}tD]C}tj |k|k}| |xx| ||z cc<DV|dkr| | || j z } |r$t| |\} }}||tjfdt|D}tjfd t|D}| ||fS) aGenerate an array with block checkerboard structure for biclustering. Read more in the :ref:`User Guide `. Parameters ---------- shape : tuple of shape (n_rows, n_cols) The shape of the result. n_clusters : int or array-like or shape (n_row_clusters, n_column_clusters) The number of row and column clusters. noise : float, default=0.0 The standard deviation of the gaussian noise. minval : float, default=10 Minimum value of a bicluster. maxval : float, default=100 Maximum value of a bicluster. shuffle : bool, default=True Shuffle the samples. random_state : int, RandomState instance or None, default=None Determines random number generation for dataset creation. Pass an int for reproducible output across multiple function calls. See :term:`Glossary `. Returns ------- X : ndarray of shape `shape` The generated array. rows : ndarray of shape (n_clusters, X.shape[0]) The indicators for cluster membership of each row. cols : ndarray of shape (n_clusters, X.shape[1]) The indicators for cluster membership of each column. See Also -------- make_biclusters : Generate an array with constant block diagonal structure for biclustering. References ---------- .. [1] Kluger, Y., Basri, R., Chang, J. T., & Gerstein, M. (2003). Spectral biclustering of microarray data: coclustering genes and conditions. Genome research, 13(4), 703-716. Examples -------- >>> from sklearn.datasets import make_checkerboard >>> data, rows, columns = make_checkerboard(shape=(300, 300), n_clusters=10, ... random_state=42) >>> data.shape (300, 300) >>> rows.shape (100, 300) >>> columns.shape (100, 300) >>> print(rows[0][:5], columns[0][:5]) [False False False True False] [False False False False False] rr>c>g|]\}}tj||SrBrSrTs r&rFz%make_checkerboard..(SSSc3  SSSr'c>g|]\}}tj||SrBrSrTs r&rFz%make_checkerboard..rir'rGrrc@g|]}tD]}|k SrBrR)rDlabelrn_col_clustersrZs r&rFz%make_checkerboard..sM   >**   %     r'c@g|]}tD]}|k SrBrl)rDrrmr]rns r&rFz%make_checkerboard.. sM   ~..   %     r')r rr^rrCrrrRrSrr_rUrr|rNr)r|r_rrOrPr rr\n_row_clustersrIrJrarbrMrar0rcrKrLrdrer]rnrZs @@@r&make_checkerboardrqsn#<00Iz9%%5)3&*44NFF%% #.??I%% #.??ISSS~1F1F -R-RSSSJSSS~1F1F -R-RSSSJXe2: . . .F > " "BB~&& B BAx aqAAH 8    1 1&& A A A     B qyy)""V\"BBB)#+FL#A#A ( ( 9     ~..     D 9     >**     D 4 r')r<r=)r)r<r<)r<)r<r)r<r)r(rG)2__doc__rrcollections.abcrrrnumpyr scipy.sparserprscipyr preprocessingrutilsr r r rZutils._param_validationr r r utils.randomrrrnrrrtuplerrrrrrrrrr r*r9r<rErNrfrqrBr'r&r|s $$$$$$""""""""//////33333333++++++KKKKKKKKKK555555    hxD@@@Ax!T&AAAB"(8QVDDDE 1d6BBBCx!T&AAABhxD@@@A!)(AtF!K!K!K L $'8D!Qv6667htQY???@[(4tI>>> dS(4D;;;\4P;'("#'%*A   !AAAA)(AHhxD@@@Ax!T&AAABhxD@@@AXh4???@8Haf===>%;+'Z((;<#8HafEEEtL"(4Af===>(4D8889; '(  #'"R   RRRR! RjhxD@@@%H;(4D888$?'(8D!Qv6667 #'   T"$T#TTTT  TnhxD@@@%H;(4D888$?'(  #'KTKKKKK\hxD@@@,Ox!T&AAABHXq$v>>> dS q$v>>> Mg;'($+  #'   d  dddd  dNhxD@@@Ax!T&AAAB(4d6:::;'(  #'H#DHHHHHVhxD@@@A(4D8889'( #' I3TIIIIIXhxD@@@A(4D8889'( #' I3TIIIIIXhxD@@@Ax!T&AAAB#8HafEEEF"(4Af===>'( #'   ]% ]%]%]%]%  ]%@hxD@@@A!(AtFCCCDx!T&AAAB$HXq$vFFFG'( #'     T&T&T&T&  T&nhxD@@@Ax!T&AAAB'( #' 9499999x(8QV<<<='(#' ,0% % % % % P(8QV<<<=(4Af5556["(4Af===>!$1V<<<= JHHH I I  ((  #' f, f,f,f,f,f,RhxD@@@A(4D8889'(  #'GCdGGGGGThxD@@@A(4D8889'( #' .4.....bt$q$v6667hxD@@@Ax!T&AAABhxD@@@A;'(#'      gggg  gT$$$$x!T&AAAB(4D88898D$Y???@8D$Y???@;'(#'      kkkk  k\x!T&AAA<P(4D88898D$Y???@8D$Y???@;'(#'        r'