Chạy chức rxLogit() có thể phải tương đương với cuộc gọi glm() và rxGlm() trên cùng một dữ liệu có thể trả lại kết quả khác nhau hơn mong đợi.
Trong ví dụ đơn giản sau hệ trở lại glm() và rxGlm() sẽ phù hợp nhưng những trả lại bởi rxLogit() có thể khác nhau.
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)
Điều quan trọng là đặt initialValues = NA trong rxLogit().
Cho rxLogit initialValues mặc định không. Từ trợ giúp về rxLogit: "giá trị ban đầu sẽ được ước tính dựa trên một hồi quy tuyến tính. Điều này có thể tăng tốc độ hội tụ đáng kể trong nhiều trường hợp. Nếu mô hình không hội tụ bằng cách sử dụng các giá trị, ước tính được tự động lại bắt đầu bằng cách sử dụng tuỳ chọn NA giá trị ban đầu." Nếu sử dụng NA "ban đầu giá trị tham số được tính bằng một bước trọng tối thiểu". Đối với rxGlm, initialValues mặc định NA.
Vì vậy, nếu kết quả từ rxLogit bất ngờ khác, nó có thể kiểu hội tụ với initialValues = NULL chức năng trả lại kết quả khác nhau. Khi chạy mẫu với initalValues = NA trong rxLogit(), tất cả kết quả phù hợp.
