# R script used in https://remcat.hatenadiary.jp/entry/20210911/gap
# URL: http://tsigeto.info/maikin/maikin-monthly.r.txt
# 2021/08/16 - 2021/09/11
# Created by TANAKA Sigeto
x <- read.delim( "maikin-monthly.dat", header=T )
# Sort by date and establishment size
x <- x[ order(x$size) , ]
x <- x[ order(x$yyyymm) , ]
x$is <- factor( paste( x$industry, x$size, sep="." ) )
x$year <- round( x$yyyymm / 100 )
x$month<- floor( x$yyyymm %% 100 )
x$worker2 <- ( x$e0 + x$e1 ) /2
x1 <- x
x2 <- x
x1$worker <- x1$e0
x2$worker <- x2$e1
x2$yyyymm <- x2$yyyymm + 0.5
temp <- rbind( x1, x2 )
x.long <- temp[ order(temp$yyyymm) , ]
x.bysize <- subset( x.long, industry=="TL" & 0