importancia de la diferencia entre dos recuentos

15

¿Hay alguna manera de determinar si una diferencia entre un recuento de accidentes de tráfico en el momento 1 es significativamente diferente de un recuento en el momento 2?

He encontrado diferentes métodos para determinar la diferencia entre grupos de observaciones en diferentes momentos (por ejemplo, comparar medias de Poisson) pero no para comparar solo dos recuentos. ¿O es inválido incluso intentarlo? Cualquier consejo o dirección sería apreciada. Estoy feliz de dar seguimiento a los conductores.

jessop
fuente

Respuestas:

11

Si está contento de asumir que cada recuento sigue una distribución de Poisson (con su propia media bajo la hipótesis alternativa; con una media común debajo de la nula) no hay problema, es solo que no puede verificar esa suposición sin réplicas. La sobredispersión puede ser bastante común con los datos de conteo.

Una prueba exacta dada recuentos y x 2 es sencilla, puesto que el total general de los recuentos de n = x 1 + x 2 es accesoria; condicionando da X 1B i n ( 1X1X2norte=X1+X2como la distribución de su estadística de prueba bajo nulo. X1siyonorte(12,norte) Es un resultado intuitivo: el recuento general, que refleja tal vez cuánto tiempo podría molestarse en observar los dos procesos de Poisson, no contiene información sobre sus tasas relativas, pero afecta el poder de su prueba; y, por lo tanto, otros recuentos generales que pueda tener son irrelevantes.

Ver prueba de hipótesis basada en verosimilitud para la prueba de Wald (una aproximación).

† Cada recuento tiene una distribución de Poisson con media λ i f X ( x i ) = λ x i iXyoλyo Reparametrizar como θ

FX(Xyo)=λyoXyomi-λyoXyo!yo=1,2
dondeθes lo que le interesa, yϕes un parámetro molesto. La función de masa articular se puede volver a escribir: f X 1 , X
θ=λ1λ1+λ2ϕ=λ1+λ2
θϕ El recuento totalnes auxiliar paraθ, teniendo una distribución de Poisson con mediaϕfN(n)
FX1,X2(X1,X2)=λ1X1λ2X2mi-(λ1+λ2)X1!X2!FX1,norte(X1,norte)=θX1(1-θ)norte-X1ϕnortemi-ϕX1!(norte-X1)!
norteθϕ
fN(n)=x1=0fX1,N(x1,n)=ϕneϕn!x1=0n!x1!(nx1)!θx1(1θ)nx1=ϕneϕn!
while the conditional distribution of X1 given n is binomial with Bernoulli probability θ & no. trials n
fX1|n(x1;n)=fX1,N(x1,n)fN(n)=θx1(1θ)nx1ϕneϕx1!(nx1)!n!ϕneϕ=n!x1!(nx1)!θx1(1θ)nx1
Scortchi - Reinstate Monica
fuente
The total number of counts is a complete sufficient statistic, isn't it? How can it be ancillary? Have I misunderstood something?
JohnK
@JohnK: The sufficient statistic is (X1,N), N being the ancillary complement to X1. Note the distribution of N doesn't depend on θ.
Scortchi - Reinstate Monica