โดยค่าเริ่มต้น คุณจะไม่ได้รับผลลัพธ์แบบเดียวกับจาก 'rxGlm' ในขณะที่คุณทำจาก 'glm'
จริง ๆ แล้วที่คุณจำเป็นต้องตั้งค่าอาร์กิวเมนต์ 'dropFirst' เป็น TRUE และ 'dropMain' เป็น FALSE ได้ด้วยเมื่อต้องการทบทวนเกิดผลลัพธ์จาก glm เนื่องจาก
RevoScaleR จะใช้ไม่ใช่ contrasts เริ่มต้นของ R contrasts SAS โดยค่าเริ่มต้น
ต่อไปนี้คือบางตัวอย่างข้อมูลและรหัสที่ใช้ในการทดสอบนี้ปัญหาที่แสดงให้เห็นถึงวิธีการเรียกฟังก์ชันสองเพื่อให้ผลลัพธ์ที่ตรงกัน:
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)