4

我正在写我的硕士论文，我在我的 R 代码中遇到了这个问题。在数学上，我想用 R 包“nleqslv”求解这个非线性方程组：

fnewton <- function(x){

y <- numeric(2)

d1 = (log(x[1]/D1)+(R+x[2]^2/2)*T)/x[2]*sqrt(T)

d2 = d1-x[2]*sqrt(T)

y1 <- SO1 - (x[1]*pnorm(d1) - exp(-R*T)*D1*pnorm(d2))

y2 <- sigmaS*SO1 - pnorm(d1)*x[2]*x[1]

y}

xstart <- c(21623379, 0.526177094846878)

nleqslv(xstart, fnewton, control=list(btol=.01), method="Newton")

我已经尝试了这个代码的几个版本，现在我得到了错误：

error: error in pnorm(q, mean, sd, lower.tail, log.p): not numerical.

Pnorm 分别是 d1 和 d2 的累积标准正态分布。我真的不知道，当我在 Teterevas 幻灯片上定位我的模型时我做错了什么（幻灯片第 5 号是她的模型代码），谁的演示文稿是谷歌搜索的第一个结果

https://www.google.de/search?q=moodys+KMV+in+R&rlz=1C1SVED_enDE401DE401&aq=f&oq=moodys+KMV+in+R&aqs=chrome.0.57.1​​3309j0&sourceid=chrome&ie=UTF-8#q=distance+到+默认+in+R

像我一样，无论多么成功，她都通过 Black-Scholes-Merton 方法计算违约距离风险度量。在这个模型中，股权的价值（通常用市值表示，->SO1）可以写成欧式看涨期权——我在上面的代码中标记为 y2，然而，之前的等式设置为 0！

其他变量是：

x[1] -> 我要导出的变量，总资产的价值

x[2] -> 我要导出的变量，总资产的波动率

D1 -> 债务账面价值（1998-2009）

R -> 无风险利率

T -> 设置为 1（时间）

sigmaS -> 估计（历史）股票波动率

已经谢谢了！！！我会很高兴，任何人都可以帮助我。卡罗

4

1 回答 1

11

我是 nleqslv 的作者，我对您使用它的方式感到非常惊讶。正如其他人所提到的，您没有返回任何明智的东西。

y1 应该是 y[1]，y2 应该是 y[2]。如果您想让我们说一些明智的话，您必须提供 D1、R、T、sigmaS 和 SO1 的数值。我试过这个：

T <- 1; D1 <- 1000; R <- 0.01; sigmaS <- .1; SO1 <- 1000

这些已在函数定义之前输入。看到这个

library(nleqslv)

T <- 1
D1 <- 1000
R <- 0.01

sigmaS <- .1
SO1 <- 1000

fnewton <- function(x){
    y <- numeric(2)
    d1 <- (log(x[1]/D1)+(R+x[2]^2/2)*T)/x[2]*sqrt(T)
    d2 <- d1-x[2]*sqrt(T)
    y[1] <- SO1 - (x[1]*pnorm(d1) - exp(-R*T)*D1*pnorm(d2))
    y[2] <- sigmaS*SO1 - pnorm(d1)*x[2]*x[1]
    y
}

xstart <- c(21623379, 0.526177094846878)

nleqslv在这种情况下找到解决方案没有问题。找到的解决方案是：c(1990.04983,0.05025)。似乎不需要设置btol参数，您可以使用 method Broyden。

于 2013-09-23T17:49:15.830 回答