健康と生活習慣に関する重回帰分析

高木 廣文

医学のあゆみ,174(3), 217-223, 1995



■リード:
多変量解析の最も基本的な手法である重回帰分析の解説を行う。とくに,生活習慣調査という実例を用いて,その手順,用語等の解説を行う。

■キーワード:
1.重回帰分析  2.生活習慣  3.変数選択  4.多変量解析


トップページ
神戸市看護大学

最終更新日: 6 April 2023
無断転載の禁止 (C)高木廣文

1.はじめに

 身長が高ければ,体重も重くなる。数学ができる人は,英語もよくできる。塩分や肉・油脂などの取りすぎは,高血圧の原因の1つになる。などなど,我々の周囲には,相互に何らかの関係がある事柄は結構多いものである。

 ゴールトン(F. Galton, 1822-1911)は,親子の身長についての研究から,「回帰」や「相関」という概念を,はじめて生物統計学の分野に導入した1,2)
 もともと回帰という概念は,大きい両親からの子供の身長も,小さな両親からの子供の身長も,全体の平均値に向かうという,平均化の傾向を示すために用いられたものである。現在では,ある変数について,他のいくつかの変数を用いて予測などを行うための統計的な手法を「回帰分析」と呼んでいる。

 回帰分析では,予測すべき変数のことを「基準変数」と呼び,予測のために使用する変数を「説明変数」と呼ぶ。説明変数が1つの場合には,「単回帰分析」と呼ばれ,説明変数が2つ以上の場合には「重回帰分析」と呼ばれる。

 今回は,実際に地域調査で得られた生活習慣調査データと血液検査値を用いて,重回帰分析を中心にその概略について解説する。


2.重回帰分析とは

 まず,重回帰分析を行う主な目的には,大きく分けて,

(1) 特定の変数の予測式を求める。
(2) 特定の変数に対する各説明変数の影響の程度を検討する。

という2通りの場合がある。

 例えば,父親と母親の身長から,子供の身長を予測する式を求める場合が(1)に該当する。また,血圧に対して,年齢,食塩摂取量,脂肪摂取量等がどのような影響を及ぼすのかを調べたいのならば,(2)の場合ということになる。実際の医学的な研究では,予測式を求めるというよりも,変数間の関係の分析を行うために,重回帰分析を用いるのが普通であろう。

 重回帰分析では,基準変数と説明変数を用いて,

基準変数の予測値= (重み1)×(説明変数1)+(重み2)×(説明変数2)
+・・・+(重みP)×(説明変数P)+定数

のような式を考える。このような式は,説明変数がP個ある場合の「重回帰式」と呼ばれる。重回帰式は各説明変数に適当な重みをかけた値の合計に,定数を加えたものである。

 さて,重回帰分析では,上記の式中の重みを求める必要がある。この重みは「偏回帰係数」と呼ばれている。説明変数の基準変数への影響の有無は,この偏回帰係数が0か否かの検定を行うことで,調べることができる。

 また,重回帰式から求められる基準変数の予測値と実際値との相関係数は,「重相関係数」と呼ばれる。重相関係数の2乗は「多重決定係数」と呼ばれ,説明変数を用いて,どの程度の予測が可能かを示す指標になる。


3.血液検査値と生活習慣データの解析

 ここでは,都市,農村,漁村に在住する主婦の314人について,我々が1991年に行った調査データをもとに,重回帰分析の例を示そう。

 長期にわたる日常の生活習慣が,ガンや循環器障害などの多くの成人病の原因になることに異論はないだろう。我々はさまざまな生活習慣を測定するために,質問紙を用いた生活習慣調査票を開発してきた3-5)。ここでは,表1に示したような生活習慣尺度を用いることにする。生活習慣を表す各尺度は,それぞれ6質問項目からなり,回答によって,0,1,2点を与え,その合計点を各尺度の得点としている。この得点が高いほど,その生活習慣の傾向があることになる。このように,質問紙を用いて目的にそった測定を行う方法は「尺度構成法」として知られている。作成した尺度については,その信頼性と妥当性についての検討が必要であるが,ここでは詳細は文献を参照してもらうことにする3-5)

表1.生活習慣尺度とその意味
生活習慣尺度尺度の意味
1)肉・油脂肉類,脂肪分の多い食事を摂る傾向
2)洋風の食事朝食がパン・牛乳などの洋風の傾向
3)高塩分漬け物など,塩分を多く摂る傾向
4)糖分菓子類,ジュース類などを多く摂る傾向
5)食事規則性食事時間の規則的な傾向
6)料理への進取性新料理への関心,健康のための料理の工夫
7)娯楽娯楽番組,観劇,ショッピングなど
8)教養教養番組,読書などの教養ある行動
9)社会奉仕ボランティアなど社会への奉仕活動
10)義理人情義理人情の重視,礼儀,恩義を重んじる
11)経済型実用本位,無駄を省いた現実的生き方
12)伝統型慣習に従い伝統的生き方を重んじる
13)清潔歯磨き等の清潔さを心がける傾向
14)妻主導型家庭内での決定権を妻が持つ傾向
15)運動実施日常生活での運動を行う傾向
16)疾病頻度病気がちの傾向
17)多愁訴身体的な訴えの多い傾向
18)情緒不安定精神状態が不安定になりやすい傾向
19)外向性心が自己だけでなく外にも向いている性格
20)共感性感受性豊かで他人の気持ちが理解できる性格
21)自発性積極的で自発的に行動する性格
22)遺伝的健康観健康は遺伝で決まると考える傾向

 生活習慣としては,この他に飲酒や喫煙についても調べているが,調査対象が女性のみであり,頻度も少ないため,ここでは取り上げない。実際の地域での調査では,対象者の同意が得られた場合には血液検査を行った。尿素窒素,クレアチニン,尿酸,血清総蛋白,アルブミン,コレステロール,GOT,GPT,γ−GTPなどであるが,ここではグリコヘモグロビン(HbAlc)について検討する。

 本来,生活習慣を調べる意味は,将来での特定の疾患の罹患や死亡が,現在の各生活習慣からどのような影響を受けているのかを調べることにある。このような点を検討するためには,現時点から今後何年も追跡調査を実施しなければならない。そのようなコホート調査を行えない場合には,特定時点での断面調査ということになる。断面調査では,因果関係の解明は行えない。本例も断面調査であり,現在の相関関係を検討している点をあらかじめ注意しておきたい。

 さて,基準変数はHbAlcとしよう。HbAlcは糖尿病の指標として有用なのではないかと考えられてきている。説明変数としては,表1に示した22個の生活習慣尺度とする。さらに,年齢と肥満度の影響は大きいものと考えられるので,これも説明変数に加える。ここで,肥満度はBMI(体重/身長2,単位はKg/m2)を用いる。

表2.重回帰分析に用いた変数の基本統計
変数名平均値 標準偏差最小最大
0)HbAlc5.80.534.79.7
1)肉・油脂4.71.950.011.0
2)洋風の食事5.32.180.012.0
3)高塩分4.21.940.09.0
4)糖分3.91.990.010.0
5)食事規則性8.41.902.012.0
6)料理への進取性6.32.480.012.0
7)娯楽4.41.950.011.0
8)教養5.62.240.012.0
9)社会奉仕5.12.950.012.0
10)義理人情8.22.163.012.0
11)経済型7.21.792.011.0
12)伝統型8.81.903.012.0
13)清潔9.81.704.012.0
14)妻主導型7.12.290.012.0
15)運動実施2.92.450.012.0
16)疾病頻度3.32.690.012.0
17)多愁訴3.92.670.012.0
18)情緒不安定4.32.710.012.0
19)外向性5.01.861.09.0
20)共感性7.82.102.012.0
21)自発性6.32.510.012.0
22)遺伝的健康観6.72.540.012.0
23)年齢48.58.8324.079.0
24)BMI22.92.9816.535.0

 表2は,重回帰分析に用いる全変数の基本的な統計量を求めたものである。 すべてのデータが完備していたのは247例であった。 生活習慣尺度は0〜12点が取り得る範囲であるが,尺度によってはそれよりも小さな範囲内におさまっていた。 年齢は24〜79歳であり,48.5歳が平均値であった。 肥満度をあらわすBMIの値も16.5〜35.0,平均22.9となっていた。

表3.HbAlcと説明変数間の相関係数(Cor)
変数名 Cor 変数名 Cor 変数名 Cor
1)肉・油脂0.0219)社会奉仕-0.128**17)多愁訴0.056
2)洋風の食事-0.125*10)義理人情0.171**18)情緒不安定0.065
3)高塩分-0.02011)経済型0.08719)外向性0.060
4)糖分-0.10512)伝統型0.06020)共感性0.061
5)食事規則性-0.02713)清潔-0.05821)自発性0.048
6)料理への進取性-0.06614)妻主導型-0.02822)遺伝的健康観0.246☆☆
7)娯楽0.02115)運動実施-0.03323)年齢0.296☆☆
8)教養-0.10016)疾病頻度0.309☆☆24)BMI0.216
* p<0.05 ** p<0.01 p<0.001 ☆☆p<0.0001

 表3はグリコヘモグロビンと,他の24変数との単相関係数(ピアソンの積率相関係数)を示したものである。また,いくつかの相関係数の横についている*印などは,相関係数のt検定によるp値(有意確率)によって,4段階で印をつけたものである。なお,検定は「母集団における相関係数(母相関係数)は0」という帰無仮説についてのもので,「無相関の検定」と呼ばれるものである。検定のp値が小さい(通常5%以下)母相関係数は0ではないと考えてよい。すなわち,変数間には何らかの関係があることになる。表3の結果は,HbAlcは年齢と肥満度とは正の相関があることを示している。また同様に,生活習慣のうち,疾病頻度,遺伝的健康観,義理人情とは正の相関が認められた。一方,洋風の食事と社会奉仕の各尺度とは負の相関が認められた。

 上記の結果の解釈としては,年齢の影響が大きいのではないかと考えられる。すなわち,洋風の食事は年齢が若いほど得点が高いし,義理人情の得点はその逆になっている。また,高齢者では病気がちになるし,そのような状態では,グリコヘモグロビンが高値になるのではないかと考えられる。身近に糖尿病患者などがいると,健康が遺伝によって決まると考える傾向があり,実際にHbAlcが高くなっているのだろうか。

 このような解釈はそれらしく思えるかもしれないが,個々の2変数間の相関係数をみて考えついた推測でしかない。実際には,他の変数との関係を考慮して統計的解析をするには,重回帰分析を行う必要がある。

表4.全ての説明変数を用いた場合の重回帰分析(基準変数 HbAlc)
変 数 偏回帰係数(標準誤差) 標準偏
回帰係数
F値 (p値) 偏相関
係 数
1)肉・油脂 0.027461( 0.018870) 0.1022 2.118(0.1470) 0.097
2)洋風の食事 -0.010671( 0.017430) -0.0443 0.375(0.5410) -0.041
3)高塩分-0.022220( 0.018296)-0.08201.475(0.2259)-0.081
4)糖分-0.042046( 0.017308)-0.15925.901(0.0159)-0.161
5)食事規則性-0.022112( 0.017923)-0.08001.522(0.2186)-0.083
6)料理への進取性0.002130( 0.015809)0.01010.018(0.8930)0.009
7)娯楽0.033040( 0.019972)0.12272.737(0.0995)0.110
8)教養-0.023875( 0.019721)-0.10191.466(0.2273)-0.081
9)社会奉仕-0.040800( 0.014739)-0.22967.663(0.0061)-0.183
10)義理人情0.026562( 0.018350)0.10962.095(0.1492)0.097
11)経済型0.001977( 0.019041)0.00680.011(0.9174)0.007
12)伝統型-0.021723( 0.019311)-0.07881.265(0.2619)-0.075
13)清潔-0.000565( 0.019977)-0.00180.001(0.9775)-0.002
14)妻主導型-0.009935( 0.014328)-0.04330.481(0.4888)-0.047
15)運動実施0.000598( 0.014554)0.00280.002(0.9673)0.003
16)疾病頻度0.042519( 0.014967)0.21818.070(0.0049)0.187
17)多愁訴-0.024898( 0.015643)-0.12682.533(0.1129)-0.106
18)情緒不安定0.007263( 0.013906)0.03750.273(0.6020)0.035
19)外向性0.000011( 0.024298)0.00000.000(0.9997)0.000
20)共感性0.013569( 0.018983)0.05430.511(0.4755)0.048
21)自発性0.011149( 0.017935)0.05330.386(0.5348)0.042
22)遺s伝的健康観0.031011( 0.013449)0.14995.317(0.0220)0.153
23)年齢0.012741( 0.003890)0.214410.730(0.0012)0.215
24)BMI0.024250( 0.011252)0.13774.645(0.0322)0.143
定数 4.7829( 0.000000)
重相関係数(二乗)0.533( 0.284) F値
(自由度)
3.672
( 24, 222)
自由度調整済重相関係数0.455( 0.207)p値0.000
自由度再調整済重相関係数0.361( 0.130)
AIC 351.7859

 表4は,HbAlcを基準変数に,他の24変数を説明変数として重回帰分析を行った結果である。表4の「偏回帰係数」とは重回帰式での,各説明変数にかける重みである。基準変数であるHbAlcの予測値は,偏回帰係数と説明変数の積の合計に,定数を加えた値によって求められる。

 表4の「標準誤差」とは,偏回帰係数の検定のときに使用するもので,偏回帰係数の推定上の誤差と考えればよい。

 ところで,基準変数に対して,どの説明変数の影響が大きいかは,偏回帰係数の大きさからは分からない。偏回帰係数は,説明変数にかける重みなので,大きな数値の変数では相対的に小さく,小さな数値の説明変数では相対的に大きくなるのが普通だからである。そのため,全ての変数を平均0,分散1に標準化した場合の重み,すなわち「標準偏回帰係数」が求められる。この標準偏回帰係数が大きい変数ほど,基準変数への影響があるものと考えてよい。表4では,標準偏回帰係数の絶対値がが 0.1 未満の説明変数は,影響が小さいと考えたので,灰色の帯をつけておいた。標準偏回帰係数でみると,社会奉仕が -0.230 で最も大きく,ついで疾病頻度の 0.218,年齢の 0.214 などの影響が大きいことが分かる。

 標準偏回帰係数で基準変数に対する影響の大きさをみることができるが,その大きさが偶然のものでないかを検定する必要がある。表4に示した上側の「F値」は,「母偏回帰係数=0」の帰無仮説を検定するためのものである。このF値が大きければ,その説明変数の基準変数への影響は大きいものと考えてよい。通常の検定のように,有意確率を「(p値)」として()内に示した。普通の検定では5%水準で,有意か有意でないかを判断することが多い。重回帰分析でもそのように考えてもよいが,後述するようにF値の数値が2以上か否かが規準として用いられることが多い。F値の大小と,前述の標準偏回帰係数の大小には,若干の相違があるものの,ほぼ対応関係がある。

 表4の最右欄の「偏相関係数」は,基準変数HbAlcと各説明変数について,他の説明変数の影響を除いた場合の相関係数を示している。したがって,基準変数と個々の説明変数単独の関係の強さは,この偏相関係数で表される。偏相関係数と標準偏回帰係数の数値はほぼ一致しているが,やはり微妙な相違がある。偏相関係数でみると,年齢がHbAlcと最も相関があり,社会奉仕は疾病頻度よりも相関が低くなっている。この関係は,F値の大小に一致する。

 表4の下側には,重相関係数 0.533 とその2乗の多重決定係数 0.284 が表示されている。その右側に,重相関係数の検定のためのF値と自由度が示されている。F値は 3.672 であり,そのp値は 0.000 と極めて高度に有意である。この結果は,求めた重回帰式が基準変数の予測に意味をもつことを示している。

 ところで,重相関係数は説明変数の増加に伴い単調に増加し,説明変数の個数が(ケース数−1)になると必ず1となる。従って,重相関係数はケース数が説明変数の個数に比べ大きくない場合,あまり適切な指標とはいえなくなる。そのような場合,ケース数と説明変数の個数を用いて,重相関係数を修正した「自由度調整済重相関係数」を求めればよい。実際にはさらに,自由度の修正を行った「自由度再調整済重相関係数」を求めることが多い。これらの2つの重相関係数は,説明変数の個数がある一定の値を超えると減少し,自由度調整済み重相関係数は負の値さえとりえる。もしも,その値が負となれば,求めた重回帰式には何の意味もない。

 最下段のAICは「赤池の情報量規準」と呼ばれるものであり,この値が小さいほどモデルのあてはまりがよいことを示すものである。これは,次に述べる説明変数の選択に関係している。


4.重回帰分析における説明変数の選択

すでに述べたように,重相関係数は説明変数の個数が増加するにつれ,その値も単調に増加し,説明変数の個数が標本数より1少ない数に達すると必ず,その値が1になる。したがって,説明変数の個数を増やして,予測の程度を見掛け上よいようにすることができる。そのようなことを避けるとともに,実際の分析で最適な変数のみを重回帰式に含めるための適当な基準が必要となる。ここでは,説明変数の選択方法と最適な説明変数の個数の決定に関する方法について簡単に述べる。

(1)総当たり法

説明変数が複数個ある場合,可能な変数の組み合わせについて重回帰式を求め,それらを比較し,最良の変数の組み合わせのものを選ぶ方法である。選択の基準には,自由度調整済み重相関係数もしくは自由度再調整済み重相関係数が最大となる変数の組み合わせを採用すればよい。しかし,単純な計算でわかるように,変数の組み合わせの個数は,変数がp個あれば,(2p−1)通りある。かりに説明変数が 10 個あるのならば 1023 通りある。このため,総当たり法は,説明変数の個数が少ない場合のみに用いられる。

(2)逐次選択法

変数選択の基本的な方法として,適当な基準により説明変数を1つずつ重回帰式に加えていく「変数増加法」と,逆にすべての説明変数を用いた分析から1つずつ変数を減らしていく「変数減少法」がある。

変数増加法では,一度重回帰式に含めた変数は,新たな変数の追加によって,回帰への寄与がほとんど無くなっても除外されることがない。そこで,一度重回帰式に含めた変数であっても,ある基準を満たさなくなった場合,重回帰式から除くようにする方法は,「変数増減法」と呼ばれる。

 同様に,変数減少法において,一度除外された変数であっても,その後他の変数との関係で回帰に寄与することがある場合,もう一度採用するという方法は,「変数減増法」と呼ばれる。

上記の4つの方法による変数選択結果は,変数選択の基準によっては,どれを用いても同様の場合もあるし,全く異なる場合もある。

 変数選択の基準として,前述の偏回帰係数の検定のF値を用いることが多い。しかし,重回帰分析での変数選択は,通常のF検定のように有意水準を5%や1%に設定するとは限らない。

 一般に,ある大きさのF値を設定し(Fin),その値と比較することで,変数を重回帰式に含めるか否かを決めるようにする。

 すなわち,ある説明変数のF値が,Finより大きければ,重回帰式に含めることにすればよい。

 ただし,重回帰式に含まれていない説明変数が複数個ある場合,最大のF値である変数を1つだけ選択するようにする。

 もし,すべての変数のF値がFinより小さければ,そこで変数の追加を終了する。 これが,F検定による変数増加法である。

 変数減少法も同様の基準を用いて行うことができる。

 全ての説明変数のF値を求め,その最小のものをFminとしよう。

 Fminの値があらかじめ設定しておいた基準値Foutと比べて小さければ,Fminをもつ変数を重回帰式から除くことにすればよい。

 求めたFminの値がFoutより大ならば,そこで変数の減少を終了する。

 変数増加法や減少法の基準であるF値をどのくらいにするかは,最終的に重回帰式に残される変数の個数に大きく影響する。

 通常は,FinやFoutの値として2を用いることが多い。

 変数増減法は,変数増加法で変数を選択していく場合,すでに重回帰式に含めた説明変数に対しても,その変数を除いた場合のF値を求め,そのうち最小の数値がFoutより小さければ,その変数を除くことにすればよい。

同様に,変数減増法は,変数減少法により変数を除いていく場合,すでに重回帰式から除かれた変数に対しても,F値を求め,その最大のものがFinより大きければ,その変数を再び重回帰式に含めるものとすればよい。

 これらの変数選択法のうち「変数減増法」が,最もよいモデルを与えると言われている。

表5.変数選択後の重回帰分析
変数偏回帰係数(標準誤差)標準偏
回帰係数
F値(p値)偏相関
係 数
4)糖分-0.043808( 0.015742)-0.16597.745(0.0058)-0.178
7)娯楽0.044704( 0.017981)0.16606.181(0.0136)0.159
8)教養-0.031136( 0.016972)-0.13293.366(0.0678)-0.118
9)社会奉仕-0.027891( 0.011813)-0.15695.574(0.0190)-0.152
16)疾病頻度0.052485( 0.013771)0.269214.525(0.0002)0.240
17)多愁訴-0.028694( 0.013229)-0.14614.705(0.0311)-0.140
22)遺伝的健康観0.033998( 0.012535)0.16437.356(0.0072)0.174
23)年齢0.011912( 0.003645)0.200410.683(0.0012)0.208
24)BMI0.023793( 0.010428)0.13515.206(0.0234)0.147
定数4.7043( 0.000000)
重相関係数(二乗)0.508( 0.258) F値
(自由度)
9.163
( 9, 237)
自由度調整済重相関係数0.480( 0.230)p値0.000
自由度再調整済重相関係数0.449( 0.202)
AIC 330.6123

 表5は変数減増法で,F値の基準であるFin,Foutともに2として,説明変数を選択したものである。説明変数として9個選択された。F値をみると,疾病頻度が 14.525 と最大であり,表4の結果とは異なっている。年齢のF値はあまり変化はない。糖分,娯楽,多愁訴,遺伝的健康観,BMIなどはいずれもF値は大きくなっている。逆に社会奉仕では,F値は小さくなっている。このように,重回帰分析の検定結果は,説明変数間の相関関係で大きく変わってくるので,その解釈も以外と簡単ではない。教養はF値が2以上なので,選択されているが,p値は0.068と5%以上になっている。

 重相関係数は減少したものの,自由度調整済み重相関係数,自由度再調整済み重相関係数のいずれも,大きくなっていることがわかる。また,変数選択のすべての過程を示していないが,赤池の情報量規準AICは最小の値となっていた。

 表5の結果から,疾病頻度,すなわち病気がちの傾向が高いとHbAlcが高いことになる。また,高年齢,肥満傾向の場合もHbAlcは高いことになる。糖分尺度で符合が負となっているのは,あまり間食等をするのはよくないとの認識があることによるものと考えられる。社会奉仕が負なのは,年齢,疾病頻度が高ければ,社会奉仕などを行えないためと考えられよう。娯楽傾向が何を意味するかは不明であるが,教養尺度が負なのは,健康情報の収集に不熱心な傾向とHbAlcが関係するのは,興味深い。また,健康は遺伝で決まるとする傾向とHbAlcが関係するのも,おもしろい結果であろう。多愁訴は身体的な訴えの多さを示すものであるが,これが負になっているのは,疾病頻度との関係によるものと考えられる。すなわち,重回帰分析では,説明変数間に極めて高い相関が認められる場合がよくある。そのような場合,一方の偏回帰係数が正であれば,他方の変数については負の値になることが多い。そのような場合には,相関の高い変数の一方を説明変数から除外したり,主成分分析をあらかじめ行い,主成分得点を用いて重回帰分析を行うなどの工夫が必要である。


5.おわりに

 実例を用いて,重回帰分析の基本的な方法について解説を行った。重回帰分析は多変量解析の中では,基本的な方法であり,その考え方を理解することは,他の手法の理解にも通じるものと考えられる。

 実際の分析には,前回と同様に統計学パッケージHALBAUを用いた6,7)。また,数式等は省略したが,興味ある読者は文献6などを参考にして欲しい。


参考文献

  1. 安藤洋美:ゴールトン(2),BASIC数学,27(12):pp.51-60,1994.
  2. 安藤洋美:ゴールトンとウェルドン,BASIC数学,28(1):pp.54-62,1995.
  3. 日野原重明,他:循環器疾患予防のための生活習慣に関する研究(第1報)生活習慣の多変量解析による分析,日本公衆衛生雑誌,29(7):309-320, 1982
  4. 佐伯圭一郎,他:LPC式生活習慣検査の作成,行動計量学,15(2),32-44,1988
  5. Takagi, H., et al.: Construction of the Life Habits Inventory, Methods and Applications in Mental Health Surveys: The Todai Health Index, (ed. by Suzuki, S. and Roberts, R. E.), University of Tokyo Press, Tokyo, 1991, 103-121 & 297-306.
  6. 柳井晴夫,高木廣文 編著:多変量解析ハンドブック,現代数学社,1989.
  7. 高木廣文,柳井晴夫 編著:HALBAUによる多変量解析の実践,現代数学社,1995.

トップページ
神戸市看護大学

最終更新日: 6 April 2023
無断転載の禁止 (C)高木廣文