日本アレルギー学会のアナフィラキシーガイドラインを読んだ。
最近プロットしてなかったので、死亡原因をプロットした。
dat <- read.delim("clipboard") anaph <- dat[, -1] yl <- c(0, 40) par(mar=c(4.5, 4.5, 2, 5), cex.lab=1.6) matplot(dat$year, anaph, ylim=yl, type="l", lty=1, lwd=3, xlab="年", ylab="アナフィラキシー原因と死亡者数") text(par()$usr[2], tail(anaph, 1), colnames(anaph), pos=4, xpd=TRUE, col=seq(ncol(anaph))) legend("topleft", legend=colnames(anaph), lty=1, col=seq(ncol(anaph)), text.col=seq(ncol(anaph)), lwd=3, bty="n")
year ハチ刺傷 食物 医薬品 血清 詳細不明 2001 26 3 17 0 12 2002 23 0 17 0 13 2003 24 3 19 1 6 2004 18 2 19 0 7 2005 26 1 31 1 14 2006 20 5 34 1 6 2007 19 5 29 1 12 2008 15 4 19 0 10 2009 13 4 26 1 7 2010 20 4 21 0 6 2011 16 5 32 0 18 2012 22 2 22 0 9 2013 24 2 37 1 13