廿TT

譬如水怙牛過窓櫺 頭角四蹄都過了 因甚麼尾巴過不得

精神疾患の軽症化

今日の川柳


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

f:id:abrahamcow:20180228220402p:plain

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

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

f:id:abrahamcow:20180228220651p:plain

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

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

f:id:abrahamcow:20180228220754p:plain

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

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

f:id:abrahamcow:20180228221255p:plain

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

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

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

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

dat1 <-estat_getStatsData(appId = myappId, statsDataId = "0000010209")
unique(dat1$`I 健康・医療`)
nyuin <-dplyr::filter(dat1,`I 健康・医療`=="#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,`I 健康・医療`=="#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,`I 健康・医療`=="#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("5 精神及び行動の障害",傷病分類_001),
                     入院期間_002!="総数",入院期間_002!="6月以上(再掲)",
                     精神疾患の有無_001 == "精神疾患(副傷病)あり") %>% 
    mutate(year=tim[i])}) %>% 
  bind_rows()
lev <-c("0~14日","15~30日","1月~3月","3月~6月","6月~1年","1年~3年","3年~5年","5年以上","不詳" )
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)