Skip to contents

Postovsky (1970) investigatged the effect of delay in oral practice at the beginning of second language learning. A control condition began oral practice with no delay, while an experimental group had a four-week delay before starting oral practice. The data consists of scores on language skills at the end of six weeks of study.

Students in this study were matched on age, education, former language training, intelligence and language aptitude.

Usage

data("oral")

Format

A data frame with 56 observations on the following 5 variables.

group

Group, a factor with levels Control Exptl

listen

Listening test, a numeric vector

speak

Speaking test, a numeric vector

read

Reading test, a numeric vector

write

Writing test, a numeric vector

Source

Timm, N. H. (1975). Multivariate Analysis with Applications in Education and Psychology. Wadsworth (Brooks/Cole), Exercise 3.12, p. 279.

References

Postovsky, V. A. (1970). Effects of delay in oral practice at the start of second language training. Unpublished doctoral dissertation, University of California, Berleley.

Examples

library(car)
library(candisc)
data(oral)

# make some boxplots
op <- par(mfrow=c(1,4), cex.lab=1.5)
clr <- c("pink", "lightblue")
Boxplot(listen ~ group, data=oral, col = clr, cex.lab = 1.5)
Boxplot(speak ~  group, data=oral, col = clr, cex.lab = 1.5)
Boxplot(read ~   group, data=oral, col = clr, cex.lab = 1.5)
#> [1] "37" "47" "51" "23" "27"
Boxplot(write ~  group, data=oral, col = clr, cex.lab = 1.5)

#> [1] "51" "27"
par(op)

# view the data ellipses
covEllipses(cbind(listen, speak, read, write) ~ group, data=oral,
    variables = 1:4,
    level = 0.40,
    pooled = FALSE,
    fill = TRUE, fill.alpha = 0.05)


oral.mod <- lm(cbind(listen, speak, read, write) ~ group, data=oral)
Anova(oral.mod)
#> 
#> Type II MANOVA Tests: Pillai test statistic
#>       Df test stat approx F num Df den Df Pr(>F)
#> group  1  0.085586   1.1933      4     51  0.325

# canonical view
oral.can <- candisc(oral.mod) |> print()
#> 
#> Canonical Discriminant Analysis for group:
#> 
#>     CanRsq Eigenvalue Difference Percent Cumulative
#> 1 0.085586   0.093596                100        100
#> 
#> Test of H0: The canonical correlations in the 
#> current row and all that follow are zero
#> 
#>   LR test stat approx F numDF denDF Pr(> F)
#> 1      0.91441   1.1933     4    51   0.325
summary(oral.can)
#> 
#> Canonical Discriminant Analysis for group:
#> 
#>     CanRsq Eigenvalue Difference Percent Cumulative
#> 1 0.085586   0.093596                100        100
#> 
#> Class means:
#> 
#> [1]  0.30042 -0.30042
#> 
#>  std coefficients:
#>     listen      speak       read      write 
#>  0.5689397 -0.7403339 -0.0060816 -0.7079937 

# reflect the structure & scores to make them positive
oral.can$structure[, "Can1"] <- -1 * oral.can$structure[, "Can1"]
oral.can$scores[, "Can1"]    <- -1 * oral.can$scores[, "Can1"]

plot(oral.can, var.lwd=2)