The R Project SVN R

#### lm, glm, aov, etc --- typically *strict* tests (no *.Rout.save)

data(mtcars)
mtcar2 <- within(mtcars, {
    mpg_c <- mpg * (1+am) + 5
    am <- factor(am)
})
fm2 <- glm(disp ~ am * mpg + mpg_c, data = mtcar2)
c2 <- coef(fm2)
V2 <- vcov(fm2)
jj <- !is.na(c2)
stopifnot(names(which(!jj)) == "am1:mpg"
    , identical(length(c2), 5L), identical(dim(V2), c(5L,5L))
    , all.equal(c2[jj],    coef(fm2, complete=FALSE))
    , all.equal(V2[jj,jj], vcov(fm2, complete=FALSE))
    , all.equal(c2[jj], c(`(Intercept)`= 626.0915, am1 = -249.4183,
                  mpg = -33.74701, mpg_c = 10.97014),
            tol = 7e-7)# 1.01e-7 [F26 Lnx 64b]
)