# file "REM" -- simple EM algorithm for an iid Exponential model

# observed data:
y = 10

# initial value:
lambda = 2

cat("i =", 0, "  lambda =", lambda, "\n")

numits = 20

for (i in 1:numits) {

    # E-step:
    z = 1/lambda

    # M-step:
    lambda = 2 / (y+z)

    # output as we go:
    cat("i =", i, "  z =", z, "  lambda =", lambda, "\n")

}