library(openxlsx) data<-read.xlsx("http://kanggc.iptime.org/book/data/rsi-jj-e.xlsx") id<-data$id name<-data$name k1986<-data$k1986 k1996<-data$k1996 k2006<-data$k2006 k2016<-data$k2016 jj1986<-data$jj1986 jj1996<-data$jj1996 jj2006<-data$jj2006 jj2016<-data$jj2016 B<-matrix(data=NA, nrow=16, ncol=4, byrow=T) for(i in 1:16) { for(j in 1:4) { B[i,1]<-abs((jj1986[i+1]/jj1986[1])-(k1986[i+1]/k1986[1])) B[i,2]<-abs((jj1996[i+1]/jj1996[1])-(k1996[i+1]/k1996[1])) B[i,3]<-abs((jj2006[i+1]/jj2006[1])-(k2006[i+1]/k2006[1])) B[i,4]<-abs((jj2016[i+1]/jj2016[1])-(k2016[i+1]/k2016[1])) } } BB<-round(B,digits=3) BB RSI<-colSums(B)/2 RSI A<-matrix(data=NA, nrow=16, ncol=8, byrow=T) for(i in 1:16) { for(j in 1:4) { A[i,1]<-(k1986[i+1]/k1986[1])*100 A[i,2]<-(k1996[i+1]/k1996[1])*100 A[i,3]<-(k2006[i+1]/k2006[1])*100 A[i,4]<-(k2016[i+1]/k2016[1])*100 } } for(i in 1:16) { for(j in 5:8) { A[i,5]<-(jj1986[i+1]/jj1986[1])*100 A[i,6]<-(jj1996[i+1]/jj1996[1])*100 A[i,7]<-(jj2006[i+1]/jj2006[1])*100 A[i,8]<-(jj2016[i+1]/jj2016[1])*100 } } AA<-round(A,digits=3) AA k1986<-A[,1] k1986_1<-sort(k1986, decreasing=T) k1986_2<-cumsum(k1986_1) k1986_3<-c(0,k1986_2) jj1986<-A[,5] jj1986_1<-sort(jj1986, decreasing=T) jj1986_2<-cumsum(jj1986_1) jj1986_3<-c(0,jj1986_2) k1996<-A[,2] k1996_1<-sort(k1996, decreasing=T) k1996_2<-cumsum(k1996_1) k1996_3<-c(0,k1996_2) jj1996<-A[,6] jj1996_1<-sort(jj1996, decreasing=T) jj1996_2<-cumsum(jj1996_1) jj1996_3<-c(0,jj1996_2) k2006<-A[,3] k2006_1<-sort(k2006, decreasing=T) k2006_2<-cumsum(k2006_1) k2006_3<-c(0,k2006_2) jj2006<-A[,7] jj2006_1<-sort(jj2006, decreasing=T) jj2006_2<-cumsum(jj2006_1) jj2006_3<-c(0,jj2006_2) k2016<-A[,4] k2016_1<-sort(k2016, decreasing=T) k2016_2<-cumsum(k2016_1) k2016_3<-c(0,k2016_2) jj2016<-A[,8] jj2016_1<-sort(jj2016, decreasing=T) jj2016_2<-cumsum(jj2016_1) jj2016_3<-c(0,jj2016_2) C<-cbind(k1986_1,jj1986_1,k1996_1,jj1996_1,k2006_1,jj2006_1,k2016_1,jj2016_1) CC<-round(C,digits=3) CC D<-cbind(k1986_3,jj1986_3,k1996_3,jj1996_3,k2006_3,jj2006_3,k2016_3,jj2016_3) DD<-round(D,digits=3) DD par(mfrow=c(1,2)) plot(id,k1986_3, type="l", xlab="", ylab="cumulative sum", col="black") lines(id,jj1986_3, col="blue") legend("bottomright", legend=c("Àü±¹","Á¦ÁÖ"), col=c("black","blue"), lty=1) plot(id,k2016_3, type="l", xlab="",ylab="cumulative sum", col="black") lines(id,jj2016_3, col="blue") legend("bottomright", legend=c("Àü±¹","Á¦ÁÖ"), col=c("black","blue"), lty=1)