crosshyou

主にクロス表(分割表)分析をしようかなと思いはじめましたが、あまりクロス表の分析はできず。R言語の練習ブログになっています。

宗教統計調査データの分析2 - R言語で各宗派ごとの団体数の推移を調べたら、キリスト教系だけが増えていた。

 

www.crosshyou.info

 の続きです。

まず、合計だけのデータセットを作成して、全体の宗教法人の数がどのようなトレンドで推移しているかを見てみましょう。R言語のfilter関数を使います。

f:id:cross_hyou:20190307190335j:plain

ggplot関数でグラフを描いてみましょう。

f:id:cross_hyou:20190307190703j:plain

f:id:cross_hyou:20190307190736j:plain


あ、宗教法人は全体としては低下傾向ですね。でも、2005年が18万2641で、2017年が18万1497と1141(率にして0.6%)ほどの低下です。

宗教のタイプ別ではどうでしょうか?合計以外のデータセットを作成して調べてみましょう。

f:id:cross_hyou:20190307192131j:plain

たしか所轄が都道府県知事と文部科学大臣があったんですよね。この二つを合計したいですね。。どうすればいいのかな?とりあえず、都道府県知事だけのデータセットと文部科学大臣だけのデータセットを作成しましょう。

f:id:cross_hyou:20190307193206j:plain

こうしたら、それぞれの所轄を削除しましょう。

f:id:cross_hyou:20190307193740j:plain

両方のNumberを都道府県知事のほうはNumber_chiji, 文部科学大臣のほうはNumber_daijinに変更しましょう。

f:id:cross_hyou:20190307194501j:plain

こうしたら、left_join関数でいけると思います。やってみましょう。

f:id:cross_hyou:20190307195353j:plain

chijiにdaijinを合体させますが、YearとTypeが同じものどうしを合体させるように、ということですね。

途中で気が付いたのですが、Typeに「都道府県知事」が混じっていたので、除外します。

f:id:cross_hyou:20190307200018j:plain

with 42 more rows となっています。さっきは with 55 more rowsでしたので削除されました。

Number_chijiとNumber_daijinを合計します。mutate関数ですね。

f:id:cross_hyou:20190307200334j:plain

これでようやく宗派別の推移をグラフにできます。ggplot関数ですね。

まずはキリスト教系から

f:id:cross_hyou:20190307201345j:plain

 

f:id:cross_hyou:20190307201358j:plain

 

キリスト教系は増加しています。

仏教系はどうでしょうか?

f:id:cross_hyou:20190307201607j:plain

f:id:cross_hyou:20190307201618j:plain

仏教系は減少です。

 

神道系はどうでしょうか?

f:id:cross_hyou:20190307201820j:plain

f:id:cross_hyou:20190307201831j:plain

神道系も減っています。

諸教はどうでしょうか?

f:id:cross_hyou:20190307202032j:plain

f:id:cross_hyou:20190307202045j:plain

諸教も減っていますね。

キリスト教系の団体だけが増えていることがわかりました。

今回は以上です。

次回は

 

www.crosshyou.info

 です。

 

最後に本日のR言語のコマンドを記しておきます。


# 2019-03-07

# 合計だけのデータセットを作成
total <- df %>%
filter(Shokatsu == "合計")
total

# YearとNumberのグラフ
ggplot(data = total, mapping = aes(x = Year, y = Number)) +
geom_line()

# 合計以外のデータセットを作成
ntotal <- df[df$Type != "合計", ]
ntotal

# 都道府県知事のデータセットを作成
chiji <- ntotal[ntotal$Shokatsu == "都道府県知事", ]
chiji

# 文部科学大臣のデータセットを作成
daijin <- ntotal[ntotal$Shokatsu == "文部科学大臣", ]
daijin

# Shokatssu列を削除
chiji <- chiji[ , -2]
chiji

daijin <- daijin[ , -2]
daijin

# 変数名の変更
colnames(chiji)[3] <- "Number_chiji"
chiji

colnames(daijin)[3] <- "Number_daijin"
daijin

# chijiとdaijinを合体
chida <- chiji %>%
left_join(daijin, by = c("Year", "Type"))
chida

# Typeに「都道府県知事合計」が混じっていたので除外
chida <- chida[chida$Type != "都道府県知事合計", ]
chida

# Number_chijiとNumber_daijinを合計
chida <- chida %>%
mutate(Number_total = Number_chiji + Number_daijin)
chida

# 宗派別グラフ:キリスト教系
chida %>%
filter(Type == "キリスト教系") %>%
ggplot(mapping = aes(x = Year, y = Number_total)) +
geom_line() +
labs(title = "キリスト教系")

 

# 宗派別グラフ:仏教系
chida %>%
filter(Type == "仏教系") %>%
ggplot(mapping = aes(x = Year, y = Number_total)) +
geom_line() +
labs(title = "仏教系")


# 宗派別グラフ:神道系
chida %>%
filter(Type == "神道系") %>%
ggplot(mapping = aes(x = Year, y = Number_total)) +
geom_line() +
labs(title = "神道系")


# 宗派別グラフ:諸教
chida %>%
filter(Type == "諸教") %>%
ggplot(mapping = aes(x = Year, y = Number_total)) +
geom_line() +
labs(title = "諸教")