【发布时间】:2018-04-14 15:48:23
【问题描述】:
我有以下代码用于运行自回归过程:
import numpy as np
def AR1(mu,phi,x0, Nt, a, b, neg = False):
AR_1 = []
x1 = np.zeros(len(x0))
#np.random.seed(0)
for j in range(Nt):
#np.random.seed(0)
for i in range(len(x0)):
x1[i] = mu[i] + phi*(x0[i] - mu[i]) + ss.truncnorm.rvs(a[i],b[i])
AR_1.append(x1)
x0 = x1.copy()
x1 = np.zeros(len(x0))
return AR_1
如您所见,我有几个职位可以放置np.random.seed(0)。我希望为每次运行选择相同的随机数集,以便我可以看到变化 phi 的效果。考虑到这一点,我应该将np.random.seed 放在代码中的哪个位置?谢谢!
【问题讨论】: