percentilsArchetypoid {Anthropometry} | R Documentation |
Helper function for computing percentiles of a certain archetypoid
Description
This helper function computes the percentiles of an archetypoid for a given variable. Once these percentile values have been calculated, they can be represented by means of a barplot.
Usage
percentilsArchetypoid(column,indiv,data,digits)
Arguments
column |
Numeric variable (column of a data frame). |
indiv |
A certain archetypoid. |
data |
Data frame that contains the columns and archetypoids to be analyzed. |
digits |
Argument of the |
Value
Numerical vector with the percentile values of an archetypoid.
Author(s)
Guillermo Vinue
References
Vinue, G., Epifanio, I., and Alemany, S., (2015). Archetypoids: a new approach to define representative archetypal data, Computational Statistics and Data Analysis 87, 102–115.
Epifanio, I., Vinue, G., and Alemany, S., (2013). Archetypal analysis: contributions for estimating boundary cases in multivariate accommodation problem, Computers & Industrial Engineering 64, 757–765.
See Also
Examples
#COCKPIT DESIGN PROBLEM:
#As a toy example, only the first 25 individuals are used.
USAFSurvey_First25 <- USAFSurvey[1:25, ]
#Variable selection:
variabl_sel <- c(48, 40, 39, 33, 34, 36)
#Changing to inches:
USAFSurvey_First25_inch <- USAFSurvey_First25[,variabl_sel] / (10 * 2.54)
#Data preprocessing:
USAFSurvey_preproc <- preprocessing(USAFSurvey_First25_inch, TRUE, 0.95, TRUE)
#For reproducing results, seed for randomness:
#suppressWarnings(RNGversion("3.5.0"))
#set.seed(2010)
#Run archetype algorithm repeatedly from 1 to numArch archetypes:
#This is a toy example. In other situation, choose numArch=10 and numRep=20.
numArch <- 5 ; numRep <- 2
lass <- stepArchetypesRawData(data = USAFSurvey_preproc$data, numArch = 1:numArch,
numRep = numRep, verbose = FALSE)
#To understand the warning messages, see the vignette of the
#archetypes package.
#screeplot(lass)
#Three archetypoids:
numArchoid <- 3
res_ns <- archetypoids(numArchoid, USAFSurvey_preproc$data, huge = 200, step = FALSE,
ArchObj = lass, nearest = "cand_ns" , sequ = TRUE)
percentilsArchetypoid(1, res_ns$archet[1], USAFSurvey_preproc$data, 0)