関数が出力されない。math domain error

Question

適合度関数のSLSCを求めるにあたり関数が必要でSiと置いたのだが、うまく出力されません。
Si（１）、Si（２）は正しく出力されるのですがSi（３）以降がエラーが出ます。
Si(3)の中を取り出して実行したらうまくいくのですが、Si（３）として出力するとうまくいきません。

以下プログラミング実行画面です。
//////////////////////////////////////////////////////////////////////////
import numpy as np
import math
import pandas as pd
df = pd.read_excel('C:/Users/yamah/OneDrive/ドキュメント/seminaaaBkanazawa1971-2020.xlsx')

def f(x):
return ((1000/24)(np.exp((-x+74)/24))(np.exp(-np.exp((-x+74)/24))))

def Si(i):
return -math.log(-math.log(f(df.iloc[i-1,0])))

print (Si(3))
/////////////////////////////////////////////////////////////////////////////

以下、結果です。
/////////////////////////////////////////////////////////////////////////////
In[12]:df.iloc[3-1,0]
Out[12]: 158.5

In[13]:runfile('C:/Users/yamah/OneDrive/Office のカスタム テンプレート/Python/Si surch.py', wdir='C:/Users/yamah/OneDrive/Office のカスタム テンプレート/Python')
Traceback (most recent call last):

File "C:\Users\yamah\OneDrive\Office のカスタム テンプレート\Python\Si surch.py", line 20, in <module>
print (Si(3))

File "C:\Users\yamah\OneDrive\Office のカスタム テンプレート\Python\Si surch.py", line 17, in Si
return -math.log(-math.log(f(df.iloc[i-1,0])))

ValueError: math domain error

Answer 1

負の数の対数を取ろうとして定義域エラー(domain error)が出ています。

対数関数の定義域(domain)は正の実数です。
対数関数のグラフと書き方【数学Ⅱ】参照

>>> import math
>>> math.log(-1)
Traceback (most recent call last):
  File "<stdin>", line 1, in <module>
ValueError: math domain error

Comments

[mayottarago]

ご回答ありがとうございます

f(df.iloc[3-1,0]) は正なので定義域の範囲だと思うのですが
Python上のエラーなのでしょうか

[ppaul]

対数を2回とっていますのでエラーになる可能性は常にあります。
>>> -math.log(10)
-2.302585092994046
>>> -math.log(-math.log(10))
Traceback (most recent call last):
File "<stdin>", line 1, in <module>
ValueError: math domain error