FindElbow
1 |
1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 | ## Not run:
tmp <- FindElbow(c(0.9, 1.1, 1.1, 1.9, 2.5, 2.8, 4.9, 8.5),
plot = TRUE) # wandering
tmp <- FindElbow(c(0.9, 1.0, 1.2, 1.3, 1.5, 1.5, 10.0, 22.0),
plot = TRUE) # late rise
tmp <- FindElbow(c(2, 4, 6, 8, 10, 12, 14, 16)^2,
plot = TRUE) # gradual, no obvious break
# Not the usual way to choose the number of PCs:
library("chemometrics")
data(glass)
pca <- prcomp(glass)
eigensum <- sum(pca$sdev * pca$sdev)
vv <- 100 * (pca$sdev * pca$sdev/eigensum)
cs <- cumsum(vv)
tmp <- FindElbow(vv, plot = TRUE)
tmp <- FindElbow(cs, plot = TRUE)
## End(Not run)
|
Add the following code to your website.
For more information on customizing the embed code, read Embedding Snippets.