LincolnPetersen模型中群体大小的置信限和置信区间.pdfVIP

中国现堑竺生型窒室笙±三塑兰查兰全笙奎 ●_-___-____-___————_-—————’————P———_———————_—————_——____-____-——___———————一一 !!!!竺!!!．．一．堕 陈家鼎 陈奇志 北京大学数学科学学院 北京大学光华管理学院 1 引言 在野生动物研究、渔业和生物学中群体数目(即该群体中个体的总数)的估计 一问题而提出的解决方法。这种方法也可以用于其它领域，例如具有某种消 费倾向的消费者的人数的估计，吸毒群体的人数估计，艾滋病病毒携带者人 数的估计等。早在18世纪P．Laplace就使用过捕获再捕获方法。时至今日，这种 方法得到蓬勃发展，针对不同的具体情形提出了多种统计模型和统计分析方 型(见Seber(1982))进行研究。这个模型可如下描述。设有一个封闭的有限总 体(例如一个池塘中的鱼)，个体总数是N(未知的正整数)。从中任取n】个 (例如捕获佗1条鱼)，将这n，个个体做上记号再放回原群体中(如将礼1条鱼 做上记号后放入原池塘)，过一段时间再从该群体中任取n2个个体(如再捕 获n2条鱼)，发现这几2个个体中有m个是带有记号的。问Ⅳ是多少?这里假定 每次从群体中抽取个体时，每个个体都有相同的概率被抽取到。n，和n2是已知 的正整数。 我tfJ 9}H道，若m≥1，则Ⅳ的最大似然估计是 砖=【丐nln-2】(fX]表示不超过x的最大整数)。 若m=0，要给出Ⅳ的合适估计值是不可能的。可以证明，Ⅳ没有无偏估 计。实际上，用f表示第二次取到的扎2个个体中带有记号的个体数，则f是服从 超几何分布的随机变量。若妒(∈)是Ⅳ的无偏估计，则对一切Ⅳ≥max(n1，，礼2)成 立E妒(∈)=N，即 n2 EI，o(m)L(mlN)=Ⅳ (1) mO +国家自然科学基金(100471007)资助项目 中囝现场统计研究会第十三届学术年会论文 2007年8月 这里L(mJⅣ)是群体数目为Ⅳ时f等于m的概率，即 L(m[N)=叼懈／四，m=0，1，…，nO． 这里伽=min(nl，竹2)，q=O(ji≥o)，锘=1．于是有 n2 ∑妒(m)叼四篇=Ⅳ嘴 m=O 此式右端是Ⅳ的钆2+1次多项式，左端是Ⅳ的至多n2次多项式，两端不可能对一 切N≥rnax(nl，佗2)成立等式。这个矛盾表明Ⅳ没有无偏估计。 在实际问题中点估计往往不能满足我们的要求，我们还希望了解到底Ⅳ的估 计值与Ⅳ的真值相差多少，本文要给出Ⅳ的置信限和置信区间及其有效的计算 方法。 2 群体数目的置信限和置信区问 我们可使用统计量方法寻求Ⅳ的置信下限和置信上限(统计量方法见陈家鼎等 (2006))。 沿用前面的记号，令 G(m，N)=R@≤m)，g(m，N)=PⅣ(∈m)，(2) 这里PⅣ(A)表示群体数目为Ⅳ时事件A的概率。 易知H(m，N)=G(m一1，Ⅳ)且 tn C(m，Ⅳ)=∑磷。四甓／嘴 (3) k--O (m=0，1，2，…，no，no=rain(n1，n2))．给定Q∈(0，{)．令 (4) (5) 我