本文整理汇总了Python中Point.Point.from_prior方法的典型用法代码示例。如果您正苦于以下问题：Python Point.from_prior方法的具体用法？Python Point.from_prior怎么用？Python Point.from_prior使用的例子？那么恭喜您, 这里精选的方法代码示例或许可以为您提供帮助。您也可以进一步了解该方法所在类Point.Point的用法示例。
在下文中一共展示了Point.from_prior方法的1个代码示例，这些例子默认根据受欢迎程度排序。您可以为喜欢或者感觉有用的代码点赞，您的评价将有助于系统推荐出更棒的Python代码示例。
示例1: __init__
# 需要导入模块: from Point import Point [as 别名]
# 或者: from Point.Point import from_prior [as 别名]
class Walker:
"""
A walker
"""
def __init__(self):
pass
def initialise(self):
self.point = Point()
self.point.from_prior()
self.num_scalars = len(self.point.scalars)
# Assign a direction for the walker to travel in
self.direction = -np.log(rng.rand(self.num_scalars))
self.direction /= self.direction.sum()
# Give a hard edge to the distribution
self.edge = -np.Inf*np.ones(self.num_scalars)
# Keep track of which choices were made
self.choices = []
self.keep = []
def advance(self, steps=1000):
self.keep.append(self.point.scalars)
choice = rng.choice(np.arange(0, self.num_scalars),\
p=self.direction)
self.choices.append(choice)
self.edge[choice] = self.point.scalars[choice]
accept = 0
for i in xrange(0, steps):
proposal = copy.deepcopy(self.point)
logH = proposal.perturb()
if logH > 0.:
logH = 0.
if np.all(proposal.scalars >= self.edge) \
and rng.rand() <= np.exp(logH):
self.point = proposal
accept += 1
print('Accepted {a}/{b}. Scalars={s}'\
.format(a=accept, b=steps, s=self.point.scalars))
if accept <= 0.05*steps:
return False
return True