http://blog.sina.com.cn/s/blog_676069b10101k1lz.html

该博主对这个问题，自己编了程序，建议参考

FCM <- function(x, K, mybeta = 2, nstart = 1, iter_max = 100, eps = 1e-06) {
  ## FCM


  ## INPUTS
  ##   x: input matrix n*d, n  d-dim samples
  ##   K: number of desired clusters
  ##   Optional :
  ##       mybeta : beta, exponent for m (defaut 2).
  ##       nstart:  how many random sets should be chosen(defaut 1)
  ##       iter_max : The maximum number of iterations allowed. (default 100)
  
  ##      
  ## OUTPUTS
  ##   u: The fuzzy membership matrix = maxtrix of size n*K;
  ##   g: matrix of size K*d of the centers of the clusters
  ##   J: objective function
  ##   histJ: all the objective function values in the iter process


  ## modified time: 2015-02-07


    FCM_onetime <- function(x, init_centers, mybeta = 2, iter_max = 100, eps = 1e-06) {
        n = dim(x)[1]
        d = dim(x)[2]
        g = init_centers
        K = dim(g)[1]
        histJ = c()
        pasfini = 1
        Jold = Inf
        D = matrix(0, n, K)
        for (j in 1:K) {
            D[, j] = rowSums(sweep(x, 2, g[j, ], "-")^2)
        }
        iter = 1
        J_old = Inf
        while (pasfini) {
            s = (1/(D + eps))^(1/(mybeta - 1))
            u = s/(s %*% matrix(1, K, K))
            t1 = t(u^mybeta) %*% x
            t2 = t(u^mybeta) %*% matrix(1, n, d)
            V = t1/t2
            g = V
            D = matrix(0, n, K)
            for (j in 1:K) {
                D[, j] = rowSums(sweep(x, 2, g[j, ], "-")^2)
            }
            J = sum(u^mybeta * D)
            pasfini = abs(J - Jold) > 0.001 && (iter < iter_max)
            Jold = J
            histJ = c(histJ, J)
            iter = iter + 1
        }
        cluster_id = apply(u, 1, which.max)
        re = list(u, J, histJ, g, cluster_id)
        names(re) = c("u", "J", "histJ", "g", "cluster_id")
        return(re)
    }
    x = as.matrix(x)
    seeds = 1:nrow(x)
    id = sample(seeds, K)
    g = as.matrix(x[id, ])
    re_best = FCM_onetime(x = x, init_centers = g, mybeta = mybeta, iter_max = iter_max, eps = eps)
    if (nstart > 1) {
        minJ = 0
        i = 2
        while (i <= nstart) {
            init_centers_id = sample(seeds, K)
            init_centers = as.matrix(x[init_centers_id, ])
            run = FCM_onetime(x, init_centers = init_centers, mybeta = mybeta, iter_max = iter_max)
            if (run$J <= re_best$J) {
                re_best = run
            }
            i = i + 1
        }
    }
    return(re_best)
}

这是一个博主自己写的小程序

恩恩，谢谢，那你会不会最优化，就是说我现在有20个自变量和1个因变量（汽油收率），我已经用svm方法做了回归预测，用前20个变量预测了最后一个汽油收率，现在我想知道的是当汽油收率最大的时候其他变量的取值范围，这个怎么做呢？还有一个问题就是：我有一个excel文本，其中有两列是文字性的，我想知道他们两个之间的关系，就是文本挖掘之类的，这个你会做吗？

cluster包里面有做FCM的函数