検量線の相関係数は濃度と応答の直線的な関係、傾きは濃度1単位あたりの応答変化、切片は濃度0で予測される応答を表します。 相関係数が高くても、低濃度側の誤差や曲がりがないとは限りません。残差、標準液の逆算濃度、独立したQC(品質管理)試料の結果を合わせて確認します。
↔ 表が見切れている場合は、左右にスクロールしてご覧ください。↔ Scroll the table horizontally to see any hidden columns.
| 確認したいこと | まず見るポイント |
|---|---|
| 相関係数はどのくらい必要? | 一律の基準はなく、適用する試験法・SOP(標準作業手順書)の基準と、残差・逆算値・QCを確認 |
| 切片がゼロでなくてもよい? | ゼロでないだけでは不適合とはいえず、ブランクや低濃度の定量への影響を確認 |
| 傾きが大きい方がよい? | 同じ測定条件では感度の指標になるが、ノイズやQCも合わせて確認 |
| 検量線から濃度を求めるには? | 回帰式 y = ax + b を求め、試料応答を x = (y − b) / a へ代入 |
| R² = 0.999は誤差0.1%という意味? | いいえ。個々の濃度の誤差は、標準液の逆算値や独立QCで確認 |
検量線は、既知濃度の標準液と装置応答の関係を表し、試料の応答から濃度を求めるために使います。分析ソフトには傾き、切片、相関係数または決定係数が表示されますが、相関係数が基準を満たしただけでは、検量線が測定目的に適しているとは判断できません。
例えば、広い濃度範囲では高濃度点の影響によって相関が非常に高く見えても、低濃度側の誤差が大きい場合があります。また、測定点が緩やかに曲がっていても、相関係数だけでは見逃すことがあります。
この記事では、一般的な直線検量線を前提に、傾き、切片、相関係数、残差、濃度範囲の見方を実務的に整理します。具体的な合否基準は、適用する試験法、規格、分析法バリデーション、試験所のSOPに従って設定してください。
検量線とは?式 y = ax + b の意味
直線検量線は、一般に次の式で表します。
y = ax + b
↔ 表が見切れている場合は、左右にスクロールしてご覧ください。↔ Scroll the table horizontally to see any hidden columns.
| 記号 | 意味 | 実務で確認すること |
|---|---|---|
| x | 標準液または試料中の分析対象濃度 | 単位、希釈倍率、内標準補正前後の定義 |
| y | 装置応答 | 発光強度、カウント、ピーク面積、応答比など |
| a | 傾き | 濃度変化に対する応答変化、日間の感度変動 |
| b | 切片 | ゼロ濃度付近の応答、ブランク、背景補正の影響 |
試料の濃度は、測定応答を式へ代入して x = (y − b) / a として求めます。そのため、傾きや切片の変化は、特に低濃度試料の換算値へ直接影響します。
検量線の求め方:標準液の測定から濃度計算まで
基本の流れは 「標準液を測る → 濃度と応答を回帰する → 妥当性を確認する → 試料濃度を逆算する」 です。横軸を既知濃度x、縦軸を装置応答yとし、試料の未知濃度を標準液の代わりに使わないようにします。
1.既知濃度と装置応答を対応させる
次は計算手順を示すための架空のデータです。実際の標準液数・濃度配置・ブランクの扱いは、適用する試験法に従います。
↔ 表が見切れている場合は、左右にスクロールしてご覧ください。↔ Scroll the table horizontally to see any hidden columns.
| 標準液濃度 x(µg/L) | 装置応答 y(任意単位) |
|---|---|
| 0 | 22 |
| 1 | 118 |
| 2 | 225 |
| 3 | 316 |
| 4 | 419 |
2.最小二乗法で傾きと切片を求める
切片を含む、重み付けなしの直線回帰では、次の式を使います。x̄は濃度の平均、ȳは応答の平均、Σは全測定点について足し合わせる記号です。
a = Σ[(x − x̄)(y − ȳ)] / Σ[(x − x̄)²]
b = ȳ − a × x̄
このデータでは x̄ = 2、ȳ = 220、分子 = 992、分母 = 10 なので、傾きaは99.2、切片bは21.6となり、検量線は y = 99.2x + 21.6 です。傾きの単位は応答単位/(µg/L)、切片の単位は応答単位です。2点だけを選んで結ぶのではなく、使用する全標準点を回帰に含めます。
最小二乗法は、実測応答と回帰式の予測応答との差の二乗和を小さくする方法です。NIST(米国国立標準技術研究所):線形最小二乗回帰
3.R²と、各点のずれを確認する
切片を含む通常の最小二乗回帰では、R² = 1 − Σ(y − ŷ)² / Σ(y − ȳ)² で計算できます。ŷは回帰式から求めた予測応答です。この例では残差平方和43.6、応答の全平方和98,450から、R² ≈ 0.99956 となります。
一方、1 µg/Lの標準液の逆算値は (118 − 21.6) / 99.2 ≈ 0.972 µg/L、設定濃度からの相対誤差は約−2.8%です。R²が高いことと、各点の誤差が十分小さいことは区別し、残差・逆算値・独立QCを確認します。この説明用データの合否を一律に決めることはできません。
4.試料応答から濃度を逆算する
試料の応答が270なら、x = (270 − 21.6) / 99.2 ≈ 2.50 µg/L です。これは測定した溶液の濃度です。試料を希釈した場合や固体試料のmg/kgへ換算する場合は、希釈倍率・試料量・定容量を別途反映します。補正済みの値へ倍率を二重に掛けないよう、装置の出力定義を確認してください。
以下では、計算に使う傾き・切片・相関係数・R²をどう評価するか、図解と実務例で説明します。
Excelで同じ検量線を求めるには?
上の架空のデータを、ExcelのA1:B6へ入力します。A1に「濃度(µg/L)」、B1に「応答」を置き、A2:A6に0、1、2、3、4、B2:B6に22、118、225、316、419を順に入力してください。
↔ 表が見切れている場合は、左右にスクロールしてご覧ください。↔ Scroll the table horizontally to see any hidden columns.
| 求める値 | 入力セル | Excelの式 | この例の結果 |
|---|---|---|---|
| 傾きa | E2 | =SLOPE(B2:B6,A2:A6) | 99.2 |
| 切片b | E3 | =INTERCEPT(B2:B6,A2:A6) | 21.6 |
| 決定係数R² | E4 | =RSQ(B2:B6,A2:A6) | 約0.99956 |
| 試料の応答 | E6 | 数値270を入力 | 270 |
| 試料の濃度(µg/L) | E7 | =(E6-E3)/E2 | 約2.50 |
関数には「応答の範囲、濃度の範囲」の順で指定します。グラフは散布図を選び、横軸を濃度、縦軸を応答にします。上の式は切片を含む、重み付けなしの直線回帰です。原点を通るように強制した回帰や重み付き回帰には、そのまま当てはめられません。
R²の大きさだけで合否は決めません。残差、標準液の逆算濃度、独立QCを確認してください。E7は測定した溶液の濃度であり、希釈倍率や元試料への換算は別に行います。
傾きは「感度」を表すが、高ければよいとは限らない
傾きは、濃度が1単位増えたときに応答がどれだけ増えるかを表します。傾きが大きいほど濃度変化に対する応答変化も大きいため、同じ測定条件内では感度の指標になります。
ただし、傾きの大小だけで分析法の良否は決められません。傾きは次の要因でも変化します。
- 装置感度、光学系、イオンレンズ、検出器状態
- ネブライザー流量、プラズマ条件、導入効率
- 内標準補正の有無と計算方法
- 標準液と試料の酸濃度・マトリックス差
- 測定波長、質量数、積分時間、希釈条件
日常測定では、傾きを単独で見るより、過去の正常な検量線と比べて急に変化していないか、内標準応答やQC試料の結果と整合しているかを確認します。傾きが変化していてもQCが良好なら直ちに不適合とは限らず、逆に傾きが通常どおりでも、低濃度点の誤差や干渉が隠れている場合があります。
傾きを数字で読む:変化量と単位
例えば y = 100x + 20 で、xがµg/L、yがカウントなら、傾き100の単位は カウント/(µg/L) です。濃度が1 µg/L増えると、予測応答は100カウント増えます。2 µg/Lから5 µg/Lへの変化なら、予測応答の差は 100 × (5 − 2) = 300カウント です。切片20は差を取ると消えるので、この計算には入りません。
同じ試料応答を比べるなら切片も確認します。切片がどちらも20で、傾きが100から80へ下がった場合、応答420の逆算濃度は (420 − 20)/100 = 4.0 µg/L から (420 − 20)/80 = 5.0 µg/L へ変わります。これは説明用の比較であり、実際の試料で傾きが変化したときは、標準液の再調製、装置状態、マトリックス、内標準、QCを調べてから結果を判断します。
↔ 表が見切れている場合は、左右にスクロールしてご覧ください。↔ Scroll the table horizontally to see any hidden columns.
| よくある疑問 | 読み方 |
|---|---|
| 傾きが2倍なら感度も2倍? | 同じ応答の定義・濃度単位・測定条件なら、濃度1単位あたりの応答変化は2倍。ただしノイズや定量下限まで2倍改善するとは限りません。 |
| 日によって傾きが違えば結果は全部誤り? | それだけでは判断できません。管理範囲、QC、ブランク、内標準応答、標準液の状態を合わせて確認します。 |
| 傾きを別の検量線と比べてよい? | 応答の単位、内標準補正、濃度単位、測定条件が一致する場合に比較します。ピーク面積と面積比の傾きを直接比べることはできません。 |
切片はブランクや背景の影響を受ける
切片は、回帰式上で濃度が0のときに予測される応答です。理想的には小さい方が分かりやすいものの、切片がゼロでないこと自体を直ちに異常とは判断できません。
切片へ影響する代表例は次のとおりです。
- 試薬、容器、環境からのブランク
- 装置バックグラウンドと背景補正
- メモリーやキャリーオーバー
- 検量線低濃度点のばらつき
- 内標準や補正計算の影響
- 標準液調製時の汚染または希釈誤差
切片の影響は濃度に対して相対的に考えます。高濃度試料では小さく見える切片でも、報告下限付近では換算濃度へ大きく影響する場合があります。
原点を通す回帰は理由なく使わない
切片を強制的にゼロへ固定すると、見た目が単純になり、点によっては相関指標も良く見えます。しかし、実際に存在するブランクや背景を無視すると、範囲全体に系統的な偏りを生じる可能性があります。
原点通過を使う場合は、分析原理、試験法、ブランク補正、バリデーション結果に基づく理由が必要です。「切片が気になるからゼロにする」という後付けの変更は避けます。
相関係数だけで良否を決めない
相関係数 r は、濃度と応答の直線的な関係の強さと向きを示します。決定係数 R² は、回帰モデルが応答の変動をどの程度説明するかを示す指標です。分析ソフトが r と R²のどちらを表示しているかを区別します。
通常の正の検量線では、rは+1に近く、R²は1に近いほど点が直線関係に近く見えます。しかし、次の理由から、これらだけを合否判定に使うのは不十分です。
- 濃度範囲を広くすると、高濃度側の大きな応答差によって高い相関が得られやすい。
- 低濃度側に大きな相対誤差があっても、全体の相関には表れにくい。
- 緩やかな曲線関係でも、範囲によっては非常に高い相関となる。
- 1点の外れ値や高濃度点が、傾きと切片を大きく動かす場合がある。
- 相関が高くても、逆算した標準液濃度やQC試料が許容範囲を外れる場合がある。
ICH(医薬品規制調和国際会議) Q2(R2)は、直線性の評価において、測定点のプロット、相関係数または決定係数、切片、傾きに加え、実測点の回帰線からのずれを評価することが有用だとしています。ICH Q2(R2):Validation of Analytical Procedures
R²を数字で読む:0.999なら誤差0.1%?
いいえ。 切片を含む通常の最小二乗回帰では、R²は「観測された応答の全変動のうち、回帰式で説明される割合」を表します。例えばR² = 0.999を 個々の濃度の誤差が0.1%以内 と読み替えることはできません。残る0.1%は応答の変動についての割合であり、各標準液の逆算濃度の相対誤差ではありません。濃度の誤差は (実測応答 − 切片)/傾き で逆算し、標準液の設定濃度と比べて確認します。
↔ 表が見切れている場合は、左右にスクロールしてご覧ください。↔ Scroll the table horizontally to see any hidden columns.
| R²で分かること | R²だけでは分からないこと |
|---|---|
| 選んだ回帰モデルが、測定範囲全体の応答変動をどの程度説明したか | 最低濃度点の逆算誤差、QCの回収、ブランク、残差の偏り、報告下限の妥当性 |
例えば0、1、10、100 µg/Lのように範囲が広い検量線では、100 µg/L付近の応答差が全体の変動を大きく占めます。1 µg/Lの点が相対的に大きく外れていても、R²は高いままになり得ます。「R²を見る → 残差を見る → 各点を逆算する → 独立QCを確認する」 の順で考えると、数字に隠れた低濃度側の問題を見つけやすくなります。
残差を見ると、直線からのずれ方が分かる
残差は、各濃度点の実測応答と、回帰式から予測される応答との差です。
残差 = 実測応答 − 予測応答
残差プロットでは、横軸を濃度、縦軸を残差として表示します。直線モデルが妥当で、誤差に強い規則性がなければ、残差はゼロの上下へおおむねランダムに分布します。
↔ 表が見切れている場合は、左右にスクロールしてご覧ください。↔ Scroll the table horizontally to see any hidden columns.
| 残差の見え方 | 考えられる状況 | 確認すること |
|---|---|---|
| ゼロの上下へランダム | 直線モデルが範囲内で概ね適合 | 個々の点の大きさ、QC、逆算濃度 |
| U字・逆U字 | 応答が曲がっている | 高濃度側の飽和、範囲短縮、別モデル |
| 濃度とともに広がる | 高濃度ほど分散が大きい | 不等分散、重み付け、反復測定 |
| 1点だけ大きい | 調製・測定・入力の個別異常 | 原液、希釈、容器、測定履歴、再調製 |
| 連続して同じ符号 | 系統的なずれ | モデル、ブランク、標準液、ドリフト |
残差が大きい点を見つけても、相関係数を上げるためだけに削除しません。調製記録、原データ、装置状態、再測定・再調製結果を調べ、事前に定めた除外基準に従います。
濃度範囲は測定目的から決める
検量線の範囲は、装置が応答を返せる範囲ではなく、必要な真度と精度で定量できる範囲として設定します。
下限付近
低濃度側では、ブランク、ノイズ、汚染、切片の影響が相対的に大きくなります。最下点を検量線へ含められたことと、その濃度を定量下限として報告できることは同じではありません。下限では、逆算濃度、真度、精度、ブランクを別途確認します。
上限付近
高濃度側では、検出器の飽和、自己吸収、イオン化影響、内標準との挙動差などにより直線性が失われることがあります。残差が高濃度側で同じ方向へずれる場合は、上限を下げる、試料を希釈する、条件を変更するなどを検討します。
範囲外の外挿
試料応答が最高濃度標準液を超えた場合、回帰式をそのまま延長して定量するのは原則として避けます。適切に希釈して検量線範囲内へ入れ、希釈操作を含む回収と計算を確認します。最低濃度より下の結果も、数値が計算できたことだけを理由に定量値として報告しません。
低濃度側の誤差と重み付け
通常の最小二乗法は、応答方向の残差平方和を最小にします。濃度とともに応答のばらつきが大きくなる場合、高濃度点が回帰へ強く影響し、低濃度側の逆算誤差が大きくなることがあります。
このような不等分散が確認された場合、1/xや1/x²などの重み付けを検討することがあります。ただし、重み付けはR²を改善するために選ぶものではありません。想定する濃度範囲で、標準液の逆算値、独立QC、残差、真度・精度を比較し、分析目的に適した方法を事前に定めます。
ICH Q2(R2)も、濃度によってデータ点の変動が異なる場合には、変換や重み付けを適用できるとしています。一方、単純な直線以外のモデルを用いる場合も、使用範囲でモデルの適合性を示す必要があります。
実務例:同じ相関係数でも確認結果は異なる
↔ 表が見切れている場合は、左右にスクロールしてご覧ください。↔ Scroll the table horizontally to see any hidden columns.
| 状況 | 相関指標 | その他の観察 | 判断の方向 |
|---|---|---|---|
| A | 非常に高い | 残差はランダム、逆算値とQCは良好 | 使用目的に対する基準を満たすか確認 |
| B | 非常に高い | 低濃度点の相対誤差が大きい | 下限、重み付け、標準液配置を再検討 |
| C | 非常に高い | 残差がU字状 | 直線範囲の短縮または妥当な別モデルを検討 |
| D | 基準付近 | 1点だけ大きく外れる | 外れた原因を調査し、規定に従って再評価 |
相関指標が同程度でも、採るべき対応は異なります。検量線の合否判定は、数値一つではなく、使用目的に対して必要な証拠を組み合わせて行います。
計算例:切片を無視すると濃度はどう変わる?
説明用の架空の直線検量線を y = 100x + 20 とします。xの単位はµg/L、yは装置応答の任意単位です。傾きは100応答単位/(µg/L)、切片は20応答単位です。
↔ 表が見切れている場合は、左右にスクロールしてご覧ください。↔ Scroll the table horizontally to see any hidden columns.
| 試料の応答 y | 切片を含めた逆算 x = (y − 20) / 100 | 同じ傾きのまま切片だけを無視した y / 100 | 正しい逆算値に対する差 |
|---|---|---|---|
| 120 | 1.00 µg/L | 1.20 µg/L | +20% |
| 1020 | 10.00 µg/L | 10.20 µg/L | +2% |
この例では、切片を無視した濃度差はどちらも+0.20 µg/Lですが、相対的な影響は低濃度側で大きくなります。実際に原点通過で再回帰すると傾きも変わり得るため、この表は原点通過回帰の結果を示したものではありません。
また、回帰切片を測定ブランクと同一視して機械的に差し引くことはできません。装置の背景補正や回帰モデルを含め、試験法で定めた計算手順を確認します。得られた測定溶液の濃度を元試料のmg/kgへ換算する場合は、ICP-MS/OES 試料濃度換算ツールで試料量・定容量・希釈倍率を確認できます。
検量線の相関係数・切片・傾きについてよくある質問
相関係数が0.999以上なら合格?
一律に合格とはいえません。試験法に相関係数の基準がある場合はそれに従いますが、高い相関でも低濃度点の逆算誤差や系統的な残差が生じることがあります。NISTのモデル評価の解説も、高いR²だけではモデルの適合を保証できず、残差の確認が重要であるとしています。
rとR²は同じ?
同じ指標ではありません。rは直線的な関係の強さと向き、R²はモデルが応答の変動を説明する割合を表します。切片を含む通常の単回帰を、重み付けなしの最小二乗法で当てはめる場合はR² = r²ですが、原点通過回帰や重み付け回帰ではソフトの定義を確認してください。
R²が前回より下がったら、すぐに測定を止める?
まず、試験法・SOPの判定基準に従います。そのうえで標準液の調製・測定順、各点の残差と逆算濃度、ブランク、独立QCを確認します。特定の点の調製ミス、低濃度側のノイズ、装置ドリフトなど原因によって対応は異なります。R²だけを戻すために点を除外したり、濃度範囲を後から変更したりしません。
R²を上げるために標準液の最高濃度を高くしてよい?
それだけを目的に最高濃度を上げるのは避けます。高濃度点が全体の変動を支配するとR²は良く見えても、低濃度側の逆算値が改善したとは限りません。必要な報告範囲と、各濃度で求める真度・精度から標準液範囲を決めます。
切片が負なら異常?
負の切片だけでは判断できません。背景補正や回帰のばらつきで負になる場合もあります。低濃度点、ブランク、逆算値、QCを確認し、理由なく切片をゼロへ固定しないことが大切です。
傾きから検出下限を判断できる?
傾きだけでは判断できません。応答のばらつきやブランクの影響も必要です。検出下限と定量下限の違い、算出後の性能確認はLOD・LOQの解説を参照してください。
日常測定で確認したい項目
- 標準液の濃度、単位、調製履歴は正しいか。
- 検量線範囲が試料濃度とQC濃度を適切に含むか。
- 傾きが過去の正常範囲から大きく変化していないか。
- 切片とブランク応答に説明できない変化がないか。
- rとR²のどちらを評価しているか明確か。
- 各点の残差と逆算濃度に系統的な偏りがないか。
- 独立QC、添加回収、内標準応答が基準を満たすか。
- 範囲外試料を外挿せず、適切に希釈・再測定したか。
- 除外した点がある場合、理由と調査結果を記録したか。
- 回帰モデルと重み付けが検証済みの手順と一致するか。
まとめ
検量線は、相関係数が基準を満たしたかだけでなく、傾き、切片、各濃度点、残差、逆算濃度、QC、濃度範囲を合わせて評価します。
- 傾きは感度の指標だが、測定条件やマトリックスでも変化する。
- 切片はブランクや背景の影響を受け、低濃度側ほど重要になる。
- 相関係数やR²が高くても、曲がりや低濃度側の誤差を見逃す場合がある。
- 残差プロットでは、ずれの大きさだけでなく規則的なパターンを見る。
- 濃度範囲は、真度・精度を確認できた範囲として設定し、外挿を避ける。
- 重み付けや別モデルは、相関指標を良くする目的ではなく、使用範囲で妥当性を示して選ぶ。
IUPAC(国際純正・応用化学連合)の検量線ガイドラインも、検量を試料の量・濃度と測定信号の関係として体系的に扱っています。IUPAC:Guidelines for calibration in analytical chemistry また、Eurachemの分析法バリデーションガイドは、方法が使用目的に適することを中心に、検量関数や残差の評価を含む考え方を示しています。Eurachem:Method validation guides
関連記事:分析法バリデーションは何のため?/ICP-MS/OESの濃度計算で起こりやすいミス
本記事は、一般的な直線検量線の読み方を説明するものです。規制分野や個別試験法では、標準液数、濃度配置、回帰モデル、重み付け、逆算値、QCについて固有の要求があります。実際の合否判定には、適用する規格、ガイドライン、承認済み分析法、試験所のSOPを使用してください。
