
    i                        d dl Z d dlmZ d dlZd dlmZ d dlmZ d dlmZ d dl	m
Z
mZmZmZmZmZ 	 	 ddeded	ed
ed   dee   defdZdededefdZ	 	 	 ddededed
ed   dee   defdZ	 	 	 ddeded
ed   dee   def
dZy)    N)Optional)Tensor)Literal)#_multiclass_confusion_matrix_update)_compute_bias_corrected_values_compute_chi_squared_drop_empty_rows_and_cols_handle_nan_in_data_nominal_input_validation&_unable_to_use_bias_correction_warningpredstargetnum_classesnan_strategy)replacedropnan_replace_valuereturnc                     | j                   dk(  r| j                  d      n| } |j                   dk(  r|j                  d      n|}t        | |||      \  } }t        | ||      S )a  Compute the bins to update the confusion matrix with for Tschuprow's T calculation.

    Args:
        preds: 1D or 2D tensor of categorical (nominal) data
        target: 1D or 2D tensor of categorical (nominal) data
        num_classes: Integer specifying the number of classes
        nan_strategy: Indication of whether to replace or drop ``NaN`` values
        nan_replace_value: Value to replace ``NaN`s when ``nan_strategy = 'replace```

    Returns:
        Non-reduced confusion matrix

          )ndimargmaxr
   r   )r   r   r   r   r   s        /Volumes/fast/ai/experiments/voice-extract-mac/.venv/lib/python3.12/site-packages/torchmetrics/functional/nominal/tschuprows.py_tschuprows_t_updater       s\    (  %zzQELLOEE!'!1V]]1vF'v|EVWME6.ufkJJ    confmatbias_correctionc                    t        |       } | j                         }t        | |      }||z  }| j                  \  }}|rt	        ||||      \  }}}	t        j                  ||	      dk(  r6t        d       t        j                  t        d      | j                        S t        j                  |t        j                  |dz
  |	dz
  z        z        }
nvt        j                  ||j                        }t        j                  ||j                        }t        j                  |t        j                  |dz
  |dz
  z        z        }
|
j                  dd      S )a  Compute Tschuprow's T statistic based on a pre-computed confusion matrix.

    Args:
        confmat: Confusion matrix for observed data
        bias_correction: Indication of whether to use bias correction.

    Returns:
        Tschuprow's T statistic

    r   zTschuprow's T)metric_namenandevice        g      ?)r	   sumr   shaper   torchminr   tensorfloatr#   sqrtclamp)r   r   cm_sumchi_squaredphi_squarednum_rowsnum_colsphi_squared_correctedrows_correctedcols_correctedtschuprows_t_valuen_rows_tensorn_cols_tensors                r   _tschuprows_t_computer8   :   s0    (0G[[]F&w@K&K Hh@^8VA
=~~ 99^^492O<<eW^^DD"ZZ(=

N]^L^cqtucuKv@w(wxXk6H6HIXk6H6HI"ZZejj-RSBSXehiXiAj6k(kl##C--r   c                     t        ||       t        t        j                  | |g      j	                               }t        | ||||      }t        ||      S )a;  Compute `Tschuprow's T`_ statistic measuring the association between two categorical (nominal) data series.

    .. math::
        T = \sqrt{\frac{\chi^2 / n}{\sqrt{(r - 1) * (k - 1)}}}

    where

    .. math::
        \chi^2 = \sum_{i,j} \ frac{\left(n_{ij} - \frac{n_{i.} n_{.j}}{n}\right)^2}{\frac{n_{i.} n_{.j}}{n}}

    where :math:`n_{ij}` denotes the number of times the values :math:`(A_i, B_j)` are observed with :math:`A_i, B_j`
    represent frequencies of values in ``preds`` and ``target``, respectively.

    Tschuprow's T is a symmetric coefficient, i.e. :math:`T(preds, target) = T(target, preds)`.

    The output values lies in [0, 1] with 1 meaning the perfect association.

    Args:
        preds: 1D or 2D tensor of categorical (nominal) data:

            - 1D shape: (batch_size,)
            - 2D shape: (batch_size, num_classes)

        target: 1D or 2D tensor of categorical (nominal) data:

            - 1D shape: (batch_size,)
            - 2D shape: (batch_size, num_classes)

        bias_correction: Indication of whether to use bias correction.
        nan_strategy: Indication of whether to replace or drop ``NaN`` values
        nan_replace_value: Value to replace ``NaN``s when ``nan_strategy = 'replace'``

    Returns:
        Tschuprow's T statistic

    Example:
        >>> from torch import randint, round
        >>> from torchmetrics.functional.nominal import tschuprows_t
        >>> preds = randint(0, 4, (100,))
        >>> target = round(preds + torch.randn(100)).clamp(0, 4)
        >>> tschuprows_t(preds, target)
        tensor(0.4930)

    )r   lenr'   catuniquer   r8   )r   r   r   r   r   r   r   s          r   tschuprows_tr=   Z   sP    f l,=>eii0779:K"5&+|M^_G /::r   matrixc                    t        ||       | j                  d   }t        j                  ||| j                        }t        j                  t        |      d      D ]l  \  }}| dd|f   | dd|f   }	}t        t        j                  ||	g      j                               }
t        ||	|
||      }t        ||      x|||f<   |||f<   n |S )a  Compute `Tschuprow's T`_ statistic between a set of multiple variables.

    This can serve as a convenient tool to compute Tschuprow's T statistic for analyses of correlation between
    categorical variables in your dataset.

    Args:
        matrix: A tensor of categorical (nominal) data, where:

            - rows represent a number of data points
            - columns represent a number of categorical (nominal) features

        bias_correction: Indication of whether to use bias correction.
        nan_strategy: Indication of whether to replace or drop ``NaN`` values
        nan_replace_value: Value to replace ``NaN``s when ``nan_strategy = 'replace'``

    Returns:
        Tschuprow's T statistic for a dataset of categorical variables

    Example:
        >>> from torch import randint
        >>> from torchmetrics.functional.nominal import tschuprows_t_matrix
        >>> matrix = randint(0, 4, (200, 5))
        >>> tschuprows_t_matrix(matrix)
        tensor([[1.0000, 0.0637, 0.0000, 0.0542, 0.1337],
                [0.0637, 1.0000, 0.0000, 0.0000, 0.0000],
                [0.0000, 0.0000, 1.0000, 0.0000, 0.0649],
                [0.0542, 0.0000, 0.0000, 1.0000, 0.1100],
                [0.1337, 0.0000, 0.0649, 0.1100, 1.0000]])

    r   r"   r   N)r   r&   r'   onesr#   	itertoolscombinationsranger:   r;   r<   r   r8   )r>   r   r   r   num_variablestschuprows_t_matrix_valueijxyr   r   s               r   tschuprows_t_matrixrJ      s    H l,=>LLOM %

=-PVP]P] ^&&u]';Q?1ad|VAqD\1%))QF+2245&q![,HYZLa_M
 	
!!Q$'*CAqD*I	 @ %$r   )r   r$   )Tr   r$   )rA   typingr   r'   r   typing_extensionsr   7torchmetrics.functional.classification.confusion_matrixr   %torchmetrics.functional.nominal.utilsr   r   r	   r
   r   r   intr*   r   boolr8   r=   rJ    r   r   <module>rR      s/       % g  09),KKK K +,	K
  K K4.6 .D .V .F !/8),6;6;6; 6; +,	6;
  6; 6;v !/8),	.%.%.% +,.%  	.%
 .%r   