估计量$$\hat{N}=x_{(k)}+x_{(k)}/k-1$$

详细计算结果表格:

观测数据 \(x=(x_{ 1},\cdots,x_{ k})\), 不放回来自于 \(1,\cdots,N\), 则最大值 \(x_{(k)}\)的分布为:
$$ P(X_{(k)}=m)=\binom{m-1}{k-1}\Big/\binom{N}{k}. $$ 利用 \(\sum\limits_{m=k}^N\binom{m}{k}=\binom{N+1}{k+1}\), 则 $$E X_{(k)}=\mu=\sum_{m=k}^Nm\binom{m-1}{k-1}\Big/\binom{N}{k}=\frac{k(N+1)}{k+1}.$$ 因此, 可以得到 $$N=\mu(1+k^{-1})-1.$$ 于是由矩估计方法可知 \(N\)的矩估计量为 $$\hat{N}=X_{(k)}(1+k^{-1})-1=m(1+k^{-1})-1.$$ 可以证明, \(\hat{N}\)是$N\)的最小方差无偏估计. 其方差为 $$Var(\hat{N})=(1+k^{-1})^2Var(X_{(k)})=\frac1k\frac{(N-k)(N+1)}{k+2}\approx\frac{N^2}{k^2} \mbox{当}k\ll N.$$ 从而标准差约为 \(N/k\).