oask <- devAskNewPage(dev.interactive(orNone=TRUE))
bp <- bpca(rock,
var.rb=TRUE,
var.rd=TRUE)
summary(bp)
# The total variance explained is satisfactory (>= .80)!
plot(bp)
# A more accurate diagnostic
bp$var.rd
# It is possible to observe that the variable 'perm'
# has not a good representation (bpca.2d)!
# Observed correlations:
cor(rock)
# Projected correlations:
bp$var.rb
# Additional diagnostic
plot(qbpca(rock,
bp))
# This variable remains as important in a dimension not contemplated
# by the biplot reduction (PC3):
bp$eigenvectors
bp1 <- bpca(rock,
d=3:4)
summary(bp1)
plot(bp1)
# The recommendation, knowing that this variable has a poor
# representation is:
# 1- Avoid to discute it;
# 2- Consider to incorporate the information with a bpca.3d
bp2 <- bpca(rock,
d=1:3,
var.rb=TRUE,
var.rd=TRUE)
summary(bp2)
plot(bp2) # Static
plot(bp2,
rgl.use=TRUE) # Dynamic
bp2$var.rd # Nice!
# Additional diagnostic
plot(qbpca(rock,
bp2))
devAskNewPage(oask)
Add the following code to your website.
For more information on customizing the embed code, read Embedding Snippets.