先贴一些写的代码library("RQuantLib")
read.csv("E:\\data\\sk.csv")

列表： data TOP.10.TRADED..BY.VOLUME. O.Q.P..CLOSE underlying strike
1 2013/9/2 C 132 2217522800
2 2013/9/3 C 123 2239523000
3 2013/9/4 C 127 2232623000
4 2013/9/5 C 137 2259823200
5 2013/9/6 C 141 2262125000
6 2013/9/9 C 84 2275123400
7 2013/9/10 C 111 2297723400
8 2013/9/11 C 111 2293723400
9 2013/9/12 C 102 2295423400
10 2013/9/13 C 138 2291523200
11 2013/9/16 C 115 2325223600
12 2013/9/17 C 81 2318123600
13 2013/9/18 C 26 2311724000
14 2013/9/19 C 96 2350323800
15 2013/9/25 C 25 2321023400
16 2013/9/26 C 41 2312523200
17 2013/9/27 C 0 2320723200
18 2013/9/30 C 680 2286023000
19 2013/10/2 C 89 2298423800
20 2013/10/3 C 200 2321423600
21 2013/10/4 C 48 2313924200
22 2013/10/7 C 120 2297423600
23 2013/10/8 C 74 2317924000
24 2013/10/9 C 94 2303423800
252013/10/10 C 162 2295123400
262013/10/11 C 89 2321823800
272013/10/15 C 164 2333723600
282013/10/16 C 130 2322823600
292013/10/17 C 48 2309523600
302013/10/18 C 165 2334023400
312013/10/21 C 49 2343823800
322013/10/22 P 77 2331623000
332013/10/23 C 21 2300023600
342013/10/24 C 34 2283623200
352013/10/25 C 30 2269823000
362013/10/28 C 27 2280723000
372013/10/29 C 27 2284723000
382013/10/30 C 0 2330423000
392013/10/31 C 383 2320623200
40 2013/11/1 C 81 2325024000
41 2013/11/4 C 98 2319023800
42 2013/11/5 C 70 2303923800
43 2013/11/6 C 104 2303723600
44 2013/11/7 C 21 2288124000
45 2013/11/8 C 78 2274423400
462013/11/11 C 152 2307023400
472013/11/12 C 108 2290123400
482013/11/13 C 41 2246423400
492013/11/14 C 85 2264923000
502013/11/15 C 108 2303223400
512013/11/18 C 154 2366024000
522013/11/19 C 90 2365824200
532013/11/20 C 75 2370124200
542013/11/21 C 30 2358024200
552013/11/22 C 45 2369624200
562013/11/25 C 31 2368424200
572013/11/26 C 30 2368124000
582013/11/27 C 69 2380623800
592013/11/28 C 0 2378924000
602013/11/29 C 475 2388123800
61 2013/12/2 C 412 2403924000
62 2013/12/3 C 337 2391024000
63 2013/12/4 C 41 2372925000
64 2013/12/5 C 174 2371324200
65 2013/12/6 C 271 2374324000
66 2013/12/9 C 33 2381125000
672013/12/10 C 571 2374424400
682013/12/11 C 77 2333824000
692013/12/12 C 74 2321823800
702013/12/13 C 131 2324623600
712013/12/16 C 22 2311524000
722013/12/17 C 101 2306923400
732013/12/18 C 142 2314423400
742013/12/19 C 33 2288923400
752013/12/20 C 84 2281223000
762013/12/23 C 92 2292223000
772013/12/24 C 47 2318023400
782013/12/27 C 31 2324323400
792013/12/31 C 110 2349024000
80 2014/1/2 C 96 2334024000
81 2014/1/3 C 81 2281723600
82 2014/1/6 C 90 2268423400
83 2014/1/7 C 165 2271323000
84 2014/1/8 C 131 2299723400
85 2014/1/9 C 90 2278723400
86 2014/1/10 C 152 2284623200
87 2014/1/13 P 80 2288922200
88 2014/1/14 C 57 2279123400
89 2014/1/15 C 126 2290223200
90 2014/1/16 C 131 2298624400
91 2014/1/17 C 55 2313323600
92 2014/1/20 C 38 2292923400
93 2014/1/21 C 24 2303323600
94 2014/1/22 C 96 2308223200
95 2014/1/23 C 49 2273423000
96 2014/1/24 C 33 2245022800
97 2014/1/27 C 8 2197622600
98 2014/1/28 C 72 2196122000
99 2014/1/29 C 0 2214222200
1002014/1/30 C 222 2203522400
risk.free.rate time.to.maturity
1 0.0087 0.1111
2 0.0087 0.1071
3 0.0087 0.1032
4 0.0087 0.0992
5 0.0087 0.0952
6 0.0087 0.0833
7 0.0087 0.0794
8 0.0087 0.0754
9 0.0087 0.0714
10 0.0087 0.0675
11 0.0087 0.0556
12 0.0087 0.0516
13 0.0087 0.0476
14 0.0087 0.0437
15 0.0087 0.0278
16 0.0087 0.0198
17 0.0087 0.0119
18 0.0087 0.3651
19 0.0087 0.1151
20 0.0087 0.1111
21 0.0087 0.1071
22 0.0087 0.0952
23 0.0087 0.0913
24 0.0087 0.0873
25 0.0087 0.0833
26 0.0087 0.0794
27 0.0087 0.0635
28 0.0087 0.0595
29 0.0087 0.0556
30 0.0087 0.0516
31 0.0087 0.0397
32 0.0087 0.0357
33 0.0087 0.0317
34 0.0087 0.0278
35 0.0087 0.0238
36 0.0087 0.0119
37 0.0087 0.0079
38 0.0087 0.0040
39 0.0087 0.1151
40 0.0686 0.1111
41 0.0686 0.0992
42 0.0686 0.0952
43 0.0686 0.0913
44 0.0686 0.0873
45 0.0686 0.0833
46 0.0686 0.0714
47 0.0686 0.0675
48 0.0686 0.0635
49 0.0686 0.0595
50 0.0686 0.0556
51 0.0686 0.0437
52 0.0686 0.0397
53 0.0686 0.0357
54 0.0686 0.0317
55 0.0686 0.0278
56 0.0686 0.0159
57 0.0686 0.0119
58 0.0686 0.0079
59 0.0686 0.0040
60 0.0686 0.1270
61 0.0750 0.1151
62 0.0750 0.1111
63 0.0750 0.1071
64 0.0750 0.1032
65 0.0750 0.0992
66 0.0750 0.0873
67 0.0750 0.4405
68 0.0750 0.0794
69 0.0750 0.0754
70 0.0750 0.0714
71 0.0750 0.0595
72 0.0750 0.0556
73 0.0750 0.0516
74 0.0750 0.0476
75 0.0750 0.0437
76 0.0750 0.0317
77 0.0750 0.0278
78 0.0750 0.0159
79 0.0750 0.1190
80 0.0619 0.1111
81 0.0619 0.1071
82 0.0619 0.0952
83 0.0619 0.0913
84 0.0619 0.0873
85 0.0619 0.0833
86 0.0619 0.0794
87 0.0619 0.0675
88 0.0619 0.0635
89 0.0619 0.0595
90 0.0619 0.2937
91 0.0619 0.0516
92 0.0619 0.0397
93 0.0619 0.0357
94 0.0619 0.0317
95 0.0619 0.0278
96 0.0619 0.0238
97 0.0619 0.0119
98 0.0619 0.0079
99 0.0619 0.0040
100 0.0006 0.1151

然后继续代码：

for(i in 1:length(d[,3]))
+ iv<-EuropeanOptionImpliedVolatility(type="call",value=d[i,3],underlying=d[i,4],strike=d[i,5],dividendYield=0,riskFreeRate=d[i,6],maturity=d[i,7],volatility=0.1)



root not bracketed: f[1e-007,4] -> [9.242558e+000,3.878750e+003]，共有16个警告 (用warnings()来显示)

然后输出结果的时候，结果只能显示第1个到第16个的，其他都为0，但是我如果将数据手动输入函数之中，例如第17个，就可以得到0.129473的，而不为0的结果。

小弟完全R语言初学+自学者，希望大神可以指导一下是怎么回事。或者往哪方面思考。