library(openxlsx)

df<-read.xlsx("http://kanggc.iptime.org/book/data/lq-e.xlsx")

industry<-df$industry
id<-df$id[2:20]
nfirm<-df$nfirm
nemp<-df$nemp
jfirm<-df$jfirm
jemp<-df$jemp

lqfirm<-rep(NA, 19)
lqemp<-rep(NA, 19)


for(i in 1:19) {
  lqfirm[i]<-(jfirm[i+1]/jfirm[1])/(nfirm[i+1]/nfirm[1])
}

for(i in 1:19) {
  lqemp[i]<-(jemp[i+1]/jemp[1])/(nemp[i+1]/nemp[1])
}

lqfirm
lqemp

plot(lqfirm,lqemp,type="n",cex=1.5,xlim=c(0,8),ylim=c(0,8))
points(lqfirm,lqemp,pch=16,cex=1.5,col="blue")
with(df,text(lqemp~lqfirm,pos=1))
abline(v=1,h=1,col="blue",lty=2)
legend("center",legend=industry[2:20])