> source("em.r") > EM.censored.poisson(5,4.2,10) 0 2.066667 -4.549411 1 1.849275 -4.357812 2 1.832689 -4.356601 3 1.831319 -4.356592 4 1.831205 -4.356592 5 1.831196 -4.356592 6 1.831195 -4.356592 7 1.831195 -4.356592 8 1.831195 -4.356592 9 1.831195 -4.356592 10 1.831195 -4.356592 11 1.831195 -4.356592 12 1.831195 -4.356592 13 1.831195 -4.356592 14 1.831195 -4.356592 15 1.831195 -4.356592 16 1.831195 -4.356592 17 1.831195 -4.356592 18 1.831195 -4.356592 19 1.831195 -4.356592 20 1.831195 -4.356592 [1] 1.831195 > EM.censored.poisson(5,4.2,1000) 0 1.015920 -319.5924 1 0.522337 -118.3403 2 0.3623037 -76.25878 3 0.2855217 -62.10686 4 0.2418962 -56.27045 5 0.2147062 -53.58567 6 0.1967718 -52.26782 7 0.1844963 -51.59381 8 0.1758799 -51.23961 9 0.1697244 -51.05002 10 0.1652714 -50.94723 11 0.1620207 -50.89099 12 0.159632 -50.86001 13 0.1578681 -50.84287 14 0.1565610 -50.83336 15 0.1555897 -50.82806 16 0.1548666 -50.82511 17 0.1543275 -50.82346 18 0.1539251 -50.82254 19 0.1536245 -50.82202 20 0.1533998 -50.82173 [1] 0.1533998