gpt4 book ai didi

r - 对数协方差到算术协方差矩阵函数?

转载 作者:行者123 更新时间:2023-12-02 14:10:30 26 4
gpt4 key购买 nike

是否有一个函数可以将使用对数返回构建的协方差矩阵转换为基于简单算术返回的协方差矩阵?

动机:我们希望使用均值-方差效用函数,其中预期 yield 和方差以算术术语指定。然而,由于对数返回的可加性,通常使用对数返回来估计返回和协方差,并且我们假设 Assets 价格遵循对数正态随机过程。

Meucci 描述了一个为Appendix page 5 上对数正态 yield 的通用/任意分布生成基于算术 yield 的协方差矩阵的过程。 .

最佳答案

这是我对公式的翻译:

linreturn <- function(mu,Sigma) {
m <- exp(mu+diag(Sigma)/2)-1
x1 <- outer(mu,mu,"+")
x2 <- outer(diag(Sigma),diag(Sigma),"+")/2
S <- exp(x1+x2)*(exp(Sigma)-1)
list(mean=m,vcov=S)
}

编辑:根据评论修复了-1问题。

尝试一个例子:

m1 <- c(1,2)
S1 <- matrix(c(1,0.2,0.2,1),nrow=2)

生成多元对数正态返回:

set.seed(1001)
r1 <- exp(MASS::mvrnorm(200000,mu=m1,Sigma=S1))-1
colMeans(r1)
## [1] 3.485976 11.214211
var(r1)
## [,1] [,2]
## [1,] 34.4021 12.4062
## [2,] 12.4062 263.7382

与公式的预期结果进行比较:

linreturn(m1,S1)

## $mean
## [1] 3.481689 11.182494

## $vcov
## [,1] [,2]
## [1,] 34.51261 12.08818
## [2,] 12.08818 255.01563

关于r - 对数协方差到算术协方差矩阵函数?,我们在Stack Overflow上找到一个类似的问题: https://stackoverflow.com/questions/7663690/

26 4 0
Copyright 2021 - 2024 cfsdn All Rights Reserved 蜀ICP备2022000587号
广告合作:1813099741@qq.com 6ren.com