if (FALSE) {
data(bearings)
# Split dataset into train (75
set.seed(3)
Bearings <- split(p = 0.75, Dataset = bearings, class = 14)
# Estimate number of components, component weights and component
# parameters for train subsets.
bearingsest <- REBMIX(model = "REBMVNORM",
Dataset = a.train(Bearings),
Preprocessing = "histogram",
cmax = 15,
Criterion = "BIC")
# Classification.
bearingscla <- RCLSMIX(model = "RCLSMVNORM",
x = list(bearingsest),
Dataset = a.test(Bearings),
Zt = a.Zt(Bearings))
bearingscla
summary(bearingscla)
}
Run the code above in your browser using DataLab