R/Fuzzy.kmeans.sim.noise.R

In mosclust: Model Order Selection for Clustering

"Fuzzy.kmeans.sim.noise" <-
function(X, c = 2, nnoisy=100, perc=0.5, s = sFM, distance="euclidean", hmethod=NULL) {

 n <- ncol(X);
 sim.vector <- numeric(nnoisy);
  
 for (i in 1:nnoisy)  {
   X1 <- perturb.by.noise(X, perc);	 	 
	 if (distance == "euclidean")
		  d <- stats::dist (t(X1))
	 else 
			d <- stats::as.dist(1 - stats::cor(X1));
	 r <- cluster::fanny(d, c, maxit = 1000);		
	 cl1 <- clusterv::Transform.vector.to.list(r$clustering);
	 M1 <- Do.boolean.membership.matrix(cl1, n, 1:n);
	 
	 X2 <- perturb.by.noise(X, perc);	 	 
	 if (distance == "euclidean")
		  d <- stats::dist (t(X2))
	 else 
			d <- stats::as.dist(1 - stats::cor(X2));
	 r <- cluster::fanny(d, c, maxit = 1000);		
	 cl2 <- clusterv::Transform.vector.to.list(r$clustering);
	 M2 <- Do.boolean.membership.matrix(cl2, n, 1:n);
	
	 sim.vector[i] <- s(M1,M2);
 }
 return(sim.vector);
}