中国mooc北京理工大学机器学习第一周（一）

    import numpy as np
    from sklearn.cluster import KMeans
       
    def loadData(filePath):#def一个读取数据的loadData
        fr = open(filePath,‘r+‘)
        lines = fr.readlines()
        retData = []
        retCityName = []
        for line in lines:
            items = line.strip().split(",")
            retCityName.append(items[0])
            retData.append([float(items[i]) for i in range(1,len(items))])
        return retData,retCityName
     　　　　　　　　　　　　　　　　　　　　　　　　　　　　　　　　　　　　　　　　　　　　　　　　　　　　　　　　　　
         　　　　　　　　　　　　　　　　　　　　　　　　　　　　　　　　　　　　　　　　　　　　　　　　　　　　 
    if __name__ == ‘__main__‘:　　　　　　　　　　　　　　　　　　　　　　　　　　　　　　　　　　　　　　　　
        data,cityName = loadData(‘city.txt‘)
        km = KMeans(n_clusters=4)#构造聚类容器，分为四类
        label = km.fit_predict(data)#进行聚类，贴标签，label=4
        expenses = np.sum(km.cluster_centers_,axis=1)#expenses是每一个聚类容器的sum
        CityCluster = [[],[],[],[]]#这里的类别大于等于n_clusters
        for i in range(len(cityName)):
            CityCluster[label[i]].append(cityName[i])
        for i in range(len(CityCluster)):
            print("Expenses:%.2f" % expenses[i])
            print(CityCluster[i])

 1   sklearn.cluster.KMeans(
 2 　　  n_clusters=8,#你想分类的种类
 3       init=‘k-means++‘,#初始簇中心的获取方法    
 4       n_init=10,#获取初始簇中心的更迭次数，为了弥补初始质心的影响，算法默认会初始10个质心，实现算法，然后返回最好的结果。 
 5       max_iter=300,#最大迭代次数 
 6       tol=0.0001,#容忍度，即kmeans运行准则收敛的条件 
 7       precompute_distances=‘auto‘, #是否需要提前计算距离
 8       verbose=0, 
 9       random_state=None, 
10       copy_x=True, 
11      n_jobs=1, #并行设置
12      algorithm=‘auto‘ #kmeans的实现算法
13      )

    import numpy as np
    import sklearn.cluster as skc
    from sklearn import metrics
    import matplotlib.pyplot as plt  
    mac2id=dict()#字典
    onlinetimes=[]#数组
    f=open(‘TestData.txt‘,encoding=‘utf-8‘)
    for line in f:
        mac=line.split(‘,‘)[2]  #A417314EEA7B
        onlinetime=int(line.split(‘,‘)[6]) #15
        starttime=int(line.split(‘,‘)[4].split(‘ ‘)[1].split(‘:‘)[0])#22
        if mac not in mac2id:
            mac2id[mac]=len(onlinetimes)#如果没有出现过，直接append
            onlinetimes.append((starttime,onlinetime))
        else:
            onlinetimes[mac2id[mac]]=[(starttime,onlinetime)]#出现过会替代（不太理解这，可能是数据的特点吧）
    real_X=np.array(onlinetimes).reshape((-1,2))#改成两列（-1的含义是任意行）

 　　X=real_X[:,0:1]#按照原形式读取第一列
    print(X)
     
    db=skc.DBSCAN(eps=0.01,min_samples=20).fit(X)#进行聚类，eps：半径，min_samples:数据的个数
    labels = db.labels_#贴标签与取标签
     
    print(‘Labels:‘)
    print(labels)
    raito=len(labels[labels[:] == -1]) / len(labels)#-1为噪声点
    print(‘Noise raito:‘,format(raito, ‘.2%‘))
     
    n_clusters_ = len(set(labels)) - (1 if -1 in labels else 0)   #set的作用是去重复，然后减去噪声
     
    print(‘Estimated number of clusters: %d‘ % n_clusters_)
    print("Silhouette Coefficient: %0.3f"% metrics.silhouette_score(X, labels))
     
    for i in range(n_clusters_):
        print(‘Cluster ‘,i,‘:‘)
        print(list(X[labels == i].flatten()))#flatten方法：高维array和matrix转化成1维。list不可用。
         
    plt.hist(X,20)
    plt.show()