Uruchomiono funkcję rxLogit(), które można spodziewać się za równoważne glm() i rxGlm() połączenia na tych samych danych mogą zwracać inne wyniki, niż oczekiwano.
W poniższym przykładzie prosty dopasuje współczynniki zwracane przez glm() i rxGlm(), ale te zwrócone przez rxLogit() może być inna.
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)
Najważniejsze jest, aby ustawić initialValues = NA w rxLogit().
RxLogit initialValues jest domyślnie na wartość NULL. Z pomocy dla rxLogit: "wartości początkowe oznaczana będzie na podstawie regresji liniowej. Można to przyspieszyć konwergencji znacznie w wielu przypadkach. Jeśli model nie zostanie uzyskana zbieżność funkcji za pomocą tych wartości, szacowania jest automatycznie ponownie uruchomiony za pomocą opcji NA wartości początkowe." Jeżeli zastosowano NA "początkowe wartości parametrów są obliczane po kroku ważonej najmniejszych kwadratów". RxGlm initialValues jest domyślnie NA.
Więc jeśli wyniki z rxLogit są nieoczekiwanie różne, może być że model jest zbieżna z initialValues = NULL i funkcja zwraca różne wyniki. Podczas uruchamiania przykładu z initalValues = NA w rxLogit(), wszystkie dopasowania wyników.
