22 августа 2023, 21:47

Тетта - поделитесь скриптом для Питона :)

Всем привет!

Если у кого есть, поделитесь плиз скриптом в Питоне для расчета тетты у коллов (по Б-Ш).

У меня почему-то результат неверный, в инете варианты кодов тоже дают ошибочные ответы…

Тетта питон опционы торговый софт

Врач-бондиатОр

Москва

786

с 4 ноября 2018

12 Комментариев

Иван-дурак
22 августа 2023, 22:24
import numpy as np
from scipy.stats import norm
def black_scholes_theta(S, K, T, r, sigma, option_type):
d1 = (np.log(S / K) + (r + 0.5 * sigma ** 2) * T) / (sigma * np.sqrt(T))
d2 = d1 — sigma * np.sqrt(T)

if option_type == 'call':
theta = — (S * norm.pdf(d1) * sigma) / (2 * np.sqrt(T)) — r * K * np.exp(-r * T) * norm.cdf(d2)
elif option_type == 'put':
theta = — (S * norm.pdf(d1) * sigma) / (2 * np.sqrt(T)) + r * K * np.exp(-r * T) * norm.cdf(-d2)
else:
raise ValueError(«Invalid option type. Use 'call' or 'put'.»)

return theta
# Параметры опциона
S = 100 # Текущая цена акции
K = 100 # Цена страйк
T = 1 # Срок до экспирации (в годах)
r = 0.05 # Безрисковая процентная ставка
sigma = 0.2 # Волатильность акции
option_type = 'call' # Тип опциона: 'call' (колл) или 'put' (пут)
theta = black_scholes_theta(S, K, T, r, sigma, option_type)
print(f«Theoretical Theta for the {option_type} option: {theta:.4f}»)

Это чатГПТ выжал на твой коммент. Но тут вставляется без отступов..
+1
- Врач-бондиатОр
  22 августа 2023, 22:35
  Иван-дурак, уже пробовал — результат не совпадает с калькуляторами
  0
3Qu
22 августа 2023, 22:49
Для теты нет, но есть вот такой:
def InsKline1m(self, sd): # thname = 'kline' """ kl = {'stream': 'btcusdt@kline_1m', 'data': {'e': 'kline', 'E': 1692194018649, 's': 'BTCUSDT', 'k': {'t': 1692193980000, 'T': 1692194039999, 's': 'BTCUSDT', 'i': '1m', 'f': 4009649234, 'L': 4009650704, 'o': '29087.10', 'c': '29094.90', 'h': '29095.20', 'l': '29087.10', 'v': '148.002', 'n': 1471, 'x': False, 'q': '4305614.76410', 'V': '88.941', 'Q': '2587451.86430', 'B': '0'}}} """ ""«if threading.get_ident() != self.threads[thname]['thread_id']: # print('aaaa ',threading.get_ident()) self.threads[thname]['thread_id'] = threading.get_ident() self.OpenDB(thname) „“» dt = self.ConvertToDTime(sd['data']['E']) var = (sd['data']['E'],sd['data']['s'],dt,sd['data']['k']['t'],sd['data']['k']['T'],sd['data']['k']['f'],sd['data']['k']['L'], sd['data']['k']['o'],sd['data']['k']['c'],sd['data']['k']['h'],sd['data']['k']['l'],sd['data']['k']['v'], sd['data']['k']['n'],sd['data']['k']['x'],sd['data']['k']['q'],sd['data']['k']['V'],sd['data']['k']['Q']) #17 sql2 = '(?,?,?,?,?,?,?,?,?,?,?,?,?,?,?,?,?)' sql1 = 'insert into Kline1m (E,Symbol,dtime,tb,tc,f,L,open,close,high,low,vol,n,x,q,V,Qp) values' + sql2 begin = time.monotonic_ns() self.cursor.execute(sql1, var) self.connect.commit() if sd['data']['k']['x']: dt = self.ConvertToDTime(sd['data']['k']['t']) print('dt ', dt) var = (sd['data']['s'],dt,sd['data']['k']['t'],sd['data']['k']['T'],sd['data']['k']['o'], sd['data']['k']['c'],sd['data']['k']['h'],sd['data']['k']['l'],sd['data']['k']['v'], sd['data']['k']['n'],sd['data']['k']['q'],sd['data']['k']['V'],sd['data']['k']['Q']) #13 sql2 = '(?,?,?,?,?,?,?,?,?,?,?,?,?)' sql1 = 'insert into Kline_1mc (Symbol,dtime,tb,tc,open,close,high,low,vol,n,q,V,Qp) values' + sql2 self.cursor.execute(sql1, var) self.connect.commit() Tq = (time.monotonic_ns() — begin)/1000 print('Time of record ', Tq)
Сорри, форматирование слетело.)
0
- Врач-бондиатОр
  22 августа 2023, 22:52
  3Qu, вах! Это прям тот самый код Грааля? )
  0
  - 3Qu
    22 августа 2023, 23:54
    Врач-бондиатОр, одна из функций кода Грааля.))
    0
_sg_
23 августа 2023, 07:10
разделите на 365 и будет Вам счастье
0
- Врач-бондиатОр
  23 августа 2023, 08:34
  _sg_, а результат какого кода делить?
  0
  - _sg_
    23 августа 2023, 08:47
    Врач-бондиатОр,
    последняя строка
    return theta / 365.0
    +1
    - Врач-бондиатОр
      23 августа 2023, 09:16
      _sg_, спасибо, попробую еще раз
      0

Читайте на SMART-LAB:

Деловая активность в услугах снижается второй месяц подряд

PMI в сфере услуг России в апреле вырос до 49,7 пункта после 49,5 в марте, второй месяц подряд оставаясь ниже нейтральной отметки 50 пунктов. Это означает, что деловая активность в секторе услуг...

Freedom Finance Global

11:37

5 идей в российских акциях. Голубые фишки снова на уровнях конца 2025

Индекс МосБиржи за неделю просел более чем на 3%. Кроме того, он обновил пятимесячный минимум. Это значит, что многие голубые фишки торгуются по лучшим ценам с конца прошлого года....

БКС Мир инвестиций

10:07

📊 «Цифра брокер» подтвердил оценку акций «МГКЛ» на уровне 4 рублей, апсайд – более 50%

Аналитики Цифра брокер подтвердили справедливую стоимость акций ПАО «МГКЛ» на уровне 4 рублей за акцию, что предполагает потенциал роста более 50% от текущих уровней. В новом...

МГКЛ

11:34

Газпром: рекордная квартальная прибыль или оставь надежду всяк сюда входящий, но сентимент разворачивается вместо цены акции?

Газпром недавно отчитался по МСФО за 2025 год и по РСБУ за 1-й квартал Рассмотрим все сразу, обновим модель + сравним прогноз-факт, как всегда Начинаем традиционно — ошибки...