← Back
Editing: test_resampling.cpython-311.pyc
� d�cd � � � d dl Zd dlZd dlmZmZmZ d dlmZm Z m Z d dlmZ d dlm Z ddlmZ d dlmZ d d lmZ d � Zej � dg d�� � ej � d g d�� � d� � � � � Zej � dg d�� � d� � � Zej � dg d�� � ej � d g d�� � ej � dddg� � d� � � � � � � Zej � dg d�� � d� � � Zdddd�Zej � de� � � � � d� � � Zdddd�Zd� Zd� Z ej � de� � � � � d� � � Z!d d!d"�Z"ej � de"� � � � � d#� � � Z#ej � dd$d%g� � ej � d d d&g� � d'� � � � � Z$ej �% d(� � ej � dg d�� � ej � d d d&g� � d)� � � � � � � Z&ej � dg d�� � d*� � � Z'ej � dg d�� � d+� � � Z(d,� Z)ej � d-d d.g� � d/� � � Z*d0� Z+ej � d1d2d3g� � d4� � � Z,ej � d5g d6�� � ej � d g d�� � d7� � � � � Z-d8� Z.ej � d g d�� � d9� � � Z/ej � dg d�� � d:� � � Z0ej j1 ej �2 d;� � d<� � � � � Z3 G d=� d>� � Z4 G d?� d@� � Z5dA� Z6ej � dBg dC�� � dD� � � Z7dS )E� N)� bootstrap�monte_carlo_test�permutation_test)�assert_allclose�assert_equal�suppress_warnings)�stats)�special� )�_resampling)�rng_integers)�rootc �� � d} t j t | �� � 5 t dt j � � d d d � � n# 1 swxY w Y d} t j t | �� � 5 t t � � t j � � d d d � � n# 1 swxY w Y d} t j t | �� � 5 t g d�dgft j � � d d d � � n# 1 swxY w Y d} t j t | �� � 5 t g d�g d�ft j d � � � d d d � � n# 1 swxY w Y d} t j t | �� � 5 t dt j d� � � d d d � � n# 1 swxY w Y d} t j t | �� � 5 t g d�ft j d�� � d d d � � n# 1 swxY w Y d} t j t | �� � 5 t g d�ft j d�� � d d d � � n# 1 swxY w Y d} t j t | �� � 5 t g d�ft j d�� � d d d � � n# 1 swxY w Y d} t j t | �� � 5 t g d�ft j d�� � d d d � � n# 1 swxY w Y d} t j t | �� � 5 t g d�ft j d�� � d d d � � n# 1 swxY w Y d} t j t | �� � 5 t g d�ft j d�� � d d d � � n# 1 swxY w Y d} t j t | �� � 5 t g d�ft j d�� � d d d � � n# 1 swxY w Y d} t j t | �� � 5 t g d�ft j d�� � d d d � � n# 1 swxY w Y d} t j t | �� � 5 t g d�ft j d �� � d d d � � n# 1 swxY w Y d!} t j t | �� � 5 t g d�ft j d"�#� � d d d � � d S # 1 swxY w Y d S )$Nz%`data` must be a sequence of samples.��match� z(`data` must contain at least one sample.z>each sample in `data` must contain two or more observations...�r r � z=When `paired is True`, all samples must have the same length �r r r � T)�paired�0`vectorized` must be `True`, `False`, or `None`.�ekki�� vectorized�`axis` must be an integer.� �?��axisz!could not convert string to float�ni)�confidence_levelz-`n_resamples` must be a non-negative integer.������n_resamples� D�@�+`batch` must be a positive integer or None.��batchz`method` must be in��methodz?`bootstrap_result` must have attribute `bootstrap_distribution'� )�bootstrap_resultz5Either `bootstrap_result.bootstrap_distribution.size`r �"'herring' cannot be used to seed a�herring��random_state)�pytest�raises� ValueErrorr �np�mean�tuple)�messages �C/usr/lib/python3/dist-packages/scipy/stats/tests/test_resampling.py�test_bootstrap_ivr9 s� � �5�G� ��z�� 1� 1� 1� � ��!�R�W����� � � � � � � � � � ���� � � � � 9�G� ��z�� 1� 1� 1� $� $��%�'�'�2�7�#�#�#�$� $� $� $� $� $� $� $� $� $� $���� $� $� $� $� O�G� ��z�� 1� 1� 1� -� -��9�9�9�q�c�"�B�G�,�,�,�-� -� -� -� -� -� -� -� -� -� -���� -� -� -� -� O�G� ��z�� 1� 1� 1� C� C��9�9�9�l�l�l�+�R�W�T�B�B�B�B�C� C� C� C� C� C� C� C� C� C� C���� C� C� C� C� A�G� ��z�� 1� 1� 1� 1� 1��!�R�W��0�0�0�0�1� 1� 1� 1� 1� 1� 1� 1� 1� 1� 1���� 1� 1� 1� 1� +�G� ��z�� 1� 1� 1� 3� 3��9�9�9�,���c�2�2�2�2�3� 3� 3� 3� 3� 3� 3� 3� 3� 3� 3���� 3� 3� 3� 3� 2�G� ��z�� 1� 1� 1� @� @��9�9�9�,���$�?�?�?�?�@� @� @� @� @� @� @� @� @� @� @���� @� @� @� @� >�G� ��z�� 1� 1� 1� <� <��9�9�9�,���U�;�;�;�;�<� <� <� <� <� <� <� <� <� <� <���� <� <� <� <� >�G� ��z�� 1� 1� 1� =� =��9�9�9�,���V�<�<�<�<�=� =� =� =� =� =� =� =� =� =� =���� =� =� =� =� <�G� ��z�� 1� 1� 1� 6� 6��9�9�9�,���u�5�5�5�5�6� 6� 6� 6� 6� 6� 6� 6� 6� 6� 6���� 6� 6� 6� 6� <�G� ��z�� 1� 1� 1� 7� 7��9�9�9�,���v�6�6�6�6�7� 7� 7� 7� 7� 7� 7� 7� 7� 7� 7���� 7� 7� 7� 7� $�G� ��z�� 1� 1� 1� 8� 8��9�9�9�,����7�7�7�7�8� 8� 8� 8� 8� 8� 8� 8� 8� 8� 8���� 8� 8� 8� 8� P�G� ��z�� 1� 1� 1� >� >��9�9�9�,���"�=�=�=�=�>� >� >� >� >� >� >� >� >� >� >���� >� >� >� >� F�G� ��z�� 1� 1� 1� 8� 8��9�9�9�,���Q�7�7�7�7�8� 8� 8� 8� 8� 8� 8� 8� 8� 8� 8���� 8� 8� 8� 8� 3�G� ��z�� 1� 1� 1� A� A��9�9�9�,���i�@�@�@�@�A� A� A� A� A� A� A� A� A� A� A� A���� A� A� A� A� A� As �A�A �A �-'B � B$�'B$� C4�4C8�;C8�#E�E�E�3F�F �#F � G0�0G4�7G4� I�I�I�, J�J�J� K,�,K0�3K0� M � M�M�( N�N�N�< O(�(O,�/O,� P<�<Q �Q �$ R�R�R�8 S%�%S)�,S)r* )�basic� percentile�BCar �r r r c �� � t j � d� � t j � ddd� � }t |ft j d | d|d�� � }t |ft j d| d|d�� � }t |j j |j j � � t |j j |j j � � t |j |j � � d S )Nr r+ � � �d )r( r* r0 r r$ )r4 �random�seed�randr r5 r �confidence_interval�low�high�standard_error)r* r �x�res1�res2s r8 �test_bootstrap_batchrL K s� � � �I�N�N�1���� � ���r�2�r�"�"�A��a�T�2�7�$�v�"#�$�C�A� A� A�D��a�T�2�7�"�V�"#�$�C�A� A� A�D� ��)�-�t�/G�/K�L�L�L���)�.��0H�0M�N�N�N���$�d�&9�:�:�:�:�:� c �� ���� t j � d� � d}t j � |� � �t j � |� � �d d��d ���fd� }t j t �� � � � }t |f|d�� � }t ��f�dd�� � }t |j |j � � t |j |j � � d S ) Nr rA ���c �: � | |z dz � |�� � S )Nr r �r5 �rI �yr s r8 �my_statisticz+test_bootstrap_paired.<locals>.my_statisticd s � ��1��q����D��)�)�)rM c �@ �� �| }�| } �||� � }|S �N� )�ir �a�b�resrT rI rS s ���r8 �my_paired_statisticz2test_bootstrap_paired.<locals>.my_paired_statisticg s+ �� � �a�D�� �a�D���l�1�a� � ��� rM r/ T)r r0 �rO ) r4 rB rC rD �arange�lenr r rE rH ) r* �nr\ rX rJ rK rT rI rS s @@@r8 �test_bootstrap_pairedra \ s� ���� � �I�N�N�1�����A� � ���q���A� � ���q���A�*� *� *� *�� � � � � � � � � �#�a�&�&���A��a�T�.�Q�?�?�?�D��a��V�\�$�Q�G�G�G�D��D�,�d�.F�G�G�G��D�'��)<�=�=�=�=�=rM r TFc �z � t j � d� � d d�}d}|| }t j � |� � }t j � |� � }t j � |� � }t |||f||| ddd�� � } | j j | j j dz k sJ �g d�} || |<