首页
学习
活动
专区
工具
TVP
发布
精选内容/技术社群/优惠产品,尽在小程序
立即前往

使用R语言随机波动模型SV处理时间序列中的随机波动率

此函数仅产生SV流程的实现,并返回svsim类的对象,该对象具有自己的print,summary和plot方法。 下面给出了使用svsim的示例代码,该模拟实例显示在图2中。...可以看出,该函数调用主MCMC采样器并将其输出转换为与coda兼容的对象。后者的完成主要是出于兼容性的考虑,并且可以直接访问在那里实施的收敛诊断检查。...图6显示了从汇率提取数据集中获得的EUR-USD汇率的示例输出。 svdraws对象的通用绘图方法将上述所有绘图合并为一个绘图。可以使用上述所有参数。请参见图7。...R> plot(res, showobs = FALSE)  为了提取标准化残差,可以在给定的svdraws对象上使用残差/残差方法。使用可选的参数类型,可以指定摘要统计的类型。...当前,类型允许为“平均值”或“中位数”,其中前者对应于默认值。此方法返回svresid类的实向量,其中包含每个时间点所请求的标准化残差的摘要统计量。

1.9K10
  • 您找到你想要的搜索结果了吗?
    是的
    没有找到

    PYTHON 用几何布朗运动模型和蒙特卡罗MONTE CARLO随机过程模拟股票价格可视化分析耐克NKE股价时间序列数据|附代码数据

    这是使用 Python 中的几个函数完成的,并使用迭代设置将后续股票价格建模为马尔可夫链,给定初始起始价格 S0。...可以很容易地观察到,收益呈正态分布,平均值约为 0.15,这是模拟中 μ 的输入值。这种收益率的正态分布也是布朗运动模型的预期结果。...两年期末的预期价格 St 预计在 98.6 左右。 两年年化收益率为 μ=0.13,这是使用雅虎财经历史记录生成的价格列表计算得出的。 两年年化波动率也从价格中获得,σ=0.05。...以下代码调用函数以通过几何布朗运动模型生成随机游走。 #使用这些价格计算回报率和波动率。...pprint(k.geical('2013-01-01', '2015-01-01')) #使用输入参数的样本值运行多个模拟dt = 0.01sim_count = 500#调用函数并运行模拟prie

    84411

    PYTHON 用几何布朗运动模型和蒙特卡罗MONTE CARLO随机过程模拟股票价格可视化分析耐克NKE股价时间序列数据|附代码数据

    这是使用 Python 中的几个函数完成的,并使用迭代设置将后续股票价格建模为马尔可夫链,给定初始起始价格 S0。...可以很容易地观察到,收益呈正态分布,平均值约为 0.15,这是模拟中 μ 的输入值。这种收益率的正态分布也是布朗运动模型的预期结果。...两年期末的预期价格 St 预计在 98.6 左右。 两年年化收益率为 μ=0.13,这是使用雅虎财经历史记录生成的价格列表计算得出的。 两年年化波动率也从价格中获得,σ=0.05。...以下代码调用函数以通过几何布朗运动模型生成随机游走。 #使用这些价格计算回报率和波动率。...pprint(k.geical('2013-01-01', '2015-01-01')) #使用输入参数的样本值运行多个模拟dt = 0.01sim_count = 500#调用函数并运行模拟prie

    1.4K30

    CVPR2020 | 显著性目标检测,多尺度信息相互融合

    另外,考虑到卷积运算固有的局部性和交叉熵函数的像素级特征,很难实现物体的均匀显著性提取。 ? 对于第一个问题,现有方法的主要解决方法是逐层整合较浅的特征。...为了进一步从提取的特征中获得丰富的尺度特异性信息,本文还设计了一个自交互模块(SIM)(图2(h))。...两个不同分辨率的交互分支被训练成从单个卷积块中学习多尺度特征, AIM和SIM有效地提高了SOD任务中处理尺度变化的能力。...在经过上采样、归一化和非线性处理后,采用FU对SIM和残差支路的双路径进行处理,将SIM集成到解码器中,使得网络在训练阶段能够自适应地处理不同样本的尺度变化。 ? ?...首先,最终预测的计算如下: ? 为了解决各种尺度引起的前/后台失衡问题,损失函数至少需要满足两个要求:1)它比背景更多地关注前景,而对象尺度上的差异并不能引起广泛的影响。

    3.5K40

    R语言乘法GARCH模型对高频交易数据进行波动性预测

    考虑的日内时间范围是09:30至16:00,即证券交易所的正式交易时间。与大多数此类关于日内数据建模的研究一样,当天的第一个收益被删除。每日数据从雅虎财经下载。...预测波动率 sigma_t是由外部的多因素风险模型得出的,但是同样可以从每日 GARCH 模型中产生这种预测。该过程的季节性部分定义为:。 ? 用残差除以波动率和日波动率,就得到了归一化残差。 ?...这是一个xts对象,也可以选择有m.sim列,这样每个独立的模拟都是基于日方差独立模拟的调整残差。下面的示例代码显示了对未来1分钟间隔的10,000个点的模拟,并说明了季节性成分的影响。...# dtime包含了数据集中唯一的区间点的集合 # (可从mcsGARCH模型的所有rugarch对象中获得) sim = sim(fit, n.sim = 10000,Var = var_sim) #...滚动的回测和风险值 ugarchroll函数对于在回测应用中测试模型的充分性非常有用,下面的代码说明了模型在数据期间的情况。

    1.5K20

    PYTHON 用几何布朗运动模型和蒙特卡罗MONTE CARLO随机过程模拟股票价格可视化分析耐克NKE股价时间序列数据|附代码数据

    这是使用 Python 中的几个函数完成的,并使用迭代设置将后续股票价格建模为马尔可夫链,给定初始起始价格 S0。...可以很容易地观察到,收益呈正态分布,平均值约为 0.15,这是模拟中 μ 的输入值。这种收益率的正态分布也是布朗运动模型的预期结果。...两年期末的预期价格 St 预计在 98.6 左右。两年年化收益率为 μ=0.13,这是使用雅虎财经历史记录生成的价格列表计算得出的。两年年化波动率也从价格中获得,σ=0.05。...以下代码调用函数以通过几何布朗运动模型生成随机游走。#使用这些价格计算回报率和波动率。...pprint(k.geical('2013-01-01', '2015-01-01'))#使用输入参数的样本值运行多个模拟dt = 0.01sim_count = 500#调用函数并运行模拟prie,

    1.2K00

    R语言GARCH模型对股市sp500收益率bootstrap、滚动估计预测VaR、拟合诊断和蒙特卡罗模拟可视化|附代码数据

    这测试了标准化残差中杠杆效应的存在(以捕捉 GARCH 模型可能的错误指定), 其中 I 是指标函数, ^t 是 GARCH 过程的估计残差。...该函数必须采用 2 个参数,即拟合对象以及用于对值进行分类的箱数。...ititin="std") > fi = grit(sc,sp5et) bot(fit, mehod = c("Pl", "Full")[1], 完整 GARCH bootstrap程序总结如下: 从估计对象中提取标准化残差...如果是具有固定参数的规范,首先使用提供的数据集进行过滤,然后从过滤后的对象中提取标准化残差。...使用 spd 或基于内核的方法从原始标准化残差中采样大小为 N 的 n.bootfit 集(原始数据集减去任何样本周期外)。

    35800

    R语言线性混合效应模型(固定效应&随机效应)和交互可视化3案例|附代码数据

    测量斑块长度 这第一个数据集是从Griffith和Sheldon(2001年,《动物行为学》61:987-993)的一篇论文中提取的,他们在两年内对瑞典哥特兰岛上的30只雄性领头鶲的白色额斑进行了测量。...重复性的解释如何改变? 从保存的lmer对象中提取参数估计值(系数)。检查随机效应的输出。随机变异的两个来源是什么?固定效应指的是什么? 在输出中,检查随机效应的标准差。...这两个来源中的哪一个对应于"(截距)",哪一个对应于 "残差"? 同时检查固定效应结果的输出。模型公式中唯一的固定效应是所有长度测量的平均值。它被称为"(截距)",但不要与随机效应的截距相混淆。...固定效应输出给了你平均值的估计值和该估计值的标准误差。注意固定效应输出是如何提供均值估计值的,而随机效应输出则提供方差(或标准差)的估计值。 从拟合模型中提取方差分量,估计各年斑块长度的可重复性*。...每条鱼的预测值和观察值之间的差异代表残差。 你在(1)中做了什么假设?创建一个残差与拟合值的图,以检查这些假设之一。 从保存的lmer对象中提取参数估计值。检查固定效应的结果。

    1.2K30

    R语言线性混合效应模型(固定效应&随机效应)和交互可视化3案例

    ---- 测量斑块长度 这第一个数据集是从Griffith和Sheldon(2001年,《动物行为学》61:987-993)的一篇论文中提取的,他们在两年内对瑞典哥特兰岛上的30只雄性领头鶲的白色额斑进行了测量...重复性的解释如何改变? 从保存的lmer对象中提取参数估计值(系数)。检查随机效应的输出。随机变异的两个来源是什么?固定效应指的是什么? 在输出中,检查随机效应的标准差。...这两个来源中的哪一个对应于"(截距)",哪一个对应于 "残差"? 同时检查固定效应结果的输出。模型公式中唯一的固定效应是所有长度测量的平均值。它被称为"(截距)",但不要与随机效应的截距相混淆。...固定效应输出给了你平均值的估计值和该估计值的标准误差。注意固定效应输出是如何提供均值估计值的,而随机效应输出则提供方差(或标准差)的估计值。 从拟合模型中提取方差分量,估计各年斑块长度的可重复性*。...每条鱼的预测值和观察值之间的差异代表残差。 你在(1)中做了什么假设?创建一个残差与拟合值的图,以检查这些假设之一。 从保存的lmer对象中提取参数估计值。检查固定效应的结果。

    8.8K61

    R语言线性混合效应模型(固定效应&随机效应)和交互可视化3案例|附代码数据

    测量斑块长度 这第一个数据集是从Griffith和Sheldon(2001年,《动物行为学》61:987-993)的一篇论文中提取的,他们在两年内对瑞典哥特兰岛上的30只雄性领头鶲的白色额斑进行了测量。...重复性的解释如何改变? 从保存的lmer对象中提取参数估计值(系数)。检查随机效应的输出。随机变异的两个来源是什么?固定效应指的是什么? 在输出中,检查随机效应的标准差。...这两个来源中的哪一个对应于"(截距)",哪一个对应于 "残差"? 同时检查固定效应结果的输出。模型公式中唯一的固定效应是所有长度测量的平均值。它被称为"(截距)",但不要与随机效应的截距相混淆。...固定效应输出给了你平均值的估计值和该估计值的标准误差。注意固定效应输出是如何提供均值估计值的,而随机效应输出则提供方差(或标准差)的估计值。 从拟合模型中提取方差分量,估计各年斑块长度的可重复性*。...每条鱼的预测值和观察值之间的差异代表残差。 你在(1)中做了什么假设?创建一个残差与拟合值的图,以检查这些假设之一。 从保存的lmer对象中提取参数估计值。检查固定效应的结果。

    1.7K00

    实例复习机器学习数学 - 2. 几种典型离散随机变量分布

    ,随机 100000 次 binom_sim = data = binom.rvs(n=10, p=0.3, size=100000) #计算本次平均值以及标准差 print ("Mean: %g" %...(binom_sim, bins=10, density=True, stacked=True) plt.show() 可以看到输出为: 可以看到平均值很接近期望,标准差也是很接近我们用公式算出来的方差开根...,随机 10000 次 poisson_sim = poisson.rvs(mu=5, size=10000) #计算本次平均值以及标准差 print ("Mean: %g" % np.mean(poisson_sim...)) print ("Std: %g" % np.std(poisson_sim, ddof=1)) #绘图,利用归一化的直方图,这样展示的很接近概率质量函数 plt.hist(poisson_sim,...,随机 100000 次 geom_rv = geom(p=0.5) geom_sim = geom_rv.rvs(size=100000) #计算本次平均值以及标准差 print ("Mean: %g

    71320

    最全的NumPy教程

    从ndarray对象提取的任何元素(通过切片)由一个数组标量类型的 Python 对象表示。下图显示了ndarray,数据类型对象(dtype)和数组标量类型之间的关系。 ?...基本的ndarray是使用 NumPy 中的数组函数创建的,如下所示: numpy.array 它从任何暴露数组接口的对象,或从返回数组的任何方法创建一个ndarray。...通过将start,stop和step参数提供给内置的slice函数来构造一个 Python slice对象。此slice对象被传递给数组来提取数组的一部分。...函数说明如下: numpy.amin() 和 numpy.amax() 这些函数从给定数组中的元素沿指定轴返回最小值和最大值。...加权平均值 = (1*4+2*3+3*2+4*1)/(4+3+2+1) 标准差 标准差是与均值的偏差的平方的平均值的平方根。

    4.2K10

    使用时空-频率模式分析从脑电数据的一些试验中提取N400成分

    因此,从EEG信号中准确提取N400波形具有重要意义。 N400的信噪比(SNR)非常差,这是由于N400的低振幅以及自发的脑电图活动和眼伪影的存在。...首先,由于N400成分具有锁相特性,可以分别从每个数据集中随机重采样获得一致波形和不一致波形,然后取平均值。从而使N400的差分波形具有更高的锁相分量信噪比。...STPA提取的N400成分反投影到每个对象的Cz电极上。恢复的N400波形具有比Cz电极上平均波形更高的SNR。...N400时间演变分析 为了分析N400的时间演变,将每个受试者的数据集分为两组。第一组由实验中的前20个试验组成,第二组由后20个试验组成。从两组中提取ERPs。...结论 清华大学高小榕团队成员在本研究提出了STPA方法,该方法可以从脑电图数据中提取N400成分。仿真和真实脑电图数据分析结果表明,STPA优于ESSPs、SIM和r-ICA。

    86410

    「R」用purrr实现迭代

    因为R是一门函数式编程语言,我们可以先将for循环包装在函数中,然后再调用函数,而不是使用for循环,因此for循环在R中不像在其他编程语言中那么重要。...接下来我们将学习和使用purrr包,它提供的函数可以替代很多常见的for循环应用。R基础包中的apply应用函数族也可以完成类似的任务,但purrr包的函数更一致,也更容易学习。...当检查多个模型时,有时候我们需要提取像R方这样的摘要统计量,要想完成这个任务,我们需要先运行summary()函数,然后提取结果中的r.squared: models %>% map(summary...这样,修改后的函数就不会抛出错误,相反,它总是返回由下面两个元素组成的列表: result - 原始结果。如果出现错误,那么它就是NULL error - 错误对象。...,一般来说,我们应该检查一下y中错误对象所对应的x值,或者使用y中的正常结果进行一些处理: is_ok = y$error %>% map_lgl(is_null) x[!

    4.8K20

    R语言中的copula GARCH模型拟合时间序列并模拟分析|附代码数据

    最近我们被客户要求撰写关于copula GARCH的研究报告,包括一些图形和统计输出。 在这个文章中,我们演示了copula GARCH方法(一般情况下) 1 模拟数据 首先,我们模拟一下创新分布。...(或使用,例如,"std") ## 使用创新模拟ARMA-GARCH模型 ## 注意: ugarchpath(): 从spec中模拟; garchpath(uspec,                 ...n.sim = n, # 模拟的路径长度                 ## 提取结果系列 X. <- fitted(X) # X_t = mu_t + eps_t (simulated process...spec(varModel, mean.model = meanModel)  ugarchfit(uspec, data = x)) 检查(标准化的)Z,即残差Z的伪观测值。...从拟合的copula 模型进行模拟。

    37200

    Python数据挖掘算法入门与实践

    一、数据挖掘简介 数据挖掘是一个通过对大量数据进行清理和处理,以发现其中隐藏的信息和模式的过程。简单来说,它是从大量数据中提取或“挖掘”知识的过程,也称为知识发现。...这个阶段的目标是根据所选择的算法,从预处理好的数据中提取出有意义的特征,并将其转换为适合特定数据挖掘算法的分析模型。 然后是数据挖掘阶段。...2.4 聚类算法 聚类分析处理的对象集合中,对象的类是未知的。它的目标是将对象集合分组为多个由类似对象组成的簇。...z-score 标准化:将数值范围缩放到0附近, 经过处理的数据符合标准正态分布。u是平均值,σ是标准差。 修正的标准z-score:修正后可以减少样本数据异常值的影响。...将z-score标准化公式中的均值改为中位数,将标准差改为绝对偏差。

    1K11

    ARIMA模型、随机游走模型RW模拟和预测时间序列趋势可视化

    = 0 模拟随机游走序列 我们现在可以通过为arima.sim 提供适当的参数来模拟 R 中的随机游走序列, 如下所示: R sim 我们可以使用该plot.ts() 函数绘制新生成的序列...为此,我们为函数提供了一个额外的参数均值/截距 。这个截距是模型的斜率。我们还可以更改模拟序列的标准差。在下面的代码中,我们提供了 1 的平均值和 5 的标准差。...> Rt sim > plot.ts 估计随机游走模型 为了拟合具有时间序列偏移的随机游走模型,我们将遵循以下步骤 取数据的一阶差分。...arima使用阶数为 的函数 将白噪声模型拟合到差分数据 c(0,0,0)。 绘制原始时间序列图。 abline通过提供通过将白噪声模型拟合为斜率得到的截距,使用该函数添加估计趋势 。 1....在我们的例子中,我们将指定白噪声模型的“a=0”和“b=intercept”。 > abline 估计的趋势线将添加到我们的图中。

    2.2K30

    机器学习 | 特征工程(数据预处理、特征抽取)

    标准化原理 其中mean为平均值,σ为标准差。 为什么要用标准化 因为标准差可以解决归一化存在的问题。标准差对异常点不敏感。 上图中红点为平均值,粉点为两个异常值。...当异常点出现时,总体数据的平均值和标准差并不会有特别大的波动。这也是在机器学习中标准化应用广泛的主要原因。...0,即平均值为0,并且标准差为1。...只需要将函数中的语句dict = DictVectorizer()改成dict = DictVectorizer(sparse=False) 结果如下: 那么问题又来了,转化为的数组是什么意思呢?...之前字典特征提取sklearnAPI中可以直接将sparse调成False,但是文本提取特征API没有该功能,也算是函数的一个bug吧。所以我们只能将最后的data转化为array形式。

    2.2K21
    领券