| #### 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] |
| ) |