R : Black-Scholes 및 이분법을 사용하여 IV 계산, 루프 작동을 거부 함

Question

CSV의 데이터가 포함 된 통화 옵션에 대한 Black-Scholes 기능과 이등분 모델이 있습니다. 그것은 관용을 초과하기 때문에 안쪽 루프에 갇혀있는 것처럼 보입니다. 내 Black-Scholes는 정확하게 계산하지 않으며 입찰가 평균을 사용하고 옵션의 실제 가격 대신 시장 가격을 요구합니다. 몇 시간 동안이 작업을하고 나면 분명히 뭔가 빠져있는 것일 수도 있습니다.

CSV에 대한 링크는 여기에 있습니다 : 당신은 true의 경우, 항상 사실이다,하는 조건에 while 루프를 사용하고

while(abs(bsmCall(as.numeric(as.character(myData[i,17])), 
        as.numeric(as.character(myData[i,1])), 
        as.numeric(as.character(myData[i,22])), 
        volMid, 
        as.numeric(as.character(myData[i,23])))-(as.numeric(as.character(myData[i,5])))) >= tol) 

: http://s000.tinyupload.com/?file_id=06213890949979926112는

######################################################################## 
#Black-Scholes-Merton Call 
bsmCall <- function(S, K, M, sig, r) { 
    yrTime=(M/252) 
    d1 <- (log(S/K)+(r+(sig^2/2))*(yrTime))/(sig*(sqrt(yrTime))) 
    d2 <- d1-sig*(sqrt(yrTime)) 
    C <- (S*(pnorm(d1)))-((pnorm(d2))*K*(exp(-r*yrTime))) 
    return(C) 
} 
######################################################################## 

myData = read.csv("09-26-16.csv", stringsAsFactors=FALSE) #DATA 
myData <- myData[,2:24] #omit first column 

####### start bisection method of CALLS and put IV in database ####### 
i <- 1 # reset counter 
tol <- 0.000001 #tolerance 

while(i <= nrow(myData)) { 
    if((myData[i,5] != 0) & (myData[i,6] != 0)) { 
    volLower <- .0001 #will need to reset with each iteration 
    volUpper <- 1   #will need to reset with each iteration 
    volMid <- (volLower + volUpper)/2 #will need to reset with each iteration 

    while(abs(bsmCall(as.numeric(as.character(myData[i,17])),as.numeric(as.character(myData[i,1])),as.numeric(as.character(myData[i,22])),volMid,as.numeric(as.character(myData[i,23])))-(as.numeric(as.character(myData[i,5])))) >= tol) { 
     if((bsmCall(as.numeric(as.character(myData[i,17])),as.numeric(as.character(myData[i,1])),as.numeric(as.character(myData[i,22])),volMid,as.numeric(as.character(myData[i,23])))-(as.numeric(as.character(myData[i,5])))) < 0) { 
     volLower <- volMid 
     volMid <- (volUpper + volMid)/2 
     } else { 
     volUpper <- volMid 
     volMid <- (volLower + volMid)/2 
     } 
    } 
    myData[i,8] <- volMid 
    } else { myData[i,8] <- 0 } 
    i=i+1 
} 

Answer 1

문제는 여기에있다 . 그것은 무한 루프입니다. 데이터의 첫 번째 행에서이 문제가 발생합니다.

이 오류를 수정하는 방법은 유스 케이스에 따라 다르지만, while을 if으로 변경하면 즉시 루프가 완료됩니다.

이분법에 대해 질문했습니다. 이 패키지에서 몇은 여기 here에서 또 다른입니다 :

bisect <- function(fn, lower, upper, tol=1.e-07, ...) { 
f.lo <- fn(lower, ...) 
f.hi <- fn(upper, ...) 
feval <- 2 

if (f.lo * f.hi > 0) stop("Root is not bracketed in the specified interval 
\n") 
chg <- upper - lower 

while (abs(chg) > tol) { 
     x.new <- (lower + upper)/2 
     f.new <- fn(x.new, ...) 
     if (abs(f.new) <= tol) break 
     if (f.lo * f.new < 0) upper <- x.new 
     if (f.hi * f.new < 0) lower <- x.new 
     chg <- upper - lower 
     feval <- feval + 1 
} 
list(x = x.new, value = f.new, fevals=feval) 
} 

# An example 
fn1 <- function(x, a) { 
exp(-x) - a*x 
} 

bisect(fn1, 0, 2, a=1) 

bisect(fn1, 0, 2, a=2) 

재귀 버전 :

이 지금까지 내 3 편집 ...

가 while 루프를 재구성 할 수 있습니다은

bisectMatt <- function(fn, lo, hi, tol = 1e-7, ...) { 

    flo <- fn(lo, ...) 
    fhi <- fn(hi, ...) 

    if(flo * fhi > 0) 
     stop("root is not bracketed by lo and hi") 

    mid <- (lo + hi)/2 
    fmid <- fn(mid, ...) 
    if(abs(fmid) <= tol || abs(hi-lo) <= tol) 
     return(mid) 


    if(fmid * fhi > 0) 
     return(bisectMatt(fn, lo, mid, tol, ...)) 

    return(bisectMatt(fn, mid, hi, tol, ...)) 
} 

Answer 2

이름 : Jeez입니다 i=1 일 때 각 반복 이후에 while 상태의 유일한 부분 인 volMid을 인쇄하십시오.

i <- 1 
volLower <- .0001 #will need to reset with each iteration 
volUpper <- 1   #will need to reset with each iteration 
volMid <- (volLower + volUpper)/2 #will need to reset with each iteration 

j <- 1 
while(abs(bsmCall(myData[i,17], myData[i,1], myData[i,22],volMid,myData[i,23])-myData[i,5]) >= tol & j < 30) { 
    if(bsmCall(myData[i,17], myData[i,1], myData[i,22],volMid,myData[i,23])-myData[i,5] < 0) { 
volLower <- volMid 
volMid <- (volUpper + volMid)/2 
    } else { 
    print("pos") 
    volUpper <- volMid 
    volMid <- (volLower + volMid)/2 
    } 
    j <- j + 1 
    print(volMid) 
} 

결과 :

#[1] 0.750025 
#[1] 0.8750125 
#[1] 0.9375062 
#[1] 0.9687531 
#[1] 0.9843766 
#[1] 0.9921883 
#[1] 0.9960941 
#[1] 0.9980471 
#[1] 0.9990235 
#[1] 0.9995118 
#[1] 0.9997559 
#[1] 0.9998779 
#[1] 0.999939 
#[1] 0.9999695 
#[1] 0.9999847 
#[1] 0.9999924 
#[1] 0.9999962 
#[1] 0.9999981 
#[1] 0.999999 
#[1] 0.9999995 
#[1] 0.9999998 
#[1] 0.9999999 
#[1] 0.9999999 
#[1] 1 
#[1] 1 
#[1] 1 
#[1] 1 
#[1] 1 
#[1] 1 

volMid 미만 30 반복 후 1으로 수렴하고, 거기에서 밖으로에,이 붙어 있습니다.