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.
