In [1]:
import numpy as np
import thinkbayes2
from thinkbayes2 import Pmf, Cdf, Suite, Beta
import thinkplot
% matplotlib inline
In [7]:
missile = Beta(2,3, label='missile')
thinkplot.Pdf(missile.MakePmf())
missile.Mean()
Out[7]:
In [9]:
blaster9000 = Beta(2,3, label='missile')
thinkplot.Pdf(missile.MakePmf())
In [10]:
blaster9000.Update((2,8))
In [13]:
(blaster9000.MakePmf()).MaximumLikelihood()
Out[13]:
In [14]:
class AlienBlaster(Suite):
def Likelihood(self, data, hypo):
x = hypo
return (x)**12 * (1-x)**12
In [17]:
blaster10000 = AlienBlaster(range(0, 101))
In [18]:
thinkplot.Pdf(blaster10000)
In [23]:
blaster10000.Update('d')
Out[23]:
In [24]:
thinkplot.Pdf(blaster10000)
In [26]:
def Likelihood(self, data, hypo):
like = (1-x)**4 + 4x**2(1-x)**2 + x**4
In [ ]: