Krippendorff's alpha

Reliability data are generated in a situation in which m ≥ 2 jointly instructed (e.g., by a code book) but independently working coders assign any one of a set of values 1,...,V to a common set of N units of analysis. In their canonical form, reliability data are tabulated in an m-by-N matrix containing N values v_ij that coder c_i has assigned to unit u_j. Define m_j as the number of values assigned to unit j across all coders c. When data are incomplete, m_j may be less than m. Reliability data require that values be pairable, i.e., m_j ≥ 2. The total number of pairable values is ${\displaystyle \sum _{j=1}^{N}m_{j}=}$ n ≤ mN.

To help clarify, here is what the canonical form looks like, in the abstract:

We denote by ${\displaystyle R}$ the set of all possible responses an observer can give. The responses of all observers for an example is called a unit (it forms a multiset). We denote a multiset with these units as the items, ${\displaystyle U}$.

Alpha is given by:

${\displaystyle \alpha =1-{\frac {D_{o}}{D_{e}}}}$

where ${\displaystyle D_{o}}$ is the disagreement observed and ${\displaystyle D_{e}}$ is the disagreement expected by chance.

${\displaystyle D_{o}={\frac {1}{n}}\sum _{c\in R}\sum _{k\in R}\delta (c,k)\sum _{u\in U}m_{u}{\frac {n_{cku}}{P(m_{u},2)}}}$

where ${\displaystyle \delta }$ is a metric function (note that this is not a metric in the mathematical sense, but often the square of a mathematical metric, see below), ${\displaystyle n}$ is the total number of pairable elements, ${\displaystyle m_{u}}$ is the number of items in a unit, ${\displaystyle n_{cku}}$ number of ${\displaystyle (c,k)}$ pairs in unit ${\displaystyle u}$, and ${\displaystyle P}$ is the permutation function. Rearranging terms, the sum can be interpreted in a conceptual way as the weighted average of the disagreements of the individual units---weighted by the number of coders assigned to unit j:

${\displaystyle D_{o}={\frac {1}{n}}\sum _{j=1}^{N}m_{j}\,\mathbb {E} (\delta _{j})}$

where ${\displaystyle \mathbb {E} (\delta _{j})}$ is the mean of the ${\displaystyle m_{j} \choose 2}$ numbers ${\displaystyle \delta (v_{ij},v_{i'j})}$ (here ${\displaystyle i>i'}$ and define pairable elements). Note that in the case ${\displaystyle m_{j}=m}$ for all ${\displaystyle j}$, ${\displaystyle D_{o}}$ is just the average all the numbers ${\displaystyle \delta (v_{ij},v_{i'j})}$ with ${\displaystyle i>i'}$. There is also an interpretation of ${\displaystyle D_{o}}$as the (weighted) average observed distance from the diagonal.

${\displaystyle D_{e}={\frac {1}{P(n,2)}}\sum _{c\in R}\sum _{k\in R}\delta (c,k)P_{ck}}$

where ${\displaystyle P_{ck}}$ is the number of ways the pair ${\displaystyle (c,k)}$ can be made. This can be seen to be the average distance from the diagonal of all possible pairs of responses that could be derived from the multiset of all observations.

${\displaystyle P_{ck}={\begin{cases}c\neq k&n_{c}n_{k}\\c=k&n_{c}(n_{c}-1)\end{cases}}}$

The above is equivalent to the usual form of ${\displaystyle \alpha }$ once it has been simplified algebraically.[10]

One interpretation of Krippendorff's alpha is: ${\displaystyle \alpha =1-{\frac {D_{{\text{within units}}={\text{in error}}}}{D_{{\text{within and between units}}={\text{in total}}}}}}$

${\displaystyle \alpha =1}$ indicates perfect reliability

${\displaystyle \alpha =0}$ indicates the complete absence of reliability. Units and the values assigned to them are statistically unrelated.

${\displaystyle \alpha <0}$ when disagreements are systematic and exceed what can be expected by chance.

In this general form, disagreements D_o and D_e may be conceptually transparent but are computationally inefficient. They can be simplified algebraically, especially when expressed in terms of the visually more instructive coincidence matrix representation of the reliability data.

A coincidence matrix cross tabulates the n pairable values from the canonical form of the reliability data into a v-by-v square matrix, where v is the number of values available in a variable. Unlike contingency matrices, familiar in association and correlation statistics, which tabulate pairs of values (cross tabulation), a coincidence matrix tabulates all pairable values. A coincidence matrix omits references to coders and is symmetrical around its diagonal, which contains all perfect matches, v_iu = v_i'u for two coders i and i' , across all units u. The matrix of observed coincidences contains frequencies:

${\displaystyle {\begin{aligned}o_{vv'}&=\sum _{u=1}^{N}{\frac {\sum _{i\neq i'}^{m}I(v_{iu}=v)\cdot I(v_{i'u}=v')}{m_{u}-1}}=o_{v'v},\\[5pt]n_{v}&=\sum _{\ell =1}^{V}o_{v\ell }=\sum _{v_{ij}}^{m,N}I(v_{ij}=v){\text{ and }}n=\sum _{\ell =1,p=1}^{V}o_{\ell p},\end{aligned}}}$

omitting unpaired values, where I(∘) = 1 if ∘ is true, and 0 otherwise.

Because a coincidence matrix tabulates all pairable values and its contents sum to the total n, when four or more coders are involved, o_ck may be fractions.

The matrix of expected coincidences contains frequencies:

${\displaystyle \ e_{vv'}={\frac {\sum _{i\neq i'}^{m}I(v_{iu}=v)\cdot I(v_{i'u}=v')}{n-1}}={\frac {1}{n-1}}\cdot \left.{\begin{cases}n_{v}(n_{v}-1)&{\text{if }}v=v'\\n_{v}n_{v'}&{\text{if }}v\neq v'\end{cases}}\right\}=e_{kc},}$

which sum to the same n_c, n_k, and n as does o_ck. In terms of these coincidences, Krippendorff's alpha becomes:

${\displaystyle \alpha =1-{\frac {D_{o}}{D_{e}}}=1-{\frac {\sum _{v=1,v'=1}^{V}o_{vv'}\delta (v,v')}{\sum _{v=1,v'=1}^{V}e_{vv'}\delta (v,v')}}.}$

Difference functions ${\displaystyle \delta (v,v')}$[11] between values v and v' reflect the metric properties (levels of measurement) of their variable.

In general:

${\displaystyle {\begin{aligned}\delta (v,v')&\geq 0\\[4pt]\delta (v,v)&=0\\[4pt]\delta (v,v')&=\delta (v',v)\end{aligned}}}$

In particular:

For nominal data ${\displaystyle \delta _{\text{nominal}}(v,v')={\begin{cases}0&{\text{if }}v=v'\\1&{\text{if }}v\neq v'\end{cases}}}$, where v and v' serve as names.

For ordinal data ${\displaystyle \delta _{\text{ordinal}}(v,v')=\left(\sum _{g=v}^{g=v'}n_{g}-{\frac {n_{v}+n_{v'}}{2}}\right)^{2}}$, where v and v′ are ranks.

For interval data ${\displaystyle \delta _{\text{interval}}(v,v')=(v-v')^{2}}$, where v and v′ are interval scale values.

For ratio data ${\displaystyle \delta _{\text{ratio}}(v,v')=\left({\frac {v-v'}{v+v'}}\right)^{2}}$, where v and v′ are absolute values.

For polar data ${\displaystyle \delta _{\text{polar}}(v,v')={\frac {(v-v')^{2}}{(v+v'-2v_{\min })(2v_{\max }-v-v')}}}$, where v_min and v_max define the end points of the polar scale.

For circular data ${\displaystyle \delta _{\text{circular}}(v,v')=\left(\sin \left[180{\frac {v-v'}{U}}\right]\right)^{2}}$, where the sine function is expressed in degrees and U is the circumference or the range of values in a circle or loop before they repeat. For equal interval circular metrics, the smallest and largest integer values of this metric are adjacent to each other and U = v_largest – v_smallest + 1.

Inasmuch as mathematical statements of the statistical distribution of alpha are always only approximations, it is preferable to obtain alpha’s distribution by bootstrapping.[12][13] Alpha's distribution gives rise to two indices:

• The confidence intervals of a computed alpha at various levels of statistical significance
• The probability that alpha fails to achieve a chosen minimum, required for data to be considered sufficiently reliable (one-tailed test). This index acknowledges that the null-hypothesis (of chance agreement) is so far removed from the range of relevant alpha coefficients that its rejection would mean little regarding how reliable given data are. To be judged reliable, data must not significantly deviate from perfect agreement.

The minimum acceptable alpha coefficient should be chosen according to the importance of the conclusions to be drawn from imperfect data. When the costs of mistaken conclusions are high, the minimum alpha needs to be set high as well. In the absence of knowledge of the risks of drawing false conclusions from unreliable data, social scientists commonly rely on data with reliabilities α ≥ 0.800, consider data with 0.800 > α ≥ 0.667 only to draw tentative conclusions, and discard data whose agreement measures α < 0.667.[14]

Let the canonical form of reliability data be a 3-coder-by-15 unit matrix with 45 cells:

Suppose “*” indicates a default category like “cannot code,” “no answer,” or “lacking an observation.” Then, * provides no information about the reliability of data in the four values that matter. Note that unit 2 and 14 contains no information and unit 1 contains only one value, which is not pairable within that unit. Thus, these reliability data consist not of mN = 45 but of n = 26 pairable values, not in N = 15 but in 12 multiply coded units.

The coincidence matrix for these data would be constructed as follows:

o₁₁ = {in u=4}: ${\textstyle {\frac {2}{2-1}}+}$ {in u=10}: ${\textstyle {\frac {2}{2-1}}+}$ {in u=11}: ${\textstyle {\frac {2}{2-1}}=6}$

o₁₃ = {in u=8}: ${\textstyle {\frac {1}{2-1}}=1=}$ o₃₁

o₂₂ = {in u=3}: ${\textstyle {\frac {2}{2-1}}+}$ {in u=9}: ${\textstyle {\frac {2}{2-1}}=4}$

o₃₃ = {in u=5}: ${\textstyle {\frac {2}{2-1}}+}$ {in u=6}: ${\textstyle {\frac {2}{3-1}}+}$ {in u=12}: ${\textstyle {\frac {2}{2-1}}+}$ {in u=13}: ${\textstyle {\frac {2}{2-1}}=7}$

o₃₄ = {in u=6}: ${\textstyle {\frac {2}{3-1}}+}$ {in u=15}: ${\textstyle {\frac {1}{2-1}}=2=}$ o₄₃

o₄₄ = {in u=7}: ${\textstyle {\frac {6}{3-1}}=3}$

In terms of the entries in this coincidence matrix, Krippendorff's alpha may be calculated from:

${\displaystyle \alpha _{\text{metric}}=1-{\frac {D_{o}}{D_{e}}}=1-{\frac {\sum _{v=1,v'=1}^{V}o_{vv'}\delta _{\text{metric}}(v,v')}{{\frac {1}{n-1}}\sum _{v=1,v'=1}^{V}n_{v}n_{v'}~\delta _{\text{metric}}(v,v')}}.}$

For convenience, because products with ${\displaystyle \delta (v,v)=0}$ and ${\displaystyle \delta (v,v')=\delta (v',v)}$, only the entries in one of the off-diagonal triangles of the coincidence matrix are listed in the following:

${\displaystyle \alpha _{\text{metric}}=1-{\frac {1\delta _{\text{metric}}(1,3)+2\delta _{\text{metric}}(3,4)}{{\frac {1}{26-1}}(4\cdot 7\delta _{\text{metric}}(1,2)+10\cdot 7\delta _{\text{metric}}(1,3)+5\cdot 7\delta _{\text{metric}}(1,4)+10\cdot 4\delta _{\text{metric}}(2,3)+5\cdot 4\delta _{\text{metric}}(2,4)+5\cdot 10\delta _{\text{metric}}(3,4))}}}$

Considering that all ${\displaystyle \delta _{\text{nominal}}(v,v')=1}$ when ${\displaystyle v{\neq }v'}$ for nominal data the above expression yields:

${\displaystyle \alpha _{\text{nominal}}=1-{\frac {1+2}{{\frac {1}{26-1}}(4\cdot 7+10\cdot 7+5\cdot 7+10\cdot 4+5\cdot 4+5\cdot 10)}}=0.691}$

With ${\displaystyle \delta _{\text{interval}}(1,2)=\delta _{\text{interval}}(2,3)=\delta _{\text{interval}}(3,4)=1^{2},\qquad \delta _{\text{interval}}(1,3)=\delta _{\text{interval}}(2,4)=2^{2},{\text{ and }}\delta _{\text{interval}}(1,4)=3^{2},}$ for interval data the above expression yields:

${\displaystyle \alpha _{\text{interval}}=1-{\frac {1\cdot 2^{2}+2\cdot 1^{2}}{{\frac {1}{26-1}}(4\cdot 7\cdot 1^{2}+10\cdot 7\cdot 2^{2}+5\cdot 7\cdot 3^{2}+10\cdot 4\cdot 1^{2}+5\cdot 4\cdot 2^{2}+5\cdot 10\cdot 1^{2})}}=0.811}$

Here, ${\displaystyle \alpha _{\text{interval}}>\alpha _{\text{nominal}}}$ because disagreements happens to occur largely among neighboring values, visualized by occurring closer to the diagonal of the coincidence matrix, a condition that ${\displaystyle \alpha _{\text{interval}}}$ takes into account but ${\displaystyle \alpha _{\text{nominal}}}$ does not. When the observed frequencies o_{v ≠ v′} are on the average proportional to the expected frequencies e_{v ≠ v'}, ${\displaystyle \alpha _{\text{interval}}=\alpha _{\text{nominal}}}$.

Comparing alpha coefficients across different metrics can provide clues to how coders conceptualize the metric of a variable.

Krippendorff's alpha brings several known statistics under a common umbrella, each of them has its own limitations but no additional virtues.

• Scott's pi[15] is an agreement coefficient for nominal data and two coders.
  $${\displaystyle \pi ={\frac {P_{o}-P_{e}}{1-P_{e}}}{\text{ where }}P_{o}=\sum _{c}{\frac {o_{cc}}{n}},{\text{ and }}P_{e}=\sum _{c}{\frac {n_{c}^{2}}{n^{2}}}.}$$

  
  When data are nominal, alpha reduces to a form resembling Scott's pi:
  $${\displaystyle _{\text{nominal}}\alpha =1-{\frac {D_{o}}{D_{e}}}={\frac {\sum _{c}o_{cc}-\sum _{c}e_{cc}}{n-\sum _{c}e_{cc}}}={\frac {\sum _{c}{\frac {O_{cc}}{n}}-\sum _{c}{\frac {n_{c}(n_{c}-1)}{n(n-1)}}}{1-\sum _{c}{\frac {n_{c}(n_{c}-1)}{n(n-1)}}}}}$$

  
  Scott's observed proportion of agreement ${\displaystyle \ P_{o}}$ appears in alpha’s numerator, exactly. Scott's expected proportion of agreement, ${\textstyle P_{e}=\sum _{c}{\frac {n_{c}^{2}}{n^{2}}}}$ is asymptotically approximated by ${\textstyle \sum _{c}{\frac {n_{c}(n_{c}-1)}{n(n-1)}}}$ when the sample size n is large, equal when infinite. It follows that Scott's pi is that special case of alpha in which two coders generate a very large sample of nominal data. For finite sample sizes: ${\displaystyle {_{\text{nominal}}\alpha }=1-{\tfrac {n-1}{n}}(1-\pi )\geq \pi }$. Evidently, ${\textstyle \lim _{n\to \infty }{_{\text{nominal}}\alpha }=\pi }$.
• Fleiss’ kappa[16] is an agreement coefficient for nominal data with very large sample sizes where a set of coders have assigned exactly m labels to all of N units without exception (but note, there may be more than m coders, and only some subset label each instance). Fleiss claimed to have extended Cohen's kappa[17] to three or more raters or coders, but generalized Scott's pi instead. This confusion is reflected in Fleiss’ choice of its name, which has been recognized by renaming it K:[18]
  $${\displaystyle K={\frac {{\bar {P}}-{\bar {P}}_{e}}{1-{\bar {P}}_{e}}}{\text{ where }}{\bar {P}}={\frac {1}{N}}\sum _{u=1}^{N}\sum _{c}{\frac {n_{cu}(n_{cu}-1)}{m(m-1)}}=\sum _{c}{\frac {o_{cc}}{mN}},{\text{ and }}{\bar {P}}_{e}=\sum _{c}{\frac {n_{c}^{2}}{(mN)^{2}}}}$$

  
  When sample sizes are finite, K can be seen to perpetrate the inconsistency of obtaining the proportion of observed agreements ${\displaystyle {\bar {P}}}$ by counting matches within the m(m − 1) possible pairs of values within u, properly excluding values paired with themselves, while the proportion ${\displaystyle {\bar {P}}_{e}}$ is obtained by counting matches within all (mN)² = n² possible pairs of values, effectively including values paired with themselves. It is the latter that introduces a bias into the coefficient. However, just as for pi, when sample sizes become very large this bias disappears and the proportion ${\textstyle \sum _{c}{\frac {n_{c}(n_{c}-1)}{n(n-1)}}}$ in _nominalα above asymptotically approximates ${\displaystyle {\bar {P}}_{e}}$ in K. Nevertheless, Fleiss' kappa, or rather K, intersects with alpha in that special situation in which a fixed number of m coders code all of N units (no data are missing), using nominal categories, and the sample size n = mN is very large, theoretically infinite.
• Spearman's rank correlation coefficient rho[19] measures the agreement between two coders’ ranking of the same set of N objects. In its original form:
  $${\displaystyle \rho =1-{\frac {6\sum D^{2}}{N(N^{2}-1)}},}$$

  
  where ${\textstyle \sum D^{2}=\sum _{u=1}^{N}{_{\text{ordinal}}\delta }_{c_{u}k_{u}}^{2}}$ is the sum of N differences between one coder's rank c and the other coder's rank k of the same object u. Whereas alpha accounts for tied ranks in terms of their frequencies for all coders, rho averages them in each individual coder's instance. In the absence of ties, ${\displaystyle \rho }$'s numerator ${\textstyle \sum D^{2}=ND_{o}}$ and ${\displaystyle \rho }$'s denominator ${\textstyle {\frac {N(N^{2}-1)}{6}}={\frac {n}{n-1}}ND_{e}}$, where n = 2N, which becomes ${\displaystyle \ ND_{e}}$ when sample sizes become large. So, Spearman's rho is that special case of alpha in which two coders rank a very large set of units. Again, ${\displaystyle {_{\text{ordinal}}\alpha }\geq \rho }$ and ${\textstyle \lim _{n\to \infty }{_{\text{ordinal}}\alpha }=\rho }$.
• Pearson's intraclass correlation coefficient r_ii is an agreement coefficient for interval data, two coders, and very large sample sizes. To obtain it, Pearson's original suggestion was to enter the observed pairs of values twice into a table, once as c − k and once as k − c, to which the traditional Pearson product-moment correlation coefficient is then applied.[20] By entering pairs of values twice, the resulting table becomes a coincidence matrix without reference to the two coders, contains n = 2N values, and is symmetrical around the diagonal, i.e., the joint linear regression line is forced into a 45° line, and references to coders are eliminated. Hence, Pearson's intraclass correlation coefficient is that special case of interval alpha for two coders and large sample sizes, ${\displaystyle {_{\text{interval}}\alpha }\geq r_{ii}}$ and ${\textstyle \lim _{n\to \infty }{_{\text{interval}}\alpha }=r_{ii}}$.
• Finally, The disagreements in the interval alpha, D_u, D_o and D_e are proper sample variances.[21] It follows that the reliability the interval alpha assesses is consistent with all variance-based analytical techniques, such as the analysis of variance. Moreover, by incorporating difference functions not just for interval data but also for nominal, ordinal, ratio, polar, and circular data, alpha extends the notion of variance to metrics that classical analytical techniques rarely address.

Krippendorff's alpha is more general than any of these special purpose coefficients. It adjusts to varying sample sizes and affords comparisons across a wide variety of reliability data, mostly ignored by the familiar measures.

Semantically, reliability is the ability to rely on something, here on coded data for subsequent analysis. When a sufficiently large number of coders agree perfectly on what they have read or observed, relying on their descriptions is a safe bet. Judgments of this kind hinge on the number of coders duplicating the process and how representative the coded units are of the population of interest. Problems of interpretation arise when agreement is less than perfect, especially when reliability is absent.

• Correlation and association coefficients. Pearson's product-moment correlation coefficient r_ij, for example, measures deviations from any linear regression line between the coordinates of i and j. Unless that regression line happens to be exactly 45° or centered, r_ij does not measure agreement. Similarly, while perfect agreement between coders also means perfect association, association statistics register any above chance pattern of relationships between variables. They do not distinguish agreement from other associations and are, hence, unsuitable as reliability measures.
• Coefficients measuring the degree to which coders are statistically dependent on each other. When the reliability of coded data is at issue, the individuality of coders can have no place in it. Coders need to be treated as interchangeable. Alpha, Scott's pi, and Pearson's original intraclass correlation accomplish this by being definable as a function of coincidences, not only of contingencies. Unlike the more familiar contingency matrices, which tabulate N pairs of values and maintain reference to the two coders, coincidence matrices tabulate the n pairable values used in coding, regardless of who contributed them, in effect treating coders as interchangeable. Cohen's kappa,[22] by contrast, defines expected agreement in terms of contingencies, as the agreement that would be expected if coders were statistically independent of each other.[23] Cohen's conception of chance fails to include disagreements between coders’ individual predilections for particular categories, punishes coders who agree on their use of categories, and rewards those who do not agree with higher kappa-values. This is the cause of other noted oddities of kappa.[24] The statistical independence of coders is only marginally related to the statistical independence of the units coded and the values assigned to them. Cohen's kappa, by ignoring crucial disagreements, can become deceptively large when the reliability of coding data is to be assessed.
• Coefficients measuring the consistency of coder judgments. In the psychometric literature,[25] reliability tends to be defined as the consistency with which several tests perform when applied to a common set of individual characteristics. Cronbach's alpha,[26] for example, is designed to assess the degree to which multiple tests produce correlated results. Perfect agreement is the ideal, of course, but Cronbach's alpha is high also when test results vary systematically. Consistency of coders’ judgments does not provide the needed assurances of data reliability. Any deviation from identical judgments – systematic or random – needs to count as disagreement and reduce the measured reliability. Cronbach's alpha is not designed to respond to absolute differences.
• Coefficients with baselines (conditions under which they measure 0) that cannot be interpreted in terms of reliability, i.e. have no dedicated value to indicate when the units and the values assigned to them are statistically unrelated. Simple %-agreement ranges from 0 = extreme disagreement to 100 = perfect agreement with chance having no definite value. As already noted, Cohen's kappa falls into this category by defining the absence of reliability as the statistical independence between two individual coders. The baseline of Bennett, Alpert, and Goldstein's S[27] is defined in terms of the number of values available for coding, which has little to do with how values are actually used. Goodman and Kruskal's lambda_r[28] is defined to vary between –1 and +1, leaving 0 without a particular reliability interpretation. Lin's reproducibility or concordance coefficient r_c[29] takes Pearson's product moment correlation r_ij as a measure of precision and adds to it a measure C_b of accuracy, ostensively to correct for r_ij's above mentioned inadequacy. It varies between –1 and +1 and the reliability interpretation of 0 is uncertain. There are more so-called reliability measures whose reliability interpretations become questionable as soon as they deviate from perfect agreement.

Naming a statistic as one of agreement, reproducibility, or reliability does not make it a valid index of whether one can rely on coded data in subsequent decisions. Its mathematical structure must fit the process of coding units into a system of analyzable terms.

• K. Krippendorff, 2013, Content Analysis: An Introduction to Its Methodology, 3rd ed. Thousand Oaks, CA, USA: Sage, PP. 221–250

• Bennett, Edward M., Alpert, R. & Goldstein, A. C. (1954). Communications through limited response questioning. Public Opinion Quarterly, 18, 303–308.
• Brennan, Robert L. & Prediger, Dale J. (1981). Coefficient kappa: Some uses, misuses, and alternatives. Educational and Psychological Measurement, 41, 687–699.
• Cohen, Jacob (1960). A coefficient of agreement for nominal scales. Educational and Psychological Measurement, 20 (1), 37–46.
• Cronbach, Lee, J. (1951). Coefficient alpha and the internal structure of tests. Psychometrika, 16 (3), 297–334.
• Fleiss, Joseph L. (1971). Measuring nominal scale agreement among many raters. Psychological Bulletin, 76, 378–382.
• Goodman, Leo A. & Kruskal, William H. (1954). Measures of association for cross classifications. Journal of the American Statistical Association, 49, 732–764.
• Hayes, Andrew F. & Krippendorff, Klaus (2007). Answering the call for a standard reliability measure for coding data. Communication Methods and Measures, 1, 77–89.
• Krippendorff, Klaus (2013). Content analysis: An introduction to its methodology, 3rd edition. Thousand Oaks, CA: Sage.
• Krippendorff, Klaus (1978). Reliability of binary attribute data. Biometrics, 34 (1), 142–144.
• Krippendorff, Klaus (1970). Estimating the reliability, systematic error, and random error of interval data. Educational and Psychological Measurement, 30 (1), 61–70.
• Lin, Lawrence I. (1989). A concordance correlation coefficient to evaluate reproducibility. Biometrics, 45, 255–268.
• Nunnally, Jum C. & Bernstein, Ira H. (1994). Psychometric Theory, 3rd ed. New York: McGraw-Hill.
• Pearson, Karl, et al. (1901). Mathematical contributions to the theory of evolution. IX: On the principle of homotyposis and its relation to heredity, to variability of the individual, and to that of race. Part I: Homotyposis in the vegetable kingdom. Philosophical Transactions of the Royal Society (London), Series A, 197, 285–379.
• Scott, William A. (1955). Reliability of content analysis: The case of nominal scale coding. Public Opinion Quarterly, 19, 321–325.
• Siegel, Sydney & Castella, N. John (1988). Nonparametric Statistics for the Behavioral Sciences, 2nd ed. Boston: McGraw-Hill.
• Tildesley, M. L. (1921). A first study of the Burmes skull. Biometrica, 13, 176–267.
• Spearman, Charles E. (1904). The proof and measurement of association between two things. American Journal of Psychology, 15, 72–101.
• Zwick, Rebecca (1988). Another look at interrater agreement. Psychological Bulletin, 103 (3), 347–387.

• Youtube video about Krippendorff’s alpha using SPSS and a macro.
• Reliability Calculator calculates Krippendorff's alpha.
• Krippendorff Alpha Javascript implementation and library
• Python implementation
• Krippendorff Alpha Ruby Gem implementation and library.
• Simpledorff Python implementation that works with Dataframes