#TRIMOWA ALGORITHM FOR ONE SIZE:
set.seed(1900)
rand <- sample(1:600,20)
dataComp <- sampleSpanishSurvey[rand, c(2, 3, 5)]
numVar <- dim(dataComp)[2]
orness <- 0.7
weightsTrimowa <- weightsMixtureUB(orness, numVar)
numClust <- 3 ; alpha <- 0.01 ; niter <- 5
ahVect <- c(28, 25, 25) ; algSteps <- 7
set.seed(1900)
Trimowa <- trimowa(dataComp, weightsTrimowa, numClust, alpha, niter,
algSteps, ahVect, verbose = FALSE)
prototypes_Trimowa <- anthrCases("anthropometry", "trimowa",
Trimowa, oneSize = TRUE)
Run the code above in your browser using DataLab