Я определил следующую функцию:
def GMM(s1, s2, s3, s4, s5, a):
"""The GMM objective function.
Arguments
---------
si: float
standard deviations of preference distribution
a: float
marginal utility of residutal income
Paramters
---------
Px: array (1,ns)
projector onto nonprice characteristic space
xk, z: arrays (J, 5) and (J, 12)
nonprice char. and instruments
invW: array (12, 12)
GMM weight matrix
Returns
-------
float."""
delta = invert(s1, s2, s3, s4, s5, a, delta0) # Invert market shares to get mean utility
bmean = np.dot(Px, delta) # Project delta onto charancteristic space
xihat = delta - np.dot(xk, bmean) # Compute implied unobservable prod. quality
temp1 = np.dot(xihat.T, z)
if np.any(np.isnan(delta)) == True:
value = 1e+10
else:
value = np.dot(np.dot(temp1, invW), temp1.T)
return np.sqrt(value)
Мой вопрос относится к переменной delta
, связанной внутри функции. Вне функции я установлю начальное значение delta0
. Теперь, в конечном счете, я скрою эту функцию. Я бы хотел, чтобы каждый раз, когда функция GMM
оценивает, delta
из предыдущей оценки используется как новый delta0
. Я попытался определить delta0
как глобальную переменную, но, похоже, это не работало... скорее всего, это была моя ошибка. Хотя, я читал здесь, что, как правило, это плохой подход. Любые предложения?