potrzebuję policzyć coś takiego:
\(\displaystyle{ E_{t-1}(logS_t - K)_{+}}\)
Gdzie wiem, że \(\displaystyle{ logS_t}\) ma warunkowy rozkład normalny:
\(\displaystyle{ E_{t-1} logS_t = \mu}\)
\(\displaystyle{ VAR_{t-1} logS_t = \sigma^2}\)
Próbuję do tego użyć
Kod: Zaznacz cały
https://en.wikipedia.org/wiki/Mills_ratio#Inverse_Mills_ratio
1. Czy dobrze to robię, tzn. czy inverse mills ratio rzeczywiście się do tego nadaje? Całość jest 'uproszczeniem' opcji call, polegającym na wykorzystaniu logarytmu z instrumentu, zamiast samego instrumentu.
2. (opcjonalnie) Czy ktoś kojarzy jakąś dobrą implementację Inverse Mills ratio w R? Teoretycznie powinno to być proste do napisania, w praktyce z uwagi na precyzję potrafi dość szybko zwracać jakieś bzdury, pomimo wykorzystania logarytmów przy podawaniu wartości gęstości / dystrybuanty rozkładu normalnego.
Dzięki za pomoc!