4.25

4.2 | 4 ratings Rate this file 30 Downloads (last 30 days) File Size: 27.1 KB File ID: #20545
image thumbnail

Clarke Error Grid Analysis

by

 

01 Jul 2008 (Updated )

Clarke EGA quantifies the accuracy of glucose estimates generated by meters as compared to reference

| Watch this File

File Information
Description

The Clarke error grid approach is used to assess the clinical significance of differences between the glucose measurement technique under test and the venous blood glucose reference measurements. The method uses a Cartesian diagram, in which the values predicted by the technique under test are displayed on the y-axis, whereas the values received from the reference method are displayed on the x-axis. The diagonal represents the perfect agreement between the two, whereas the points below and above the line indicate, respectively, overestimation and underestimation of the actual values. Zone A (acceptable) represents the glucose values that deviate from the reference values by ±20% or are in the hypoglycemic range (<70 mg/dl), when the reference is also within the hypoglycemic range. The values within this range are clinically exact and are thus characterized by correct clinical treatment. Zone B (benign errors) is located above and below zone A; this zone represents those values that deviate from the reference values, which are incremented by 20. The values that fall within zones A and B are clinically acceptable, whereas the values included in areas C-E are potentially dangerous, and there is a possibility of making clinically significant mistakes. [1-4]

 Syntax:
 
 [total, percentage] = clarke(y,yp)
 
 Inputs:
 y = reference values (mg/dl)
 yp = predicted/estimtated values (mg/dl)
 
 Outputs:
 total = total points per zone:
 total(1) = zone A,
 total(2) = zone B, and so on

 percentage = percentage of data which fell in certain region:
 percentage(1) = zone A,
 percentage(2) = zone B, and so on.
 
 Example:
 load example_data.mat
 [tot, per] = clarke(y,yp)

 References:
 [1] A. Maran et al., “Continuous subcutaneous glucose monitoring in diabetic patients: a multicenter analysis,” Diabetes Care, vol. 25, no. 2, pp. 347–352, Feb. 2002.
 [2] B. P. Kovatchev et al. “Evaluating the accuracy of continuous glucose-monitoring sensors: continuous glucose-error grid analysis illustrated by TheraSense Freestyle Navigator data,” Diabetes Care, vol. 27, no. 8, pp. 1922–1928, Aug. 2004.
 [3] E. Guevara and F. J. Gonzalez, “Prediction of Glucose Concentration by Impedance Phase Measurements,” in MEDICAL PHYSICS: Tenth Mexican Symposium on Medical Physics, Mexico City (Mexico), 2008, vol. 1032, pp. 259–261.
 [4] E. Guevara and F. J. Gonzalez, “Joint optical-electrical technique for noninvasive glucose monitoring,” REVISTA MEXICANA DE FISICA, vol. 56, no. 5, pp. 430–434, Sep. 2010.
 
© Edgar Guevara Codina
codina@REMOVETHIScactus.iico.uaslp.mx
File Version 1.2
March 29 2013
 
Ver. 1.2 Statistics verified, fixed some errors in the display; thanks to Tim Ruchti from Hospira Inc. for the corrections
Ver. 1.1 corrected upper B-C boundary, lower B-C boundary slope ok; thanks to Steven Keith from BD Technologies for the corrections!
MATLAB ver. 7.10.0.499 (R2010a)

MATLAB release MATLAB 7.10 (R2010a)
Tags for This File   Please login to tag files.
Please login to add a comment or rating.
Comments and Ratings (9)
09 Feb 2014 Wei

excellent!

26 Sep 2013 Edgar Guevara

@Jan Simon: Thanks for the kind words!

26 Sep 2013 Jan Simon

The purpose of this function is not my field of science, therefore I cannot rate the actual calculations. But this code appeared in discussion and I'm impressed by the clean and complete help section.

01 Apr 2013 Edgar Guevara

Ver. 1.2 Statistics verified, fixed some errors in the display; thanks to Tim Ruchti from Hospira Inc. for the corrections

09 Jan 2011 Peter

Thank you for the interesting code! However, I suspect that there is an inconsistency in the interpretation of zone A and B. The totals and the plots are not compatible e.g. zoom in on clarke(200,243) and clarke(200, 160).

22 Mar 2010 Peter Jacobs

Thank you very much for the useful code. Well done!

04 Nov 2008 Edgar Guevara

Steven,

First of all, thank you for the suggestions, I think I might have made some mistakes on the regions boundaries, since I measured all boundaries from a printed plot (on the referred articles) so nothing was very exact. I hope you could provide me with the original work from Clarke, to make the necessary corrections.

09 Sep 2008 Steven Keith

Overall, the logic is excellent. There are 2 suggestions: for the lower C region, the slope of the boundary should be 7/5 rather than 6/5. For the upper B-C boundary, there are two schools of though displayed in error grids being published. I think Clarke still uses a line parallel to 45deg rather than a line parallel to the +20% line. Very nice code.

03 Sep 2008 Guillermo Quintas

A really useful job!.

Updates
18 Nov 2008

Ver. 1.1 corrected upper B-C boundary, lower B-C boundary slope ok, thanks to Steven Keith from BD Technologies for the corrections!

30 May 2012

Added some more references

29 Mar 2013

Ver. 1.2 Statistics verified, fixed some errors in the display; thanks to Tim Ruchti from Hospira Inc. for the corrections.

01 Apr 2013

Ver. 1.2 Statistics verified, fixed some errors in the display; thanks to Tim Ruchti from Hospira Inc. for the corrections

Contact us