ó
ù`]c           @  s=  d  Z  d d l m Z m Z 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 d d l m Z d d	 l m Z d d
 l m Z e j e ƒ Z e d e f d „  ƒ  Yƒ Z e d k r9d d l Z d d l Z d d l m Z e j ƒ  Z e j  d d d d d d d d ƒe j  d d d d d d d d ƒe j  d d d d d d  d g  d d! ƒe j  d" d# d d$ d d  d g  d d% ƒe j  d& d' d d( d d) ƒe j  d* d+ d d, d d- d d. ƒe j  d/ d0 d d1 d d2 d d3 ƒe j  d4 d5 d d6 d d7 d d8 ƒe j  d9 d: d d; d e! d d< ƒe j  d= d> d d? d e" d d@ d dA ƒe j# ƒ  \ Z$ Z% e$ j& rŸe j' ƒ  e( ƒ  n  e j) dB dC dD e* e$ j+ ƒ ƒ g  Z, e- e$ j& dE ƒ â Z. xØ e. D]Ð Z/ e/ j0 e$ j1 ƒ Z2 e2 dF e3 e2 dG d !ƒ e4 e2 d j5 ƒ  j0 e$ j6 ƒ ƒ f \ Z7 Z8 Z9 e$ j: e$ j; k p–e< e$ j: ƒ dF k ore7 e$ j: k p–e< e$ j; ƒ dF k o–e7 e$ j; k r²e, j= e7 e8 e9 f ƒ n  qâWWd QXe$ j> rìe e, e? e e$ j ƒ e$ j> ƒ ƒ Z@ n e e, e? e e$ j ƒ ƒ Z@ e$ jA rn eB e? e@ e$ jC ƒ ƒ  ƒ e jD ƒ  n  d S(H   uN	  
Implementations of inter-annotator agreement coefficients surveyed by Artstein
and Poesio (2007), Inter-Coder Agreement for Computational Linguistics.

An agreement coefficient calculates the amount that annotators agreed on label
assignments beyond what is expected by chance.

In defining the AnnotationTask class, we use naming conventions similar to the
paper's terminology.  There are three types of objects in an annotation task:

    the coders (variables "c" and "C")
    the items to be annotated (variables "i" and "I")
    the potential categories to be assigned (variables "k" and "K")

Additionally, it is often the case that we don't want to treat two different
labels as complete disagreement, and so the AnnotationTask constructor can also
take a distance metric as a final argument.  Distance metrics are simply
functions that take two arguments, and return a value between 0.0 and 1.0
indicating the distance between them.  If not supplied, the default is binary
comparison between the arguments.

The simplest way to initialize an AnnotationTask is with a list of triples,
each containing a coder's assignment for one object in the task:

    task = AnnotationTask(data=[('c1', '1', 'v1'),('c2', '1', 'v1'),...])

Note that the data list needs to contain the same number of triples for each
individual coder, containing category values for the same set of items.

Alpha (Krippendorff 1980)
Kappa (Cohen 1960)
S (Bennet, Albert and Goldstein 1954)
Pi (Scott 1955)


TODO: Describe handling of multiple coders and missing data

Expected results from the Artstein and Poesio survey paper:

    >>> from nltk.metrics.agreement import AnnotationTask
    >>> import os.path
    >>> t = AnnotationTask(data=[x.split() for x in open(os.path.join(os.path.dirname(__file__), "artstein_poesio_example.txt"))])
    >>> t.avg_Ao()
    0.88
    >>> t.pi()
    0.7995322418977615...
    >>> t.S()
    0.8199999999999998...

    This would have returned a wrong value (0.0) in @785fb79 as coders are in
    the wrong order. Subsequently, all values for pi(), S(), and kappa() would
    have been wrong as they are computed with avg_Ao().
    >>> t2 = AnnotationTask(data=[('b','1','stat'),('a','1','stat')])
    >>> t2.avg_Ao()
    1.0

    The following, of course, also works.
    >>> t3 = AnnotationTask(data=[('a','1','othr'),('b','1','othr')])
    >>> t3.avg_Ao()
    1.0

iÿÿÿÿ(   t   print_functiont   unicode_literalst   divisionN(   t   groupby(   t
   itemgetter(   t	   iteritems(   t   FreqDistt   ConditionalFreqDist(   t
   deprecated(   t   python_2_unicode_compatible(   t   binary_distancet   AnnotationTaskc           B  s  e  Z d  Z d e d „ Z d „  Z d „  Z d d „ Z d „  Z	 d „  Z
 d „  Z e d ƒ d d d d	 „ ƒ Z d d
 „ Z d „  Z d „  Z d „  Z d d „ Z d d „ Z d „  Z d „  Z d „  Z d „  Z d „  Z d „  Z d „  Z d „  Z d d „ Z d d „ Z RS(   u/  Represents an annotation task, i.e. people assign labels to items.

    Notation tries to match notation in Artstein and Poesio (2007).

    In general, coders and items can be represented as any hashable object.
    Integers, for example, are fine, though strings are more readable.
    Labels must support the distance functions applied to them, so e.g.
    a string-edit-distance makes no sense if your labels are integers,
    whereas interval distance needs numeric values.  A notable case of this
    is the MASI metric, which requires Python sets.
    c         C  sV   | |  _  t ƒ  |  _ t ƒ  |  _ t ƒ  |  _ g  |  _ | d k	 rR |  j | ƒ n  d S(   u¯  Initialize an annotation task.

        The data argument can be None (to create an empty annotation task) or a sequence of 3-tuples,
        each representing a coder's labeling of an item:
            (coder,item,label)

        The distance argument is a function taking two arguments (labels) and producing a numerical distance.
        The distance from a label to itself should be zero:
            distance(l,l) = 0
        N(   t   distancet   sett   It   Kt   Ct   datat   Nonet
   load_array(   t   selfR   R   (    (    s5   lib/python2.7/site-packages/nltk/metrics/agreement.pyt   __init__f   s    		c         C  s   d j  t d „  |  j ƒ ƒ S(   Nu   
c         S  s2   d |  d |  d j  d d ƒ d j |  d ƒ f S(   Nu   %s	%s	%su   coderu   itemu   _u   	u   ,u   labels(   t   replacet   join(   t   x(    (    s5   lib/python2.7/site-packages/nltk/metrics/agreement.pyt   <lambda>|   s   (   R   t   mapR   (   R   (    (    s5   lib/python2.7/site-packages/nltk/metrics/agreement.pyt   __str__y   s    c         C  ss   xl | D]d \ } } } |  j  j | ƒ |  j j | ƒ |  j j | ƒ |  j j i | d 6| d 6| d 6ƒ q Wd S(   u×   Load an sequence of annotation results, appending to any data already loaded.

        The argument is a sequence of 3-tuples, each representing a coder's labeling of an item:
            (coder,item,label)
        u   coderu   labelsu   itemN(   R   t   addR   R   R   t   append(   R   t   arrayt   codert   itemt   labels(    (    s5   lib/python2.7/site-packages/nltk/metrics/agreement.pyR   ‚   s
    c           så   | p |  j  } t ‡  ‡ ‡ f d †  | Dƒ ƒ } | d ˆ  k rc t ‡ ‡ f d †  | Dƒ ƒ } n t ‡  ‡ f d †  | Dƒ ƒ } d t |  j | d | d ƒ ƒ } t j d ˆ  ˆ ˆ | ƒ t j d | d | d d | ƒ | S(	   u6   Agreement between two coders on a given item

        c         3  s;   |  ]1 } | d  ˆ  ˆ f k r | d ˆ k r | Vq d S(   u   coderu   itemN(    (   t   .0R   (   t   cAt   cBt   i(    s5   lib/python2.7/site-packages/nltk/metrics/agreement.pys	   <genexpr>–   s    u   coderc         3  s5   |  ]+ } | d  ˆ  k r | d ˆ k r | Vq d S(   u   coderu   itemN(    (   R"   R   (   R$   R%   (    s5   lib/python2.7/site-packages/nltk/metrics/agreement.pys	   <genexpr>˜   s    c         3  s5   |  ]+ } | d  ˆ  k r | d ˆ k r | Vq d S(   u   coderu   itemN(    (   R"   R   (   R#   R%   (    s5   lib/python2.7/site-packages/nltk/metrics/agreement.pys	   <genexpr>š   s    g      ð?u   labelsu.   Observed agreement between %s and %s on %s: %fu"   Distance between "%r" and "%r": %f(   R   t   nextt   floatR   t   logt   debug(   R   R#   R$   R%   R   t   k1t   k2t   ret(    (   R#   R$   R%   s5   lib/python2.7/site-packages/nltk/metrics/agreement.pyt   agrŽ   s    ""$c           s#   t  t ‡  f d †  |  j Dƒ ƒ ƒ S(   Nc         3  s%   |  ] } | d  ˆ  k r d Vq d S(   u   labelsi   N(    (   R"   R   (   t   k(    s5   lib/python2.7/site-packages/nltk/metrics/agreement.pys	   <genexpr>§   s    (   R'   t   sumR   (   R   R.   (    (   R.   s5   lib/python2.7/site-packages/nltk/metrics/agreement.pyt   Nk¦   s    c           s&   t  t ‡  ‡ f d †  |  j Dƒ ƒ ƒ S(   Nc         3  s5   |  ]+ } | d  ˆ  k r | d ˆ k r d Vq d S(   u   itemu   labelsi   N(    (   R"   R   (   R%   R.   (    s5   lib/python2.7/site-packages/nltk/metrics/agreement.pys	   <genexpr>ª   s    (   R'   R/   R   (   R   R%   R.   (    (   R%   R.   s5   lib/python2.7/site-packages/nltk/metrics/agreement.pyt   Nik©   s    c           s&   t  t ‡  ‡ f d †  |  j Dƒ ƒ ƒ S(   Nc         3  s5   |  ]+ } | d  ˆ  k r | d ˆ k r d Vq d S(   u   coderu   labelsi   N(    (   R"   R   (   t   cR.   (    s5   lib/python2.7/site-packages/nltk/metrics/agreement.pys	   <genexpr>­   s    (   R'   R/   R   (   R   R2   R.   (    (   R2   R.   s5   lib/python2.7/site-packages/nltk/metrics/agreement.pyt   Nck¬   s    u   Use Nk, Nik or Nck insteadc         C  sÞ   | d k	 r6 | d k r6 | d k r6 |  j | ƒ } n‹ | d k	 ro | d k	 ro | d k ro |  j | | ƒ } nR | d k	 r¨ | d k	 r¨ | d k r¨ |  j | | ƒ } n t d | | | f ƒ ‚ t j d | | | | ƒ | S(   uH   Implements the "n-notation" used in Artstein and Poesio (2007)

        u7   You must pass either i or c, not both! (k=%r,i=%r,c=%r)u   Count on N[%s,%s,%s]: %dN(   R   R0   R1   R3   t
   ValueErrorR(   R)   (   R   R.   R%   R2   R,   (    (    s5   lib/python2.7/site-packages/nltk/metrics/agreement.pyt   N¯   s    $$$c         C  s4   | p |  j  } t t | d t | ƒ ƒt | ƒ ƒ S(   Nt   key(   R   R   t   sortedR   (   R   t   fieldR   (    (    s5   lib/python2.7/site-packages/nltk/metrics/agreement.pyt   _grouped_dataÁ   s    c           sq   ˆ j  d ‡  ‡ f d †  ˆ j Dƒ ƒ } t ‡  ‡ ‡ f d †  | Dƒ ƒ t ˆ j ƒ } t j d ˆ  ˆ | ƒ | S(   u=   Observed agreement between two coders on all items.

        u   itemc         3  s+   |  ]! } | d  ˆ  ˆ f k r | Vq d S(   u   coderN(    (   R"   R   (   R#   R$   (    s5   lib/python2.7/site-packages/nltk/metrics/agreement.pys	   <genexpr>Ê   s    c         3  s-   |  ]# \ } } ˆ j  ˆ  ˆ | | ƒ Vq d  S(   N(   R-   (   R"   R    t	   item_data(   R#   R$   R   (    s5   lib/python2.7/site-packages/nltk/metrics/agreement.pys	   <genexpr>Ì   s    u(   Observed agreement between %s and %s: %f(   R9   R   R/   t   lenR   R(   R)   (   R   R#   R$   R   R,   (    (   R#   R$   R   s5   lib/python2.7/site-packages/nltk/metrics/agreement.pyt   AoÅ   s    ""c         C  sx   d } d } |  j  j ƒ  } xL |  j  D]A } | j | ƒ x+ | D]# } | | | | ƒ 7} | d 7} q? Wq% W| | } | S(   uP   
        Calculates the average of function results for each coder pair
        i    i   (   R   t   copyt   remove(   R   t   functiont   totalt   nt   sR#   R$   R,   (    (    s5   lib/python2.7/site-packages/nltk/metrics/agreement.pyt   _pairwise_averageÒ   s    
c         C  s&   |  j  |  j ƒ } t j d | ƒ | S(   uA   Average observed agreement across all coders and items.

        u   Average observed agreement: %f(   RC   R<   R(   R)   (   R   R,   (    (    s5   lib/python2.7/site-packages/nltk/metrics/agreement.pyt   avg_Aoá   s    g      ð?c   	        s    d } ‡  ‡ f d †  |  j  Dƒ } xJ |  j d | ƒ D]6 \ } } | |  j t | ƒ d t | ƒ d ƒ 7} q5 W| t |  j ƒ | } t j d ˆ  ˆ | ƒ | S(   uG   The observed disagreement for the weighted kappa coefficient.

        g        c         3  s+   |  ]! } | d  ˆ  ˆ f k r | Vq d S(   u   coderN(    (   R"   R   (   R#   R$   (    s5   lib/python2.7/site-packages/nltk/metrics/agreement.pys	   <genexpr>î   s    u   itemu   labelsu+   Observed disagreement between %s and %s: %f(   R   R9   R   R&   R;   R   R(   R)   (	   R   R#   R$   t   max_distanceR@   R   R%   t   itemdataR,   (    (   R#   R$   s5   lib/python2.7/site-packages/nltk/metrics/agreement.pyt   Do_Kw_pairwiseé   s    .c           s/   ˆ j  ‡  ‡ f d †  ƒ } t j d | ƒ | S(   u$   Averaged over all labelers

        c           s   ˆ j  |  | ˆ  ƒ S(   N(   RG   (   R#   R$   (   RE   R   (    s5   lib/python2.7/site-packages/nltk/metrics/agreement.pyR   ü   t    u   Observed disagreement: %f(   RC   R(   R)   (   R   RE   R,   (    (   RE   R   s5   lib/python2.7/site-packages/nltk/metrics/agreement.pyt   Do_Kw÷   s    c         C  s/   d t  |  j ƒ } |  j ƒ  | d | } | S(   u,   Bennett, Albert and Goldstein 1954

        g      ð?(   R;   R   RD   (   R   t   AeR,   (    (    s5   lib/python2.7/site-packages/nltk/metrics/agreement.pyt   S  s    c         C  s„   d } t  d „  |  j Dƒ ƒ } x( t | ƒ D] \ } } | | d 7} q, W| t |  j ƒ t |  j ƒ d } |  j ƒ  | d | S(   u_   Scott 1955; here, multi-pi.
        Equivalent to K from Siegel and Castellan (1988).

        g        c         s  s   |  ] } | d  Vq d S(   u   labelsN(    (   R"   R   (    (    s5   lib/python2.7/site-packages/nltk/metrics/agreement.pys	   <genexpr>  s    i   i   (   R   R   R   R;   R   R   RD   (   R   R@   t   label_freqsR.   t   fRJ   (    (    s5   lib/python2.7/site-packages/nltk/metrics/agreement.pyt   pi
  s    $c         C  su   d } t  t |  j ƒ ƒ } t d „  |  j Dƒ ƒ } x: | j ƒ  D], } | | | | | | | | | 7} qA W| S(   Ng        c         s  s#   |  ] } | d  | d f Vq d S(   u   labelsu   coderN(    (   R"   R   (    (    s5   lib/python2.7/site-packages/nltk/metrics/agreement.pys	   <genexpr>  s    (   R'   R;   R   R   R   t
   conditions(   R   R#   R$   RJ   t   nitemsRL   R.   (    (    s5   lib/python2.7/site-packages/nltk/metrics/agreement.pyt   Ae_kappa  s    *c         C  sJ   |  j  | | ƒ } |  j | | ƒ | d | } t j d | | | ƒ | S(   u
   

        g      ð?u(   Expected agreement between %s and %s: %f(   RQ   R<   R(   R)   (   R   R#   R$   RJ   R,   (    (    s5   lib/python2.7/site-packages/nltk/metrics/agreement.pyt   kappa_pairwise  s    c         C  s   |  j  |  j ƒ S(   uN   Cohen 1960
        Averages naively over kappas for each coder pair.

        (   RC   RR   (   R   (    (    s5   lib/python2.7/site-packages/nltk/metrics/agreement.pyt   kappa'  s    c         C  s(   |  j  |  j ƒ } |  j ƒ  | d | S(   ul   Davies and Fleiss 1982
        Averages over observed and expected agreements for each coder pair.

        g      ð?(   RC   RQ   RD   (   R   RJ   (    (    s5   lib/python2.7/site-packages/nltk/metrics/agreement.pyt   multi_kappa.  s    c         C  sŠ   t  | j ƒ  ƒ } d } x[ t | ƒ D]M \ } } x> t | ƒ D]0 \ } } | t | | ƒ |  j | | ƒ 7} q> Wq% Wd | | | d S(   Ng        g      ð?i   (   R/   t   valuesR   R'   R   (   R   RL   t   total_labelst   pairst   jt   njt   lt   nl(    (    s5   lib/python2.7/site-packages/nltk/metrics/agreement.pyt   Disagreement6  s    ,c         C  sS  t  |  j ƒ d k r$ t d ƒ ‚ n  t  |  j ƒ d k rJ t j d ƒ d St  |  j ƒ d k rƒ t  |  j ƒ d k rƒ t d ƒ ‚ n  d } d } t g  ƒ } d } xx |  j d ƒ D]g \ } } t d „  | Dƒ ƒ } t	 | j
 ƒ  ƒ } | d	 k  r÷ q± n  | | 7} | |  j | ƒ | 7} q± W| t	 | j
 ƒ  ƒ }	 |  j | ƒ }
 d
 |	 |
 } | S(   u   Krippendorff 1980

        i    u(   Cannot calculate alpha, no data present!i   u.   Only one annotation value, allpha returning 1.u8   Cannot calculate alpha, only one coder and item present!g        u   itemc         s  s   |  ] } | d  Vq d S(   u   labelsN(    (   R"   R   (    (    s5   lib/python2.7/site-packages/nltk/metrics/agreement.pys	   <genexpr>Q  s    i   g      ð?(   R;   R   R4   R(   R)   R   R   R   R9   R/   RU   R\   (   R   t   total_disagreementt   total_ratingst   all_valid_labels_freqt   total_doR%   RF   RL   t   labels_countt   dot   det   k_alpha(    (    s5   lib/python2.7/site-packages/nltk/metrics/agreement.pyt   alpha>  s,    *
c           sØ   d } t  ‡  ‡ f d †  |  j Dƒ ƒ } xS |  j D]H } x? |  j D]4 } | | ˆ  | | ˆ | |  j | | ƒ 7} qB Wq2 W| | t t |  j ƒ d ƒ } t j d ˆ  ˆ | ƒ |  j	 ˆ  ˆ ƒ }	 d |	 | }
 |
 S(   u   Cohen 1968

        g        c         3  s9   |  ]/ } | d  ˆ  ˆ f k r | d  | d f Vq d S(   u   coderu   labelsN(    (   R"   R   (   R#   R$   (    s5   lib/python2.7/site-packages/nltk/metrics/agreement.pys	   <genexpr>f  s    i   u+   Expected disagreement between %s and %s: %fg      ð?(
   R   R   R   R   t   powR;   R   R(   R)   RG   (   R   R#   R$   RE   R@   RL   RX   RZ   t   Det   DoR,   (    (   R#   R$   s5   lib/python2.7/site-packages/nltk/metrics/agreement.pyt   weighted_kappa_pairwise`  s    6 c           s   ˆ j  ‡  ‡ f d †  ƒ S(   u   Cohen 1968

        c           s   ˆ j  |  | ˆ  ƒ S(   N(   Ri   (   R#   R$   (   RE   R   (    s5   lib/python2.7/site-packages/nltk/metrics/agreement.pyR   v  RH   (   RC   (   R   RE   (    (   RE   R   s5   lib/python2.7/site-packages/nltk/metrics/agreement.pyt   weighted_kappaq  s    N(   t   __name__t
   __module__t   __doc__R   R
   R   R   R   R-   R0   R1   R3   R   R5   R9   R<   RC   RD   RG   RI   RK   RN   RQ   RR   RS   RT   R\   Re   Ri   Rj   (    (    (    s5   lib/python2.7/site-packages/nltk/metrics/agreement.pyR   X   s4   																			"u   __main__(   R   u   -du
   --distancet   destu   distancet   defaultu   binary_distancet   helpu   distance metric to useu   -au   --agreementu	   agreementu   kappau"   agreement coefficient to calculateu   -eu	   --excludeu   excludet   actionu   appendu8   coder names to exclude (may be specified multiple times)u   -iu	   --includeu   includeu.   coder names to include, same format as excludeu   -fu   --fileu   fileuP   file to read labelings from, each line with three columns: 'labeler item labels'u   -vu	   --verboseu   verboseu   0u+   how much debugging to print on stderr (0-4)u   -cu   --columnsepu	   columnsepu   	uI   char/string that separates the three columns in the file, defaults to tabu   -lu
   --labelsepu   labelsepu   ,u[   char/string that separates labels (if labelers can assign more than one), defaults to commau   -pu
   --presenceu   presenceu=   convert each labeling into 1 or 0, based on presence of LABELu   -Tu
   --thoroughu   thoroughu
   store_trueu6   calculate agreement for every subset of the annotatorst   leveli2   i
   u   ri    i   (E   Rm   t
   __future__R    R   R   t   loggingt	   itertoolsR   t   operatorR   t   sixR   t   nltk.probabilityR   R   t   nltk.internalsR   t   nltk.compatR	   t   nltk.metrics.distanceR
   t	   getLoggerRk   R(   t   objectR   t   ret   optparset   nltk.metricsR   t   OptionParsert   parsert
   add_optionR   t   Falset
   parse_argst   optionst	   remaindert   filet
   print_helpt   exitt   basicConfigt   intt   verboseR   t   opent   infileRZ   t   splitt	   columnsept   tokst   strt	   frozensett   stript   labelsepR   t   object_R!   t   includet   excludeR;   R   t   presencet   getattrt   taskt   thorought   printt	   agreementt   shutdown(    (    (    s5   lib/python2.7/site-packages/nltk/metrics/agreement.pyt   <module>F   sÐ   ÿ "


!.$$#	$	