def __init__(self, value=None, mean=0.0, sd=1.0, proposal_sd=1.0, **kwargs): Stochastic.__init__(self, value=value, **kwargs) self.mean = mean self.sd = sd self.proposal_sd = proposal_sd if value is None: self.set_value(ilogit(norm.rvs(loc=mean, scale=sd)))
def propose(self): ret = copy(self) ret.value = ilogit(norm.rvs(loc=logit(self.value), scale=self.proposal_sd)) return ret, 0.0 # symmetric
def propose(self): ret = copy(self) ret.value = ilogit( norm.rvs(loc=logit(self.value), scale=self.proposal_sd)) return ret, 0.0 # symmetric