12つの測定方法の差と平均値をあらわす変数の作成

グラフに使用する変数の作成

ブランド-アルトマンプロットは、横軸に2つの測定方法の差、縦軸に2つの測定方法の平均値をあらわします。グラフ作成の準備として、2つの測定方法の差の変数「DIFF」と2つの測定方法の平均値の変数「MMEAN」をデータセットに追加します。新しい変数の作成は、変数の計算メニューを利用します。

以上で、2つの測定方法の差をあらわす変数「DIFF」と、2つの測定方法の平均値をあらわす変数「MMEAN」がデータセットに追加されます。

誤差の許容範囲(limits of agreement, LOA)を計算するために、「記述統計」でY軸（DIFF）変数の平均値と標準偏差を確認しておきます。

2つの測定方法の差の平均値が0.688、標準偏差が0.14083であることが確認できました。差の平均値が0に近いほど、2つの測定方法間のバイアスが小さいと解釈されます。例えば、2つの測定方法による値が完全に一致している場合、差の平均値は0になります。この値を中心に、差の平均値±1.96*標準偏差の範囲が誤差の許容範囲となります。

2ブランド-アルトマンプロットの作成

散布図の作成

次に、散布図の機能を利用してブランド-アルトマンプロットを作成します。ここでは「レガシーダイアログ」メニューを使用した手順をご紹介します。

32つの測定方法の差（Y軸）の原点（0）を示す参照線の追加

Y軸の原点

固定誤差の有無を確認するために、Y軸の原点（0）をあらわす参照線を追加します。

42つの測定方法の差（Y軸）の平均値を示す参照線の追加

Y軸の平均値

次に、原点（0）線と区別してグラフを見やすくするために線の種類を変更しておきます。

5許容範囲を示す参照線の追加

下限の参照線

許容範囲の下限は「差の平均値-1.96*差の標準偏差」で計算します。この例では、「-0.20723」（=0.0688-1.96*0.14083）となります。

上限の参照線

許容範囲の上限の計算は「差の平均値+1.96*差の標準偏差」で計算します。この例では「0.344827」（=0.0688+1.96*0.14083）となります。

以上の手順で、ブランド-アルトマンプロットが作成されました。X軸は2つの測定方法の平均値、Y軸は2つの測定方法の差をあらわしています。つまり、X軸が被験者間の測定平均値のばらつき、Y軸が2つの測定方法の差のばらつきとして解釈することができます。

プロットは、Y軸の0よりやや上側にばらついているように見えます。つまり、測定方法Aよりも測定方法Bの方がやや小さく測定されており、固定誤差の存在が示唆されます。プロットはランダムにばらついており、比例誤差はなさそうです。

プロットは誤差の許容範囲（LOA）内でランダムにばらついていますので、2つの測定方法の一致性は高そうです。ただし、固定誤差が存在しますので、新しい測定方法Bを調整して偏りを補正することで、従来の測定方法Aと同等の結果を得ることができそうです。

6ｔ検定による固定誤差の評価方法

固定誤差の評価

固定誤差は、t検定を利用して確認することができます。2つの測定方法の差の平均が0であるかどうかを、有意確率や差の95%信頼区間（Confidence Interval, CI）に基づいて評価します。

1サンプルの検定テーブルに、t検定の結果が出力されます。有意確率はp=0.004（p < 0.05）で有意差が認められます。差の95%信頼区間は0.0237～0.1138であり、区間内に0を含みません。区間内に0を含まない場合、2つの測定値の差が正または負の一方向に分布していることを示すため、固定誤差の存在が示唆されます。なお、元の測定値の変数を使用したt検定は、「記述統計」メニュー >「平均の比較」> 対応のあるサンプルのt検定」から実行可能です。

7回帰分析による比例誤差の評価手順

比例誤差の評価

比例誤差は、回帰分析や相関分析を利用して確認することができます。比例誤差が存在する場合、ブランド-アルトマンプロットはX軸の増減に比例してY軸も増減する分布が示されますので、2つの変数間に有意な関係性が認められるかを確認します。ここでは、回帰分析を利用して確認してみます。

係数テーブルに、回帰係数の有意性検定の結果が出力されます。有意確率はp=0.723で有意差は認められません。つまり、回帰式に含まれる回帰係数が0であることを否定できませんので、比例誤差は存在するとはいえないと解釈することができます。

8比例誤差が存在する場合の一致性の評価

視覚的な評価

誤差の許容範囲（LOA）は測定範囲全体で測定値の差が変化しないことを想定しています。たとえば、測定値によって測定値の差の大きさが変化するような（比例誤差が含まれる）場合は、測定値の大小によってLOAが必要以上に広くなったり狭くなったりしてしまいます。

このような場合は、元の変数を対数変換（ln）し、再度ブランド-アルトマンプロットを作成する方法があります。ただし、対数変換を行う際は、変換後のデータに基づいて分析結果を解釈する必要があり、元の尺度に戻して結果を解釈する際には注意が必要です。対数変換が常に解決策であるわけではないこと、また変換を行う際にはデータの特性をよく理解し、適切な解釈が必要であることを明示することが望ましいです。