Rで何かをしたり、読書をするブログ

政府統計の総合窓口のデータや、OECDやUCIやのデータを使って、Rの練習をしています。ときどき、読書記録も載せています。

鉱工業出荷内訳表の分析2 - データの推移をグラフにして観る。(plot関数)

今回は、鉱工業出荷内訳表のデータをグラフにしようと思います。とはいっても、データの種類が48もあるので、前回の分析で判明した、平均値が一番大きかった、はん用.国内、平均値が一番小さかった、情報通.国内、標準偏差が一番大きかった、情報通.輸出、標準偏差が一番小さかった、化学工.出荷の4つのデータ系列についてグラフを作成しようと思います。

まずは、read.csv関数でデータを読み込み、summary関数でデータの要約を表示します。

f:id:cross_hyou:20180808163001j:plain

f:id:cross_hyou:20180808163024j:plain

それではグラフを作成します。plot関数を使います。type = "l" とオプション引数を追加してライングラフにします。

まずは、平均値が一番大きいデータ、はん用.国内です。

f:id:cross_hyou:20180808164249j:plain

f:id:cross_hyou:20180808164301j:plain

ylim=c(0,150)というオプション引数でY軸の範囲を0から150に設定しています。はじめのほうに急落して、だんだんと回復という動きです。

次は、平均値が一番小さいデータ、情報通.国内です。

f:id:cross_hyou:20180808170025j:plain

f:id:cross_hyou:20180808170040j:plain

はじめのころに大きく落ち込み、大きく回復した後、減少傾向が続いています。

 

次は、標準偏差が一番大きいデータ、情報通.輸出です。

f:id:cross_hyou:20180808170338j:plain

f:id:cross_hyou:20180808170350j:plain

はじめのころに大きく落ち込んだ後、回復しないで減少傾向です。

 

次は、標準偏差が一番小さいデータ、化学工.出荷です。

f:id:cross_hyou:20180808170719j:plain

標準偏差が一番小さいだけあって、動きが今までの3つとは違ってかなり小さいですね。

それでは最後に、この4つのデータを一つの画面にグラフ化します。

f:id:cross_hyou:20180808172447j:plain

まず、はん用.国内のグラフを表示します。このとき、main=""でグラフのタイトルを空白に、col="black"でグラフの色を黒に、ylab=""でY軸のラベルを空白にしています。

そして、そのあとに、par(new=TRUE)でグラフを重ねるように指示しています。

f:id:cross_hyou:20180808172716j:plain

このように、データによってはグラフの形状がかなり違うことがわかります。