玄関雑学の部屋雑学コーナー統計学入門

第13章 用量反応解析

この章では生物検定と薬効評価の違い、用量反応直線またはプロビット変換を利用した平行線検定法、勾配比検定法、そしてD50の推定法について解説します。

13.1 用量反応直線

(1) 生物検定法と薬効評価

薬学分野では生物——ヒト、動物、植物、微生物等——の反応を利用して物質の生物活性の質や量を調べることがあります。 例えばある化学物質の毒性を調べるために、マウスを使ってLD50(対象とする母集団の50%の個体が死亡する量、50%致死量中央致死量、median lethal dose)を推定するのがその代表例です。 これは生物を秤として生物活性を定性・定量する試験法であり、生物(学的)検定法(bioassay、biological assay)と呼ばれます。

それに対して生物の反応を利用して物質の生物活性自体の薬理的な評価を行うこともあります。 例えばある化学物質の薬理的な効果つまり薬効を調べるために、その化学物資を一定量投与した患者と、投与しない患者の反応を比較するのがその代表例です。 これは生物を秤として薬効を調べる試験法であり、薬効評価(pharmacometrics)と呼ばれます。 大雑把に言えば反応を固定して用量の変化を検討するのが生物検定法であり、用量を固定して反応の変化を検討するのが薬効評価と考えれば良いでしょう。

以前は自然界の物質を薬物として利用していて、生物活性を持つ化学物質が明確に同定されていないものが多かったので生物検定法が盛んに利用されました。 しかし近年になり生物活性を持つ化学物質が明確に同定されるようになると、生物検定法の重要性は低くなり、生物活性自体の薬理的な評価の方が重要になってきました。 そのため最近では薬効評価という用語は生物検定法も含めた広義の意味で使われることが多くなり、生物検定法という用語が使われることは少なくなっているようです。

(2) 用量反応直線

生物検定法では薬物の用量と反応の関係を調べ、それに基づいて特定の反応を起こす時の用量を推定します。 そのため薬物の用量反応関係(dose-response relationship)を分析するのに適したデータが得られるように試験をデザインします。 最も単純な試験デザインとして、例えば40匹のマウスを無作為に10匹ずつ4群に分け、それぞれの群にある薬物を1mg/kg(マウスの体重1kgあたり1mgの用量)、2mg/kg、4mg/kg、8mg/kgを投与するというものが考えられます。 このデザインの試験を実施して、マウスの反応を適当な方法で数値化したところ表13.1.1のようになったとします。

表13.1.1 用量反応試験のデータ
群内No.1mg投与群2mg投与群4mg投与群8mg投与群全体
1774848925919 
2776875917969
3816858912924
4808823912902
5773858949880
6707891908956
7774843908964
8819870989948
9852872931983
10816881909985
791586199260943035224
平均値791.5861.9926.0943.0880.6
標準誤差12.66.48.111.210.7

多くの薬物は用量を対数変換した対数用量と反応が正比例するので、薬理学分野では用量のデータとして対数用量を用います。 そのため用量反応試験では表13.1.1のように用量を対数スケールで等間隔になるようにデザインするのが普通です。 そして対数用量と反応の関係を直線で近似して用量反応関係を分析します。 その直線のことを用量反応直線(Dose-Response Line)といい、次のような式で表されます。 (→第14章 薬物動態学)

y = α + βx
y:反応  x=log(D):用量Dのメタメター(metameter)   α:切片  β:回帰係数(傾き)

この式は実際には対数用量反応直線であり、対数用量を元の用量に戻した場合は用量反応曲線(Dose-Response Curve)になります。 しかし普通は用量を広義に解釈して、この式で表される直線のことを用量反応直線と呼んでいます。

一般に対数用量と反応の関係を直線で近似するには回帰分析を用います。 例えば表13.1.1のデータでは4群の平均値に回帰直線をあてはめます。 しかし表13.1.1のデータは一元配置型のデータですから、一元配置分散分析と回帰分析を組み合わせた一元配置型回帰分析を適用することができます。 表13.1.1にその手法を適用すると次のようになります。 (注1)

表13.1.2 用量反応試験データの一元配置型回帰分析表
要因平方和SS自由度φ平均平方和Ms(分散V)分散比F
回帰1344731134473138.057
ズレ7961.2223980.614.08671
用量142434347478.148.7435
残差35065.436974.039 
全体17750039 
○回帰
用量反応直線:y = 802.81 + 172.275x … 図13.1.1の赤色の直線
寄与率:r2 = 0.94411(94.4%)
Fβ = 138.057(p=7.09615×10-14) > F(1,36,0.05) = 4.113 … 有意水準5%で有意
○ズレ(LOF:Lack Of Fit)
FLOF = 4.08671(p=0.0252) > F(2,36,0.05) = 3.259 … 有意水準5%で有意
○用量
FD = 48.7435(p=9.25371×10-13) > F(3,36,0.05) = 2.866 … 有意水準5%で有意
図13.1.1 平均値と用量反応直線

表13.1.2の「用量」は4群の平均値の変動を表し、「回帰」はその変動の中で用量反応直線で説明できる変動を表し、「ズレ」は用量反応直線で説明できない変動つまり非直線的な変動を表します。 そのため回帰の平方和とズレの平方和を合計したものが用量の平方和になり、回帰の自由度とズレの自由度を合計したものが用量の自由度になります。 また「用量」以下は通常の一元配置分散分析表に相当します。

このデータでは回帰の検定結果もズレの検定結果も用量の検定結果も有意になっています。 このことから反応の平均値は用量によって変動する、つまり用量依存性——用量反応関係——があり、その用量反応関係は直線で表される部分と直線では表せない非直線的な部分とがあるということがわかります。 図13.1.1を見ると反応の平均値は1mg投与群から4mg投与群まではほぼ直線的に増加し、8mg投与群だけが少し増加分が減っている、つまり用量が8mgになると反応が少し鈍くなっていることがわかります。 そこで8mg投与群を除外し、あらためて一元配置型回帰分析を適用すると次のようになります。

表13.1.3 用量反応試験データの一元配置型回帰分析表
(3群の場合)
要因平方和SS自由度φ平均平方和Ms(分散V)分散比F
回帰90451.3190451.3102.598
ズレ66.15166.150.075
用量90517.4245258.751.3366
残差23803.427881.607 
全体11432129 
○回帰
用量反応直線:y = 792.55 + 223.4x … 図13.1.1の青色の直線
寄与率:r2 = 0.99927(99.9%)
Fβ = 102.598(p=1.07754×10-10) > F(1,27,0.05) = 4.210 … 有意水準5%で有意
○ズレ(Lack Of Fit:LOF)
FLOF = 0.0751(p=0.7862) < F(1,27,0.05) = 4.210 … 有意水準5%で有意ではない
○用量
FD = 51.3366(p=6.30864×10-10) > F(2,27,0.05) = 3.354 … 有意水準5%で有意

この場合は回帰の検定結果も用量の検定結果も有意ですが、ズレの検定結果が有意ではなくなっています。 このことから用量反応関係は少なくとも1mg〜4mgの範囲では直線的であり、それ以上の用量になると反応が少し鈍くなる可能性があることがわかります。

ちなみに用量を対数変換せず、実用量のまま一元配置型回帰分析を適用すると次のようになります。

表13.1.4 実用量の一元配置型回帰分析表(4群の場合)
要因平方和SS自由度φ平均平方和Ms(分散V)分散比F
回帰1068691106869109.717
ズレ35565.2217782.618.2565
用量142434347478.148.7435
残差35065.436974.039 
全体17750039 
○回帰
用量反応直線:y = 808.3 + 19.28x … 図13.1.2の赤色の直線
寄与率:r2 = 0.75030(75.0%)
Fβ = 109.717(p=1.77361×10-12) > F(1,36,0.05) = 4.113 … 有意水準5%で有意
○ズレ(Lack Of Fit:LOF)
FLOF = 18.2565(p=3.35699×10-6) > F(2,36,0.05) = 3.259 … 有意水準5%で有意
○用量
FD = 48.7435(p=9.25371×10-13) > F(3,36,0.05) = 2.866 … 有意水準5%で有意
表13.1.5 実用量の一元配置型回帰分析表(3群の場合)
要因平方和SS自由度φ平均平方和Ms(分散V)分散比F
回帰86315.3186315.397.9068
ズレ4202.0614202.064.76637
用量90517.4245258.751.3366
残差23803.427881.607 
全体11432129 
○回帰
用量反応直線:y = 759.45 + 43.0071x … 図13.1.2の青色の直線
寄与率:r2 = 0.95358(95.4%)
Fβ = 97.9068(p=1.78066×10-10) > F(1,27,0.05) = 4.210 … 有意水準5%で有意
○ズレ(Lack Of Fit:LOF)
FLOF = 4.76637(p=0.0379) > F(1,27,0.05) = 4.210 … 有意水準5%で有意
○用量
FD = 51.3366(p=6.30864×10-10) > F(2,27,0.05) = 3.354 … 有意水準5%で有意
図13.1.2 平均値と実数用量反応直線

上記の結果と図13.1.2から実用量を用いた時は用量反応直線の当てはまり具合が少し悪く、しかも4群を用いた時も3群を用いた時もズレが大きいことがわかります。 このことから、この薬物の用量反応関係を分析する時は対数用量を用いた方が良いことがわかります。 このように一元元配置型回帰分析は用量反応関係と用量反応直線の検定だけでなくズレの検定もできるので、用量反応直線の妥当性を詳細に検討することができます。

ただし普通の一元配置分散分析と同様に、この手法の検定もたいていは有意性検定になります。 そのため検定結果よりも用量反応直線そのものや寄与率を科学的に検討する方が有意義です。 このデータの場合はどの用量反応直線も寄与率が高く、回帰係数も大きな値のため、各用量反応直線を比較して最も妥当なものを検討をすることに意味があると考えられます。

(3) 用量の逆推定

用量反応直線を利用すると特定の用量の時の反応を予測することができます。 しかし生物検定法では特定の反応の時の用量を用量反応直線を利用して逆推定します。 その方法は標準物質濃度と吸光度の検量線を利用して、未知物質の濃度を逆推定する一般的な方法と同じです。 つまり生物の反応を吸光度のように考え、用量反応直線を検量線にして薬物の生物活性の量を測定するわけです。 例えばある薬物をマウスに投与して、その反応が850だったとします。 その時の用量とその95%信頼区間を、3群を用いた用量反応直線を利用して求めると次のようになります。 (注2)

用量反応直線:y = 792.55 + 223.4x
y=850の時の用量 = 1.80785mg/kg (x=log(用量)=0.257162)
95%信頼区間:下限=0.93767mg/kg (xL=-0.0279502)  上限=3.45558mg/kg (xU=0.538521)

また例えば反応の数値には±50程度の測定誤差があるので、反応の数値が50以上変化しなければ薬理学的には意義がないことが判明しているとします。 この時、用量反応直線を利用すれば、用量が何倍になれば薬理学的に意義のある反応の変化が得られるかを推定することができます。 上の用量反応直線の回帰係数223.4は、x=log(用量)が「1」変化した時つまり用量が10倍になった時の反応yの変化量を表します。 このことから反応yが50変化する時の用量の倍率を推定すると次のようになります。

50 = 223.4x
x= 50 223.4 =0.2238138
∴100.2238138 = 1.674225

この計算結果から用量が約1.7倍になると反応が50以上増加し、薬理学的に意義のある変化になることがわかります。 したがって表13.1.1のように用量を2倍にした4群の反応は、それぞれ薬理学的に有意義に増加していることになります。

(4) 多重比較と最大対比法

用量反応関係を検討する時、一元配置型回帰分析を用いずに一元配置分散分析だけを用いたり、各用量群の平均値を多重比較したりすることがあります。 例えば表13.1.1のデータに一元配置分散分析とダネット型多重比較を適用すると次のようになります。 (注3) (→4.1 多標本の計量値)

表13.1.6 用量反応試験データの一元配置分散分析表
要因平方和SS自由度φ平均平方和Ms(分散V)分散比F
用量142434347478.148.7435
残差35065.436974.039 
全体17750039 
○用量
FD = 48.7435(p=9.25371×10-13) > F(3,36,0.05) = 2.866 … 有意水準5%で有意
○ダネット型多重比較:1mg投与群を対照とした場合
1mg投与群対2mg投与群:d12 = 5.044(p=3.8459×10-5) > d(4,36,0.05) = 2.452 … 有意水準5%で有意
1mg投与群対4mg投与群:d14 = 9.636(p=3.37143×10-9) > d(4,36,0.05) = 2.452 … 有意水準5%で有意
1mg投与群対8mg投与群:d18 = 10.854(p=2.22045×10-16) > d(4,36,0.05) = 2.452 … 有意水準5%で有意

上記のように一元配置分散分析表は表13.1.2の「用量」以下と同じになり、用量の検定結果も同じになります。 そして多重比較では1mg投与群の平均値と他の3群の平均値の差は全て有意水準5%で有意になります。

また用量反応関係を検討する場合、普通のダネット型多重比較ではなく、逐次的なダネット型多重比較であるウィリアムズの多重比較(Williams's multiple comparison)という手法を用いることもあります。 この手法は2mg投与群〜8mg投与群の平均値は必ず増加または横ばいになり、低下することはないという前提で、用量の多い群から順番にダネット型多重比較を行っていきます。 そして有意ではなくなったところで、それ以下の用量は全て有意ではないという結論を採用して検定を終了する手法です。 検定を途中で終了するので普通のダネット型多重比較よりも検定回数が少なくなり、その分だけ検定効率が高くなります。

ちなみに表13.1.1のデータは普通のダネット型多重比較でも全ての比較が有意になっているので、ウィリアムズの多重比較を適用しても結果は変わりません。

ダネット型の多重比較ではなく、シェッフェ型の多重比較の中の一般対比という手法を応用して、用量反応関係が特定のパターンかどうかを検討する最大対比法(maximum contrast method)という手法もあります。 例えば表13.1.1のデータについて1mg投与群から8mg投与群まで直線的に増加すると仮定した直線的増加パターンと、1mg投与群から4mg投与群までは直線的に増加し、8mg投与群は4mg投与群と同じと仮定した頭打ちパターンの2種類について検定すると次のようになります。 (注4)

(1) 直線的増加パターン
係数:1mg投与群 = -3 2mg投与群 = -1 4mg投与群 = +1 8mg投与群 = +3
寄与率:r2 = 0.94411(94.4%)
Fo = 46.019(p=2.10214×10-12) > F(3,36,0.05) = 2.866 … 有意水準5%で有意
(2) 頭打ちパターン
係数:1mg投与群 = -5 2mg投与群 = -1 4mg投与群 = +3 8mg投与群 = +3
寄与率:r2 = 0.98979(99.0%)
Fo = 48.246(p=1.07175×10-12) > F(3,36,0.05) = 2.866 … 有意水準5%で有意
図13.1.3 2種類の用量反応パターン

上の結果中の寄与率は一元配置型回帰分析の寄与率と同様の値であり、平均値の変動のうち係数で表されるパターンによって説明できる割合つまりパターンの当てはまり具合を表します。 上記の結果を見るとどちらの検定結果も有意ですが、直線的増加パターンよりも頭打ちパターンの方が寄与率が高く、用量反応関係はどちらかと言えば後者のパターンの方に近いことがわかります。

多重比較法や最大対比法は用量反応関係を定性的なパターンとして検討するための手法であり、用量と反応の計量的な関係を検討するための手法ではありません。 実際、これらの手法は用量の値を計算に用いず、用量の順序という情報しか用いていないので、対数用量でも実用量でも全く同じ結果になります。 そのため用量反応直線のように特定の反応の時の用量を逆推定したり、薬理学的に有意義な反応をする時の用量の倍率を求めたりすることはできません。

(5) ノンパラメトリック手法

多重比較法や最大対比法は用量を順序尺度扱いしますが、用量だけでなく反応も順序尺度扱いするヨンクヒール・タプストラ(Jonckheere-Terpstra)検定という手法もあります。 この手法はマン・ホイットニィのU検定を利用したノンパラメトリック手法であり、用量が多くなるほど反応の順位平均値が増加するかどうかを検定することによって用量反応関係を定性的に検討します。

しかし用量も反応も順序尺度扱いするのなら、むしろスピアマンの順位相関係数を利用して用量の順位と反応の順位の間に相関性があるかどうかを検討する方がわかりやすいでしょう。 例えば表13.1.1のデータにヨンクヒール・タプストラ検定とスピアマンの順位相関係数を適用すると次のようになります。 (注5) (→3.4 2標本の計数値5.3 計数値の相関分析と回帰分析)

○ヨンクヒール・タプストラ検定:JT = 560(93.3333%)  Tmax = 600
J=zo = 6.277(p=3.45172×10-10) > t(∞,0.05) = 1.96 … 有意水準5%で有意
○順位相関係数:rs = 0.887399  寄与率:rs2 = 0.787476
95%信頼区間 下限:ρsL = 0.795903 上限:ρsL = 0.939264
検定:zo = 5.5418(p=2.99374×10-8) > t(∞,0.05) = 1.96 … 有意水準5%で有意

ヨンクヒール・タプストラ検定の統計量JTはマン・ホイットニィのU検定のU統計量を拡張した指標です。 そしてJTの理論的最大値がTmaxであり、JTをTmaxで割って100を掛けた値がU検定の勝率に相当します。 この勝率が50%よりも大きければ、用量が多くなるほど反応の順位平均値が増加することを表します。 このデータの場合、JTの割合が約93%とかなり大きな値です。 しかしこの値から用量反応関係の大きさを具体的に解釈するのは難しいでしょう。

それに対してスピアマンの順位相関係数は解釈が容易です。 このデータの場合、順位相関係数が約0.9もあるので、用量の順序と反応の順序の間に高い相関性があることがすぐわかります。 そして上記のように、ヨンクヒール・タプストラ検定とスピアマンの順位相関係数の検定結果はよく似ています。 このことから、ノンパラメトリックな用量反応関係を検討するにはヨンクヒール・タプストラ検定よりもスピアマンの順位相関係数の方が便利なことがわかると思います。

ただし生物検定法では用量反応関係が直線的な部分を検量線のように利用して生物活性を定量し、薬剤の常用量も用量反応関係が直線的な範囲内で設定するのが普通です。 そのため多重比較法や最大対比法やノンパラ手法を用いて用量反応関係を検討しても、最終的には用量反応直線を求める必要があります。 したがって用量反応解析としては一元配置型回帰分析によって用量反応直線を求めるのが合理的であり、多重比較法や最大対比法やノンパラ手法を用いるのはお勧めできません。


(注1) 表13.1.1を一般化すると次のようになります。

表13.1.7 一元配置型回帰分析の
一般的データ
群内No.対数用量全体
x1xixa
1y11yi1ya1
::::
jy1jyijyaj
::::
riy1・r1yi・riya・ra
T1.Ti.Ta.TT
平均値m1.mi.ma.mT

一元配置型回帰分析ではデータyijを次のように2通りに分解して考えます。

基本式:
:直線回帰式による推定値   xi = log(Di):用量Diのメタメター

この基本式に対応する平方和と自由度と分散、そして直線回帰式を求めると次のようになります。

総例数:
全体:  φT = N - 1   
用量:  φD = a - 1   
回帰:   
    
 φβ = 1        a = mT - b mx
ズレ:   φLOF = φD - φβ=a - 2   
残差:   φR = φT - φD=N - a   
寄与率:

上式からわかるように、回帰の平方和は個々のデータyiiを使わずに用量xiと用量群ごとの平均値mi.に例数riかけて計算することができます。 つまり一元配置型回帰分析では用量xiを説明変数に、用量群の平均値mi.を目的変数にし、用量群の例数の平方根√(ri)を重みとした重み付け最小2乗法によって直線回帰式を求めているわけです。

通常、平均値を用いる重み付け最小2乗法では平均値の標準誤差の逆数(1/SE)を重みにします。 しかし一元配置分散分析は各群の分散は等しいと仮定して計算するので、一元配置型回帰分析でも各用量群の分散は等しいと仮定して計算します。 そのため(1/SE)を重みとした重み付け最小2乗法は、√nを重みとした重み付け最小2乗法同じ結果になるのです。

これらの平方和の関係をグラフで表すと図13.1.4のようになります。 そして平方和と自由度と分散を表13.1.8のような一元配置型回帰分析表にまとめます。

図13.1.4 一元配置型回帰分析のグラフ的解釈
表13.1.8 一元配置型回帰分析表
要因平方和SS自由度φ平均平方和Ms(分散V)分散比F
回帰SβφβVβFβ=Vβ/VR
ズレSLOFφLOFVLOFFLOF=VLOF/VR
用量SDφDVDFD=VD/VR
残差SRφRVR 
全体STφT 
○回帰の検定
帰無仮説 H0:用量反応直線の傾きが0である (β=0)
Fβ > F(φβR,α)の時、有意水準100α%で有意
○ズレの検定
帰無仮説 H0:用量反応直線と平均値の変動が一致している
FLOF > F(φLOFR,α)の時、有意水準100α%で有意
○用量の検定
帰無仮説 H0:平均値が用量によって変動しない (m1.=…=ma.=mT)
FD > F(φDR,α)の時、有意水準100α%で有意

表13.1.1のデータについて実際に計算してみましょう。

N = 40
ST = 7742 + 7762 + … + 9852 - 40×880.62 = 31195754 - 31018254 = 177500   φT = 40 - 1 = 39
SD = 10×791.52 + 10×861.92 + 10×9262 + 10×9432 - 40×880.62 = 31160688 - 31018254 = 142434
φD = 4 - 1 = 3  

Sxx = 10×{log(1)2 + log(2)2 + log(4)2 + log(8)2} - 40×0.4515452 = 12.68667 - 8.155715 = 4.53095
Sxy = 10×{log(1)×791.5 + log(2)×861.9 + log(4)×926 + log(8)×943} - 40×0.451545×880.6 = 16685.8 - 15905.2 = 780.571
  φβ = 1   Vβ = Sβ = 134473
  a = 880.6 - 172.275×0.451545 = 802.81
SLOF = 142434 - 134473=7961  φLOF = 3 - 1 = 2   
SR = 177500 - 142434 = 35066  φR = 40 - 4 = 36   

これらの統計量を用いて表13.1.2の一元配置型回帰分析表を作成し、用量反応直線を求めて回帰、ズレ、用量の検定を行うことができます。

実際の薬物の用量反応関係はシグモイド曲線(S字状曲線)になるのが普通です。 そこでシグモイド曲線の中でも数学的な取り扱いが比較的容易なロジスティック曲線で用量反応関係を近似することが考えられます。 事実、第4節で説明するプロビット分析はロジスティック曲線を用いたロジット分析で近似することがあります。 (→9.6 ロジスティック曲線 (注1))

ロジスティック関数:
:y0のロジット=初期ロジット
y:時点tにおける人口  M:人口最大値  c:比例定数  y0:t=0における人口初期値
図9.6.2 ロジスティック曲線

このロジスティック関数においてyを反応、tを対数用量x、Mを反応の最大値、cを比例対数としたものが用量反応曲線になります。 そして反応yをロジット変換すると用量ロジット直線になります。

用量反応曲線:
用量ロジット直線:ly = ly0 + c・x
反応yのロジット:
y:反応  M:反応最大値  x=log(D):用量Dのメタメター  c:比例定数
ly0:x=0におけるロジット初期値  y0:x=0における反応初期値

この用量ロジット直線を利用すると、用量反応直線と同様の一元配置型回帰分析による用量反応解析を行うことができます。 ただしyをロジット変換するには何らかの方法でMを推定する必要があります。 例えば非線形最小2乗法によってロジスティック関数のパラメータM、c、y0を求め、そのMを用いるという方法が考えられます。

そこで表13.1.1のデータに非線形最小2乗法を適用してロジスティック関数のパラメータを求めたところ、M=1007.99という推定値が得られました。 このMを用いて用量ロジット直線を利用した用量反応解析を行うと次のような結果になります。

表13.1.9 用量ロジット直線の一元配置型回帰分析表
要因平方和SS自由度φ平均平方和Ms(分散V)分散比F
回帰13.8126113.812672.8475
ズレ0.28459220.1422960.750466
用量14.097234.6990724.7828
残差6.82596360.18961 
全体20.923239 
○回帰
用量ロジット直線:ly = 1.31684 + 1.746x
用量反応曲線:
M = 1007.99  ly0 = 1.31648  c = 1.746  y0 = 794.959
寄与率:r2 = 0.979812(98.0%)
Fβ = 72.8475(p=3.57443×10-10) > F(1,36,0.05) = 4.113 … 有意水準5%で有意
○ズレ(LOF:Lack Of Fit)
FLOF = 0.750466(p=0.479389) < F(2,36,0.05) = 3.259 … 有意水準5%で有意ではない
○用量
FD = 24.7828(p=7.11738×10-9) > F(3,36,0.05) = 2.866 … 有意水準5%で有意
○用量の逆推定
y=850の時の用量 = 1.62008mg/kg(x=log(D)=0.209536)
95%信頼区間(g=0とした近似値):下限=0.494582mg/kg(xL=-0.305762) 上限=5.30683(xU=0.724835)
図13.1.5 用量反応直線と用量反応曲線 図13.1.6 用量ロジット直線

上の結果と図13.1.5を見ると、用量反応直線よりも用量反応曲線の方が平均値により近似していて寄与率が高く、用量の逆推定値もより正確な値であることがわかります。 このロジスティック曲線を利用した用量反応解析は用量群が4個未満だったり、最大用量群の反応平均値が反応最大値Mよりもかなり小さい時(だいたいM×0.8未満)はうまく計算できません。 しかし実際の用量反応関係を用量反応直線よりもうまく近似できるので、用量群の数が多い時には非常に有用な方法です。

(注2) 用量反応直線を利用して特定の反応の時の用量を逆推定する方法は次のとおりです。 この方法は標準物質濃度と吸光度の検量線を利用して、未知物質の吸光度からその濃度を逆推定する方法と同じです。

用量反応直線:y = a + bx
反応y0の時のxをx0とすると:y0 = a + bx0
x=log(D)の時:D0 = 10x0

上記のように、逆推定した用量x0は(y0-a)とbの比になります。 この場合、aにもbにも誤差があるのでx0の信頼区間は少し複雑になります。 一般に2つの確率変数の比の100(1-α)%信頼区間は、フィーラー(Fieller)の式と呼ばれる次のような式で求めます。

確率変数xとyの比をとすると

ただし
  Vxx = V(x)   Vyy = V(y)  Vxy = C(x,y):xとyの共分散
t=t(φ,α):自由度φのt分布における100α%点

このフィーラーの式を利用してx0の100(1-α)%信頼区間を求めると次のようになります。 (→5.5 各手法の相互関係 (注3))

a = mT - bmx より
x = x0 - mx、y = y0 - mT と置き直すと

r0回の観測で平均値y0を得た時はV(y0)=VR/r0になる。 上式はr0=1の時のもの
  C(y,b)=0 より


ただし    t=t(φR,α)
※g<0.1の時はg=0とした次のような近似式も用いられる。

用量反応直線に関する反応yの100(1-α)%予測限界は次のようになります。 この予測限界曲線とy=y0の交点つまりy=y0の時のxを求めてみましょう。

yの100(1-α)%予測限界:

図5.5.4 信頼区間・許容限界・予測限界

このように2本の予測限界曲線上でy=y0になる時のxは、g=0とした時のx0の100(1-α)%信頼区間に相当します。 つまり反応y0の時の用量x0はy=y0と用量反応直線の交点の用量になり、その信頼区間はy=y0と2本の予測限界曲線の交点付近の用量になるわけです。

表13.1.1の1mg〜4mgのデータだけを使った時の用量反応直線を利用して、反応が850の時の用量とその95%信頼区間を計算してみましょう。

用量反応直線:y = 792.55 + 223.4x
→ D0 = 100.257162 = 1.80785
VR=881.607  t(27,0.05)=2.052  mT=859.8   mx=0.301030  Sxx=1.81238
  
x0 - gmx = 0.257162 - 0.0410338×0.30103=0.24481
  


D0LU = 10-0.02795〜100.53852 = 0.93767〜3.45558
※g=0とした時の近似解は次のようになる。

D0LU = 10-0.020195〜100.534519 = 0.954564〜3.42388

(注3) 表13.1.1のデータに一元配置型分散分析を適用した結果は、(注1)で計算した一元配置型回帰分析の用量以下の分散分析表と用量の検定結果と同じにものになります。 この結果を用いて1mg投与群を対照としたダネット型多重比較を計算してみましょう。 (→4.1 多標本の計量値 (注4))

○1mg投与群対2mg投与群

|d12|=5.044(p=0.00004)>d(4,36,0.05)=2.452 … 有意水準5%で有意
○1mg投与群対4mg投与群

|d14|=9.363(p=3.845×10-9)>d(4,36,0.05)=2.452 … 有意水準5%で有意
○1mg投与群対8mg投与群

|d18|=10.854(p=4.011×10-10)>d(4,36,0.05)=2.452 … 有意水準5%で有意

ウィリアムズの多重比較では、最初に用量が最も多い群と対照群についてダネット型の多重比較と同じ検定統計量dを計算します。 それをダネット型数表を少し変更したウィリアムズ型数表の100α%点の値w(a,φR,α)と比較し、d<w(a,φR,α)なら全ての用量が有意水準αで有意ではないとして検定を終了します。

d>w(a,φR,α)なら用量が最も多い群と2番目に用量が多い群を一緒にして例数と平均値を計算し、これと対照群について検定統計量dを計算します。 それをw(a-1,φR,α)と比較し、d≦w(a-1,φR,α)なら2番目に用量が多い群以下の用量は全て有意水準αで有意ではないとして検定を終了します。 d>w(a-1,φR,α)なら次は3番目に用量が多い群を一緒にして同様の検定を行い、これを対照群以外の群について繰り返します。

この方法は途中で検定を中止するため、普通のダネット型検定よりも効率的です。 しかし「用量が多くなるほど反応が増加するか横ばいである」という前提が必要なため、ダネット型よりも適用範囲が狭くなります。 またダネット型多重比較もウィリアムズの多重比較も単に用量反応があるかどうかを定性的に検定するための手法であり、用量と反応の計量的な関係を検討するための手法ではありません。 そのため用量反応解析では使い道はあまりありません。

(注4) 最大対比法は一般対比を利用して次のように計算します。 (→4.1 多標本の計量値 (注4))

Ci:用量反応パターンを表す群の係数、ただし
一般対比の検定:>F(φDR,α)の時、有意水準100α%で有意
※1種類の用量反応パターンだけを検定する時は、次のようにフィッシャー型の検定を行う。
>F(1,φR,α)の時、有意水準100α%で有意

表13.1.7において用量xiの代わりに係数Ciを用いて一元配置型回帰分析を行うと、回帰の検定が最大対比法のフィッシャー型検定になり、寄与率が最大対比法の用量反応パターンの寄与率になります。 表13.1.1のデータを用いて、直線的増加パターンと頭打ちパターンについて実際に計算してみましょう。

(1) 直線的増加パターン
C1=-3  C2=-1  C3=+1   C4=+3

Fo=46.019>F(3,36,0.05)=2.866
Sxx=10×{(-3)2 + (-1)2 + 12+32} - 40×02=200
Sxy=10×{(-3)×791.5 + (-1)×861.9 + 1×926 + 3×943} - 40×0×880.6=5186
  
※xiが等間隔なので、回帰平方和と寄与率は用量反応直線と同じ値になる。
(2) 頭打ちパターン
C1=-5  C2=-1  C3=+3   C4=+3

Fo=48.246>F(3,36,0.05)=2.866
Sxx=10×{(-5)2 + (-1)2 + 32 + 32} - 40×02=440
Sxy=10×{(-5)×791.5 + (-1)×861.9 + 3×926 + 3×943} - 40×0×880.6=7876
  

最大対比法は用量反応関係を定性的なパターンとして検討するための手法であり、用量と反応の計量的な関係を検討するための手法ではありません。 そのため用量反応解析では補助的な手法になります。

(注5) ヨンクヒール・タプストラ検定は、U検定のU統計量を利用してノンパラメトリックな用量反応関係を検討する手法です。 U統計量は群1のデータと群2のデータを総当たりで比較し、大きい方を勝ち(Upper)とした時の群2の勝ち数を表します。 そして試合数つまりUの最大値は(群1のデータ数×群2のデータ数)であり、Uをその最大値で割った値は勝率に相当します。 (→3.4 2標本の計数値 (注3))

そのため表13.1.7のデータについて用量の小さい方を群1、用量の大きい方を群2としてU検定を適用すると、用量が増えると反応が増加する(または減少する)時は群2の勝率は50%よりも大きく(または小さく)なります。 そこで表13.1.7のデータについてx1群 vs x2群、x1群 vs x3群、…、x1群 vs xa群、x2群 vs x3群、…、x2群 vs xa群、…、xa-1群 vs xa群についてU検定を適用し、それらのU統計量を合計して勝率を求め、それが50%よりも大きいかどうかを検定すれば用量が増えるほど反応が増加するかどうかを定性的に検定できるはずです。

この原理に基づいて、ヨンクヒール・タプストラ検定は次のように計算します。

xi群とxj群のU統計量(xi<xj):Uij      Uij max:Uijの最大値
     Tmax:JTの最大値


  
gk:同位のデータ数。 全体データをこみにして順位を付け、同位のデータがt組あった時のk番目のデータ数
  |zo|>t(∞,α)の時、有意水準100α%で有意
JTの割合(勝率):   勝率>0.5なら用量が増えると反応は増加、勝率<0なら用量が増えると反応は減少