Polynomial Chaos Expansion Approach to Interest Rate Models

<div>Absolute errors of the two quantiles <svg height="15.4737pt" id="M293" 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="M294" 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 GBM at time <svg height="8.85534pt" id="M295" 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>, whose parameters are <svg height="9.20766pt" id="M296" style="vertical-align:-0.5195694pt" version="1.1" viewbox="-0.0498162 -8.68809 37.606 9.20766" width="37.606pt" 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="M397 380C406 395 404 411 396 425S369 451 350 451C302 451 239 372 192 294H189L199 338C214 407 207 451 180 451C152 451 83 405 30 345L48 318C87 354 117 372 125 372S135 362 127 324L55 -5L61 -12C87 -5 117 3 139 5C154 87 168 162 179 207C198 250 240 310 260 332C281 355 297 366 307 366C321 366 333 360 347 346C351 342 360 342 369 348C379 355 389 366 397 380Z" id="g50-115"></path><glyph.data ascent="3443" descent="-2856" horiz-adv-x="427" vert-adv-y="427"></glyph.data></g><g transform="matrix(.013,0,0,-0.013,9.229,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,20.847,0)"><path d="M290 377C321 398 342 415 358 430C378 450 389 473 389 502C389 578 329 635 238 635H237C184 635 137 610 109 578L64 515L88 493C112 529 154 573 208 573S303 542 303 482C303 409 233 370 141 341L149 308C165 313 190 319 215 319C272 319 341 283 341 193C342 98 292 43 222 43C163 43 122 72 96 94C88 101 79 100 70 94C61 87 47 73 46 60C44 47 48 37 62 23C76 10 118 -12 165 -12C238 -12 430 62 430 223C430 297 379 359 290 375V377Z" id="g50-52"></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,27.178,0)"><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>, <svg height="9.20766pt" id="M297" style="vertical-align:-0.5195694pt" version="1.1" viewbox="-0.0498162 -8.68809 45.905 9.20766" width="45.905pt" 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="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,29.115,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><g transform="matrix(.013,0,0,-0.013,35.493,0)"><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 starting value <svg height="11.9661pt" id="M298" style="vertical-align:-3.16057pt" version="1.1" viewbox="-0.0498162 -8.80553 45.6058 11.9661" width="45.6058pt" 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="M463 633C453 636 427 644 417 648C388 660 349 667 316 667C183 667 100 592 100 483C100 398 176 335 223 305L259 282C321 242 358 203 358 149C358 66 301 22 234 22C110 22 70 131 57 208L24 205L28 48C51 23 128 -17 215 -17C342 -17 446 60 446 174C446 262 394 310 317 361L284 383C229 420 183 455 183 515C183 569 222 628 298 628C384 628 418 563 428 481L461 487C459 527 458 593 463 633Z" id="g50-84"></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,6.149,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,14.771,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,26.389,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,32.72,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,39.051,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 8