下面程序不知哪有问题,哪位大侠帮忙看看?total.callvix是数据。报错信息
Error in s * (pnorm((log(s/K) + (x + 0.5 * sigma^2) * t)/(sigma * sqrt(t)))) : 's' is missing
for (i in 1:9946){
fun<-function (x,s,t,sigma,K,C) {s*(pnorm((log(s/K) + (x +0.5 * sigma^2) * t)/(sigma * sqrt(t)))) -K*exp(-x * t) * (pnorm((log(s/K) + (x - 0.5 * sigma^2) *t)/(sigma * sqrt(t))))-C}
total.callvix[i,"r"]<-uniroot(fun(x),c(-10,10),s=total.callvix[i,"Close"],t=as.numeric(total.callvix[i,"t"])/365,sigma=total.callvix[i,"Voal"],K=0.001*total.callvix[i,"K"],C=(total.callvix[i,"High"]+total.callvix[i,"Low"])/2)$root
}