xmp13.01 {Devore5}R Documentation

data from Example 13.1

Description

The xmp13.01 data frame has 14 rows and 2 columns.

Format

This data frame contains the following columns:

rate
a numeric vector
emission
a numeric vector

Source

Devore, J. L. (2000) Probability and Statistics for Engineering and the Sciences (5th ed), Duxbury

Examples

data(xmp13.01)
plot(emission ~ rate, data = xmp13.01,
     xlab = "Burner area liberation rate",
     ylab = expression(plain("NO")["x"]*
                 plain("emissions")), las = 1,
     main = "Data from Example 13.1, page 545")
fm1 <- lm(emission ~ rate, data = xmp13.01)
abline(fm1)                   # plot 1, Figure 13.1
if (require(MASS)) {
   sres <- stdres(fm1)
   plot(sres ~ fitted(fm1), ylab = "Standardized residuals",
       xlab = "Fitted values")  # plot 2, Figure 13.1
   abline(h = 0, lty = 2, lwd = 0)   # horizontal reference
}
plot(resid(fm1) ~ fitted(fm1), ylab = "Residuals",
     xlab = "Fitted values")  # alternative plot 2
abline(h = 0, lty = 2, lwd = 0)   # horizontal reference
plot(fitted(fm1) ~ emission, data = xmp13.01)
abline(0, 1)                  # plot 3, Figure 13.1
if (require(MASS)) {
   plot(sres ~ rate, data = xmp13.01) # plot 4
   abline(h = 0, lty = 2, lwd = 0)
   qqnorm(sres)               # plot 5
} else {
   plot(resid(fm1) ~ rate, data = xmp13.01) # plot 4
   qqnorm(resid(fm1))         # plot 5
}
## The residuals versus fitted plot and the normal
## probability plot of the standardized residuals 
plot(fm1, which = 1:2)

[Package Contents]