糖質制限ダイエットを始めると
たった1日で2kgも痩せた件

@Med_KU

20170701
StanとRでベイズ統計モデリング読書会
(Osaka.Stan#5)
関西学院大学西宮上ケ原キャンパス

今日の話

実家に帰ったとき

糖質制限ダイエットを始めたら1日で2kgも痩せたの!!

計画的にデータを取ってから言えよ…

本当に??

本当に糖質制限ダイエットで
2kg/day 痩せるのかもしれない(?)

糖質制限ダイエット

low-carbohydrate diet (LC)

厳格には，一日20g 以下の糖質摂取で，総エネルギー量は減らさなくてよい(1500kcal/day)．

15kg/年 (BMI 33.4) (Nutrition. 2016 Sep;32(9):1033-6.)
9.5kg/年 (BMI 34.6) (J Intern Med. 2016 Oct;280(4):388-97.)
8.6-9.6kg/半年(Obesity (Silver Spring). 2016 Jan;24(1):79-86.)
5.3kg/年 (Nutrients. 2015 Sep 17;7(9):7978-94.)
総説では，けっこう痩せる(Br J Sports Med. 2017 Jan;51(2):133-139.)
メタアナリシスでは，痩せる(Lancet Diabetes Endocrinol. 2015 Dec;3(12):968-79.)

長期的には 糖質制限ダイエットで痩せそう．

セミナーに行ってきた話

日本行動計量学会第19回春の合宿セミナー

Rstan とベイズモデリングの実習

Rstan Advent Calender 2016

時系列データを用いて1日あたりの体重変動を推定しよう
(※このデータは糖質制限ダイエットではありません !)

状態空間モデル

時系列データ\(w_t\) とその観測\(W_t\)

欠測値があっても推定可能

Rstan でのモデル

parameters{
  real<lower=0, upper=200> w[Day]; # 真の体重
  real<lower=0> b;                 # 隠れ層の分散
  real<lower=0> s;                 # 観測層の分散
}
model{
  for(i in 3:Day){
    w[i] ~ normal(2*w[i-1] - w[i-2], b); # 二階差分
  }
  for(i in 1:Na){                        # 欠測値を回避
    W[ na[i] ] ~ normal(w[ na[i] ], s);  # 観測値が正規分布から得られるモデル    #1
    W[ na[i] ] ~ cauchy(w[ na[i] ], s);  # 観測値がcauchy 分布から得られるモデル #2
  }
  b ~ cauchy(0, 0.5);                    # 適当にコーシー分布
  s ~ cauchy(0, 0.5);                    # 適当にコーシー分布
}

分布の選択

何も考えなければ，正規分布

裾を広くしたければ，t分布もしくはcauchy分布

結果
糖質制限ダイエットは
たった1日で2kg 痩せるのか?

長期に頑張れば，普通は痩せる

7ヶ月で8kg

推定された1日あたりの減量効果量は

norm モデルで 0.0366 kg/day
cauchy モデルで 0.0385 kg/day

2kg/day 痩せる減量効果量は ?

1日で2kg 痩せる95% 確率点を実現する減量効果量は

norm モデルで1.2729 kg/day
cauchy モデルで0.7652 kg/day

1ヶ月も頑張れば
10 kg 以上痩せる !?!?

重篤な疾患が隠れている可能性があるので

病院に行きましょう

糖質制限ダイエットを始めるとたった1日で2kgも痩せた件