a ߙfb @sPddlmZddlmZmZddlmZmZddZddZdd Z d d Z d S) )assert_allclose)bayesian_info_criterionbayesian_info_criterion_lsq)akaike_info_criterionakaike_info_criterion_lsqcCsNd}d}d}d}t|d|d|}t|d|d|}t|||dddS) N)g fge)dg(\@rg?Zatol)rr)lnLn_params n_samplesanswerZbic_gZbic_trClib/python3.9/site-packages/astropy/stats/tests/test_info_theory.pytest_bayesian_info_criterionsrcCsNd}d}d}d}t|d|d|}t|d|d|}t|||dddS) Ny)gRQ gGz)gffffff?rr {Gz?r )rr)rr r rZaic_1Zaic_2rrrtest_akaike_info_criterionsrcCszd}d}d}d}t|dt|d|d|ddt|dt|d|d|ddt|d t|d |d |dddS) Nr )rr)g9@g:@g;@)gQF`gQ`gQ+_rrr r r)rrrr Zssrrrrrtest_akaike_info_criterion_lsqs rcCszd}d}d}d}t|dt|d|d|ddt|dt|d|d|ddt|d t|d |d |ddd S) a This test is from: http://www.statoek.wiso.uni-goettingen.de/veranstaltungen/non_semi_models/ AkaikeLsg.pdf Note that in there, they compute a "normalized BIC". Therefore, the answers presented here are recalculated versions based on their values. )r rr )i?ii\)g;Oh@g;O6g@gQKg@rrr r rN)rrrrrr test_bayesian_info_criterion_lsq0s,rN) Z numpy.testingrZastropy.stats.info_theoryrrrrrrrrrrrrs