Ejecuta una función rxLogit() quizás tenga que ser equivalente a un glm() y rxGlm() de llamada en los mismos datos pueden devolver resultados diferentes de lo esperado.

En el siguiente ejemplo, coincidirán con los coeficientes devueltos por glm() y rxGlm() pero los devueltos por rxLogit() pueden ser diferentes.

 myFormula <- as.formula("y ~ x1 + x2 + x3 ")model <- rxLogit( myFormula, data = sampleData)
modelGLM <- glm(myFormula,family=binomial(logit),data=sampleData)
modelrxGLM <- rxGlm(myFormula,family=binomial(logit), data=sampleData)

La clave es establecer initialValues = NA en rxLogit().

Para rxLogit, initialValues el valor predeterminado es NULL. De la Ayuda de rxLogit: "los valores iniciales se determinará basándose en una regresión lineal. Esto puede acelerar la convergencia significativamente en muchos casos. Si el modelo no converge con estos valores, la estimación es volver a iniciados automáticamente mediante la opción de NA para los valores iniciales". Si se utiliza NA, "valores iniciales de los parámetros se calculan por un paso de mínimos cuadrados ponderados". Para rxGlm, initialValues el valor predeterminado es NA.

Así que si los resultados de rxLogit son diferentes inesperadamente, es posible que el modelo converge con initialValues = NULL y la función devuelve resultados diferentes. Cuando se ejecuta el ejemplo con initalValues = NA en el rxLogit(), la coincidencia de todos los resultados.

¿Necesita más ayuda?

Ampliar sus conocimientos
Explorar los cursos
Obtener nuevas características primero
Unirse a Microsoft Insider

¿Le ha sido útil esta información?

¿Cómo de satisfecho está con la calidad de la traducción?
¿Qué ha afectado a tu experiencia?

¡Gracias por sus comentarios!

×