Esecuzione di una funzione di rxLogit() che si può prevedere a da un glm() e rxGlm() chiamata sugli stessi dati può restituire risultati diversi del previsto.
In questo semplice esempio, corrisponderanno a coefficienti restituiti da glm() e rxGlm(), ma quelli restituiti da rxLogit() potrebbe essere diversi.
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)
È necessario impostare initialValues = NA in rxLogit().
Per rxLogit, initialValues predefinito è NULL. La Guida per rxLogit: "i valori iniziali deve essere valutati in base alla regressione lineare. Ciò può velocizzare la convergenza in modo significativo in molti casi. Se il modello non riesce a convergere utilizzando questi valori, la stima viene automaticamente riavviata utilizzando l'opzione NA per i valori iniziali." Se viene utilizzato NA, "valori iniziali dei parametri vengono calcolati da un passaggio di Media ponderata dei minimi quadrati". Per rxGlm, initialValues per impostazione predefinita NA.
Se i risultati di rxLogit sono diversi in modo imprevisto, potrebbe essere che il modello di convergere con initialValues = NULL e la funzione restituisce risultati diversi. Quando si esegue l'esempio con initalValues = NA nel rxLogit(), tutti risultati corrispondenza.