Изпълнява rxLogit() функция, която се очаква да е равен на glm() и rxGlm() обаждане на същите данни може да върне различни резултати от очакваното.
В примера просто коефициенти, върната от glm() и rxGlm() ще се търси, но тези върнати от rxLogit() може да са различни.
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)
Ключът е да зададете initialValues = NA rxLogit().
За rxLogit initialValues по подразбиране е празен. От помощта за rxLogit: "първоначални стойности се оценява въз основа на линейна регресия. Това може да се ускори конвергенция значително в много случаи. Ако моделът не съвпадат с помощта на тези стойности, оценката е автоматично повторно започва с помощта на възможност за началните стойности." Ако се използва NA "първоначалните стойностите на параметрите се изчисляват от стъпка претеглена най-малко квадратчета". За rxGlm initialValues по подразбиране NA.
Така че ако неочаквано различни резултати от rxLogit, е че моделът съвпадат с initialValues = NULL и връща различни резултати. Когато използвате примера с initalValues = NA rxLogit(), всички резултати match.