みなさん，こんにちは．
シンノユウキ（y_stadio）です．
今回はExcelを用いてエネルギー調整を行う方法を紹介します．栄養素摂取量を解釈する場合，エネルギー調整はほぼ必須になります．以下でその方法を紹介していきます．
では行きましょう！
なぜエネルギー調整をしなければならないのか
はじめに，なぜエネルギー調整を行う必要があるのかを簡単に解説します．
その答えの１つはエネルギーの摂取量が他の多くの要素に影響を与えているためです．皆さんは，エネルギー摂取量の多い人と言われると，一般的にどのような人を想像されるでしょうか．やはり男性で体格が大きく，また普段の活動量が多い人を想像されるのではないでしょうか．また，年齢的に若い，成長期に該当する年齢の人の方を思い浮かべる方もいらっしゃるかもしれません．
また，エネルギー摂取量が多いと，それに伴って３大栄養素やビタミン，ミネラルの摂取量も増える傾向にあります．３大栄養素は当然ですよね．それらの栄養素からエネルギーが生み出されるので，３大栄養素が比例して多くなるのは当然です．また，エネルギー摂取量が多い人は，食事量が多い傾向にありますので，ビタミンやミネラルも必然的に多くなる場合があります．
すなわち，エネルギー摂取量は，性別や年齢，体格，運動量，多くの栄養素摂取量に直接的，もしくは間接的に影響を与えているという事になります．一般的に，食品や栄養素と疾病との関連を検討する栄養疫学研究においては，エネルギー摂取量が影響を与えるこれらの要素を，無視して考えることはできません．
たとえば，研究を行った場合に，飽和脂肪酸の摂取量が増えると，脳血管障害が増えるという結果が得られたとしましょう．そうした場合に，エネルギー摂取量を調整していないと，エネルギー摂取量が多いから脳血管障害が増えたという可能性や，体格の大きい人に脳血管障害が多いのかもしれないという可能性を捨て去ることができなくなります．そのため，栄養疫学研究においては，エネルギー調整は必須だと言えるでしょう．
密度法
密度法とは？
密度法というのは，栄養素の粗摂取量をエネルギー摂取量で割ったものです．たとえば，エネルギーの摂取量が1500kcal，ビタミンAの摂取量が300μg RAEだとすると，以下のような式になります．
エネルギー調整値(μg RAE/kcal) = 300(μg RAE) ÷ 1500(kcal)
密度法におけるエネルギー調整は，主に2つのパターンに分かれます．それは，
- 1000kcalあたりの栄養素摂取量で表す
- PFCバランスとして表す
のパターンです．
1000kcalあたりの栄養素摂取量で表すパターンでは，
エネルギー調整値(μg RAE/kcal) = 300(μg RAE) ÷ 1500(kcal) × 1000
となります．
PFCバランスとして示す場合はもっと単純ですね．たんぱく質を例として紹介しますと，たんぱく質60g，エネルギー1500kcalとすると，以下のような式になります．
P比(%) = 60 × 4 ÷ 1500 × 100
Excelで行う方法
では，これをExcelで行う場合を見ていきましょう．
まずは1000kcalあたりのエネルギー調整値の算出法から．
実際の計算式などは以下のようになります．
また，PFCバランスの場合は以下のようになります．
これは比較的理解しやすいのではないかと思います．
では，次は残差法について紹介していきます．
残差法
残差法とは？
残差法は，総エネルギー摂取量の影響を完全に取り除くことのできる方法です．実は，密度法では，総エネルギー摂取量の影響を完全に取り除くことができません．少し影響は弱くはなりますが，すべてを取り除くことはできないのですね．しかし，残差法を用いると，完全に取り除くことができるのです．
残差法による計算式は以下のようになっています．
エネルギー調整値 = 集団の平均エネルギー摂取量 ＋ 残差
ここでいう残差というのは，その対象者の，実際の栄養素摂取量と，回帰式から予測される予測値との差のことを指します．残差はプラスの値もマイナスの値も取りうることになります．以下に残差を算出するための計算式を示します．
残差 = 栄養素摂取量の実測値 - 栄養素摂取量の予測値
予測値は，先ほども示しましたように，回帰式から算出します．栄養素摂取量を従属変数，エネルギー摂取量を独立変数とし，回帰分析を行います．これで予測値を算出することができます．
Excelで行う場合
では，Excelで残差法によるエネルギー調整を行う方法を示したいと思います．まずは，以下からサンプルファイルをダウンロードしてください．以降ではこれを使用して説明を行います．
では以下に手順を示します：
１：回帰式の「傾き」を求める
まずは残差を求めていく必要があります．たんぱく質摂取量を従属変数（y=ax+bにおけるy）に，エネルギー摂取量を独立変数（y=ax+bにおけるx）とした回帰式における，傾きaと切片bを求めます．
ここでは，まず回帰式の「傾き」を求めましょう．
これらはExcelの関数で求めることができ，それぞれ，
・切片：=INTERCEPT()
が用意されています．まずは傾きから求めていきましょう．セルH1に以下の計算式を入力してください．
ここで，最初の引数がたんぱく質摂取量を，2番めの引数がエネルギー摂取量をそれぞれ示しています．
２：残差を求める（切片）
次に，回帰式の切片を求めましょう．切片はINTERCEPT関数で求めることができます．
以下の計算式をセルH2に入力してください：
引数は先ほどと同じですね．
３：栄養素摂取量の「予測値」を求める
次に，それぞれの対象者における，たんぱく質摂取量の予測値を計算していきます．
y=ax+bの一次回帰式に，先ほど算出した傾きと切片を当てはめ，対象者のエネルギー摂取量から，その集団におけるそのくらいのエネルギーを摂取している場合に予測される，たんぱく質摂取量を算出します．
以下の計算式を，セルC2に入力してください．
H1が傾き，A2がエネルギー摂取量，H2が切片です．y=ax+bに当てはめただけですね．後ほど下にオートフィルでコピーするので，ドルマークは忘れないようにしてくださいね．
４：オートフィルでコピー
では，セルC2を，下方向にオートフィルを使用してコピーしてください．これで，各対象者の，たんぱく質摂取量の予測値を算出することができましたね．
５：実測値と予測値から残差を求める
次は，各対象者毎に残差を求めます．以下の計算式を，セルD2に入力してください．
簡単ですよね．実測値-予測値で残差を算出することができます．これを，予測値を算出した時の同じ様に，下方向にコピーしてください．これで残差が求まります．
６：わかりやすくするために平均値を足す
残差が求まれば，それが残差法によるエネルギー調整値となりますが，残差は，プラスにもマイナスの値にもなりうるので，実際にはわかりにくくなります．そのため，集団の栄養素摂取量を加えた値をエネルギー調整栄養素摂取量として示すのが一般的です．以下の計算式を，セルE2に入力し，下方向にオートフィルでコピペしてください．
AVERAGE($B$2:$B$20)は栄養素摂取量の平均値ですね．これに個々人の残差を足しています．
これで，エネルギー調整栄養素摂取量を算出することができました！
なお，完成版のExcelファイルは以下からダウンロードできます．確認の際にご利用ください．
まとめ
今回は，栄養疫学で必要となるエネルギー調整について，密度法による方法と，残差法による方法の2つをお示ししました．2つの方法は，その使い方や解釈に違いがあります．それらをしっかりと理解した上で使用するようにしてください．
＊残差法を簡単に行えるExcelシートを公開しました．くわしくは以下よりご覧ください．
参考文献
- 田中平三, 横山徹爾. "栄養疫学における総エネルギー摂取量に対する解釈と取り扱い方." 日本栄養・食糧学会誌 50.4 (1997): 316-320.
- 伊達ちぐさ, 徳留裕子, 吉池信男. "食事調査マニュアルはじめの一歩から実践・応用まで, 改定." (2009): 4-12.
