Nothing
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)
Any scripts or data that you put into this service are public.
Add the following code to your website.
For more information on customizing the embed code, read Embedding Snippets.