library(Rcoclust);
#load data
data(data_news4);
envrdata=get_envrdata(A_ijx,lbs,name,0);
#retrieve matrix size and number of classes
n=envrdata$n;
d=envrdata$d;
g=length(unique(envrdata$lbs));
#coclus
bestresu=NULL;
for (m in 1:100) {
zi_coclus=as.integer( sample(x = 1:g,size = n,replace = TRUE) );
wj_coclus=as.integer( sample(x = 1:g,size = d,replace = TRUE) );
resu=Rcoclust::cc_coclus(g,envrdata,zi_coclus,wj_coclus,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