Percentil y bondad de ajuste de binomial negativo inflado cero (ZINB) y binomial negativo (NB)Hola a todos,
(1) Me gustaría obtener la bondad de ajuste de la distribución binomial negativa inflada cero (ZINB) y binomial negativa (NB) del siguiente conjunto de datos.
(2) También quiero calcular el percentil 95 de cada distribución (ZINB y NB) con JMP 14 o JMP 14 Pro
(3) ¿Es posible trazar las funciones de densidad de ZINB y NB frente a la densidad empírica de los datos simulados? (Y: probabilidad acumulada, X: cuentas)
Gracias de antemano.: )
| número de muestra | Datos |
| 1 | 0 |
| 2 | 0 |
| 3 | 0 |
| 4 | 0 |
| 5 | 0 |
| 6 | 0 |
| 7 | 0 |
| 8 | 0 |
| 9 | 0 |
| 10 | 0 |
| 11 | 0 |
| 12 | 0 |
| 13 | 0 |
| 14 | 0 |
| 15 | 2 |
| dieciséis | 0 |
| 17 | 0 |
| 18 | 0 |
| 19 | 0 |
| 20 | 0 |
| 21 | 0 |
| 22 | 0 |
| 23 | 0 |
| 24 | 0 |
| 25 | 0 |
| 26 | 1 |
| 27 | 0 |
| 28 | 0 |
| 29 | 0 |
| 30 | 0 |
| 31 | 0 |
| 32 | 0 |
| 33 | 0 |
| 34 | 0 |
| 35 | 0 |
| 36 | 0 |
| 37 | 0 |
| 38 | 0 |
| 39 | 0 |
| 40 | 0 |
| 41 | 0 |
| 42 | 0 |
| 43 | 0 |
| 44 | 0 |
| 45 | 0 |
| 46 | 0 |
| 47 | 0 |
| 48 | 0 |
| 49 | 0 |
| 50 | 0 |
| 51 | 0 |
| 52 | 0 |
| 53 | 0 |
| 54 | 4 |
| 55 | 0 |
| 56 | 0 |
| 57 | 0 |
| 58 | 0 |
| 59 | 0 |
| 60 | 0 |
| 61 | 0 |
| 62 | 0 |
| 63 | 0 |
| 64 | 0 |
| sesenta y cinco | 0 |
| 66 | 0 |
| 67 | 0 |
| 68 | 0 |
| 69 | 0 |
| 70 | 2 |
| 71 | 0 |
| 72 | 0 |
| 73 | 0 |
| 74 | 0 |
| 75 | 0 |
| 76 | 0 |
| 77 | 0 |
| 78 | 0 |
| 79 | 0 |
| 80 | 0 |
| 81 | 0 |
| 82 | 0 |
| 83 | 0 |
| 84 | 0 |
| 85 | 0 |
| 86 | 0 |
| 87 | 0 |
| 88 | 0 |
| 89 | 0 |
| 90 | 0 |
| 91 | 0 |
| 92 | 0 |
| 93 | 0 |
| 94 | 0 |
| 95 | 0 |
| 96 | 0 |
| 97 | 0 |
| 98 | 0 |
| 99 | 0 |
| 100 | 0 |
| 101 | 0 |
| 102 | 0 |
| 103 | 0 |
| 104 | 0 |
| 105 | 0 |
| 106 | 0 |
| 107 | 0 |
| 108 | 0 |
| 109 | 0 |
| 110 | 0 |
| 111 | 0 |
| 112 | 0 |
| 113 | 0 |
| 114 | 0 |
| 115 | 0 |
| 116 | 0 |
| 117 | 0 |
| 118 | 0 |
| 119 | 0 |
| 120 | 0 |
| 121 | 1 |
| 122 | 0 |
| 123 | 1 |
| 124 | 0 |
| 125 | 0 |
| 126 | 0 |
| 127 | 0 |
| 128 | 0 |
| 129 | 0 |
| 130 | 0 |
| 131 | 0 |
| 132 | 0 |
| 133 | 0 |
| 134 | 0 |
| 135 | 0 |
| 136 | 0 |
| 137 | 0 |
| 138 | 0 |
| 139 | 0 |
| 140 | 0 |
| 141 | 1 |
| 142 | 0 |
| 143 | 0 |
| 144 | 0 |
| 145 | 0 |
| 146 | 0 |
| 147 | 0 |
| 148 | 0 |
| 149 | 0 |
| 150 | 0 |
| 151 | 0 |
| 152 | 0 |
| 153 | 0 |
| 154 | 0 |
| 155 | 0 |
| 156 | 0 |
| 157 | 0 |
| 158 | 0 |
| 159 | 2 |
| 160 | 0 |
| 161 | 0 |
| 162 | 0 |
| 163 | 0 |
| 164 | 0 |
| 165 | 0 |
| 166 | 0 |
| 167 | 0 |
| 168 | 0 |
| 169 | 0 |
| 170 | 0 |
| 171 | 0 |
| 172 | 0 |
| 173 | 0 |
| 174 | 0 |
| 175 | 0 |
| 176 | 0 |
| 177 | 0 |
| 178 | 0 |
| 179 | 0 |
| 180 | 0 |
| 181 | 0 |
| 182 | 0 |
| 183 | 0 |
| 184 | 0 |
| 185 | 0 |
| 186 | 0 |
| 187 | 0 |
| 188 | 0 |
| 189 | 0 |
| 190 | 0 |
| 191 | 0 |
| 192 | 0 |
| 193 | 0 |
| 194 | 0 |
| 195 | 0 |
| 196 | 0 |
| 197 | 1 |
Re: Percentil y bondad de ajuste de binomial negativo inflado cero (ZINB) y binomial negativo (NB)Gracias por la clara respuesta
Re: Percentil y bondad de ajuste de binomial negativo inflado cero (ZINB) y binomial negativo (NB)1.) Utilice las estimaciones de parámetros para el modelo de distribución ZINB o NB para estimar el percentil del 95 % a partir de la función de distribución acumulada inversa.
2.) AIC y BIC son criterios para la selección de modelos. No evalúan la bondad de ajuste. JMP no tiene un gráfico QQ para estas distribuciones que yo sepa. Sin embargo, existen en la plataforma de distribución para gamma Poisson.
Re: Percentil y bondad de ajuste de binomial negativo inflado cero (ZINB) y binomial negativo (NB)¡Gracias Marcos!
Y... ¿Puedo obtener más de su consejo? 
Si el otro conjunto de datos coincide con ZINB o NB, necesito calcular el percentil 95 y el percentil 99 de esa distribución según la guía.
(1) ¿Es posible obtener estos percentiles de la plataforma de regresión generalizada JMP Pro o de otra forma?
(2) Si el conjunto de datos no se ajusta bien con Poisson o NB (tal vez pueda verificarlo con GammaPoisson), debo verificar la bondad del ajuste para ZINB. ¿Hay alguna forma de verificarlo visualmente (gráfico QQ) o con otro método que no sea AIC o BIC?
¡Muchas gracias!Re: Percentil y bondad de ajuste de binomial negativo inflado cero (ZINB) y binomial negativo (NB)Normalmente usaría la plataforma de distribución para ajustar dichos modelos, pero la distribución binomial negativa no es una de las opciones. Aquí modelé sus datos con la distribución gamma de Poisson, obtuve las estadísticas de la prueba de bondad de ajuste y estimé el cuantil del 95 %:

Parece que va bien.
Los modelos de distribución binomial negativa y binomial negativa inflada con ceros están disponibles en la plataforma de regresión generalizada JMP Pro. Usé un predictor lineal de solo intercepción para obtener estos resultados:


El AICc sugiere que la distribución gamma de Poisson es una de las que mejor se ajusta a los datos.