1.参考：田中孝文[2008]第二章に新たな景気の谷を追加
| x <- read.csv("gdp_ex.csv",skip=3) gdp <- x[,2] ex <- x[,3] ts_g <- ts(gdp,start=c(1980,1),frequency=4);ts_g ns1 <- (1993-1980)*4+4#1993q4 ne1 <- ns1+8 #2年後=8q先のデータ x1 <- ts_g[ns1:ne1]/ts_g[ns1]*100 ns2 <- (1999-1980)*4+1#1999q1 ne2 <- ns2+8 #2年後=8q先のデータ x2 <- ts_g[ns2:ne2]/ts_g[ns2]*100 ns3 <- (2002-1980)*4+1#2002q1 ne3 <- ns3+8 #2年後=8q先のデータ x3 <- ts_g[ns3:ne3]/ts_g[ns3]*100 ns4 <- (2009-1980)*4+1#2009q1 ne4 <- ns4+8 #2年後=8q先のデータ x4 <- ts_g[ns4:ne4]/ts_g[ns4]*100 #x1 plot,x2-3 lines z <- 0:8 yl <- "景気の谷=100" xl <- "景気の谷からの期間（単位：四半期）" ttl= "実質・原系列GDPの景気の谷からの回復パターン" plot(z,x1,type="o",ylim=c(y1,y2),ylab=yl,xlab=xl,main=ttl) lines(z,x2,type="o",col=2,lty=2) lines(z,x3,type="o",col=3,lty=3) lines(z,x4,type="o",col=4,lty=4) abline=(h=100) kikan <- c("1993q4-","1999q1-","2002q1-","2009q1-") legend(locator(1),legend=kikan,pch=1,lty=c(1:4),col=c(1:4)) |
2.季節調済で確認
| x <- read.csv("r_gpd_sa.csv",skip=1) gdp <- x[,2] ts_g <- ts(gdp,start=c(1980,1),frequency=4);ts_g ns1 <- (1993-1980)*4+4#1993q4 ne1 <- ns1+8 #2年後=8q先のデータ x1 <- ts_g[ns1:ne1]/ts_g[ns1]*100 ns2 <- (1999-1980)*4+1#1999q1 ne2 <- ns2+8 #2年後=8q先のデータ x2 <- ts_g[ns2:ne2]/ts_g[ns2]*100 ns3 <- (2002-1980)*4+1#2002q1 ne3 <- ns3+8 #2年後=8q先のデータ x3 <- ts_g[ns3:ne3]/ts_g[ns3]*100 ns4 <- (2009-1980)*4+1#2009q1 ne4 <- ns4+8 #2年後=8q先のデータ x4 <- ts_g[ns4:ne4]/ts_g[ns4]*100 #x1 plot,x2-3 lines y1 <- min(x1,x2,x3,x4[1:5]) y2 <- max(x1,x2,x3,x4[1:5]) z <- 0:8 yl <- "景気の谷=100" xl <- "景気の谷からの期間（単位：四半期）" ttl= "実質・季調済みGDPの景気の谷からの回復パターン" plot(z,x1,type="o",ylim=c(y1,y2),ylab=yl,xlab=xl,main=ttl) lines(z,x2,type="o",col=2,lty=2) lines(z,x3,type="o",col=3,lty=3) lines(z,x4,type="o",col=4,lty=4) abline=(h=100) kikan <- c("1993q4-","1999q1-","2002q1-","2009q1-") legend(locator(1),legend=kikan,pch=1,lty=c(1:4),col=c(1:4)) |
|
#3.純輸出 x<- read.csv("r_gpd_sa.csv",skip=1) |
4.輸出
| x <- read.csv("r_gpd_sa.csv",skip=1) ex <- x[,4] ts_ex <- ts(ex,start=c(1980,1),frequency=4);ts_ex ns1 <- (1993-1980)*4+4#1993q4 ne1 <- ns1+8 #2年後=8q先のデータ x1 <- ts_ex[ns1:ne1]/ts_ex[ns1]*100 ns2 <- (1999-1980)*4+1#1999q1 ne2 <- ns2+8 #2年後=8q先のデータ x2 <- ts_ex[ns2:ne2]/ts_ex[ns2]*100 ns3 <- (2002-1980)*4+1#2002q1 ne3 <- ns3+8 #2年後=8q先のデータ x3 <- ts_ex[ns3:ne3]/ts_ex[ns3]*100 ns4 <- (2009-1980)*4+1#2009q1 ne4 <- ns4+8 #2年後=8q先のデータ x4 <- ts_ex[ns4:ne4]/ts_ex[ns4]*100 #x1 plot,x2-3 lines y1 <- min(x1,x2,x3,x4[1:5]) y2 <- max(x1,x2,x3,x4[1:5]) z <- 0:8 yl <- "景気の谷=100" xl <- "景気の谷からの期間（単位：四半期）" ttl= "実質・季調済み輸出の景気の谷からの回復パターン" plot(z,x1,type="o",ylim=c(y1,y2),ylab=yl,xlab=xl,main=ttl) lines(z,x2,type="o",col=2,lty=2) lines(z,x3,type="o",col=3,lty=3) lines(z,x4,type="o",col=4,lty=4) abline=(h=100) kikan <- c("1993q4-","1999q1-","2002q1-","2009q1-") legend(locator(1),legend=kikan,pch=1,lty=c(1:4),col=c(1:4)) |
5.増減率
| x<- read.csv("r_gpd_sa.csv",skip=1) gdp <- x[,2] ex <- x[,4] dd <- x[,5] FD <- ts(ex,start=c(1980,1),frequency=4) DD <- ts(dd,start=c(1980,1),frequency=4) GDP <- ts(gdp,start=c(1980,1),frequency=4) WD <- DD/GDP; WF <- FD/GDP rGDP <- diff(log(GDP))*100 rDD <- diff(log(DD))*100 rFD <- diff(log(FD))*100 cDD <- rDD*WD; cFD <- rFD*WF year <- 1981:2004 kekka <- data.frame(year,rGDP,cDD,cFD);kekka ts.plot(rGDP,cDD,cFD,type="l",lty=c(1:3),col=c(1:3)) abline(h=0) youin <- c("季調済GDP前年比（対数）","国内需要前年比（対数）","輸出前年比（対数）") legend(locator(1),legend=youin,pch=1,lty=c(1:3),col=c(1:3)) |
| x<- read.csv("r_gpd_sa.csv",skip=1) gdp <- x[,2] ex <- x[,4] dd <- x[,5] FD <- ts(ex,start=c(1980,1),frequency=4) DD <- ts(dd,start=c(1980,1),frequency=4) GDP <- ts(gdp,start=c(1980,1),frequency=4) WD <- DD/GDP; WF <- FD/GDP X <- rbind(DD,FD) X1 <- X*(X>=0);X2 <- X*(X<0) rownames(X1) <- c("DD","FD") ttt <- "FD and DD" yl <- c(-4,4) mp <- barplot(X1,main=ttl,ylim=yl,ylab="%") barplot(X2,add=TRUE) abline(h=0) lines(mp,GDP,type="o") hanrei <- c("DD","FD");clr <- c(gray(0.2),1) legend(mp[7],-2,legend=hanrei,pch=c(15,0),col=clr) |
