Proč není rxGlm vrací stejné výsledky jako glm?

Ve výchozím nastavení nebude získat stejné výsledky z "rxGlm" stejně jako z "glm".

Ve skutečnosti je třeba nastavit argumenty "dropFirst" na hodnotu TRUE a "dropMain" na hodnotu FALSE a reprodukovat výsledky z glm, protože
RevoScaleR použije SAS kontrastů ve výchozím nastavení než výchozí kontrastů R.

 

Zde je některé ukázkových dat a kód použitý v tomto problému, který ukazuje, jak získat dvě funkce odpovídající výsledky testování:

basictestdata <- data.frame(  Factor1 = as.factor(c(1,1,1,1,2,2,2,2)), 
 Factor2 = as.factor(c(1,1,2,2,1,1,2,2)), 
 Discount = c(1,2,1,2,1,2,1,2), 
 Exposure = c(24000, 40000, 7000, 14000, 7500, 15000, 2000, 5600), 
 PurePrem = c(46,32,73,58,48,25,220,30))

GLM.1 <- glm(PurePrem ~ Factor1 * Factor2 - 1, 
 family = tweedie(var.power = 1.5, link.power = 0), 
 data = basictestdata, weights = Exposure 
 , offset = log(Discount))

rxGlm.1 <- rxGlm(PurePrem ~ Factor1 * Factor2 - 1 + offset(log(Discount)), 
 family = rxTweedie(var.power = 1.5, link.power = 0), 
 data = basictestdata, fweights = "Exposure", dropFirst = TRUE, dropMain = FALSE)

coef(GLM.1) 
coef(rxGlm.1)

Potřebujete další pomoc?

Rozšiřte své znalosti a dovednosti
Projít školení
Získejte nové funkce jako první
Připojit se k programu Microsoft Insider

Byly tyto informace užitečné?

Děkujeme za váš názor!

Děkujeme vám za váš názor! Pravděpodobně bude užitečné, když vás spojíme s některým z našich agentů podpory Office.

×