Stima del numero di palline selezionando successivamente una pallina e segnandola


9

Diciamo che ho N palline in una borsa. Al mio primo pareggio, segnare la palla e sostituirla nella borsa. Al mio secondo pareggio, se raccolgo una palla segnata la rimetto nella borsa. Se, tuttavia, raccolgo una palla non segnata, la segnerò e la rimetto nella borsa. Lo continuo per qualsiasi numero di pareggi. Qual è il numero atteso di palloni nella borsa dato un numero di pareggi e la storia marcata / non marcata dei disegni?


1
Possibilmente correlato: hai esaminato il metodo di cattura-riconquista per stimare l'abbondanza della popolazione? en.wikipedia.org/wiki/Mark_and_recapture
a.arfe,

"Numero previsto" non può essere inteso nel suo senso usuale tecnica di un valore atteso, perché non c'è distribuzione di probabilità per N . Sembra che si sta chiedendo uno stimatore di N .
whuber

Risposte:


2

Ecco un'idea. Lasciate essere un sottoinsieme finito di numeri naturali che servirà come i possibili valori per . Supponiamo di avere una distribuzione precedente su . Fissare un non casuale intero positivo . Sia la variabile casuale che indica il numero di volte che segniamo una palla in pesca dal sacchetto. L'obiettivo è trovare . Questa sarà la funzione di e del precedente. N I M k M E ( N | k ) M , kINIMkME(N|k)M,k

Per regola di Bayes abbiamo

P(N=j|k)=P(k|N=j)P(N=j)P(k)=P(k|N=j)P(N=j)rIP(k|N=r)P(N=r)

Il calcolo è un calcolo noto che è una variante del problema dei collezionisti di coupon. è la probabilità che osserviamo distinti coupon in disegni quando ci sono coupon in totale. Vedi qui per un argomento perP ( k | N = j ) k M jP(k|N=j)P(k|N=j)kMj

P(k|N=j)=(jk)k!S(M,k)jM

dove indica un numero stirling del secondo tipo . Possiamo quindi calcolareS

E(N|k)=jIjP(N=j|k)

Questi sono alcuni calcoli di diverse e . In ogni caso usiamo un'uniforme prima diM [ k , 10 k ]kM[k,10k]

MkE(N)1057.991555.60151023.69301520.00302039.53
Utilizzando il nostro sito, riconosci di aver letto e compreso le nostre Informativa sui cookie e Informativa sulla privacy.
Licensed under cc by-sa 3.0 with attribution required.