python 实现马科维茨_Python实现量化选股

什么是选股?

选股(stock selection)是一种主动性投资策略,先按照某种规则或算法分析单只股票的前景,然后构建一个投资组合,长期持有。一般情况下要求组合的股票具有低相关性,这样才能对冲系统性风险,否则在大盘走弱的时候投资组合也会面临巨大的下跌风险。

运用什么模型?

关于如何选股,学术界提出过很多不同的模型,最经典的莫过于马科维茨投资组合理论。这里我们使用MM趋势模型(Mark Minervini’s Trend Template),这是国外一位传奇投资大师提出的技术面选股方法,核心思想是通过技术指标来度量股票动能,从中筛选最有潜力的股票,买入并持有。

MM趋势模型

  1. 股票价格高于150天均线和200天均线
  2. 150日均线高于200日均线
  3. 200日均线上升至少1个月
  4. 50日均线高于150日均线和200日均线
  5. 股票价格高于50日均线
  6. 股票价格比52周低点高30%
  7. 股票价格在52周高点的25%以内
  8. 相对强弱指数(RS)大于等于70,这里的相对强弱指的是股票与大盘对比,RS = 股票1年收益率 / 基准指数1年收益率

关于Mark Minervini

全美最富盛名的交易员之一,曾经获得30000%的收益率,在34岁前称为亿万富翁,详情见<金融怪杰>一书。

选股面临的技术性难题?

  1. 从哪里获取大量股票的历史数据?
  2. 当股票数量很多时,如何提高计算性能?

本文将用Python实现MM模型的量化选股,并解决上述提出的两个技术难题。

import os
import datetime as dt
import time
from typing import Any, Dict, Optional, List

import requests
import pickle
import numpy as np
import pandas as pd
import matplotlib.pyplot as plt
import seaborn as sns
import talib
import multiprocessing as mp
from requests.exceptions import ConnectionError, Timeout

%matplotlib inline
plt.style.use("fivethirtyeight")

1. 从蜂鸟数据获取历史数据

蜂鸟数据是新兴的金融数据提供商,提供包括股票,外汇,商品期货和数字货币的实时报价和历史数据,并提供API接口,是所有金融从业者获取免费数据的便捷渠道。

## 撰写自定义函数,通过API获取数据

def fetch_trochil(url: str,
                  params: Dict[str, str],
                  attempt: int = 3,
                  timeout: int = 3) -> Dict[str, Any]:
    """装饰requests.get函数"""
    for i in range(attempt):
        try:
            resp = requests.get(url, params, timeout=timeout)
            resp.raise_for_status()
            data = resp.json()["data"]
            if not data:
                raise Exception("empty dataset")
            return data
        except (ConnectionError, Timeout) as e:
            print(e)
            i += 1
            time.sleep(i * 0.5)


def fetch_cnstocks(apikey: str) -> pd.DataFrame:
    """从蜂鸟数据获取A股产品列表"""
    url = "https://api.trochil.cn/v1/cnstock/markets"
    params = {
    "apikey": apikey}

    res = fetch_trochil(url, params)

    return pd.DataFrame.from_records(res)


def fetch_daily_ohlc(symbol: str,
                     date_from: dt.datetime,
                     date_to: dt.datetime,
                     apikey: str) -> pd.DataFrame:
    """从蜂鸟数据获取A股日图历史K线"""
    url = "https://api.trochil.cn/v1/cnstock/history"
    params = {
    
        "symbol": symbol,
        "start_date": date_
实现马科维茨有效边界图,需要以下步骤: 1. 安装必要的Python库,如numpy、pandas、matplotlib等。 2. 创建一个投资组合的收益率和波动率数据集。可以使用历史数据或模拟数据。 3. 计算投资组合的预期收益率、预期波动率和协方差矩阵。 4. 使用numpy中的linspace函数生成各个预期收益率水平下的最小方差投资组合的权重组合。 5. 计算每个最小方差投资组合的预期收益率和波动率。 6. 使用matplotlib中的scatter函数绘制最小方差投资组合的预期收益率和波动率散点图。 7. 使用scipy.optimize库中的minimize函数计算每个预期收益率水平下的最优投资组合。 8. 计算每个最优投资组合的预期收益率和波动率。 9. 使用matplotlib中的plot函数绘制最优投资组合的有效边界。 下面是示例代码: ```python import numpy as np import pandas as pd import matplotlib.pyplot as plt from scipy.optimize import minimize # 创建数据集 returns = np.random.normal(0.05, 0.1, 100) volatility = np.random.normal(0.1, 0.05, 100) data = pd.DataFrame({'returns': returns, 'volatility': volatility}) # 计算预期收益率、预期波动率和协方差矩阵 mu = data['returns'].mean() sigma = data['volatility'].std() cov_matrix = np.cov(data['returns'], data['volatility']) # 生成权重组合 num_portfolios = 1000 results = np.zeros((3 + len(data.columns), num_portfolios)) for i in range(num_portfolios): weights = np.random.random(len(data.columns)) weights /= np.sum(weights) portfolio_return = np.dot(weights, data['returns']) portfolio_volatility = np.sqrt(np.dot(weights.T, np.dot(cov_matrix, weights))) results[0, i] = portfolio_return results[1, i] = portfolio_volatility results[2, i] = portfolio_return / portfolio_volatility for j in range(len(weights)): results[j+3, i] = weights[j] # 绘制最小方差投资组合散点图 min_volatility_idx = np.argmin(results[1]) min_volatility_return = results[0, min_volatility_idx] min_volatility_volatility = results[1, min_volatility_idx] plt.scatter(results[1,:], results[0,:], c=results[2,:], cmap='YlOrRd', alpha=0.5) plt.scatter(min_volatility_volatility, min_volatility_return, c='b', marker='*', s=500) plt.xlabel('Volatility') plt.ylabel('Returns') plt.colorbar() # 计算最优投资组合 def objective(weights): portfolio_return = np.dot(weights, data['returns']) portfolio_volatility = np.sqrt(np.dot(weights.T, np.dot(cov_matrix, weights))) return -portfolio_return / portfolio_volatility constraints = ({'type': 'eq', 'fun': lambda x: np.sum(x) - 1}) bounds = [(0, 1)] * len(data.columns) initial_weights = np.ones(len(data.columns)) / len(data.columns) optimal_weights = minimize(objective, initial_weights, method='SLSQP', bounds=bounds, constraints=constraints) optimal_weights = optimal_weights.x # 计算最优投资组合的预期收益率和波动率 optimal_return = np.dot(optimal_weights, data['returns']) optimal_volatility = np.sqrt(np.dot(optimal_weights.T, np.dot(cov_matrix, optimal_weights))) # 绘制最优投资组合的有效边界 frontier_returns = [] frontier_volatilities = [] for r in np.linspace(min(results[0,:]), max(results[0,:]), num=100): constraints = ({'type': 'eq', 'fun': lambda x: np.sum(x) - 1}, {'type': 'eq', 'fun': lambda x, r=r: np.dot(x, data['returns']) - r}) bounds = [(0, 1)] * len(data.columns) initial_weights = np.ones(len(data.columns)) / len(data.columns) optimal_weights = minimize(objective, initial_weights, method='SLSQP', bounds=bounds, constraints=constraints) optimal_weights = optimal_weights.x frontier_returns.append(r) frontier_volatilities.append(np.sqrt(np.dot(optimal_weights.T, np.dot(cov_matrix, optimal_weights)))) plt.plot(frontier_volatilities, frontier_returns, 'g--', linewidth=2) plt.show() ``` 此代码将生成一个包含最小方差投资组合散点图和有效边界图的matplotlib图形。
评论
添加红包

请填写红包祝福语或标题

红包个数最小为10个

红包金额最低5元

当前余额3.43前往充值 >
需支付:10.00
成就一亿技术人!
领取后你会自动成为博主和红包主的粉丝 规则
hope_wisdom
发出的红包
实付
使用余额支付
点击重新获取
扫码支付
钱包余额 0

抵扣说明:

1.余额是钱包充值的虚拟货币,按照1:1的比例进行支付金额的抵扣。
2.余额无法直接购买下载,可以购买VIP、付费专栏及课程。

余额充值