Modelos lineales generalizados sin intercepciónHola,
Necesito tu ayuda. Quiero hacer un análisis con GLM, pero sin intercepción. No puedo elegir la opción "No interceptar".
Quiero obtener los mismos resultados que allí: http://support.sas.com/kb/24/992.html
No tengo idea de cómo obtener la opción "No interceptar" en GLM, distribución = Binomial, Enlace = Logit
Estaré muy agradecido por tu ayuda.
Re: Modelos lineales generalizados sin intersecciónHola @mhans ,
Es posible en JMP estándar, pero requiere una ligera reestructuración de su tabla. Lo que necesita es una variable de resultado nominal de Ganar o Perder, y lo que tiene es una tabla resumida de las frecuencias de esas ocurrencias. Adjunté una tabla ampliada aquí, con el script logístico nominal guardado (que se analiza más adelante). Observe la forma de estas tablas, en la imagen a continuación está la tabla original con la primera fila seleccionada, en la parte inferior está la tabla expandida, con las 13 filas seleccionadas que corresponden a esa fila resumida. Tenga en cuenta que hay 13, que son los 13 concursos entre Mil y Det, y hay 7 victorias y 6 derrotas, que reflejan la victoria = 7 en la tabla superior. 
El modelo es sencillo de configurar: Analizar > Ajustar modelo, con Resultado como Y, todos los equipos como efectos del modelo y No Interceptar marcado. 
Esto devuelve el siguiente resultado, que verá que coincide:
Fit Model(
Y( :Outcome ),
Effects( :mil, :det, :tor, :new, :bos, :cle, :bal ),
No Intercept( 1 ),
Choose High Target( 1 ),
Personality( "Nominal Logistic" ),
Run( Likelihood Ratio Tests( 1 ), Wald Tests( 0 ) ))
El valor agregado de la regresión generalizada aquí es que puede especificar directamente la estructura ganadora/total al definir su variable y, que maneja la expansión en segundo plano.
@Julian
Re: Modelos lineales generalizados sin intersecciónHola @Julian ,
¡guau, gracias! Ya veo que es posible en JMP. Sólo tengo la versión "normal" de JMP. ¿Es posible hacer este modelo con la versión normal (como regresión logística sin intercepción o GLM)?
En caso negativo, ¿pueden enviarme el resultado en formato pdf?
Me quito el sombrero ante esta solución.
Re: Modelos lineales generalizados sin intersecciónHola @mhans ,
Ahora entiendo mejor lo que intentas hacer. Para hacer coincidir esos resultados (modelo Bradley-Terry), deberá utilizar la plataforma de regresión generalizada, que está disponible solo en JMP Pro. Adjunté una versión JMP de la tabla de datos utilizada en ese ejemplo de SAS y guardé el script en la tabla. Estos son los pasos para reproducir esto:
1. Inicie Analizar > Ajustar modelo.
2. Especifique Regresión Generalizada como Personalidad y Distribución como Binomial.
3. Ingrese Ganancia y Total como Y, Columnas, e ingrese todos los predictores como efectos del modelo.
4. Marque "No interceptar" en la parte inferior.
Su cuadro de diálogo Ajustar modelo ahora debería verse así:
5. Haga clic en Ejecutar, luego, en el cuadro de diálogo resultante, especifique Regresión logística como método de estimación y haga clic en Ir:
Esto devolverá el siguiente resultado, que puede confirmar que coincide con los resultados de SAS.
Aquí está el JSL para este análisis, que también se guarda en la tabla adjunta.
Fit Model(
Y( :win, :total ),
Effects( :mil, :det, :tor, :new, :bos, :cle, :bal ),
No Intercept( 1 ),
Personality( "Generalized Regression" ),
Generalized Distribution( "Binomial" ),
Run(
Fit( Estimation Method( Logistic Regression ), Validation Method( None ) )
)
);
¡Espero que esto ayude!
@Julian
Re: Modelos lineales generalizados sin intersecciónHola @Julian ,
gracias por su respuesta. Me cansé de utilizar la personalidad logística nominal. La variable de respuesta y debería ser nominal, mi variable y era métrica (entre 0 y 1, un porcentaje de victorias). Hice una nueva estructura para mis datos con la variable y: 1 o -1 y copié los datos en el formulario
1 -1 0 0 0 0 0 7 1
1 0 -1 0 0 0 0 9 1
(...)
0 0 0 0 0 1 -1 6 1
(new part of data)
-1 1 0 0 0 0 0 6 -1
-1 0 1 0 0 0 0 4 -1
(...)
0 0 0 0 0 -1 1 7 -1
Los resultados no son los mismos que en SAS.No tengo idea de cómo repararlo...
Re: Modelos lineales generalizados sin intersecciónHola @mhans ,
A diferencia del Modelo lineal generalizado personalidad, la Logística nominal La personalidad en Ajustar modelo se puede ejecutar sin intercepción marcando la casilla "Sin intercepción" en la ventana de diálogo Ajustar modelo. Esta personalidad devolverá resultados equivalentes a la personalidad del modelo lineal generalizado con una distribución de error binomial y una función de enlace logit, suponiendo que no haya marcado las estimaciones ajustadas por sesgo de Firth (con esta casilla marcada, las estimaciones del modelo diferirán). ¿Quizás esto funcione para su situación?
@Julian
