如何用python利用Gamma-Normal方法计算VaR。 假设股票收益率样本内外数据量比为1/2，用样本内数据估计模型参数，再进行滚动窗口预测，在α=5%的前提下计算VaR，并画出VaR的N个预测值的时间序列图。

首先，我们需要安装scipy和numpy库。

然后，我们可以使用以下代码计算Gamma-Normal方法的VaR：

import numpy as np
from scipy.stats import gamma, norm

def gamma_normal_var(alpha, mu, sigma, gamma_shape, gamma_scale, data):
    """
    计算Gamma-Normal方法的VaR
    :param alpha: 置信水平
    :param mu: 正态分布的均值
    :param sigma: 正态分布的标准差
    :param gamma_shape: Gamma分布的形状参数
    :param gamma_scale: Gamma分布的尺度参数
    :param data: 数据
    :return: VaR
    """
    z_alpha = norm.ppf(alpha)
    q_alpha = gamma.ppf(alpha, gamma_shape, scale=gamma_scale)
    var = mu + sigma * z_alpha + q_alpha * (np.exp(mu + sigma * z_alpha) - 1)
    return var

接下来，我们可以使用以下代码来进行滚动窗口预测，并且计算VaR，并画出时间序列图：

import pandas as pd

data = pd.read_csv("data.csv")  # 读取数据
data = data["returns"]  # 只使用收益率数据

n = len(data)  # 数据量
window_size = n // 3  # 窗口大小

# 使用前两个窗口的数据估计模型参数
mu = data[:window_size * 2].mean()
sigma = data[:window_size * 2].std()
gamma_shape, _, gamma_scale = gamma.fit(data[:window_size * 2], floc=0)

var_list = []  # 保存VaR的列表
for i in range(window_size * 2, n):
    # 计算当前窗口的VaR
    var = gamma_normal_var(0.05, mu, sigma, gamma_shape, gamma_scale, data[i - window_size * 2:i])
    var_list.append(var)
    
    # 更新模型参数
    mu = data[i - window_size:i].mean()
    sigma = data[i - window_size:i].std()
    gamma_shape, _, gamma_scale = gamma.fit(data[i - window_size:i], floc=0)

# 画出时间序列图
import matplotlib.pyplot as plt

plt.plot(data[window_size * 2:], label="returns")
plt.plot(var_list, label="VaR")
plt.legend()
plt.show()

希望这个代码能够帮助你计算Gamma-Normal方法的VaR，并画出时间序列图。