# source("http://minato.sip21c.org/ldaR/PEI.R") library(fmsb) NMbyPREF2015M <- c(-3993, -3008, -1636, 386, -1931, -1691, 825, -3648, -1631, 895, 9468, 2514, 37966, 6610, -2875, -163, 265, -1000, -1351, -1368, -3053, -2672, 7105, -2092, -1240, -1265, -1684, -4727, -2406, -1879, -769, -512, -1184, -1429, -2034, -947, -47, -1645, -938, -191, -1339, -2783, -2031, -943, -1549, -2152, -228) NMbyPREF2015F <- c(-4423, -3585, -2657, -175, -2543, -2338, -2892, -4279, -2091, -1410, 8609, 5525, 46265, 10666, -3612, -874, -635, -1192, -1435, -1876, -3520, -3717, 3413, -2484, -861, 627, 2590, -2639, -1550, -1938, -762, -892, -1204, -3005, -2257, -1239, -523, -2178, -1400, 1204, -1383, -3483, -2087, -1557, -1782, -2557, 136) URL <- "https://www.e-stat.go.jp/stat-search/file-download?statInfId=000031473212&fileKind=1" x <- read.csv(URL, skip=12, header=FALSE, as.is=TRUE) y <- subset(x, V3%%1000==0&V3!=0) yM <- subset(y, V2==201) zM <- sapply(yM[, 7:117], as.integer) TotalM <- rowSums(zM, na.rm=TRUE) NMRbyPREF2015M <- NMbyPREF2015M / TotalM PEIM <- apply(zM, 1, PEI) yF <- subset(y, V2==301) zF <- sapply(yF[, 7:117], as.numeric) TotalF <- rowSums(zF, na.rm=TRUE) NMRbyPREF2015F <- NMbyPREF2015F / TotalF PEIF <- apply(zF, 1, PEI) layout(t(1:2)) plot(PEIM, NMRbyPREF2015M, type="p", pch=16, xlab="ふくらみ指数(改)", ylab="純移動率", main="都道府県別ふくらみ指数(改)と純移動率の関係\n(2015年日本,男性)") plot(PEIF, NMRbyPREF2015F, type="p", pch=16, xlab="ふくらみ指数(改)", ylab="純移動率", main="都道府県別ふくらみ指数(改)と純移動率の関係\n(2015年日本,女性)")