-
Notifications
You must be signed in to change notification settings - Fork 315
/
reg-tests-3.R
35 lines (29 loc) · 995 Bytes
/
reg-tests-3.R
1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
29
30
31
32
33
34
35
### Regression tests for which the printed output is the issue
### May fail, e.g. by needing Recommended packages
# postscript("reg-tests-3.ps")
## str() for character & factors with NA (levels), and for Surv objects:
ff <- factor(c(2:1, NA), exclude = NULL)
str(levels(ff))
str(ff)
str(ordered(ff, exclude=NULL))
if(require(survival)) {
data(aml)
(sa <- Surv(aml$time, aml$status))
str(sa)
detach("package:survival")
}
## were different, the last one failed in 1.6.2 (at least)
## lm.influence where hat[1] == 1
if(require(MASS)) {
fit <- lm(formula = 1000/MPG.city ~ Weight + Cylinders + Type + EngineSize + DriveTrain, data = Cars93)
print(lm.influence(fit))
## row 57 should have hat = 1 and resid=0.
summary(influence.measures(fit))
}
## only last two cols in row 57 should be influential
## PR#6640 Zero weights in plot.lm
if(require(MASS)) {
fm1 <- lm(time~dist, data=hills, weights=c(0,0,rep(1,33)))
plot(fm1)
}
## gave warnings in 1.8.1