Anirban DasGupta writes:

Well done to Yixin Wang (Columbia University) and Vivian Meng (McGill University), who provided correct solutions to the various parts of the problem.

Yixin Wang Vivian Meng

The problem asked was the following: suppose X1,X2,,Xn are iid N(μ,σ2) where the mean μ is some unknown positive integer, and σ is a completely unknown standard deviation. We are to find the unique MLE of μ and σ2, and show that the MLE of μ converges to the true μ exponentially fast and that the MLE of σ2 is consistent.

Denote the likelihood function by l(μ,σ). Then, directly, whenever x¯12, l(μ+1,σ)l(μ,σ)μ[x¯12], where [z] stands for the integer part of z. Thus, for x¯12, the unique MLE of μ is 1+[x¯12]=[x¯+12]; this is the same as the integer closest to the unrestricted MLE X¯, a very intuitive result. If X¯<12, l(μ,σ) is monotone decreasing in μ over the set of positive integers, and in that case the MLE of μ is 1. By a standard calculus argument, once we have found the MLE μ^ of μ, the MLE of σ2 is σ2^=1ni=1n(Xiμ^)2. To complete the solution, if the true μ>1, Pμ(μ^=μ)=Pμ(μ12X¯μ+12).

The complementary probability, Pμ(μ^μ) therefore is 2[1Φ(n2σ)], which is of the order of en8σ2n.

This is a (somewhat faster than) exponential rate. The argument for the case μ=1 is exactly similar; the expression just needs a very small modification. Obviously, therefore, by the Borel-Cantelli lemma, with probability one, for all large n,μ^=μ. This means, the MLE of σ2 is in fact
even strongly consistent, because 1ni=1n(Xiμ)2 converges a.s. to σ2 by the usual SLLN.