xmp13.01 {Devore6} | 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. (2003) Probability and Statistics for Engineering and the Sciences (6th 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]