Polynomial Chaos Expansion Approach to Interest Rate Models

<div>Absolute errors of the two quantiles <svg height="15.4737pt" id="M406" style="vertical-align:-3.3799pt" version="1.1" viewbox="-0.0498162 -12.0938 26.108 15.4737" width="26.108pt" xmlns="http://www.w3.org/2000/svg" xmlns:xlink="http://www.w3.org/1999/xlink"><g transform="matrix(.013,0,0,-0.013,.861,-3.027)"><path d="M686 564C571 601 458 635 356 679H331C229 635 115 600 0 564L11 533L343 603L675 533L686 564Z" id="g56-105"></path><glyph.data ascent="3443" descent="-2856" horiz-adv-x="686" vert-adv-y="686"></glyph.data></g><g transform="matrix(.013,0,0,-0.013,0,0)"><path d="M750 355C750 552 603 667 428 667C199 667 24 497 24 286C24 118 146 -4 316 -16L336 -29C480 -122 544 -145 582 -154C606 -159 666 -170 714 -172L722 -146C652 -117 562 -75 472 -15L454 -3C623 38 750 176 750 355ZM646 353C646 210 556 50 414 25L385 41L304 27C193 50 128 152 128 282C128 474 254 626 418 626C571 626 646 519 646 353Z" id="g50-82"></path><glyph.data ascent="3443" descent="-2856" horiz-adv-x="774" vert-adv-y="774"></glyph.data></g><g transform="matrix(.0091,0,0,-0.0091,9.399,3.001)"><path d="M247 635C116 635 39 518 39 421C39 315 113 237 221 237C239 237 258 241 281 253L349 289C315 140 204 42 60 19L66 -15C90 -15 153 -5 207 16C344 70 447 199 447 385C447 521 373 635 247 635ZM231 598C329 598 356 478 356 390C356 371 355 347 352 325C330 309 297 297 262 297C179 297 127 369 127 456C127 515 156 598 231 598Z" id="g50-58"></path><glyph.data ascent="3443" descent="-2856" horiz-adv-x="487" vert-adv-y="487"></glyph.data></g><g transform="matrix(.0091,0,0,-0.0091,13.832,3.001)"><path d="M247 635C116 635 39 518 39 421C39 315 113 237 221 237C239 237 258 241 281 253L349 289C315 140 204 42 60 19L66 -15C90 -15 153 -5 207 16C344 70 447 199 447 385C447 521 373 635 247 635ZM231 598C329 598 356 478 356 390C356 371 355 347 352 325C330 309 297 297 262 297C179 297 127 369 127 456C127 515 156 598 231 598Z" id="g50-58"></path><glyph.data ascent="3443" descent="-2856" horiz-adv-x="487" vert-adv-y="487"></glyph.data></g><g transform="matrix(.0091,0,0,-0.0091,18.296,3.001)"><path d="M611 626L575 651L170 -9L204 -35L611 626ZM200 635C93 635 50 532 50 449C50 365 94 264 199 264C306 264 349 365 349 450C349 534 304 635 200 635ZM199 601C264 601 273 501 273 450S264 298 200 298C139 298 126 398 126 450S139 601 199 601ZM590 359C485 359 440 256 440 173C440 89 485 -12 589 -12C694 -12 738 89 738 173C738 256 694 359 590 359ZM589 325C653 325 662 225 662 173C662 124 653 22 590 22C529 22 516 122 516 174C516 227 529 325 589 325Z" id="g58-35"></path><glyph.data ascent="1007" descent="-369" horiz-adv-x="786" vert-adv-y="786"></glyph.data></g></svg> and <svg height="15.4737pt" id="M407" style="vertical-align:-3.3799pt" version="1.1" viewbox="-0.0498162 -12.0938 32.721 15.4737" width="32.721pt" xmlns="http://www.w3.org/2000/svg" xmlns:xlink="http://www.w3.org/1999/xlink"><g transform="matrix(.013,0,0,-0.013,.861,-3.027)"><path d="M686 564C571 601 458 635 356 679H331C229 635 115 600 0 564L11 533L343 603L675 533L686 564Z" id="g56-105"></path><glyph.data ascent="3443" descent="-2856" horiz-adv-x="686" vert-adv-y="686"></glyph.data></g><g transform="matrix(.013,0,0,-0.013,0,0)"><path d="M750 355C750 552 603 667 428 667C199 667 24 497 24 286C24 118 146 -4 316 -16L336 -29C480 -122 544 -145 582 -154C606 -159 666 -170 714 -172L722 -146C652 -117 562 -75 472 -15L454 -3C623 38 750 176 750 355ZM646 353C646 210 556 50 414 25L385 41L304 27C193 50 128 152 128 282C128 474 254 626 418 626C571 626 646 519 646 353Z" id="g50-82"></path><glyph.data ascent="3443" descent="-2856" horiz-adv-x="774" vert-adv-y="774"></glyph.data></g><g transform="matrix(.0091,0,0,-0.0091,9.399,3.001)"><path d="M247 635C116 635 39 518 39 421C39 315 113 237 221 237C239 237 258 241 281 253L349 289C315 140 204 42 60 19L66 -15C90 -15 153 -5 207 16C344 70 447 199 447 385C447 521 373 635 247 635ZM231 598C329 598 356 478 356 390C356 371 355 347 352 325C330 309 297 297 262 297C179 297 127 369 127 456C127 515 156 598 231 598Z" id="g50-58"></path><glyph.data ascent="3443" descent="-2856" horiz-adv-x="487" vert-adv-y="487"></glyph.data></g><g transform="matrix(.0091,0,0,-0.0091,13.832,3.001)"><path d="M247 635C116 635 39 518 39 421C39 315 113 237 221 237C239 237 258 241 281 253L349 289C315 140 204 42 60 19L66 -15C90 -15 153 -5 207 16C344 70 447 199 447 385C447 521 373 635 247 635ZM231 598C329 598 356 478 356 390C356 371 355 347 352 325C330 309 297 297 262 297C179 297 127 369 127 456C127 515 156 598 231 598Z" id="g50-58"></path><glyph.data ascent="3443" descent="-2856" horiz-adv-x="487" vert-adv-y="487"></glyph.data></g><g transform="matrix(.0091,0,0,-0.0091,18.264,3.001)"><path d="M117 -12C153 -12 178 12 178 50C178 84 153 110 118 110C84 110 59 84 59 50C59 12 84 -12 117 -12Z" id="g50-47"></path><glyph.data ascent="3443" descent="-2856" horiz-adv-x="237" vert-adv-y="237"></glyph.data></g><g transform="matrix(.0091,0,0,-0.0091,20.421,3.001)"><path d="M247 635C116 635 39 518 39 421C39 315 113 237 221 237C239 237 258 241 281 253L349 289C315 140 204 42 60 19L66 -15C90 -15 153 -5 207 16C344 70 447 199 447 385C447 521 373 635 247 635ZM231 598C329 598 356 478 356 390C356 371 355 347 352 325C330 309 297 297 262 297C179 297 127 369 127 456C127 515 156 598 231 598Z" id="g50-58"></path><glyph.data ascent="3443" descent="-2856" horiz-adv-x="487" vert-adv-y="487"></glyph.data></g><g transform="matrix(.0091,0,0,-0.0091,24.909,3.001)"><path d="M611 626L575 651L170 -9L204 -35L611 626ZM200 635C93 635 50 532 50 449C50 365 94 264 199 264C306 264 349 365 349 450C349 534 304 635 200 635ZM199 601C264 601 273 501 273 450S264 298 200 298C139 298 126 398 126 450S139 601 199 601ZM590 359C485 359 440 256 440 173C440 89 485 -12 589 -12C694 -12 738 89 738 173C738 256 694 359 590 359ZM589 325C653 325 662 225 662 173C662 124 653 22 590 22C529 22 516 122 516 174C516 227 529 325 589 325Z" id="g58-35"></path><glyph.data ascent="1007" descent="-369" horiz-adv-x="786" vert-adv-y="786"></glyph.data></g></svg> of the PCE approximation of Vasicek interest rate model at time <svg height="8.85534pt" id="M408" style="vertical-align:-0.04981041pt" version="1.1" viewbox="-0.0498162 -8.80553 30.1943 8.85534" width="30.1943pt" xmlns="http://www.w3.org/2000/svg" xmlns:xlink="http://www.w3.org/1999/xlink"><g transform="matrix(.013,0,0,-0.013,0,0)"><path d="M655 676H628C608 656 600 650 571 650H156C118 650 114 651 102 676H77C63 615 45 548 24 485L60 486C81 531 97 562 114 580C134 603 148 611 252 611H308L217 127C201 41 196 38 105 32L98 0H384L393 32C300 38 288 41 304 127L393 611H466C553 611 571 604 582 580C592 559 597 530 595 485L630 489C636 548 645 624 655 676Z" id="g50-85"></path><glyph.data ascent="3443" descent="-2856" horiz-adv-x="673" vert-adv-y="673"></glyph.data></g><g transform="matrix(.013,0,0,-0.013,12.064,0)"><path d="M556 331V384H55V331H556ZM556 140V193H55V140H556Z" id="g54-34"></path><glyph.data ascent="3443" descent="-2856" horiz-adv-x="611" vert-adv-y="611"></glyph.data></g><g transform="matrix(.013,0,0,-0.013,23.683,0)"><path d="M389 0V32C297 38 291 46 291 118V635C234 613 175 595 109 583V556L161 554C203 552 207 547 207 497V118C207 46 201 38 110 32V0H389Z" id="g50-50"></path><glyph.data ascent="3443" descent="-2856" horiz-adv-x="487" vert-adv-y="487"></glyph.data></g></svg>. The parameters are set as <svg height="8.54217pt" id="M409" style="vertical-align:-0.2063894pt" version="1.1" viewbox="-0.0498162 -8.33578 38.8195 8.54217" width="38.8195pt" xmlns="http://www.w3.org/2000/svg" xmlns:xlink="http://www.w3.org/1999/xlink"><g transform="matrix(.013,0,0,-0.013,0,0)"><path d="M558 106L533 129C501 84 476 67 464 67C446 67 435 116 414 239C457 297 508 367 550 440L539 451L484 438C458 390 429 338 404 300H402C375 407 351 451 285 451C168 451 24 312 24 153C24 56 66 -12 132 -12C206 -12 284 64 345 156H347C367 24 386 -12 422 -12C457 -12 503 7 558 106ZM338 206C272 99 215 56 173 56C140 56 116 96 116 170C116 303 194 405 255 405C304 405 324 299 338 206Z" id="g50-223"></path><glyph.data ascent="3443" descent="-2856" horiz-adv-x="582" vert-adv-y="582"></glyph.data></g><g transform="matrix(.013,0,0,-0.013,11.242,0)"><path d="M556 331V384H55V331H556ZM556 140V193H55V140H556Z" id="g54-34"></path><glyph.data ascent="3443" descent="-2856" horiz-adv-x="611" vert-adv-y="611"></glyph.data></g><g transform="matrix(.013,0,0,-0.013,22.861,0)"><path d="M245 635C92 635 37 457 37 312C37 149 91 -12 244 -12C395 -12 449 166 449 312C449 469 395 635 245 635ZM243 598C332 598 358 454 358 312C358 173 334 26 245 26C158 26 128 174 128 313S152 598 243 598Z" id="g50-49"></path><glyph.data ascent="3443" descent="-2856" horiz-adv-x="487" vert-adv-y="487"></glyph.data></g><g transform="matrix(.013,0,0,-0.013,29.192,0)"><path d="M117 -12C153 -12 178 12 178 50C178 84 153 110 118 110C84 110 59 84 59 50C59 12 84 -12 117 -12Z" id="g50-47"></path><glyph.data ascent="3443" descent="-2856" horiz-adv-x="237" vert-adv-y="237"></glyph.data></g><g transform="matrix(.013,0,0,-0.013,32.273,0)"><path d="M389 0V32C297 38 291 46 291 118V635C234 613 175 595 109 583V556L161 554C203 552 207 547 207 497V118C207 46 201 38 110 32V0H389Z" id="g50-50"></path><glyph.data ascent="3443" descent="-2856" horiz-adv-x="487" vert-adv-y="487"></glyph.data></g></svg>, <svg height="12.6917pt" id="M410" style="vertical-align:-3.40337pt" version="1.1" viewbox="-0.0498162 -9.28833 38.8848 12.6917" width="38.8848pt" xmlns="http://www.w3.org/2000/svg" xmlns:xlink="http://www.w3.org/1999/xlink"><g transform="matrix(.013,0,0,-0.013,0,0)"><path d="M563 583C563 664 501 710 433 710C381 710 334 692 285 648C212 582 178 506 149 350L74 -56C49 -190 36 -222 24 -226L30 -257C52 -256 95 -247 121 -229C131 -222 134 -206 165 -21L219 303C254 510 296 665 397 665C447 665 481 626 481 572C481 512 446 457 390 420C367 405 351 400 323 397L302 349C398 338 466 285 466 197C466 117 419 45 306 45C268 45 224 59 198 74L186 50C197 20 228 -14 274 -14C318 -14 359 3 401 26C478 68 554 144 554 225C554 322 493 367 412 396C490 446 563 500 563 583Z" id="g50-224"></path><glyph.data ascent="3443" descent="-2856" horiz-adv-x="587" vert-adv-y="587"></glyph.data></g><g transform="matrix(.013,0,0,-0.013,11.308,0)"><path d="M556 331V384H55V331H556ZM556 140V193H55V140H556Z" id="g54-34"></path><glyph.data ascent="3443" descent="-2856" horiz-adv-x="611" vert-adv-y="611"></glyph.data></g><g transform="matrix(.013,0,0,-0.013,22.926,0)"><path d="M245 635C92 635 37 457 37 312C37 149 91 -12 244 -12C395 -12 449 166 449 312C449 469 395 635 245 635ZM243 598C332 598 358 454 358 312C358 173 334 26 245 26C158 26 128 174 128 313S152 598 243 598Z" id="g50-49"></path><glyph.data ascent="3443" descent="-2856" horiz-adv-x="487" vert-adv-y="487"></glyph.data></g><g transform="matrix(.013,0,0,-0.013,29.257,0)"><path d="M117 -12C153 -12 178 12 178 50C178 84 153 110 118 110C84 110 59 84 59 50C59 12 84 -12 117 -12Z" id="g50-47"></path><glyph.data ascent="3443" descent="-2856" horiz-adv-x="237" vert-adv-y="237"></glyph.data></g><g transform="matrix(.013,0,0,-0.013,32.338,0)"><path d="M414 144C384 79 371 75 317 75H135L276 221C367 316 408 376 408 465C408 570 327 635 237 635C179 635 131 609 100 575L42 494L67 471C94 510 138 565 205 565C277 565 321 517 321 435C321 348 258 270 195 195C146 137 88 81 33 26V0H411C423 44 433 88 446 135L414 144Z" id="g50-51"></path><glyph.data ascent="3443" descent="-2856" horiz-adv-x="487" vert-adv-y="487"></glyph.data></g></svg>, <svg height="8.54217pt" id="M411" style="vertical-align:-0.2063894pt" version="1.1" viewbox="-0.0498162 -8.33578 45.096 8.54217" width="45.096pt" xmlns="http://www.w3.org/2000/svg" xmlns:xlink="http://www.w3.org/1999/xlink"><g transform="matrix(.013,0,0,-0.013,0,0)"><path d="M558 454L532 470C512 451 497 442 465 442C421 442 376 450 318 450C146 450 24 313 24 164C24 46 91 -12 183 -12C317 -12 439 105 439 247C439 306 428 346 382 385L383 389C414 384 450 379 485 379C517 379 539 412 558 454ZM357 273C357 191 313 28 205 28C149 28 111 79 111 160C111 266 175 396 280 396C298 396 328 386 341 357C353 331 357 315 357 273Z" id="g50-240"></path><glyph.data ascent="3443" descent="-2856" horiz-adv-x="576" vert-adv-y="576"></glyph.data></g><g transform="matrix(.013,0,0,-0.013,11.166,0)"><path d="M556 331V384H55V331H556ZM556 140V193H55V140H556Z" id="g54-34"></path><glyph.data ascent="3443" descent="-2856" horiz-adv-x="611" vert-adv-y="611"></glyph.data></g><g transform="matrix(.013,0,0,-0.013,22.784,0)"><path d="M245 635C92 635 37 457 37 312C37 149 91 -12 244 -12C395 -12 449 166 449 312C449 469 395 635 245 635ZM243 598C332 598 358 454 358 312C358 173 334 26 245 26C158 26 128 174 128 313S152 598 243 598Z" id="g50-49"></path><glyph.data ascent="3443" descent="-2856" horiz-adv-x="487" vert-adv-y="487"></glyph.data></g><g transform="matrix(.013,0,0,-0.013,29.115,0)"><path d="M117 -12C153 -12 178 12 178 50C178 84 153 110 118 110C84 110 59 84 59 50C59 12 84 -12 117 -12Z" id="g50-47"></path><glyph.data ascent="3443" descent="-2856" horiz-adv-x="237" vert-adv-y="237"></glyph.data></g><g transform="matrix(.013,0,0,-0.013,32.196,0)"><path d="M389 0V32C297 38 291 46 291 118V635C234 613 175 595 109 583V556L161 554C203 552 207 547 207 497V118C207 46 201 38 110 32V0H389Z" id="g50-50"></path><glyph.data ascent="3443" descent="-2856" horiz-adv-x="487" vert-adv-y="487"></glyph.data></g><g transform="matrix(.013,0,0,-0.013,38.527,0)"><path d="M158 548H390L417 615L410 623H122L83 318C105 326 143 337 185 337C296 337 350 275 350 188C350 116 308 42 225 42C164 42 122 74 100 93C90 101 82 99 72 92C60 82 51 68 50 59C48 46 52 38 66 24C82 9 125 -12 172 -12C225 -11 292 15 346 59C408 108 437 166 437 226C437 309 371 397 242 397C214 397 170 382 133 369L158 548Z" id="g50-54"></path><glyph.data ascent="3443" descent="-2856" horiz-adv-x="487" vert-adv-y="487"></glyph.data></g></svg>, and <svg height="11.7574pt" id="M412" style="vertical-align:-3.16066pt" version="1.1" viewbox="-0.0498162 -8.59674 47.6675 11.7574" width="47.6675pt" xmlns="http://www.w3.org/2000/svg" xmlns:xlink="http://www.w3.org/1999/xlink"><g transform="matrix(.013,0,0,-0.013,0,0)"><path d="M636 22C611 29 589 38 560 74C536 104 518 134 495 173L430 300C518 328 592 392 592 492C592 555 563 594 528 618C491 643 435 650 377 650H140L131 618C222 610 223 601 208 518L137 136C121 43 116 40 27 32L18 0H295L300 32C212 40 206 42 223 138L248 279H296C334 279 350 272 369 230C393 174 413 140 441 94C482 24 523 1 626 -8L636 22ZM497 480C497 422 468 374 417 343C387 324 353 315 304 315H255L298 556C303 583 309 598 319 603C331 610 348 613 367 613C448 613 497 573 497 480Z" id="g50-83"></path><glyph.data ascent="3443" descent="-2856" horiz-adv-x="654" vert-adv-y="654"></glyph.data></g><g transform="matrix(.0091,0,0,-0.0091,8.203,3.001)"><path d="M245 635C92 635 37 457 37 312C37 149 91 -12 244 -12C395 -12 449 166 449 312C449 469 395 635 245 635ZM243 598C332 598 358 454 358 312C358 173 334 26 245 26C158 26 128 174 128 313S152 598 243 598Z" id="g50-49"></path><glyph.data ascent="3443" descent="-2856" horiz-adv-x="487" vert-adv-y="487"></glyph.data></g><g transform="matrix(.013,0,0,-0.013,16.825,0)"><path d="M556 331V384H55V331H556ZM556 140V193H55V140H556Z" id="g54-34"></path><glyph.data ascent="3443" descent="-2856" horiz-adv-x="611" vert-adv-y="611"></glyph.data></g><g transform="matrix(.013,0,0,-0.013,28.443,0)"><path d="M389 0V32C297 38 291 46 291 118V635C234 613 175 595 109 583V556L161 554C203 552 207 547 207 497V118C207 46 201 38 110 32V0H389Z" id="g50-50"></path><glyph.data ascent="3443" descent="-2856" horiz-adv-x="487" vert-adv-y="487"></glyph.data></g><g transform="matrix(.013,0,0,-0.013,34.774,0)"><path d="M389 0V32C297 38 291 46 291 118V635C234 613 175 595 109 583V556L161 554C203 552 207 547 207 497V118C207 46 201 38 110 32V0H389Z" id="g50-50"></path><glyph.data ascent="3443" descent="-2856" horiz-adv-x="487" vert-adv-y="487"></glyph.data></g><g transform="matrix(.013,0,0,-0.013,41.105,0)"><path d="M245 635C92 635 37 457 37 312C37 149 91 -12 244 -12C395 -12 449 166 449 312C449 469 395 635 245 635ZM243 598C332 598 358 454 358 312C358 173 334 26 245 26C158 26 128 174 128 313S152 598 243 598Z" id="g50-49"></path><glyph.data ascent="3443" descent="-2856" horiz-adv-x="487" vert-adv-y="487"></glyph.data></g></svg>.</div>

Journal of Probability and Statistics

Polynomial Chaos Expansion Approach to Interest Rate Models

Figure 12