측정 불확도 반영 4단계 — CRM 포함 한계값 계산

측정 불확도 반영 관련 이미지

측정 불확도 반영은 CRM 인증값의 불확실성 u_CRM을 QC 표준편차와 합성해 관리한계를 다시 계산하는 작업이다. 인증서의 확장불확도를 포함인자로 나눠 u_CRM을 구하고, s_total = √(s_QC² + u_CRM²)로 경보·관리한계를 재설정한다. u_CRM이 s_QC의 0.3배 미만이면 한계 확장률이 5% 이하로 미미하지만, 0.5배를 넘으면 한계를 넓히기 전에 CRM 등급 재검토가 먼저다.

측정 불확도 반영 관련 이미지

측정 불확도 반영이 필요한 이유

정도관리 관리도의 중심선을 CRM 인증값에 두면서 한계값은 QC 측정값의 표준편차만으로 계산하는 관행이 흔하다. 이 방식은 중심선 자체가 흔들릴 수 있다는 사실을 무시한 것이다.

인증값에는 균질성, 특성화, 안정성 기여분이 합성된 불확도가 따라붙는다. 즉 기준선은 하나의 점이 아니라 구간이며, 이 구간 폭을 빼놓으면 관리한계가 실제보다 좁게 잡힌다.

한계가 과도하게 좁으면 거짓 경보가 늘고, 원인 없는 OOL 조사에 인력이 소모된다. 반대로 무작정 넓히면 실제 부적합을 놓치므로, 근거 있는 폭을 계산해 두는 것이 관건이다.

ISO/IEC 17025:2017은 7.6에서 유의미한 모든 불확도 기여분을 고려하도록 요구하고, 7.7에서 결과 유효성 확보 수단으로 표준물질 사용을 명시한다. 관리한계는 두 조항이 만나는 지점이므로 측정 불확도 반영 여부가 심사에서 확인 대상이 된다.

1단계 — 인증서에서 u_CRM 확보

첫 단계는 CRM 인증서에서 확장불확도 U와 포함인자 k를 그대로 읽어오는 것이다. 대부분의 인증서는 k=2, 약 95% 신뢰수준 기준으로 U를 표기한다.

표준불확도는 u_CRM = U / k로 환산한다. 인증값 5.00 mg/L, U = 0.06 mg/L(k=2)이면 u_CRM은 0.03 mg/L이다.

k가 명시되지 않은 인증서는 그대로 쓰지 않는다. 공급업체에 포함인자와 유효자유도를 문의해 기록으로 남기고, 확인이 불가능하면 k=2로 가정한 근거를 절차서에 적는다.

인증서에 유효기간, 보관조건, 개봉 후 안정성 기여분이 별도로 표기돼 있으면 이를 u_CRM에 포함할지 판단한다. 측정 불확도 반영 계산은 입력값의 출처가 문서로 추적될 때만 감사에서 인정된다.

인증서에서 읽을 값

항목 예시 처리
인증값 5.00 mg/L 관리도 중심선
확장불확도 U 0.06 mg/L k로 나눔
포함인자 k 2 미표기 시 문의
u_CRM 0.03 mg/L 합성 입력값

2단계 — 방법 표준편차 s_QC 산출

두 번째는 자체 측정 변동을 대표하는 표준편차 s_QC를 구하는 단계다. 동일 CRM을 정상 운영 조건에서 반복 측정한 QC 데이터가 근거가 된다.

표본수는 최소 20점 이상을 권장하며, 서로 다른 날·분석가·시약 로트를 포함해 일 간 재현성이 담기도록 한다. 하루에 몰아 찍은 20점은 반복성만 반영해 s_QC를 과소평가한다.

이상점은 제거하기 전에 원인을 조사한다. 기록된 원인 없이 뺀 점은 s_QC를 인위적으로 줄여, 이후의 측정 불확도 반영 효과 판단까지 왜곡한다.

앞 예시에서 s_QC를 0.10 mg/L로 얻었다고 하자. 이 값과 u_CRM 0.03 mg/L이 다음 단계의 두 입력이 된다.

3단계 — 합성 표준편차로 한계값 재계산

두 성분은 서로 독립으로 보고 제곱합의 제곱근으로 합성한다. s_total = √(s_QC² + u_CRM²)이며, GUM이 규정한 합성표준불확도 산출 방식과 같은 형태다.

예시 값을 넣으면 s_total = √(0.10² + 0.03²) ≈ 0.1044 mg/L가 된다.

경보한계는 인증값 ± 2·s_total, 관리한계는 인증값 ± 3·s_total로 둔다. Nordtest 내부 정도관리 지침이 제시하는 ±2σ 경보·±3σ 관리 구조를 그대로 쓰되, σ 자리에 s_total을 넣는 것이 측정 불확도 반영의 실제 계산이다.

앞 예시의 경보한계는 5.00 ± 0.209 mg/L로, 불확도를 무시했을 때의 ±0.200 mg/L보다 넓다. 관리한계는 5.00 ± 0.313 mg/L가 된다.

한계값 합성 계산 순서

  1. u_CRM 환산
    U 0.06을 k=2로 나눠 0.03 mg/L
  2. s_QC 확인
    20점 이상 QC로 0.10 mg/L 산출
  3. 제곱합 합성
    √(0.10²+0.03²) ≈ 0.1044 mg/L
  4. 한계 산출
    경보 ±2s_total, 관리 ±3s_total

u_CRM 비중이 한계 폭에 미치는 영향

측정 불확도 반영의 실효는 u_CRM과 s_QC의 비율이 결정한다. 한계 확장률은 √(1+(u_CRM/s_QC)²)−1로 계산된다.

비율이 0.1이면 확장률은 0.5%에 그친다. 이 구간에서는 반영해도 한계선이 사실상 움직이지 않으므로, 계산은 수행하되 결과 영향이 없다는 판단을 기록으로 남기는 편이 낫다.

비율 0.3에서 확장률은 4.4%, 0.5에서 11.8%로 커진다. 이 지점부터는 반영하지 않은 한계를 그대로 쓸 경우 거짓 경보가 체감될 만큼 늘어난다.

비율이 1.0이면 확장률이 41.4%에 이른다. 기준값 불확실성이 자체 측정 변동만큼 크다는 뜻이므로, 한계를 넓히기보다 더 낮은 U를 가진 CRM으로 교체하는 쪽이 정도관리 목적에 맞는다.

u_CRM/s_QC 비율별 판정

비율 한계 확장률 판정
0.1 0.5% 반영 생략 가능
0.3 4.4% 반영 권장
0.5 11.8% CRM 등급 재검토
1.0 41.4% CRM 교체 검토

※ 확장률은 √(1+비율²)−1로 계산한 값이다.

4단계 — 측정 불확도 반영 결과 검증과 문서화

재계산한 한계를 과거 QC 데이터에 되짚어 적용해 경보·OOL 발생률이 어떻게 변하는지 확인한다. 발생률이 목표 수준에서 크게 벗어나면 s_QC 산출 구간부터 다시 본다.

한계가 넓어진 만큼 탐지력은 떨어진다. 넓힌 한계가 고객 요구 허용폭이나 규제 기준을 침범하지 않는지 반드시 대조하고, 침범하면 방법 개선으로 s_QC를 줄이는 방향을 택한다.

문서에는 인증서 번호와 로트, U와 k, u_CRM, s_QC, 표본수, s_total, 계산일, 승인자를 남긴다. 계산 시트는 원본 QC 데이터와 연결 가능한 형태로 보존한다.

CRM 로트가 바뀌면 인증값과 U가 모두 바뀌므로 재계산이 필수다. 로트 전환 시점은 관리도에 표시해 이전 구간과 시각적으로 구분한다.

관련해서 관리한계 CRM 불확도 반영 4단계 — 기준값 불확실성 포함하기 글도 참고할 만하다.

정리와 체크포인트

순서는 u_CRM 확보, s_QC 산출, 제곱합 합성, 검증·문서화의 네 단계다. 계산 자체는 단순하지만 입력값의 출처와 수집 조건이 결과를 좌우한다.

가장 흔한 실수는 인증서의 U를 k로 나누지 않고 그대로 합성에 넣는 것이다. 이 경우 한계가 필요 이상으로 넓어져 실제 부적합을 놓친다.

두 번째 실수는 반복성만 담긴 s_QC를 쓰는 것이다. 분모가 작으면 u_CRM 비중이 과대평가돼 한계가 부풀려진다.

측정 불확도 반영은 한 번 하고 끝나는 작업이 아니다. 로트 변경, 기기 교체, 분석법 개정 시마다 재계산 대상인지 판정하고 그 판단 근거를 기록으로 남긴다.

반영 후 점검 항목

  • 인증서 U를 포함인자 k로 나눴는가
  • s_QC에 일 간 재현성이 포함됐는가
  • 과거 QC 데이터로 역적용해 봤는가
  • 넓힌 한계가 고객 허용폭 이내인가
  • 로트 전환 시점을 관리도에 표시했는가

참고 자료

답글 남기기

이메일 주소는 공개되지 않습니다. 필수 필드는 *로 표시됩니다