R

21世紀の資本グラフ3

Quandl、xts、lattice、gridExtra、plotrix、RColorBrewer パッケージ

「21世紀の資本」のデータが公開されてるのでRを使ってグラフ化してみます。
（グラフもすでに公開されているのであまり意味はありません。）
（参考）
Piketty Codes
『21世紀の資本』日本語版サポートページ

積み上げ折れ線グラフ・・純外国資本の「マイナス」が表現できていません。「21世紀の資本グラフ4」では色の透過率を使って表現してみました。

Chapter 3: The Metamorphoses of Capital

Capital in Britain, 1700-2010
1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
29
30
31
32
library(xts)
library(Quandl)
T3_1<-Quandl("PIKETTY/T3_1")
#並べ替え
sortlist <- order(T3_1[,1])
dat <- T3_1[sortlist,]
T3_1<-dat
rownames(T3_1) <- c(1:nrow(T3_1))
#save("T3_1", file="T3_1.dat")
#load("T3_1.dat")
##### T3_1 #####
#1 Year 年
#2 National capital Wn 国民資本
#3 incl. Land 土地
#4 incl. Housing 住宅
#5 incl. Other domestic capital assets その他の国内資本資産
#6 incl. Net foreign capital 純外国資本
#7 Public capital Wg 公的資本
#8 incl. Public assets 公的資産
#9 incl. Public debt 公的債務
#10 Private capital W 民間資本
library("plotrix")
library(RColorBrewer)
#png("Piketty3_1A.png",width=1000,height=800)
stackpoly(T3_1[,3:6],ylim=c(0,800),axis4=F,main="イギリスの資本 1700-2010年",col=brewer.pal(4,"Pastel1"),
xaxlab=c("1700","1750","1810","1850","1880","1910","1920","1950","1970","1990","2010"),border="gray36",staxx=F,stack=TRUE)
#boxed.labels(3,180,"土地",border=F,bg=NA,col="gray20")
#boxed.labels(4,300,"住宅",border=F,bg=NA,col="gray20")
#boxed.labels(4,500,"その他国内資本",border=F,bg=NA,col="gray20")
#boxed.labels(5.4,620,"純外国資本",border=F,bg=NA,col="gray20")
boxed.labels(c(3,4,4,5.4),c(180,300,500,620),c("土地","住宅","その他国内資本","純外国資本"),border=F,bg=NA,col="gray20",cex=2)
#dev.off()

1
2
library(knitr)
kable(T3_1[,c(1,8,9)]);kable(T3_1[,c(1,2,7,10)])
Year incl. Public assets incl. Public debt
1700-12-31 30.00000 23.00000
1750-12-31 50.00000 107.00000
1810-12-31 70.00000 187.00000
1850-12-31 52.82894 108.50827
1880-12-31 34.25263 57.62411
1910-12-31 46.19908 26.91437
1920-12-31 93.26282 142.75725
1950-12-31 110.00000 229.89309
1970-12-31 100.48936 76.36428
1990-12-31 90.31268 59.29002
2010-12-31 92.46243 91.78362
Year National capital Wn Public capital Wg Private capital W
1700-12-31 703.4483 7.00000 696.4483
1750-12-31 681.1568 -57.00000 738.1568
1810-12-31 687.3726 -117.00000 804.3726
1850-12-31 694.3233 -55.67933 750.0026
1880-12-31 672.0366 -23.37148 695.4081
1910-12-31 679.3364 19.28471 660.0517
1920-12-31 287.7173 -49.49443 337.2117
1950-12-31 234.6592 -119.89309 354.5523
1970-12-31 333.0994 24.12508 308.9744
1990-12-31 462.2802 31.02266 431.2576
2010-12-31 522.5548 0.67881 521.8760 
1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
library(lattice)
library(gridExtra)
T3_1.xts <- as.xts(read.zoo(T3_1))
##### T3_1.xts #####
#1 National capital Wn 国民資本
#2 incl. Land 土地
#3 incl. Housing 住宅
#4 incl. Other domestic capital assets その他の国内資本資産
#5 incl. Net foreign capital 純外国資本
#6 Public capital Wg 公的資本
#7 incl. Public assets 公的資産
#8 incl. Public debt 公的債務
#9 Private capital W 民間資本
#png("Piketty3_1B.png",width=1200,height=800)
p1<-xyplot(T3_1.xts[,c(7,8)],superpose=TRUE,xlab="year",ylab="",main="イギリスの公的資本 1700-2010年",lwd=2,grid = TRUE, scales = list( y = list( rot = 0 )))
p2<-xyplot(T3_1.xts[,c(1,6,9)],superpose=TRUE,xlab="year",ylab="",main="イギリスの民間資本と公的資本 1700-2010年",lwd=2,grid = TRUE, scales = list( y = list( rot = 0 )))
grid.arrange(p1,p2,ncol=2)
#dev.off()

Capital in the France, 1700-2010
1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
29
30
31
32
33
34
library(xts)
library(Quandl)
T3_2<-Quandl("PIKETTY/T3_2")
#並べ替え
sortlist <- order(T3_2[,1])
dat <- T3_2[sortlist,]
T3_2<-dat
rownames(T3_2) <- c(1:nrow(T3_2))
#save("T3_2", file="T3_2.dat")
#load("T3_2.dat")
##### T3_2 #####
#1 Year 年
#2 National capital Wn 国民資本
#3 incl. Land 土地
#4 incl. Housing 住宅
#5 incl. Other domestic capital assets その他の国内資本資産
#6 incl. Net foreign capital 純外国資本
#7 Public capital Wg 公的資本
#8 incl. Public assets 公的資産
#9 incl. Public debt 公的債務
#10 Private capital W 民間資本
#11 incl. private assets 民間資産
#12 incl. private debt 民間負債
library("plotrix")
library(RColorBrewer)
#png("Piketty3_2A.png",width=1000,height=800)
stackpoly(T3_2[,3:6],ylim=c(0,800),axis4=F,main="フランスの資本 1700-2010年",col=brewer.pal(4,"Pastel1"),
xaxlab=c("1700","1750","1780","1810","1850","1880","1910","1920","1950","1970","1990","2000","2010"),border="gray36",staxx=F,stack=TRUE)
#boxed.labels(3,180,"土地",border=F,bg=NA,col="gray20")
#boxed.labels(4,400,"住宅",border=F,bg=NA,col="gray20")
#boxed.labels(4,600,"その他国内資本",border=F,bg=NA,col="gray20")
#boxed.labels(6,670,"純外国資本",border=F,bg=NA,col="gray20")
boxed.labels(c(3,4,4,6),c(180,400,600,670),c("土地","住宅","その他国内資本","純外国資本"),border=F,bg=NA,col="gray20",cex=2)
#dev.off()

1
2
library(knitr)
kable(T3_2[,c(1,8,9)]);kable(T3_2[,c(1,2,7,10)])
Year incl. Public assets incl. Public debt
1700-12-31 40 50
1750-12-31 40 70
1780-12-31 45 90
1810-12-31 45 14
1850-12-31 88 47
1880-12-31 89 92
1910-12-31 62 73
1920-12-31 70 73
1950-12-31 110 36
1970-12-31 90 40
1990-12-31 97 71
2000-12-31 119 91
2010-12-31 145 114
Year National capital Wn Public capital Wg Private capital W
1700-12-31 717 -10 727
1750-12-31 705 -30 735
1780-12-31 706 -45 751
1810-12-31 734 31 703
1850-12-31 720 41 679
1880-12-31 719 -3 722
1910-12-31 681 -10 692
1920-12-31 291 -3 293
1950-12-31 278 74 204
1970-12-31 363 50 311
1990-12-31 368 26 341
2000-12-31 503 28 474
2010-12-31 605 31 575 
1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
library(lattice)
library(gridExtra)
T3_2.xts <- as.xts(read.zoo(T3_2))
##### T3_2.xts #####
#1 National capital Wn 国民資本
#2 incl. Land 土地
#3 incl. Housing 住宅
#4 incl. Other domestic capital assets その他の国内資本資産
#5 incl. Net foreign capital 純外国資本
#6 Public capital Wg 公的資本
#7 incl. Public assets 公的資産
#8 incl. Public debt 公的債務
#9 Private capital W 民間資本
#10 incl. private assets 民間資産
#11 incl. private debt 民間負債
#png("Piketty3_2B.png",width=1200,height=800)
p1<-xyplot(T3_2.xts[,c(7,8)],superpose=TRUE,xlab="year",ylab="",main="フランスの公的資本 1700-2010年",lwd=2,grid = TRUE, scales = list( y = list( rot = 0 )))
p2<-xyplot(T3_2.xts[,c(1,6,9)],superpose=TRUE,xlab="year",ylab="",main="フランスの民間資本と公的資本 1700-2010年",lwd=2,grid = TRUE, scales = list( y = list( rot = 0 )))
grid.arrange(p1,p2,ncol=2)
#dev.off()