From: bstudent <marc.ruetten_at_gmx.de>

Date: Sat, 23 Apr 2011 11:57:49 -0700 (PDT)

*+
*

+ BS <- function(x) {

*+
*

*+ vE <- citin[i,2]
*

*+ sigE <- citin[i,3]
*

*+ D <- citin[i,4]
*

*+ R <- citin[i,5]
*

*+
*

*+ T = 1
*

*+
*

*+ f <- rep(NA, length(x))
*

*+
*

*+ f[1] <- (x[1] * pnorm(log (x[1]/D) + (R + ( (x[1]^2) / 2) ) * T ) / ( x[2]
*

* sqrt(T)) - exp(-R*T) * D * pnorm(log (x[1]/D) + (R - ( (x[1]^2) / 2) ) * T ) / ( x[2] * sqrt(T))) - vE

*+
*

+ f[2] <- ((x[1] * exp(-T) * pnorm(log (x[1]/D) + (R + ( (x[1]^2) / 2) ) * T

) / ( x[2] * sqrt(T)) * x[2]) / vE) - sigE

*+
*

*+ f
*

*+ }
*

*+
*

*+ p0 <- c( vE + D, sigE * (vE / (vE + D)) )
*

*+
*

*+ ans <- dfsane(par=p0, fn=BS)
*

*+
*

+ print(as.matrix(ans$par))

*+
*

*+ }
*

Date: Sat, 23 Apr 2011 11:57:49 -0700 (PDT)

Hello,

for my diploma thesis I need to program a solver for Merton´s respectively Black´s and Scholes´ Option pricing formula, which should be achieved for several dates.

What I want to do is to estimate the value of a firm´s assets "vA" (x[2] denotes vA) and the option-implied volatility of firm´s assets "sigA" (x[1] denotes sigA) by solving it simultaneous using the Black and Scholes formulas. This solution should be computed for several dates, at which the data input changes with the next date (it should be computed row by row of the data input matrix).

What I did so far is to program a for-loop in which I wrote the solver:

># Dataset: > > head(citin) Q vE sigE D R 1 1 39.81095 0.08957312 51.64004 0.00930000 2 2 39.76028 0.09127646 54.98504 0.01072000 3 3 37.00382 0.08177820 60.29025 0.01489545 4 4 35.17477 0.09221447 60.41061 0.020178795 5 36.54418 0.07852334 55.00648 0.02553438 6 6 37.20026 0.07949604 62.59768 0.02908462

> > > ># Loop with Solver for 6 dates: > > citinbs <- for (i in 1:6) {

+ BS <- function(x) {

* sqrt(T)) - exp(-R*T) * D * pnorm(log (x[1]/D) + (R - ( (x[1]^2) / 2) ) * T ) / ( x[2] * sqrt(T))) - vE

+ f[2] <- ((x[1] * exp(-T) * pnorm(log (x[1]/D) + (R + ( (x[1]^2) / 2) ) * T

) / ( x[2] * sqrt(T)) * x[2]) / vE) - sigE

+ print(as.matrix(ans$par))

The problem is, that the results I get aren´t really plausible. The next thing is, that I need the relevant Values of Output (ans$par, which includes two values - x[1]=vA and x[2]=sigA - per date) as matrix. This Matrix should look like this:

vA sigA

a x b y c z . . . . . .

The entry I wrote ("print(as.matrix(ans$par))") doesn´t achieve that.

I hope you understood what I´m trying to do. I´m an absolute beginner in programming in R, so these are some of my first steps. Please be patient ;) Also I hope my English isn´t to bad.

Thank you very much for helping me out!!!

Kind regards,

bstudent.

-- View this message in context: http://r.789695.n4.nabble.com/Loop-and-Solver-with-Black-Scholes-Formula-tp3470488p3470488.html Sent from the R help mailing list archive at Nabble.com. ______________________________________________ R-help_at_r-project.org mailing list https://stat.ethz.ch/mailman/listinfo/r-help PLEASE do read the posting guide http://www.R-project.org/posting-guide.html and provide commented, minimal, self-contained, reproducible code.Received on Sat 23 Apr 2011 - 20:50:24 GMT

Archive maintained by Robert King, hosted by
the discipline of
statistics at the
University of Newcastle,
Australia.

Archive generated by hypermail 2.2.0, at Sat 23 Apr 2011 - 21:10:33 GMT.

*
Mailing list information is available at https://stat.ethz.ch/mailman/listinfo/r-help.
Please read the posting
guide before posting to the list.
*