library(openxlsx)

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

k<-sample1$Korea
k2010<-sample1$K2010
k2017<-sample1$K2017
jj<-sample1$Jeju
jj2010<-sample1$J2010
jj2017<-sample1$J2017

A<-matrix(data=NA, nrow=17, ncol=9, byrow=T)
for(i in 1:17) {
 for(j in 1:7) {
    A[i,1]<-((jj2017[i]-jj2010[i])/jj2010[i])*100
    A[i,2]<-((k2017[i]-k2010[i])/k2010[i])*100
    A[i,3]<-((jj2010[i]/jj2010[1]))*100
    A[i,4]<-((k2010[i]/k2010[1]))*100
    A[i,5]<-((A[i,1]*A[i,3]))/100 
    A[i,6]<-((A[i,2]*A[i,4]))/100 
    A[i,7]<-(A[i,5]-A[i,6])
     }
   }
for(i in 1:16) {
  for(j in 8:9) {
    
    A[i+1,8]<-(A[i+1,3]-A[i+1,4])*A[i+1,2]/100
    A[i+1,9]<-(A[i+1,1]-A[i+1,2])*A[i+1,3]/100
    }
  }
CS<-colSums(A[,8:9], na.rm=T)
CS
A[1,8]<-CS[1]
A[1,9]<-CS[2]
AA<-round(A,digits=3)
AA