library(dplyr, warn.conflicts = FALSE)
nested_pca <- alta_lake_geochem %>%
nested_data(
qualifiers = c(depth, zone),
key = param,
value = value,
trans = scale
) %>%
nested_prcomp()
# get variance info
nested_pca %>% unnested_data(variance)
# get loadings info
nested_pca %>% unnested_data(loadings)
# scores, requalified
nested_pca %>% unnested_data(c(qualifiers, scores))
Run the code above in your browser using DataLab