前面介绍过《基于Python的指数基金量化投资-通过市盈率和市净率对指数估值》，通过估值来进行指数基金的投资。
当估值百分位为90%-100%的时候，指数处于严重高估状态，可以考虑全仓卖出；
当估值百分位为80%-90%的时候，指数处于高估状态，可以考虑大份额卖出；
当估值百分位为60%-80%的时候，指数处于正常偏高状态，可以考虑小份额卖出；
当估值百分位为40%-60%的时候，指数处于正常状态，可以考虑持有；
当估值百分位为20%-40%的时候，指数处于正常偏低状态，可以考虑小份额定投；
当估值百分位为10%-20%的时候，指数处于低估状态，可以考虑大份额定投；
当估值百分位为0%-10%的时候，指数处于严重低估状态，可以考虑大份额定投+增量买入；
有了这个策略，就需要观察每个指数的估值百分位情况，如果一个一个查看会很麻烦，最简单的办法就是把各个指数直观的通过图形的方式表现出来，如下图所示：
图中按照严重低估、低估、正常偏低、正常、正常偏高、高估和严重高估划分了对应的区域，并用不同的颜色进行了区分，然后分别计算每一个指数的估值百分位，然后把计算出来的结果画在相应的区域，这样就可以非常直观的看到所有指数的一个估值情况。
源码
# 源码中用到的估值数据文件g_*.csv可以联系小将获取。
import pandas as pd
import numpy as np
import matplotlib.pyplot as plt
index_name_csv = ['g_hs300.csv', #沪深300 - 0
'g_zz500.csv', #中证500 - 1
'g_zz100.csv', #中证100 - 2
'g_shz50.csv', #上证50 - 3
'g_hsyy300.csv', #沪深医药300 - 4
'g_zzyh.csv', #中证银行 -5
'g_zzxf.csv', #中证消费 -6
'g_zzbj.csv', #中证白酒 -7
'g_db500.csv', # 500低波动-8
'g_jz300.csv', # 300价值 -9
'g_yy100.csv', #医药100 -10
'g_zzyyao.csv', #中证医药 -11
'g_jbm50.csv', #基本面50 -12
'g_shzhl.csv', #上证红利 -13
'g_zzhl.csv', #中证红利 -14
'g_zzjg.csv', #中证军工 -15
'g_spyl.csv', #食品饮料 -16
'g_zqgs.csv', #证券公司 -17
'g_ylcy.csv', #养老产业 -18
'g_szhl.csv', #深证红利 -19
'g_zzhb.csv', #中证环保 -20
'g_cyb.csv', #创业板 -21
'g_hszs.csv', #恒生指数 -22
'g_hsgqzs.csv', #恒生国企指数 -23
'g_zghl50.csv', #中国互联50 -24
'g_xgdp.csv', #香港大盘 -25
'g_xgzx.csv'] #香港中小 -26
index_info = np.zeros([len(index_name_csv),1])
for i in range(len(index_name_csv)):
index_data = pd.read_csv('./importfile/indexSeries/indexValuation/g/' +index_name_csv[i])
if index_name_csv[i] == 'g_zzyh.csv' or index_name_csv[i] =='g_zqgs.csv':
data_calc = index_data['pb']
else:
data_calc = index_data['pe']
xx = np.where(data_calc < data_calc.values[-1])
data_percentage = len(xx[0]) / data_calc.shape[0]
if index_name_csv[i] == 'g_zzyh.csv':
data_calc = index_data['pe']
xx = np.where(data_calc < data_calc.values[-1])
data_percentage_add = len(xx[0]) / data_calc.shape[0]
index_info[i] = 0.8*data_percentage + 0.2*data_percentage_add
else:
index_info[i] = data_percentage
plt.figure(3)
plt.rcParams['font.sans-serif'] =['Microsoft YaHei']
plt.plot(index_info, '--ow', ms=12,linewidth=2)
plt.axhspan(0, 0.1, facecolor='#4D8070')
plt.axhspan(0.1, 0.2, facecolor='#21DB85')
plt.axhspan(0.2, 0.4, facecolor='#B0D91C')
plt.axhspan(0.4, 0.6, facecolor='#E3BD00')
plt.axhspan(0.6, 0.8, facecolor='#FF9421')
plt.axhspan(0.8, 0.9, facecolor='#F56600')
plt.axhspan(0.9, 1.4, facecolor='#B22222')
plt.xlim(-1,len(index_info))
plt.ylim(-0.15, 1.4)
scale_y = 0.07
offset_x = 0.25
offset_y = 0.062
offset_x_num = 0.1
offset_x_alpha = 0.15
font_size = 22
index_name = ['沪深300',
'中证500',
'中证100',
'上证50',
'沪深医药300',
'中证银行',
'中证消费',
'中证白酒',
'500低波动',
'300价值',
'医药100',
'中证医药',
'基本面50',
'上证红利',
'中证红利',
'中证军工',
'食品饮料',
'证券公司',
'养老产业',
'深证红利',
'中证环保',
'创业板',
'恒生指数',
'H股指数',
'中国互联50',
'香港大盘',
'香港中小']
font = {'size': font_size, 'color': 'w','weight': 'black'}
for i in range(len(index_name)):
index_name_word = index_name[i]
for j in range(len(index_name_word)):
if index_name_word[j].isdigit():
plt.text(i - offset_x + offset_x_num, index_info[i] - j * scale_y +len(index_name_word) * offset_y,
index_name_word[j],fontdict=font)
else:
plt.text(i - offset_x, index_info[i] - j * scale_y +len(index_name_word) * offset_y, index_name_word[j],
fontdict=font)
time_end = '2020/12/21'
plt.plot([-1,len(index_info)],[-0.025,-0.025],color='#CCCCCC',linewidth=2,linestyle='--')
shift_x = 4.23
text_base_x = -0.8
text_base_shift_x = 0.55
plt.gca().add_patch(plt.Rectangle((-1+0*shift_x,-0.15),2.6,0.1,color='#4D8070'))
plt.text(text_base_x+0*shift_x,-0.123,'严重低估',fontdict=font)
plt.gca().add_patch(plt.Rectangle((-1+1*shift_x,-0.15),2.6,0.1,color='#21DB85'))
plt.text(text_base_x+text_base_shift_x+1*shift_x,-0.123,'低估',fontdict=font)
plt.gca().add_patch(plt.Rectangle((-1+2*shift_x,-0.15),2.6,0.1,color='#B0D91C'))
plt.text(text_base_x+2*shift_x,-0.123,'正常偏低',fontdict=font)
plt.gca().add_patch(plt.Rectangle((-1+3*shift_x,-0.15),2.6,0.1,color='#E3BD00'))
plt.text(text_base_x+text_base_shift_x+3*shift_x,-0.123,'正常',fontdict=font)
plt.gca().add_patch(plt.Rectangle((-1+4*shift_x,-0.15),2.6,0.1,color='#FF9421'))
plt.text(text_base_x+4*shift_x,-0.123,'正常偏高',fontdict=font)
plt.gca().add_patch(plt.Rectangle((-1+5*shift_x,-0.15),2.6,0.1,color='#F56600'))
plt.text(text_base_x+text_base_shift_x+5*shift_x,-0.123,'高估',fontdict=font)
plt.gca().add_patch(plt.Rectangle((-1+6*shift_x,-0.15),2.6,0.1,color='#B22222'))
plt.text(text_base_x+6*shift_x,-0.123,'严重高估',fontdict=font)
plt.title(time_end,size=28)
plt.axis('off')
g_globalMarket =pd.read_csv('./importfile/indexSeries/indexValuation/g_globalMarket/g_globalMarket.csv')
index_data_pe =g_globalMarket['pe'].values[1:len(g_globalMarket['pe']):1]
index_data_pb = g_globalMarket['pb'].values[1:len(g_globalMarket['pb']):1]
size_title = 28
size_label = 23
size_text = 35
size_line = 3
size_rotation = 20
plt_gap = 10
plt.show()
程序中用到的数据如果有问题，大家可以留言获取，欢迎大家一起交流沟通^_^
课程参考：网易云课堂 基于Python的量化指数基金投资
