dataTrimowa <- sampleSpanishSurvey
numVar <- dim(dataTrimowa)[2]
bust <- dataTrimowa$bust
bustSizes <- bustSizesStandard(seq(74, 102, 4), seq(107, 131, 6))
orness <- 0.7
weightsTrimowa <- weightsMixtureUB(orness,numVar)
numClust <- 3 ; alpha <- 0.01 ; niter <- 10 ; algSteps <- 7
ah <- c(23, 28, 20, 25, 25)
set.seed(2014)
res_trimowa <- list()
for (i in 1 : (bustSizes$nsizes - 1)){
data = dataTrimowa[(bust >= bustSizes$bustCirc[i]) & (bust < bustSizes$bustCirc[i + 1]), ]
res_trimowa[[i]] <- trimowa(data, weightsTrimowa, numClust, alpha, niter,
algSteps,ah,verbose=FALSE)
}
Run the code above in your browser using DataLab