太った日本人は増えているのか?

少々、荒い解析ですが、「太った日本人は増えているのか」 という疑問があり、R でグラフを描き、軸ラベルは後から入れ替えました。

男女、年齢階級別に、年齢内の%を下記 3 群で、グラフに示します。

  • BMI<20: やせている人
  • 25<= BMI <30: 太っている人 (日本では「肥満」とされる)
  • 30<= BMI : より太った肥満の人 (世界的に「肥満」)

データは厚生労働省の、主な健康指標の経年変化 から入手して、およそ5年毎に抜き出しました。

男性は、太っている人が増えている考えるのが良さそうですね。

折れ線グラフを重ねていますが、R では、下記のようにして描きました。for ループの中に、"Par(new)=T" を入れると描けるようです。 ただし、軸ラベルは9回重ね描きになりますので、PDFで保存して、後から、イラストレーターで描きなおしています。 イラストレーターで加工時は、クリッピングマスクを解除すると、加工しやすくなります。

bmi2014<-read.csv("~xxxxxx.csv",header=T)
mbmi2014<-bmi2014[bmi2014$sex=="Males",]
x<-mbmi2014$less20
for (i in 1:9) {
par(new=T)
plot(mbmi2014$year[mbmi2014$group==i],x[mbmi2014$group==i],type="o", ylim=c(0,max(x)*1.2), col=i,pch=i, lty=i)
}
labels <- c("15-19", "20-29", "30-39", "40-49", "50-59", "60-69", "70<=", "20<=", "Total")
legend("topleft", legend = labels, col = c(1:9), pch = c(1:9), lty = c(1:9))

読み込みデータの一部は、下記のような感じです。 

year	sex	group	name	less20	over25	over30
1975	Males	9	Total	25.9	14.3	1
1975	Males	1	15〜19	51.3	5	0.8
1975	Males	2	20〜29	30.5	8.2	0.7
1975	Males	3	30〜39	17.8	20	0.7
1975	Males	4	40〜49	16.5	19.6	1.2
1975	Males	5	50〜59	21.8	14	1.3
1975	Males	6	60〜69	29.5	12.8	1.7
1975	Males	7	70 and over	34.1	12.2	0.9
1975	Males	8	20 and over	23.1	15.3	1
1975	Females	9	Total	25.7	16.2	2.2
1975	Females	1	15〜19	38	3.7	0.2
1975	Females	2	20〜29	40.7	6	0.6
1975	Females	3	30〜39	23.9	11.7	1.8
1975	Females	4	40〜49	16.2	21	3.1
1975	Females	5	50〜59	18.8	26.1	3
1975	Females	6	60〜69	23.9	24.2	4.3
1975	Females	7	70 and over	28.7	20.3	2.2
1975	Females	8	20 and over	24.8	17.1	2.4
1980	Males	9	Total	23.4	15.9	0.8
1980	Males	1	15〜19	42.7	5.6	0.7
以下略
B! LINE