Post your answers to the April 27 Let's Experiment problem here.

Best, Chris

Mu: 100 = 0.00336062777778 Mu: 1000 = 0.00144251111111 Mu: 10000 = 0.000346444444444

If star does not appear in the document, with mu: 1000 = 000109177777778

a) for mu = 100 M(t) = .00336 for mu = 1000 M(t) = .00144 for mu = 10000 M(t) = .000346

b) for mu = 1000 and frequency of 0 M(t) = .00011

term_freq = 2

length_doc = 500

mus = [1000, 10000, 100000]

McT = float(73695) / 450000000

for mu in mus:

nom =  (term_freq + mu * McT)

denom = length_doc + mu

print nom / denom

Mc(t) = 73695 / 450000000 = 0.00016376666

For μ = 100 Mμd(t) = (2 + 100 * 0.00016376666) / (500 + 100) = 0.00336062777

For μ = 1000 Mμd(t) = (2 + 1000 * 0.00016376666) / (500 + 1000) = 0.0014425111

For μ = 10000 Mμd(t) = (2 + 10000 * 0.00016376666) / (500 + 10000) = 0.00034644443

If the frequency of term "star" is 0, Mμd(t) = (1000 * 0.00016376666) / (500 + 1000) = 0.00010917777

For μ = 100,

M_μ_d(t) = (2 + (100 * (73695/450000000)))/(500 + 100) = 0.00336

For μ = 1000,

M_μ_d(t) = (2 + (1000 * (73695/450000000)))/(500 + 1000) = 0.00144

For μ = 10000,

M_μ_d(t) = (2 + (10000 * (73695/450000000)))/(500 + 10000) = 0.00034

If term "star" does not appear in d,

M_μ_d(t) = (0 + (1000 * (73695/450000000)))/(500 + 1000) = 0.000109

MC(t)=lt / lC = 73695 / 450000000 = 0.00016376666

Mμd(t) for μ=100, Mμd(t) = (ft,d+μ⋅MC(t)) / (ld+μ) = (2 + 100 * 0.00016376666) / (500 + 100) = 0.00336062777

Mμd(t) for μ=1000, Mμd(t) = (ft,d+μ⋅MC(t)) / (ld+μ) = (2 + 1000 * 0.00016376666)/(500 + 1000) = 0.0014425111

Mμd(t) for μ=10000, Mμd(t) = (ft,d+μ⋅MC(t)) / (ld+μ) =(2 + 10000 * 0.00016376666)/(500 + 10000) = 0.00034644443

If star does not appear, then,

Mμd(t) for μ=1000, Mμd(t) = (ft,d+μ⋅MC(t)) / (ld+μ) = (0 + 1000 * 0.00016376666)/(500 + 1000) = 0.00109177773