Análisis probit de dosis-respuesta<meta http-equiv="Content-Type" content="text/html; charset=UTF-8" />hola grupo,
Estoy tratando de reproducir un análisis simple de respuesta a la dosis que había realizado en SPSS. Entonces, en JMP estoy usando el módulo "Ajustar modelo", seleccionando modelos lineales generalizados, distribución binomial y enlace probit.
Tengo tres variables, número de ensayos realizados (10 ensayos), número de ensayos con resultado positivo y concentración (de 0 a 36)
0 10 2
1 10 4
2 10 9
6 10 16
8 10 25
9 10 36
0 10 0
Selecciono el número de ensayos realizados como variable "Y", el número de ensayos con resultado positivo como "Frecuencia" y la concentración en registros como variable para modelar los efectos. Pero cuando presiono ejecutar, aparece este mensaje en los resultados:
"Advertencia: Hessian no definido positivo. Es probable que la singularidad se deba a un ajuste deficiente del modelo o a dependencias lineales entre las covariables del modelo. El ajuste y el resultado tienen un valor cuestionable"
Y no se realiza el análisis.
¿Alguien puede compartir algunos comentarios?
Gracias,
dave
Re: Análisis probit de dosis-respuestaDe hecho, ese es el orden correcto para alimentar las variables. Gracias por señalarlo, Michael.
dave
Re: Análisis probit de dosis-respuestaLa última fila 0 10 0 de los datos originales, por alguna razón, evita que el análisis necesite el orden correcto. Si lo dejas fuera, el orden sí importa.
Miguel
Re: Análisis probit de dosis-respuestaHe estado siguiendo este hilo con interés.
De los documentos JMP:
Si sus datos se resumen como frecuencias de eventos y ensayos, especifique dos
columnas en este orden : una cuenta del número de éxitos, y una cuenta del número de
juicios Alternativamente, puede especificar el número de fallas en lugar de éxitos.
(mi énfasis).
Entonces, si cambia el orden en que se ingresan las 2 variables y, creo que verá más de lo que esperaba, especialmente si cambia el eje x para que sea logarítmico.
Dicho esto, no puedo explicar por qué el orden no parece importar con sus datos originales.
Miguel
Re: Análisis probit de dosis-respuestaEstimado DanO,
Disculpa que te vuelva a molestar, pero cambiando los números a la demo, parece que me dan resultados erróneos. Eche un vistazo a este conjunto de datos, con las mismas variables pero con números diferentes. Para las tres concentraciones más altas obtengo 100% de positivos
Número de resultados positivos Número de ensayos realizados Concentración
0 1140 0.01
241 1140 0.1
590 1140 0,5
963 1140 1
1134 1140 5
1140 1140 10
1140 1140 50
1140 1140 100
Ejecuto el Probit como sugeriste
Después de ejecutar, obtengo un resultado extraño que no tiene sentido para mí. En lugar de hacer porcentajes de correcciones en la Y, parece estar haciendo algo diferente. ¿Qué estoy haciendo mal?
Muchas gracias por tu ayuda
dave
Re: Análisis probit de dosis-respuestaMuchas gracias DanO, esa es una respuesta muy clara. Es interesante que JMP sepa que están vinculados.
Con respecto a la transformación del registro, sí, lo tendré en cuenta (de hecho, se mencionó en mi publicación original, ¡aunque debo admitir que es bastante vago!)
Mejor,
dave
Re: Análisis probit de dosis-respuestadavid,
Su "número de ensayos con resultado positivo" también debe ingresar como una variable "Y". El rol de frecuencia es para indicar cuántas veces ocurrió realmente una fila en la tabla de datos. En este caso, el modelo GLM espera que el número de ensayos realizados y el número de resultados positivos sean valores Y "separados". El análisis los vinculará.
Una cosa que no mencionaste que se suele hacer para este tipo de trabajo es llevar la bitácora de la concentración. Hice esto en el cuadro de diálogo del modelo que se muestra arriba. Resalte el efecto del modelo de concentración, haga clic en el botón de transformación y elija registro. Los resultados se publican aquí.
Dan
