Python MHSampler.compute_posteriorの例

コード例 #1

ファイル: TabooMCMC.py プロジェクト: jthurst3/LOTlib

 def compute_posterior(self, h, data):
     """
             Wrap the posterior with a penalty for how often we've seen h. Computes the penalty on the prior
     """
     mypenalty = self.seen[h] * self.penalty
     np, nl = MHSampler.compute_posterior(self, h, data)
     return np+mypenalty, nl

コード例 #2

ファイル: TabooMCMC.py プロジェクト: sa-/LOTlib

 def compute_posterior(self, h, data):
     """
             Wrap the posterior with a penalty for how often we've seen h. Computes the penalty on the prior
     """
     mypenalty = self.seen[h] * self.penalty
     np, nl = MHSampler.compute_posterior(self, h, data)
     return np + mypenalty, nl

コード例 #3

ファイル: MemoizedMHSampler.py プロジェクト: ebigelow/LOTlib

 def compute_posterior(self, h, data):
     if h in self.mem:
         ret = self.mem[h]
         h.posterior_score = ret # set this because it may not be set
         return ret
     else:
         ret = MHSampler.compute_posterior(self, h, data)
         self.mem[h] = ret
         return ret

コード例 #4

 def compute_posterior(self, h, data):
     if h in self.mem:
         ret = self.mem[h]
         h.posterior_score = sum(ret) # set this because it may not be set
         return ret
     else:
         ret = MHSampler.compute_posterior(self, h, data)
         self.mem[h] = ret
         return ret

コード例 #5

ファイル: MemoizedMHSampler.py プロジェクト: joshrule/LOTlib

    def compute_posterior(self, h, data, shortcut=-Infinity):

        if h in self.mem:
            ret = self.mem[h]
            h.posterior_score = ret # set this because it may not be set
            return ret
        else:
            ret = MHSampler.compute_posterior(self, h, data, shortcut=-Infinity) # calls update to posterior counter
            self.mem[h] = ret
            return ret

コード例 #6

    def compute_posterior(self, h, data, shortcut=-Infinity):

        if h in self.mem:
            ret = self.mem[h]
            h.posterior_score = ret  # set this because it may not be set
            return ret
        else:
            ret = MHSampler.compute_posterior(
                self, h, data,
                shortcut=-Infinity)  # calls update to posterior counter
            self.mem[h] = ret
            return ret