library(Rcoclust);
#load data
data(data_news4);
envrdata=get_envrdata(A_ijx,lbs,name,1);
#retrieve matrix size and number of classes
n=envrdata$n;
d=envrdata$d;
g=length(unique(envrdata$lbs));
#ddkm
bestresu=NULL;
for (m in 1:100) {
zi_ddkm=as.integer( sample(x = 1:g-1,size = n,replace = TRUE) );
wj_ddkm=as.integer( sample(x = 1:g-1,size = d,replace = TRUE) );
resu=Rcoclust::cc_ddkm(g,envrdata,zi_ddkm,wj_ddkm,-1,3,80,0);
if (m==1) bestresu=resu;
if (m>1) {
if (resu$obj[length(resu$obj)]<bestresu$obj[length(bestresu$obj)])
bestresu=resu;
}
}
print(table(envrdata$lbs,bestresu$zi));
Run the code above in your browser using DataLab