相关代码:
- library("fOptions")
- ##计算正股历史最近一年波动率
- wzg<- window(zg, start=(as.Date(Sys.timeDate())-365),end =as.Date(Sys.timeDate()))
- 年波幅 <- colSds(returns(wzg))*dim(wzg)[1]^0.5
- print(c("股票年波幅", 年波幅))
- ####两叉树(网)和B-S定价计算
- 到期时间 <- as.numeric(到期日期-as.Date(Sys.timeDate()))/365
- 股票波幅 <- 年波幅
- 股价<- zg[dim(zg)[1]]
- 权证现价<- qz[dim(qz)[1]]
- ## Plot CRR Option Tree:(两叉树法)
- CRRTree = BinomialTreeOption(TypeFlag = "ca", S = 股价, X = 行权价,
- Time = 到期时间, r = 无风险利率, b = 年化持有成本利率, sigma = 股票波幅, n = 12)
- BinomialTreePlot(CRRTree, dy = 1, cex = 0.8, ylim = c(-6, 7),
- xlab = "n", ylab = "权证理论价值")
- title(main = c(paste(colnames(qz),"权证树"),paste("(行权比例",行权比例,"/1)")))
- ## The Generalized Black Scholes Option Formula(B-S公式)
- bs<- GBSCharacteristics(TypeFlag = "c", S = 股价, X = 行权价,
- Time = 到期时间, r = 无风险利率, b = 年化持有成本利率, sigma = 股票波幅)
- 隐含波动率<- GBSVolatility(price=权证现价*行权比例 , TypeFlag = "c", S = 股价, X = 行权价,
- Time = 到期时间, r = 无风险利率, b = 年化持有成本利率)
- bs
- print(c(colnames(qz),"权证理论价值=",bs$premium/行权比例,"元"))



雷达卡





京公网安备 11010802022788号







