U mdz#@s~dddgZddlZddlZddlmZddlmZddlm Z m Z ddZ d d Z d dZ d d ZGdddeZddZdS)balanced demo_data LookupFactorN) PatsyError)C) no_picklingassert_no_picklingc s|dd}g}t|}|D].|}|fddtd|dDqttj|}i}t||D]\}t|||<qh|S)aTbalanced(factor_name=num_levels, [factor_name=num_levels, ..., repeat=1]) Create simple balanced factorial designs for testing. Given some factor names and the number of desired levels for each, generates a balanced factorial design in the form of a data dictionary. For example: .. ipython:: In [1]: balanced(a=2, b=3) Out[1]: {'a': ['a1', 'a1', 'a1', 'a2', 'a2', 'a2'], 'b': ['b1', 'b2', 'b3', 'b1', 'b2', 'b3']} By default it produces exactly one instance of each combination of levels, but if you want multiple replicates this can be accomplished via the `repeat` argument: .. ipython:: In [2]: balanced(a=2, b=2, repeat=2) Out[2]: {'a': ['a1', 'a1', 'a2', 'a2', 'a1', 'a1', 'a2', 'a2'], 'b': ['b1', 'b2', 'b1', 'b2', 'b1', 'b2', 'b1', 'b2']} repeatcsg|]}d|fqS)z%s%s).0inamer H/home/sam/Atlas/atlas_env/lib/python3.8/site-packages/patsy/user_util.py 1szbalanced..)popsortedappendrangezip itertoolsproductlist)kwargsr levelsnamesZ level_countvaluesdatavaluer rrrs $c Cstddd}|dddddddgks(t|ddd d dd d gksDttdddd }|dddddddddddddg kszt|ddd d dd d dd d dd d g kstdS) N)abr"a1a2r#b1b2Zb3)r"r#r )rAssertionError)rr r r test_balanced9s(  r)c Os|dd}|dd}|r*td|ft}i}|D]@}|ddkrR|||<q8|ddkrj||q8td |fq8tjt|t d }t t |d |}||} t fd |i|} tj d} t|D]}| j| d | |<q| S)a?demo_data(*names, nlevels=2, min_rows=5) Create simple categorical/numerical demo data. Pass in a set of variable names, and this function will return a simple data set using those variable names. Names whose first letter falls in the range "a" through "m" will be made categorical (with `nlevels` levels). Those that start with a "p" through "z" are numerical. We attempt to produce a balanced design on the categorical variables, repeating as necessary to generate at least `min_rows` data points. Categorical variables are returned as a list of strings. Numerical data is generated by sampling from a normal distribution. A fixed random seed is used, so that identical calls to demo_data() will produce identical results. Numerical data is returned in a numpy array. Example: .. ipython: In [1]: patsy.demo_data("a", "b", "x", "y") Out[1]: {'a': ['a1', 'a1', 'a2', 'a2', 'a1', 'a1', 'a2', 'a2'], 'b': ['b1', 'b2', 'b1', 'b2', 'b1', 'b2', 'b1', 'b2'], 'x': array([ 1.76405235, 0.40015721, 0.97873798, 2.2408932 , 1.86755799, -0.97727788, 0.95008842, -0.15135721]), 'y': array([-0.10321885, 0.4105985 , 0.14404357, 1.45427351, 0.76103773, 0.12167502, 0.44386323, 0.33367433])} nlevelsr min_rowszunexpected keyword arguments %rrZabcdefghijklmnZ pqrstuvwxyzz bad name %r)dtypeg?r )size)r TypeErrorsetaddrnpprodrrintceilrrandomZ RandomStaternormal) rrr*r+Z numericalZ categoricalrZbalanced_design_sizer Znum_rowsrrr r rrCs(!        c Csptddd}t|dddgks&t|dddddddddgksFt|dddddddddgksft|djttks~t|djdksttdd }t|dd gkstt|dt|d krd ksntttdd d dd kstttdddd d dd kstttdddd ddddks>tddl }| t tddd|j t tdddddS)Nr"r#xr$r%r&r')yr, )r+ r!)r+r*rZ__123{)Zasdfasdf) rrkeysr(r-r2floatshapelenpytestraisesrr/)Zd1Zd2rDr r rtest_demo_datazs    *"$rFc@sfeZdZdZdddZddZdd Zd d Zd d ZddZ ddZ ddZ ddZ ddZ eZdS)ra8A simple factor class that simply looks up a named entry in the given data. Useful for programatically constructing formulas, and as a simple example of the factor protocol. For details see :ref:`expert-model-specification`. Example:: dmatrix(ModelDesc([], [Term([LookupFactor("x")])]), {"x": [1, 2, 3]}) :arg varname: The name of this variable; used as a lookup key in the passed in data dictionary/DataFrame/whatever. :arg force_categorical: If True, then treat this factor as categorical. (Equivalent to using :func:`C` in a regular formula, but of course you can't do that with a :class:`LookupFactor`. :arg contrast: If given, the contrast to use; see :func:`C`. (Requires ``force_categorical=True``.) :arg levels: If given, the categorical levels; see :func:`C`. (Requires ``force_categorical=True``.) :arg origin: Either ``None``, or the :class:`Origin` of this factor for use in error reporting. .. versionadded:: 0.2.0 The ``force_categorical`` and related arguments. FNcCsH||_||_||_||_||_|jsD|dk r4td|dk rDtddS)Nz)contrast= requires force_categorical=Truez'levels= requires force_categorical=True)_varname_force_categorical _contrast_levelsorigin ValueError)selfvarnameforce_categoricalcontrastrrKr r r__init__szLookupFactor.__init__cCs|jSN)rGrMr r rrszLookupFactor.namecCsd|jj|jfS)Nz%s(%r)) __class____name__rGrSr r r__repr__szLookupFactor.__repr__cCs:t|to8|j|jko8|j|jko8|j|jko8|j|jkSrR) isinstancerrGrHrIrJrMotherr r r__eq__s     zLookupFactor.__eq__cCs ||k SrRr rXr r r__ne__szLookupFactor.__ne__cCstt|j|j|j|jfSrR)hashrrGrHrIrJrSr r r__hash__s zLookupFactor.__hash__cCsdS)Nrr )rMstateZeval_envr r rmemorize_passes_neededsz#LookupFactor.memorize_passes_neededcCs dstdSNFr()rMr^ which_passrr r rmemorize_chunkszLookupFactor.memorize_chunkcCs dstdSr`ra)rMr^rbr r rmemorize_finishszLookupFactor.memorize_finishcCs&||j}|jr"t||j|jd}|S)N)rPr)rGrHrrIrJ)rMZmemorize_staterrr r revals zLookupFactor.eval)FNNN)rU __module__ __qualname____doc__rQrrVrZr[r]r_rcrdrer __getstate__r r r rrs  cCshtd}|dkst|tdks(t|tdks8tt|ttdksPtt|ttdksht|iddidkst|iddidkstt|dkst|jdksttddd}|jdksttdd d d d }|iddddgi}|jdddgks t|jd kst|j d ks,td dl }|j t tdd d|j t tdd dt tddS)Nr"r#r r zLookupFactor('a')Zasdf)rKcTZCONTRAST)r r )rOrPrrnc)rP)r)rrr(r\rereprrKrrPrrDrErLr)Zl_aZ l_with_originZl_cboxrDr r rtest_LookupFactors0 rn)__all__rnumpyr2ZpatsyrZpatsy.categoricalrZ patsy.utilrrrr)rrFobjectrrnr r r r s   ( 7N