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
def compute_prior(self): return norm.logpdf(logit(self.value), loc=self.mean, scale=self.sd)