Execută o funcţie rxLogit() care pot rezulta echivalente la un apel glm() și rxGlm() pe același datele poate returna rezultate diferite decât cea așteptată.
În următorul exemplu simplă, coeficienții returnat de glm() și rxGlm() se va potrivi, dar cele returnat de rxLogit() pot fi diferite.
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)
Cheia este să stabiliți initialValues = NA în rxLogit().
Pentru rxLogit, initialValues implicit NULL. De la centrul de ajutor pentru rxLogit: "valorile inițială va fi estimat bazate pe o regresie liniară. Acest lucru poate accelera convergență semnificativ în multe cazuri. Dacă modelul nu converge utilizând aceste valori, estimarea este automat din nou a început utilizând opțiunea de NA pentru valorile inițiale." Dacă se utilizează NA, "inițială valorile parametrilor sunt calculate un pas ponderată mai mici pătrate." Pentru rxGlm, initialValues implicite la NA.
Astfel, dacă rezultatele din rxLogit sunt diferite în mod neașteptat, este posibil ca modelul converge cu initialValues = NULL și funcția returnează rezultate diferite. Când executați exemplul cu initalValues = NA în rxLogit(), toate rezultatele corespondență.
