城投债（WIND口径）结构分析
载入数据：
setwd("C:/Users/shihchosen/Desktop")
library(ggplot2)
library(scales)
bond<-read.csv("城投债20150804.csv")
names(bond)
## [1] "证券代码" "证券简称"
## [3] "发行人中文名称" "债券期限.年...单位..年"
## [5] "起息日期" "到期日期"
## [7] "发行总额..单位..亿元" "息票品种"
## [9] "利率类型" "计息方式"
## [11] "票面利率.发行时." "税率..单位..."
## [13] "上市地点" "是否城投债"
## [15] "城投行政级别" "发行时债项评级"
## [17] "发行时主体评级"
head(bond)
## 证券代码 证券简称 发行人中文名称 债券期限.年...单位..年
## 1 7110.IB 02苏交通 江苏交通控股有限公司 15
## 2 120204.SH 02苏交通 江苏交通控股有限公司 15
## 3 120301.SH 03沪轨道 上海久事公司 15
## 4 038003.IB 03沪轨道债 上海久事公司 15
## 5 058009.IB 05粤交通债 广东省交通集团有限公司 15
## 6 111026.SZ 05粤交通 广东省交通集团有限公司 15
## 起息日期 到期日期 发行总额..单位..亿元 息票品种 利率类型 计息方式
## 1 2002-12-12 2017-12-12 15 附息 固定利率 单利
## 2 2002-12-12 2017-12-12 15 附息 固定利率 单利
## 3 2003-02-19 2018-02-19 40 附息 固定利率 单利
## 4 2003-02-19 2018-02-19 40 附息 固定利率 单利
## 5 2005-06-29 2020-06-29 15 附息 固定利率 单利
## 6 2005-06-29 2020-06-29 15 附息 固定利率 单利
## 票面利率.发行时. 税率..单位... 上市地点 是否城投债 城投行政级别
## 1 4.51 20 银行间 否
## 2 4.51 20 上海 否
## 3 4.51 20 上海 否
## 4 4.51 20 银行间 否
## 5 5.30 20 银行间 否
## 6 5.30 20 深圳 否
## 发行时债项评级 发行时主体评级
## 1 AAA AAA
## 2 AAA AAA
## 3 AAA AAA
## 4 AAA AAA
## 5 AAA AAA
## 6 AAA AAA
定义函数
# function isobar: 输入待分析列，标题，作出柱形图
isobar <- function(bond,xlab,ylab="频数",title="柱形图"){
df <- data.frame(x=factor(bond))
dfTab <- as.data.frame(table(df))
colnames(dfTab)[1] <- "x"
ggplot(df) +
geom_bar(aes(x,fill=x), color="black") +
geom_text(data=dfTab,aes(x=x,y=Freq,label=Freq),vjust=0) +
xlab(xlab)+ylab(ylab) +
ggtitle(title) +
theme(legend.position="none")
}
# function isobar: 输入待分析列，标题，作出累积柱形图#
cumbar<-function(bond,xlab,ylab="频数",title="柱形累积图"){
df <- data.frame(x=factor(bond))
dfTab <- as.data.frame(table(df))
colnames(dfTab)[1] <- "x"
ggplot(dfTab,aes(x="", y=Freq, fill=x)) + geom_bar(stat = "identity", color="black")+
geom_text(data=dfTab,aes(label=Freq),hjust = 0.5, vjust = 3,position = "stack")+
xlab(xlab)+ylab(ylab) +
ggtitle(title)
}
# function pieChart:画饼状图
pieChart <- function(bond,xlab,ylab="频数",title="饼状图"){
df <- data.frame(x=factor(bond))
dfTab <- as.data.frame(table(df))
colnames(dfTab)[1] <- "x"
ggplot(dfTab,aes(x="", y=Freq, fill=x)) + geom_bar(stat = "identity", color="black")+
xlab(xlab)+ylab(ylab) +
ggtitle(title)+
coord_polar("y", start=0)
}
债券期限
summary(bond[,4])
## Min. 1st Qu. Median Mean 3rd Qu. Max. NA's
## 0.2459 5.0000 7.0000 5.7970 7.0000 20.0000 2
p <- ggplot(bond, aes(factor(0), bond[,4]))
p + geom_boxplot()+ xlab("债券期限")+ ylab("单位（年）")+title("债券期限箱形图")
## Error in title("债券期限箱形图"): plot.new has not been called yet
发行总额
summary(bond[,7])
## Min. 1st Qu. Median Mean 3rd Qu. Max. NA's
## 0.035 6.000 10.000 11.340 15.000 100.000 6
summary(bond[,7],na.rm = TRUE)
## Min. 1st Qu. Median Mean 3rd Qu. Max. NA's
## 0.035 6.000 10.000 11.340 15.000 100.000 6
p <- ggplot(bond, aes(factor(0), bond[,7]))
p + geom_boxplot()+ xlab("发行总额")+ ylab("单位（亿元）")+ title("发行总额箱形图")
## Error in title("发行总额箱形图"): plot.new has not been called yet
息票品种
summary(bond[,8])
## 到期一次还本付息 附息
## 2 453 5056
isobar(bond[,8],xlab="息票品种", ylab="频率", title="息票品种柱形图")
pieChart(bond[,8],xlab="息票品种", ylab="频率", title="息票品种柱形图")
利率类型
summary(bond[,9])
## 浮动利率 固定利率 累进利率
## 2 33 4717 759
isobar(bond[,9],xlab="利率类型", ylab="频率", title="利率类型柱形图")
pieChart(bond[,9],xlab="利率类型", ylab="频率", title="利率类型柱形图")
计息方式
summary(bond[,10])
## 单利
## 2 5509
isobar(bond[,10],xlab="计息方式", ylab="频率", title="计息方式柱形图")
pieChart(bond[,10],xlab="计息方式", ylab="频率", title="计息方式柱形图")
票面利率
summary(bond[,11])
## Min. 1st Qu. Median Mean 3rd Qu. Max. NA's
## 2.980 5.800 6.610 6.579 7.310 12.000 7
p <- ggplot(bond, aes(factor(0), bond[,4]))
p + geom_boxplot()+ xlab("票面利率")+ ylab("单位（%）")+title("票面利箱形图")
## Error in title("票面利箱形图"): plot.new has not been called yet
税率
城投债参考公司债的税率20%
summary(bond[,12])
## Min. 1st Qu. Median Mean 3rd Qu. Max. NA's
## 20 20 20 20 20 20 2
上市地点
summary(bond[,13])
## 其他 上海 深圳
## 2 1 1445 41
## 天津股权交易所 银行柜台债券 银行间 浙江股权交易中心
## 10 18 3971 23
isobar(bond[,13],xlab="上市地点", ylab="频率", title="上市地点柱形图")
pieChart(bond[,13],xlab="上市地点", ylab="频率", title="上市地点饼状图")
是否城投债
summary(bond[,14])
## 否 是
## 2 460 5049
isobar(bond[,14],xlab="是否城投债", ylab="频率", title="是否城投债柱形图")
pieChart(bond[,14],xlab="是否城投债", ylab="频率", title="是否城投债饼状图")
城投行政级别
summary(bond[,15])
## 地级市 省及省会(单列市) 县及县级市
## 1477 1524 1982 528
isobar(bond[,15],xlab="城投行政级别", ylab="频率", title="城投行政级别柱形图")
pieChart(bond[,15],xlab="城投行政级别", ylab="频率", title="城投行政级别债饼状图")
发行时债项评级
summary(bond[,16])
## - A-1 A+ AA AA- AA+ AAA
## 2 1023 297 2 2300 58 1263 566
isobar(bond[,16],xlab="债项评级", ylab="频率", title="发行时债项评级柱形图")
pieChart(bond[,16],xlab="债项评级", ylab="频率", title="发行时债项评级饼状图")
发行时主体评级
summary(bond[,17])
## A A- A+ AA AA- AA+ AAA BBB BBB-
## 123 4 5 12 2807 407 1265 885 1 2
isobar(bond[,17],xlab="主体评级", ylab="频率", title="发行时主体评级柱形图")
pieChart(bond[,17],xlab="主体评级", ylab="频率", title="发行时主体评级饼状图")
没有评论:
发表评论