Description Usage Format Details Source Examples
This is data form mardia, Kent and Bibby on 88 students who took examinations in 5 subjects. Some where with open book and other with closed book.
1 |
A data frame with 88 observations on the following 5 variables.
mechanics, closed book note
vectors, closed book note
algebra, open book note
analysis, open book note
statistics, open book note
The book uses this for bootstrap in principal component analysis.
Efron, B. and Tibshirani, R. (1993) An Introduction to the Bootstrap. Chapman and Hall, New York, London.
1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 | str(scor)
if(interactive())par(ask=TRUE)
plot(scor)
# The parameter of interest (theta) is the fraction of variance explained
# by the first principal component.
# For principal components analysis svd is better numerically than
# eigen-decomposistion, but for bootstrapping the latter is _much_ faster.
theta <- function(ind) {
vals <- eigen(var(scor[ind,]), symmetric=TRUE, only.values=TRUE)$values
vals[1] / sum(vals) }
scor.boot <- bootstrap(1:88, 500, theta)
sd(scor.boot$thetastar) # bootstrap standard error
hist(scor.boot$thetastar)
abline(v=theta(1:88), col="red2")
abline(v=mean(scor.boot$thetastar), col="blue")
|
Add the following code to your website.
For more information on customizing the embed code, read Embedding Snippets.