Displaying 1 result from an estimated 1 matches for "pr1d".
Did you mean:
  pr1
  
2002 Jul 26
0
predict() still not yet "safe".. (PR#1840)
...ars)
cars1  <- lm(dist ~      speed,              data = cars)# the `same'
all.equal(predict(cars.1), predict(cars1))# the fitted values are ok
cars.4 <- lm(dist ~ poly(speed, degree = 4), data = cars)
newd <- data.frame(speed = d <- seq(0, 25, length = 51)) # shorter than in ?cars
pr1d <- predict(cars.1, newd)
pr4d <- predict(cars.4, newd)
plot(cars, xlab = "Speed (mph)", ylab = "Stopping distance (ft)",
     las = 1, xlim = c(0, 25))
abline(cars1, col = "blue")# the correct one
## in R 1.5.0 and 1.5.1, these look wrong (particularly to the le...