精神疾患の軽症化

精神疾患の軽症化が進んでいる。統合失調症の新規発症が激減し、入院を要するレベルの気分障害がほとんどいなくなり、病棟経営の支えである慢性の統合失調症患者は急速に高齢化が進んでいる。地域移行を進めるまでもなく精神科病床数は今後10年で確実に激減する。

— 斎藤環 (@pentaxxx) November 14, 2017

精神疾患の軽症化というのは精神科医の実感としては常識レベルの話なのだけど、ついているコメントを見ると、懐疑的な反応が多いのが逆に意外。早期受診とか薬剤の副作用が少なくなったとか理由は複合的だと思うけど。 https://t.co/JldXcdghyp

— 風のハルキゲニア (@hkazano) November 15, 2017

精神科に新規で入院する患者数は増加傾向にあるようです。でもちょっと頭打ちになってきてそう。

社会・人口統計体系
都道府県データ / 社会生活統計指標より

一方で、平均在院日数（用語の説明）は減少傾向にあります。

社会・人口統計体系
都道府県データ / 社会生活統計指標より

精神病床数も減少傾向です。

社会・人口統計体系
都道府県データ / 社会生活統計指標より

入院期間ごとに近年の患者数をみても、減少傾向が伺えます。

患者調査 / 平成26年患者調査/ 閲覧（報告書非掲載表）
患者調査 / 平成23年患者調査/ 閲覧（報告書非掲載表）
患者調査 / 平成20年患者調査/ 閲覧（報告書非掲載表）より

こういうデータの出典ってこの書き方で大丈夫ですか？

以下にRのコードを載せます。

library(estatapi)library(tidyverse)
myappId <-"ここには自分のアプリケーションIDを入れる"

dat1 <-estat_getStatsData(appId = myappId, statsDataId ="0000010209")
unique(dat1$`Ｉ　健康・医療`)
nyuin <-dplyr::filter(dat1,`Ｉ　健康・医療`=="#I04103_精神科病院年間新入院患者数（人口10万人当たり）",地域=="全国")%>% 
  mutate(year=as.integer(substr(調査年,1,4)))

theme_set(theme_bw(14,"Osaka")+
            theme(axis.text = element_text(colour="black")))

p_nyuin <-ggplot(nyuin,aes(x=year,y=value))+
  geom_line()+
  labs(y="人口10万人当たり",x="年")+
  ggtitle("精神科病院年間新入院患者数（全国）")

ggsave("~/Desktop/p_nyuin.png",p_nyuin)

zaiin <-dplyr::filter(dat1,`Ｉ　健康・医療`=="#I10205_精神科病院平均在院日数",地域=="全国")%>% 
  mutate(year=as.integer(substr(調査年,1,4)))

p_zaiin <-ggplot(zaiin,aes(x=year,y=value))+
  geom_line()+
  labs(y="人口10万人当たり",x="年")+
  ggtitle("精神科病院平均在院日数（全国）")

ggsave("~/Desktop/p_zaiin.png",p_zaiin)

byosho <-dplyr::filter(dat1,`Ｉ　健康・医療`=="#I0910205_精神病床数（人口10万人当たり）",地域=="全国")%>% 
  mutate(year=as.integer(substr(調査年,1,4)))

p_byosho <-ggplot(byosho,aes(x=year,y=value))+
  geom_line()+
  labs(y="人口10万人当たり",x="年")+
  ggtitle("精神病床数（全国）")

ggsave("~/Desktop/p_byosho.png",p_byosho)

datId <-c("0003128737","0003071461","0003027589")
tim <- seq(2008,by=3,length.out =3)
dat2 <- lapply(datId,function(x){estat_getStatsData(appId = myappId,statsDataId = x)})

df <-sapply(1:3,function(i){dplyr::filter(dat2[[i]],grepl("５　精神及び行動の障害",傷病分類_001),入院期間_002!="総数",入院期間_002!="６月以上（再掲）",精神疾患の有無_001 =="精神疾患（副傷病）あり")%>% 
    mutate(year=tim[i])})%>% 
  bind_rows()
lev <-c("０～１４日","１５～３０日","１月～３月","３月～６月","６月～１年","１年～３年","３年～５年","５年以上","不詳")
df <-df %>% 
  mutate(term=factor(入院期間_002,levels = lev))

p1 <-ggplot(df,aes(x=term,y=value,fill=factor(year)))+
  geom_col(position ="dodge")+
  labs(fill="調査年",y="千人",x="入院期間")+
  facet_wrap(~傷病分類_001,scales ="free_y",nrow=4)
print(p1)
ggsave("~/Desktop/p1.png",p1,width =10)