蒙特卡羅

使用蒙特卡羅和最小二乘回歸估計條件期望

  • November 24, 2021

我希望了解用於百慕大期權估值的最小二乘蒙地卡羅問題,但要從更簡單的背景來看。

說我有隨機變數 $ X $ 和 $ Y $ 是統一的

$$ 0,1 $$並且獨立。定義 $ Z=X^2+Y^2+XY $ . 假設我想評估期望 $ E(X|Z=a) $ 使用蒙地卡羅。在這種情況下,最小二乘 MC 可以提供幫助嗎?如果是這樣,任何人都可以概述該過程嗎? 我試圖理解算法的核心,而無需在閱讀論文和其他解釋 MC 最小二乘法的文章時必須經過繁瑣的符號。

您評論中的連結提到了 Kenneth Judd 在經濟學中的數值方法第 11.6 節。我建議也閱讀一下。它只有幾頁。下面的一些程式碼為您給出的問題實現了最小二乘蒙地卡羅:

set.seed(42)
fun <- function(x, y) x^2 + y^2 + x * y
N <- 1e3L
X <- runif(N)
Y <- runif(N)
Z <- fun(X, Y)
plot(Z, X)

# We can create a function of Z that gives an estimate of X:
model <- lm(X ~ Z + sqrt(Z))
print(model)

# Call:
# lm(formula = X ~ Z + sqrt(Z))
#
# Coefficients:
# (Intercept)            Z      sqrt(Z)  
#     0.02860      0.07183      0.44900  

a <- seq(0, 3, by = 0.01)
x_hat <- 
 model$coefficients[[1L]] + 
a * model$coefficients[[2L]] + sqrt(a) * model$coefficients[[3L]]
lines(a, x_hat, col = 'blue')

一個反對 x_hat 顯示合理的適合

引用自:https://quant.stackexchange.com/questions/68881