#Loading the data to apply the trimowa algorithm:
dataDef <- dataDemo
num.variables <- dim(dataDef)[2]
bust <- dataDef$bust
orness <- 0.7
w <- WeightsMixtureUB(orness,num.variables)
bustCirc_4 <- seq(74,102,4) ; bustCirc_6 <- seq(107,131,6) ; bustCirc <- c(bustCirc_4,bustCirc_6)
nsizes <- length(bustCirc)
K <- 3 ; alpha <- 0.01 ; niter <- 10 ; Ksteps <- 7
ahVect <- c(23,28,20,25,25)
res_trimowa <- list()
for (i in 1 : (nsizes-1)){
data = dataDef[(bust >= bustCirc[i]) & (bust < bustCirc[i + 1]), ]
res_trimowa[[i]] <- trimowa(data,w,K,alpha,niter,Ksteps,ahVect=ahVect)
}
Run the code above in your browser using DataLab