xmp11.06 {Devore5}R Documentation

data from Example 11.6

Description

The xmp11.06 data frame has 24 rows and 3 columns.

Format

This data frame contains the following columns:

Resp
a numeric vector of the mean number of responses emitted by each subject during single and compound stimuli presentations over a 4-day period.
Stimulus
a numeric vector of stimulus levels. These levels correspond to L1 (moderate intensity light), L2 (low intensity light), T (tone), L1+L2, L1+T, and L2+T.
Subject
a numeric vector identifying the subject (rat).

Source

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

(1971), ``Compounding of discriminative stimuli from the same and different sensory modalities'', J. Experimental Analysis and Behavior, 337-342

Examples

data(xmp11.06)
xmp11.06$Stimulus <- factor(xmp11.06$Stimulus,
  labels = c("L1", "L2", "T", "L1 + L2", "L1 + T", "L2 + T"))
xmp11.06$Subject <- factor(xmp11.06$Subject)
plot(Resp ~ Stimulus, data = xmp11.06, col = "lightgray",
     main = "Data from Example 11.6",
     ylab = "Mean number of responses")
for (i in seq(along = levels(xmp11.06$Subject))) {
    attach(xmp11.06[ xmp11.06$Subject == i, ])
    lines(Resp ~ as.integer(Stimulus), col = i+1, type = "b")
}
legend(0.8, 95, paste("Subject", levels(xmp11.06$Subject)),
       col = 1 + seq(along = levels(xmp11.06$Subject)),
       lty = 1)
fm1 <- lm(Resp ~ Stimulus + Subject, data = xmp11.06)
anova(fm1)     # compare to Table 11.5, page 443
attach(xmp11.06)
means <- sort(tapply(Resp, Stimulus, mean))
means
diff(means)    # successive differences
qtukey(0.95, nmeans = 6, df = 15)  #for Tukey comparisons
detach()