コード例 #1
0
 def corrupt_result(self, a=100, b=200):
     #tau = gaussian(u=random_sequence_mc(1), mu = 0, sigma = np.exp(np.array([self.c])))
     #tau = np.random.normal(scale=np.exp(self.c))
     add_time = a+random_sequence_mc(1)*(b-a)
     self.add_time = add_time.flatten()
     self.corrupted_times = self.stop_times +self.add_time
コード例 #2
0
ファイル: compare_gaussians.py プロジェクト: buchhola/ABC
def brownian_self(steps):
    u = random_sequence_mc(steps)
    brownian_increments = np.cumsum(gaussian(u))
    return brownian_increments