模糊聚类R语言实现代码

    首次在这上面写Blog，格式什么的先这样了。

    对于模糊聚类均值的公式及其推到，大致如下：

                                             

    主要代码参见下面：（其中使用kmeans作比较。然后通过svm分类测验训练）

# 设置伪随机种子
set.seed(100)

# 生产数据样本
simple.data = function (n=200, nclass=2)
{
    require(clusterGeneration)
    require(mvtnorm)
    # Center of Gaussians
    xpos = seq(-nclass*2, nclass*2, length=nclass)
    ypos = runif(nclass, min=-2*nclass, max=2*nclass)
    
    func = function(i,xpos,ypos,n) {
        # Create a random covariance matrix
        cov = genPositiveDefMat(2, covMethod="eigen",
                            rangeVar=c(1, 10), lambdaLow=1, ratioLambda=10)
        # 保存随机数据
        data = rmvnorm(n=n, mean=c(xpos[i], ypos[i]), sigma=cov$Sigma)
        # 保存每一次的结果
        list(means=cbind(xpos[i], ypos[i]), covars=cov$Sigma, data=data,class=rep(i*1.0, n))
    }
    # do call 合并列表 为 数据框
    strL=do.call(rbind,lapply(1:nclass,func,xpos,ypos,n))
    data=list()
    data$means=do.call(rbind,strL[,1])
    data$covars =