STATISTICAL METHODS FOR ASSESSING AGREEMENT BETWEEN TWO METHODS OF CLINICAL MEASUREMENT J Martin Bland. Douglas G. Altman Department of Clinical Epidemiology and Social Medicine, St. George's Hospital Medical School. London SW17 ORE: and Division of Medical Statistics MRC Clinical Research Centre, Northwick Park Hospital, Harrow, Middlese SUMMARY In clinical measurement comparison of a new measurement technique with an established one is often needed to see whether they agree sufficiently for the new to replace the old. Such investigations are often analysed inappropriately, notably by using correlation coefficients The use of correlation is misleading. An alternative approach, based on graphical techniques and simple calculations, is described, together with the relation between this analysis and the assessment of repeatability ( Lancet,1986;i:307-310) INTRODUCTION Clinicians often wish to have data on, for example, cardiac stroke volume or blood pressure where direct measurement without adverse effects is difficult or impossible. The true values remain unknown. Instead indirect methods are used. and a new method has to be evaluated by comparison with an established technique rather than with the true quantity. If the new method agrees sufficiently well with the old, the old may be replaced. This is very different from calibration, where known quantities are measured by a new method and the result compared with the true value or with measurements made by a highly accurate method When two methods are compared neither provides an unequivocally correct measurement, so e try to assess the degree of agreement. But how? The correct statistical approach is not obvious. Many studies give the product-moment correlation coefficient(r) between the results of the two measurement methods as an indicator of agreement. It is no such thing. In a statistical journal we have proposed an altenative analysis, and clinical colleagues have suggested that we describe it for a medical readership Most of the analysis will be illustrated by a set of data(table 1)collected to compare two methods of measuring peak expiratory flow rate(PeFr) SAMPLE DATA The sample comprised colleagues and family of J M.B. chosen to give a wide range of PEFR but in no way representative of any defined population. Two measurements were made with a Wright peak flow meter and two with a mini Wright meter, in random order. All measurements were taken by J M. B, using the same two instruments. (These data were collected to demonstrate the statistical method and provide no evidence on the comparabilit of these two instruments. )We did not repeat suspect readings and took a single reading as our measurement of PEFR. Only the first measurement by each method is used to illustrate the comparison of methods, the second measurement being used in the study of repeatability PLOTTING DATA The first step is to plot the data and draw the line of equality on which all points would lie if the two meters gave exactly the same reading every time(fig 1). This helps the eye in gauging the degree of agreement between measurements, though, as we shall show, another type of plot is more informative
1 STATISTICAL METHODS FOR ASSESSING AGREEMENT BETWEEN TWO METHODS OF CLINICAL MEASUREMENT J. Martin Bland, Douglas G. Altman Department of Clinical Epidemiology and Social Medicine, St. George's Hospital Medical School, London SW17 ORE; and Division of Medical Statistics, MRC Clinical Research Centre, Northwick Park Hospital, Harrow, Middlesex SUMMARY In clinical measurement comparison of a new measurement technique with an established one is often needed to see whether they agree sufficiently for the new to replace the old. Such investigations are often analysed inappropriately, notably by using correlation coefficients. The use of correlation is misleading. An alternative approach, based on graphical techniques and simple calculations, is described, together with the relation between this analysis and the assessment of repeatability. (Lancet, 1986; i: 307-310) INTRODUCTION Clinicians often wish to have data on, for example, cardiac stroke volume or blood pressure where direct measurement without adverse effects is difficult or impossible. The true values remain unknown. Instead indirect methods are used, and a new method has to be evaluated by comparison with an established technique rather than with the true quantity. If the new method agrees sufficiently well with the old, the old may be replaced. This is very different from calibration, where known quantities are measured by a new method and the result compared with the true value or with measurements made by a highly accurate method. When two methods are compared neither provides an unequivocally correct measurement, so we try to assess the degree of agreement. But how? The correct statistical approach is not obvious. Many studies give the product-moment correlation coefficient (r) between the results of the two measurement methods as an indicator of agreement. It is no such thing. In a statistical journal we have proposed an alternative analysis, 1 and clinical colleagues have suggested that we describe it for a medical readership. Most of the analysis will be illustrated by a set of data (Table 1) collected to compare two methods of measuring peak expiratory flow rate (PEFR). SAMPLE DATA The sample comprised colleagues and family of J.M.B. chosen to give a wide range of PEFR but in no way representative of any defined population. Two measurements were made with a Wright peak flow meter and two with a mini Wright meter, in random order. All measurements were taken by J.M.B., using the same two instruments. (These data were collected to demonstrate the statistical method and provide no evidence on the comparability of these two instruments.) We did not repeat suspect readings and took a single reading as our measurement of PEFR. Only the first measurement by each method is used to illustrate the comparison of methods, the second measurement being used in the study of repeatability. PLOTTING DATA The first step is to plot the data and draw the line of equality on which all points would lie if the two meters gave exactly the same reading every time (fig 1). This helps the eye in gauging the degree of agreement between measurements, though, as we shall show, another type of plot is more informative
PEFR MEASURED WITH WRIGH PEAK FLOW AND MINI WRIG:T PEAK LOW METER Wright peak flow meter Mini Wright peak flow meter First pefr Second pefr first pefr Second pefr ubject (1/min) (1/min) (1/min) 525 395 397 430 415 516 512 508 476 500 6 611 600 625 7 413 415 364 431 650 638 658 642 433 445 432 11 417 420 432 420 12 656 633 626 275 478 477 178 165 259 268 372 370 17 421 451 443 已00 5 400 200 0 50o R by large meter (l/min) Fig 1. PEFR measured with large Wright peak flow meter and mini Wright peak flow meter with line of equality INAPPROPRIATE USE OF CORRELATION COEFFICIENT The second step is usually to calculate the correlation coefficient(r)between the two methods. For the data in fig 1, r=0.94(p<0.001). The null hypothesis here is that the measurements by the two methods are not linearly related. The probability is very small and we can safely conclude that PEFR measurements by the mini and large meters are related However this high correlation does not mean that the two methods agree (1)r measures the strength of a relation between two variables, not the agreement between them. We have perfect agreement only if the points in fig 1 lie along the line of equality, but we will have perfect correlation if the points lie along any straight line
2 PEFR MEASURED WITH WRIGHT PEAK FLOW AND MINI WRIGHT PEAK FLOW METER Wright peak flow meter Mini Wright peak flow meter First PEFR Second PEFR First PEFR Second PEFR Subject (l/min) (l/mi) (l/min) (l/min) 1 494 490 512 525 2 395 397 430 415 3 516 512 520 508 4 434 401 428 444 5 476 470 500 500 6 557 611 600 625 7 413 415 364 460 8 442 431 380 390 9 650 638 658 642 10 433 429 445 432 11 417 420 432 420 12 656 633 626 605 13 267 275 260 227 14 478 492 477 467 15 178 165 259 268 16 423 372 350 370 17 427 421 451 443 Fig 1. PEFR measured with large Wright peak flow meter and mini Wright peak flow meter, with line of equality. INAPPROPRIATE USE OF CORRELATION COEFFICIENT The second step is usually to calculate the correlation coefficient (r) between the two methods. For the data in fig 1, r = 0.94 (p < 0.001). The null hypothesis here is that the measurements by the two methods are not linearly related. The probability is very small and we can safely conclude that PEFR measurements by the mini and large meters are related. However, this high correlation does not mean that the two methods agree: (1) r measures the strength of a relation between two variables, not the agreement between them. We have perfect agreement only if the points in fig 1 lie along the line of equality, but we will have perfect correlation if the points lie along any straight line
(2)A change in scale of measurement does not affect the correlation, but it certainly affec the agreement. For example, we can measure subcutaneous fat by skinfold calipers. The calipers will measure two thicknesses of fat. If we were to plot calipers measurement against half-calipers measurement, in the style of fig 1, we should get a perfect straight line with slope 20. The correlation would be 1.0, but the two measurements would not agree-we could not mix fat thicknesses obtained by the two methods, since one is twice the other ()Correlation depends on the range of the true quantity in the sample. If this is wide, the correlation will be greater than if it is narrow. For those subjects whose PeFr (by peak flow meter)is less than 500 I/min, r is 0.88 while for those with greater PEFRs r is 0.90. Both are less than the overall correlation of 0.94, but it would be absurd to argue that agreement is worse below 500 I/min and worse above 500 l/min than it is for everybody. Since investigators usually try to compare two methods over the whole range of values typically encountered, a high correlation is almost guaranteed (4) The test of significance may show that the two methods are related, but it would be amazing if two methods designed to measure the same quantity were not related. The test of significance is irrelevant to the question of agreement (5)Data which seem to be in poor agreement can produce quite high correlations. For example, Serfontein and Jaroszewicz compared two methods of measuring gestational age Babies with a gestational age of 35 weeks by one method had gestations between 34 and 39.5 weeks by the other, but r was high(0. 85). On the other hand, Oldham et al. compared the mini and large Wright peak flow meters and found a correlation of 0.992. They then connected the meters in series so that both measured the same flow and obtained a"material improvement"(0.996). If a correlation coefficient of 0.99 can be materially improved upon, we need to rethink our ideas of what a high correlation is in this context. As we show below the high correlation of 0.94 for our own data conceals considerable lack of agreement between the two instruments MEASURING AGREEMENT It is most unlikely that different methods will agree exactly, by giving the identical result for all individuals. We want to know by how much the new method is likely to differ from the old: if this is not enough to cause problems in clinical interpretation we can replace the old give readings which differed by more than, say, 10 I/min, we could replace the large meter by the mini meter because so small a difference would not affect decisions on patient management. On the other hand, if the meters could differ by 100 l/ min, the mini meter would be unlikely to be satisfactory. How far apart measurements can be without causing difficulties will be a question of judgment. Ideally, it should be defined in advance to help in the interpretation of the method comparison and to choose the sample size The first step is to examine the data. A simple plot of the results of one method against those of the other(fig 1)though without a regression line is a useful start but usually the data points ill be clustered near the line and it will be difficult to assess between-method differences. A plot of the difference between the methods against their mean may be more informative. Fig 2 displays considerable lack of agreement between the large and mini meters, with discrepancies of up to 80 I/min, these differences are not obvious from fig 1. The plot of difference against mean also allows us to investigate any possible relationship between the measurement error and the true value. We do not know the true value. and the mean of the two measurements is the best estimate we have. It would be a mistake to plot the difference against either value separately because the difference will be related to each, a well-known statistical artefact. 4
3 (2) A change in scale of measurement does not affect the correlation, but it certainly affects the agreement. For example, we can measure subcutaneous fat by skinfold calipers. The calipers will measure two thicknesses of fat. If we were to plot calipers measurement against half-calipers measurement, in the style of fig 1, we should get a perfect straight line with slope 2.0. The correlation would be 1.0, but the two measurements would not agree — we could not mix fat thicknesses obtained by the two methods, since one is twice the other. (3) Correlation depends on the range of the true quantity in the sample. If this is wide, the correlation will be greater than if it is narrow. For those subjects whose PEFR (by peak flow meter) is less than 500 l/min, r is 0.88 while for those with greater PEFRs r is 0.90. Both are less than the overall correlation of 0.94, but it would be absurd to argue that agreement is worse below 500 l/min and worse above 500 l/min than it is for everybody. Since investigators usually try to compare two methods over the whole range of values typically encountered, a high correlation is almost guaranteed. (4) The test of significance may show that the two methods are related, but it would be amazing if two methods designed to measure the same quantity were not related. The test of significance is irrelevant to the question of agreement. (5) Data which seem to be in poor agreement can produce quite high correlations. For example, Serfontein and Jaroszewicz 2 compared two methods of measuring gestational age. Babies with a gestational age of 35 weeks by one method had gestations between 34 and 39.5 weeks by the other, but r was high (0.85). On the other hand, Oldham et al. 3 compared the mini and large Wright peak flow meters and found a correlation of 0.992. They then connected the meters in series, so that both measured the same flow, and obtained a "material improvement" (0.996). If a correlation coefficient of 0.99 can be materially improved upon, we need to rethink our ideas of what a high correlation is in this context. As we show below, the high correlation of 0.94 for our own data conceals considerable lack of agreement between the two instruments. MEASURING AGREEMENT It is most unlikely that different methods will agree exactly, by giving the identical result for all individuals. We want to know by how much the new method is likely to differ from the old: if this is not enough to cause problems in clinical interpretation we can replace the old method by the new or use the two interchangeably. If the two PEFR meters were unlikely to give readings which differed by more than, say, 10 l/min, we could replace the large meter by the mini meter because so small a difference would not affect decisions on patient management. On the other hand, if the meters could differ by 100 l/min, the mini meter would be unlikely to be satisfactory. How far apart measurements can be without causing difficulties will be a question of judgment. Ideally, it should be defined in advance to help in the interpretation of the method comparison and to choose the sample size. The first step is to examine the data. A simple plot of the results of one method against those of the other (fig 1) though without a regression line is a useful start but usually the data points will be clustered near the line and it will be difficult to assess between-method differences. A plot of the difference between the methods against their mean may be more informative. Fig 2 displays considerable lack of agreement between the large and mini meters, with discrepancies of up to 80 l/min, these differences are not obvious from fig 1. The plot of difference against mean also allows us to investigate any possible relationship between the measurement error and the true value. We do not know the true value, and the mean of the two measurements is the best estimate we have. It would be a mistake to plot the difference against either value separately because the difference will be related to each, a well-known statistical artefact. 4
- Mean + 2SD 吕0 巴sa Average PEFR by two meters (I/min) Fig 2. Difference against mean for peFr data For the pefr data, there is no obvious relation between the difference and the mean. Under these circumstances we can summarise the lack of agreement by calculating the bias estimated by the mean difference d and the standard deviation of the differences(s ). If there is a consistent bias we can adjust for it by subtracting d from the new method. For the PefR data the mean difference(large meter minus small meter)is -2 1 I/min and s is 38 I/min. We would expect most of the differences to lie between d-2s and d+ 2s(fig 2) If the differences are Normally distributed( Gaussian), 95% of differences will lie between these limits(or, more precisely, between d-196s and d+1.96s). Such differences are likely to follow a Normal distribution because we have removed a lot of the variation between subjects and are left with the measurement error. The measurements themselves do not have to follow a Normal distribution, and often they will not. We can check the distribution of the differences by drawing a histogram. If this is skewed or has very long tails the assumption of Normality may not be valid(see below) Provided differences within d 2s would not be clinically important, we could use the two measurement methods interchangeably. We shall refer to these as the"limits of agreement For the pefr data we get d-2s=-2.1-(2×388)=-79.7lmin d+2=-21+(2×38.8)=-75.5lmin Thus, the mini meter may be 80 I/min below or 76 I min above the large meter, which wo be unacceptable for clinical purposes. This lack of agreement is by no means obvious in PRECISION OF ESTIMATED LIMITS OF AGREEMENT The limits of agreement are only estimates of the values which apply to the whole population A second sample would give different limits. We might sometimes wish to use standard errors and confidence intervals to see how precise our estimates are, provided the differences follow a distribution which is approximately Normal. The standard error of d is vs2/n where n is the sample size, and the standard error of d-2s and d+2s is about v382/n 95% confidence intervals can be calculated by finding the appropriate point of the t distribution with n-1 degrees of freedom. on most tables the columns marked 5% or 0.0
4 Fig 2. Difference against mean for PEFR data. For the PEFR data, there is no obvious relation between the difference and the mean. Under these circumstances we can summarise the lack of agreement by calculating the bias, estimated by the mean difference d and the standard deviation of the differences (s ). If there is a consistent bias we can adjust for it by subtracting d from the new method. For the PEFR data the mean difference (large meter minus small meter) is - 2.1 l/min and s is 38.8 l/min. We would expect most of the differences to lie between d - 2s and d + 2s (fig 2). If the differences are Normally distributed (Gaussian), 95% of differences will lie between these limits (or, more precisely, between d -1.96s and d +1.96s). Such differences are likely to follow a Normal distribution because we have removed a lot of the variation between subjects and are left with the measurement error. The measurements themselves do not have to follow a Normal distribution, and often they will not. We can check the distribution of the differences by drawing a histogram. If this is skewed or has very long tails the assumption of Normality may not be valid (see below). Provided differences within d ± 2s would not be clinically important, we could use the two measurement methods interchangeably. We shall refer to these as the "limits of agreement". For the PEFR data we get: d - 2s = -2.1- (2 ´ 38.8) = -79.7 l/min d + 2s = -2.1 + (2 ´ 38.8) = -75.5 l/min Thus, the mini meter may be 80 l/min below or 76 l/min above the large meter, which would be unacceptable for clinical purposes. This lack of agreement is by no means obvious in fig 1. PRECISION OF ESTIMATED LIMITS OF AGREEMENT The limits of agreement are only estimates of the values which apply to the whole population. A second sample would give different limits. We might sometimes wish to use standard errors and confidence intervals to see how precise our estimates are, provided the differences follow a distribution which is approximately Normal. The standard error of d is s / n 2 , where n is the sample size, and the standard error of d - 2s and d + 2s is about 3s / n 2 . 95% confidence intervals can be calculated by finding the appropriate point of the t distribution with n -1 degrees of freedom, on most tables the columns marked 5% or 0.05
and then the confidence interval will be from the observed value minus t standard errors to the observed value plus t standard errors For the pefr data s=38.8. The standard error of d is thus 9 4. For the 95% confidence interval we have 16 degrees of freedom and t=2. 12. Hence the 95% confidence interval for the bias is -2.1-(2.12 x9.4)to-2.1+(2.12 x9.4), giving -22.0 to 17.8 I/min. The standard error of the limit d-2s is 16.3 1/min. The 95% confidence interval for the lower limit of agreement is-797-(2.12 x 16.3) to-797+(2. 12 x 16.3), giving -114.3 to 45. 1 I/min. For the upper limit of agreement the 95% confidence interval is 40.9 to 110.1 /min. These intervals are wide, reflecting the small sample size and the great variation of the differences. They show, however, that even on the most optimistic interpretation there can be considerable discrepancies between the two meters and that the degree of agreement is not EXAMPLE SHOWING GOOD AGREEMENT Fig 3 shows a comparison of oxygen saturation measured by an oxygen saturation monitor 0.42 percentage points with 95% confidence interval 0. 13 to 0.70. Thus pulsed oximeter< and pulsed oximeter saturation, a new non-invasive technique. Here the mean differenc saturation tends to give a lower reading, by between 0.13 and 0. 70. Despite this, the limits of agreement(-2.0 and 2.8)are small enough for us to be confident that the new method can be used in place of the old for clinical purposes Mean 2SD 0 Mean 2SD 100 Average of osm and pos Fig 3 Oxygen saturation monitor and pulsed saturation oximeter RELATION BETWEEN DIFFERENCE AND MEAN In the preceding analysis it was assumed that the differences did not vary in any systematic way over the range of measurement. This may not be so. Fig 4 compares the measurement of mean velocity of circumferential fibre shortening(VCF) by the long axis and short axis in M- mode echocardiography. The scatter of the differences increases as the VCF increases. W could ignore this, but the limits of agreement would be wider apart than necessary for smal VCF and narrower than they should be for large VCF. If the differences are proportional to the mean, a logarithmic transformation should yield a picture more like that of figs 2 and 4 and we can then apply the analysis described above to the transformed data
5 and then the confidence interval will be from the observed value minus t standard errors to the observed value plus t standard errors. For the PEFR data s = 38.8. The standard error of d is thus 9.4, For the 95% confidence interval we have 16 degrees of freedom and t = 2.12. Hence the 95% confidence interval for the bias is - 2.1- (2.12 ´ 9.4) to - 2.1+ (2.12 ´ 9.4), giving - 22.0 to 17.8 l/min. The standard error of the limit d - 2s is 16.3 l/min. The 95% confidence interval for the lower limit of agreement is - 79.7 - (2.12 ´16.3) to - 79.7 + (2.12 ´16.3) , giving -114.3 to - 45.1 l/min. For the upper limit of agreement the 95% confidence interval is 40.9 to 110.1 l/min. These intervals are wide, reflecting the small sample size and the great variation of the differences. They show, however, that even on the most optimistic interpretation there can be considerable discrepancies between the two meters and that the degree of agreement is not acceptable. EXAMPLE SHOWING GOOD AGREEMENT Fig 3 shows a comparison of oxygen saturation measured by an oxygen saturation monitor and pulsed oximeter saturation, a new non-invasive technique. 5 Here the mean difference is 0.42 percentage points with 95% confidence interval 0.13 to 0.70. Thus pulsed oximeter saturation tends to give a lower reading, by between 0.13 and 0.70. Despite this, the limits of agreement ( - 2.0 and 2.8) are small enough for us to be confident that the new method can be used in place of the old for clinical purposes. Fig 3. Oxygen saturation monitor and pulsed saturation oximeter RELATION BETWEEN DIFFERENCE AND MEAN In the preceding analysis it was assumed that the differences did not vary in any systematic way over the range of measurement. This may not be so. Fig 4 compares the measurement of mean velocity of circumferential fibre shortening (VCF) by the long axis and short axis in Mmode echocardiography. 6 The scatter of the differences increases as the VCF increases. We could ignore this, but the limits of agreement would be wider apart than necessary for small VCF and narrower than they should be for large VCF. If the differences are proportional to the mean, a logarithmic transformation should yield a picture more like that of figs 2 and 4, and we can then apply the analysis described above to the transformed data
Mean -----. Mean Mean 2SD ,8 Average vcr by long and short axes Fig 4. Mean VCf by long and short axis measurements Fig 5 shows the log-transformed data of fig 4. This still shows a relation between the difference and the mean VCF, but there is some improvement. The mean difference is 0.003 on the log scale and the limits of agreement are -0.098 and 0.106.However, although there is only negligible bias, the limits of agreement have somehow to be related to the original scale of measurement. If we take the antilogs of these limits, we get 0. 80 and 1.27 However, the antilog of the difference between two values on a log scale is a dimensionless ratio. The limits tell us that for about 95% of cases the short axis measurement of vcf will be between 0.80 and 1. 27 times the long axis VCF. Thus the short axis measurement may differ from the long axis measurement by 20% below to 27% above. (The log transformation is the only transformation giving back-transformed differences which are easy to interpret and we do not recommend the use of any other in this context. Mean 2sD 0 ---Mea∩ g 8 一.1 Mean 25D Average log vcf by long and short axes Fig 5. Data of fig 4 after logarithmic transformation. i These numbers were incorrectly printed as 0.008, -0.226, and 0.243 in the Lancet. This mistake arose because when revising the paper we dithered over whether to use natural logs or logs to base 10 and got hopelessly confused
6 Fig 4. Mean VCF by long and short axis measurements. Fig 5 shows the log-transformed data of fig 4. This still shows a relation between the difference and the mean VCF, but there is some improvement. The mean difference is 0.003 † on the log scale and the limits of agreement are - 0.098 † and 0.106.† However, although there is only negligible bias, the limits of agreement have somehow to be related to the original scale of measurement. If we take the antilogs of these limits, we get 0.80 and 1.27. However, the antilog of the difference between two values on a log scale is a dimensionless ratio. The limits tell us that for about 95% of cases the short axis measurement of VCF will be between 0.80 and 1.27 times the long axis VCF. Thus the short axis measurement may differ from the long axis measurement by 20% below to 27% above. (The log transformation is the only transformation giving back-transformed differences which are easy to interpret, and we do not recommend the use of any other in this context.) Fig 5. Data of fig 4 after logarithmic transformation. † These numbers were incorrectly printed as 0.008, -0.226, and 0.243 in the Lancet. This mistake arose because when revising the paper we dithered over whether to use natural logs or logs to base 10 and got hopelessly confused
Sometimes the relation between difference and mean is more complex than that shown in fi 4 and log transformation does not work. Here a plot in the style of fig 2 is very helpful in comparing the methods. Formal analysis, as described above, will tend to give limits of agreement which are too far apart rather than too close, and so should not lead to the acceptance of poor methods of measurement REPEATABILITY Repeatability is relevant to the study of method comparison because the repeatabilities of the two methods of measurement limit the amount of agreement which is possible. If one method has poor repeatability -i.e. there is considerable variation in repeated measurements on the me subject- the agreement between the two methods is bound to be poor too. When the old method is the more variable one, even a new method which is perfect will not agree with it. If both methods have poor repeatability, the problem is even worse The best way to examine repeatability is to take repeated measurements on a series of subjects. The table shows paired data for PEFR. We can then plot a figure similar to fig 2, showing differences against mean for each subject. If the differences are related to the mean, we can apply a log transformation. We then calculate the mean and standard deviation of the differences as before. The mean difference should here be zero since the same method was used.(If the mean difference is significantly different from zero, we will not be able to use the second or the process of measurement is altering the quantity )We expect 901on, Tecting the data to assess repeatability because either knowledge of the first measurement is affe differences to be less than two standard deviations. This is the definition of a repeatability coefficient adopted by the British Standards Institution. If we can assume the main difference to be zero this coefficient is very simple to estimate: we square all the differences, add them up, divide by n, and take the square root, to get the standard deviation of the differences 20 40 g-80 100 200 800 Average PEFR by mini meter (I/min) Fig 6. Repeated measures of PEFR using mini Wright peak flow meter Fig 6 shows the plot for pairs of measurements made with the mini wright peak flow meter There does not appear to be any relation between the difference and the size of the PEFr. There is, however, a clear outlier. We have retained this measurement for the analysis although we suspect that it was technically unsatisfactory. (In practice, one could omit this subject )The sum of the differences squared is 13479 so the standard deviation of differences between the 17 pairs of repeated measurements is 28. 2 I/ min. The coefficient of repeatability
7 Sometimes the relation between difference and mean is more complex than that shown in fig 4 and log transformation does not work. Here a plot in the style of fig 2 is very helpful in comparing the methods. Formal analysis, as described above, will tend to give limits of agreement which are too far apart rather than too close, and so should not lead to the acceptance of poor methods of measurement. REPEATABILITY Repeatability is relevant to the study of method comparison because the repeatabilities of the two methods of measurement limit the amount of agreement which is possible. If one method has poor repeatability — i.e. there is considerable variation in repeated measurements on the same subject — the agreement between the two methods is bound to be poor too. When the old method is the more variable one, even a new method which is perfect will not agree with it. If both methods have poor repeatability, the problem is even worse. The best way to examine repeatability is to take repeated measurements on a series of subjects. The table shows paired data for PEFR. We can then plot a figure similar to fig 2, showing differences against mean for each subject. If the differences are related to the mean, we can apply a log transformation. We then calculate the mean and standard deviation of the differences as before. The mean difference should here be zero since the same method was used. (If the mean difference is significantly different from zero, we will not be able to use the data to assess repeatability because either knowledge of the first measurement is affecting the second or the process of measurement is altering the quantity.) We expect 95% of differences to be less than two standard deviations. This is the definition of a repeatability coefficient adopted by the British Standards Institution. 7 If we can assume the main difference to be zero this coefficient is very simple to estimate: we square all the differences, add them up, divide by n , and take the square root, to get the standard deviation of the differences. Fig 6. Repeated measures of PEFR using mini Wright peak flow meter. Fig 6 shows the plot for pairs of measurements made with the mini Wright peak flow meter. There does not appear to be any relation between the difference and the size of the PEFR. There is, however, a clear outlier. We have retained this measurement for the analysis, although we suspect that it was technically unsatisfactory. (In practice, one could omit this subject.) The sum of the differences squared is 13479 so the standard deviation of differences between the 17 pairs of repeated measurements is 28.2 l/min. The coefficient of repeatability
is twice this. or 56.4 1/min for the mini meter. For the large meter the coefficient is 43.2 I/min If we have more than two repeated measurements the calculations are more complex. w plot the standard deviation of the several measurements for that subject against their mean and then use one-way analysis of variance, which is beyond the scope of this article MEASURING AGREEMENT USING REPEATED MEASUREMENTS If we have repeated measurements by each of the two methods on the same subjects we can calculate the mean for each method on each subject and use these pairs of means to compare the two methods using the analysis for assessing agreement described above. The estimate of bias will be unaffected but the estimate of the standard deviation of the differences will be too small, because some of the effect of repeated measurement error has been removed. We can correct for this. Suppose we have two measurements obtained by each method, as in the table. We find the standard deviation of differences between repeated measurements for each method separately, S, and S,, and the standard deviation of the differences between the means for each method, Sp. The corrected standard deviation of differences, s 152+4142. his is approximately v2s,, but if there are differences between the two methods not explicable by repeatability errors alone (i.e. interaction between subject and measurement method) this approximation may produce an overestimate. For the Pefr, we hc32=2152321 min.sis thus32×126+4×282 or 37.7 I/ min. Compare this with the estimate 38.8 I/min which was obtained using a single measurement. On the other hand, the approximation v2sD gives an overestimate(47.0 I/min) DISCUSSION In the analysis of measurement method comparison data, neither the correlation coefficient (as we show here)nor techniques such as regression analysis are appropriate. We suggest replacing these misleading analyses by a method that is simple both to do and to interpret Further, the same method may be use to analyse the repeatability of a single measurement method or to compare mea by two observers Why has a totally inappropriate method, the correlation coefficient, become almost universally used for this purpose?. Two processes may be at work here -- namely, pattern recognition and imitation. A likely first step in the analysis of such data is to plot a scatter diagram(fig 1). a glance through almost any statistical textbook for a similar picture will lead to the correlation coefficient as a method of analysis of such a plot, together with a test of the null hypothesis of no relationship. Some texts even use pairs of measurements by two different methods to illustrate the calculation of r. Once the correlation approach has been published, others will read of a statistical problem similar to their own being solved in this way and will use the same technique with their own data. Medical statisticians who ask"why did you use this statistical method? "will often be told"because this published paper used it Journals could help to rectify this error by returning for reanalysis papers which use incorrect statistical techniques. This may be a slow process. Referees, inspecting papers in which two methods of measurement have been compared, sometimes complain if no correlation coefficients are provided, even when the reasons for not doing so are given This was incorrectly printed as s in the Lancet and in Biochimica Clinica
8 is twice this, or 56.4 l/min for the mini meter. For the large meter the coefficient is 43.2 l/min. If we have more than two repeated measurements the calculations are more complex. We plot the standard deviation of the several measurements for that subject against their mean and then use one-way analysis of variance, 8 which is beyond the scope of this article. MEASURING AGREEMENT USING REPEATED MEASUREMENTS If we have repeated measurements by each of the two methods on the same subjects we can calculate the mean for each method on each subject and use these pairs of means to compare the two methods using the analysis for assessing agreement described above. The estimate of bias will be unaffected, but the estimate of the standard deviation of the differences will be too small, because some of the effect of repeated measurement error has been removed. We can correct for this. Suppose we have two measurements obtained by each method, as in the table. We find the standard deviation of differences between repeated measurements for each method separately, 1 s and 2 s , and the standard deviation of the differences between the means for each method, D s . The corrected standard deviation of differences, c s , is 2 2 2 1 2 4 1 4 1 s s s D + + . This is approximately 2 2 D s , but if there are differences between the two methods not explicable by repeatability errors alone (i.e. interaction between subject and measurement method) this approximation may produce an overestimate. For the PEFR, we have D s = 33.2, 1 s = 21.6, 2 s = 28.2 ‡ l/min. c s is thus 2 2 2 28.2 4 21.6 1 4 33.2 + 1 ´ + ´ or 37.7 l/min. Compare this with the estimate 38.8 l/min which was obtained using a single measurement. On the other hand, the approximation 2 2 D s gives an overestimate (47.0 l/min). DISCUSSION In the analysis of measurement method comparison data, neither the correlation coefficient (as we show here) nor techniques such as regression analysis 1 are appropriate. We suggest replacing these misleading analyses by a method that is simple both to do and to interpret. Further, the same method may be use to analyse the repeatability of a single measurement method or to compare measurements by two observers. Why has a totally inappropriate method, the correlation coefficient, become almost universally used for this purpose?. Two processes may be at work here --- namely, pattern recognition and imitation. A likely first step in the analysis of such data is to plot a scatter diagram (fig 1). A glance through almost any statistical textbook for a similar picture will lead to the correlation coefficient as a method of analysis of such a plot, together with a test of the null hypothesis of no relationship. Some texts even use pairs of measurements by two different methods to illustrate the calculation of r . Once the correlation approach has been published, others will read of a statistical problem similar to their own being solved in this way and will use the same technique with their own data. Medical statisticians who ask "why did you use this statistical method?" will often be told "because this published paper used it". Journals could help to rectify this error by returning for reanalysis papers which use incorrect statistical techniques. This may be a slow process. Referees, inspecting papers in which two methods of measurement have been compared, sometimes complain if no correlation coefficients are provided, even when the reasons for not doing so are given. ‡ This was incorrectly printed as c s in the Lancet and in Biochimica Clinica
We thank of our colleagues for their interest and assistance, including Dr David robson whe first brought the problem to us, Dr P D Arbela and Dr H. Seeley for the use of their data; and Mrs s Stevens for typing the manuscript REFERENCES 1. Altman DG, Bland JM. Measurement in medicine: the analysis of method comparison studies. The Statistician 1983: 32.307-317 2. Serfontein GL, Jaroszewicz AM. Estimation of gestational age at birth: comparison of two methods. Arch Dis Child 1978: 53: 509-11 3. Oldham HG, Bevan MM, McDermott M. Comparison of the new miniature Wright peak flow meter with the standard Wright peak flow meter. Thorax 1979: 34: 807-08 4. Gill JS, Zezulka AV, Beevers DG, Davies P. Relationship between initial blood pressure and its fall with treatment. Lancet 1985: i: 567-69 5. Tytler JA, Seeley HF. The Nellcor N-101 pulse oximeter-a clinical-evaluation in anesthesia and intensive-care. Anaesthesia 1986: 41: 302-305 6. D Arbela PG, Silayan ZM, Bland JM. Comparability of M-mode echocardiographic long axis and short axis left ventricular function derivatives. British heart Journal 1986: 56: 445- 7. British Standards Institution. Precision of test methods /: Guide for the determination and reproducibility for a standard test method(BS 597, Part 1). London: BSI (1975) 8. Armitage P. Statistical methods in medical research. Oxford: Blackwell Scientific Publications. 1971: chap 7 Reproduced by kind permission of the lancet
9 We thank of our colleagues for their interest and assistance, including Dr David Robson who first brought the problem to us, Dr P. D'Arbela and Dr H. Seeley for the use of their data; and Mrs S Stevens for typing the manuscript. REFERENCES 1. Altman DG, Bland JM. Measurement in medicine: the analysis of method comparison studies. The Statistician 1983; 32, 307-317. 2. Serfontein GL, Jaroszewicz AM. Estimation of gestational age at birth: comparison of two methods. Arch Dis Child 1978; 53: 509-11. 3. Oldham HG, Bevan MM, McDermott M. Comparison of the new miniature Wright peak flow meter with the standard Wright peak flow meter. Thorax 1979; 34: 807-08. 4. Gill JS, Zezulka AV, Beevers DG, Davies P. Relationship between initial blood pressure and its fall with treatment. Lancet 1985; i: 567-69. 5. Tytler JA, Seeley HF. The Nellcor N-101 pulse oximeter - a clinical-evaluation in anesthesia and intensive-care. Anaesthesia 1986; 41: 302-305. 6. D'Arbela PG, Silayan ZM, Bland JM. Comparability of M-mode echocardiographic long axis and short axis left ventricular function derivatives. British Heart Journal 1986; 56: 445- 9. 7. British Standards Institution. Precision of test methods 1: Guide for the determination and reproducibility for a standard test method (BS 597, Part 1). London: BSI (1975). 8. Armitage P. Statistical methods in medical research. Oxford: Blackwell Scientific Publications, 1971: chap 7. Reproduced by kind permission of the Lancet
