## A quantitative example with the famous geyser data set
data(geyser)
## add 10 missing values
x = geyser;
x[round(runif(5,1,nrow(geyser))), 1] <- NA
x[round(runif(5,1,nrow(geyser))), 2] <- NA
## use graphics functions
set.seed(2)
model <- clusterGamma( data=x, nbCluster=2:3
, models="gamma_pk_ajk_bjk"
, strategy = clusterFastStrategy())
if (FALSE) {
plot(model)
}
## get summary
summary(model)
## print model
if (FALSE) {
print(model)
}
## get estimated missing values
missingValues(model)
Run the code above in your browser using DataLab