99999久久久久久亚洲,欧美人与禽猛交狂配,高清日韩av在线影院,一个人在线高清免费观看,啦啦啦在线视频免费观看www

熱線電話:13121318867

登錄
首頁精彩閱讀k-means聚類”—數(shù)據(jù)分析、數(shù)據(jù)挖掘
k-means聚類”—數(shù)據(jù)分析、數(shù)據(jù)挖掘
2016-05-04
收藏

k-means聚類”—數(shù)據(jù)分析、數(shù)據(jù)挖掘

一、概要

    分類作為一種監(jiān)督學習方法,要求必須事先明確知道各個類別的信息,并且斷言所有待分類項都有一個類別與之對應。但是很多時候上述條件得不到滿足,尤其是在處理海量數(shù)據(jù)的時候,如果通過預處理使得數(shù)據(jù)滿足分類算法的要求,則代價非常大,這時候可以考慮使用聚類算法。聚類屬于無監(jiān)督學習,相比于分類,聚類不依賴預定義的類和類標號的訓練實例。本文介紹一種常見的聚類算法——k 均值和k 中心點聚類,最后會舉一個實例:應用聚類方法試圖解決一個在體育界大家頗具爭議的問題——中國男足近幾年在亞洲到底處于幾流水平。

二、聚類問題
    所謂聚類問題,就是給定一個元素集合D,其中每個元素具有n 個可觀察屬性,使用某種算法將D 劃分成k 個子集,要求每個子集內部的元素之間相異度盡可能低,而不同子集的元素相異度盡可能高。其中每個子集叫做一個簇。
    與分類不同,分類是示例式學習,要求分類前明確各個類別,并斷言每個元素映射到一個類別,而聚類是觀察式學習,在聚類前可以不知道類別甚至不給定類別數(shù)量,是無監(jiān)督學習的一種。目前聚類廣泛應用于統(tǒng)計學、生物學、數(shù)據(jù)庫技術和市場營銷等領域,相應的算法也非常的多。本文僅介紹一種最簡單的聚類算法——k 均值(k-means)算法。

三、概念介紹

區(qū)分兩個概念:

hard clustering:一個文檔要么屬于類w,要么不屬于類w,即文檔對確定的類w是二值的1或0。

soft clustering:一個文檔可以屬于類w1,同時也可以屬于w2,而且文檔屬于一個類的值不是0或1,可以是0.3這樣的小數(shù)。

K-Means就是一種hard clustering,所謂K-means里的K就是我們要事先指定分類的個數(shù),即K個。

k-means算法的流程如下:

1)從N個文檔隨機選取K個文檔作為初始質心

2)對剩余的每個文檔測量其到每個質心的距離,并把它歸到最近的質心的類

3)重新計算已經得到的各個類的質心

4)迭代2~3步直至滿足既定的條件,算法結束

在K-means算法里所有的文檔都必須向量化,n個文檔的質心可以認為是這n個向量的中心,計算方法如下:

                   

這里加入一個方差RSS的概念:

RSSk的值是類k中每個文檔到質心的距離,RSS是所有k個類的RSS值的和。

算法結束條件:

1)給定一個迭代次數(shù),達到這個次數(shù)就停止,這好像不是一個好建議。

2)k個質心應該達到收斂,即第n次計算出的n個質心在第n+1次迭代時候位置不變。

3)n個文檔達到收斂,即第n次計算出的n個文檔分類和在第n+1次迭代時候文檔分類結果相同。

4)RSS值小于一個閥值,實際中往往把這個條件結合條件1使用

回過頭用RSS討論質心的計算方法是否合理

   

為了取得RSS的極小值,RSS對質心求偏導數(shù)應該為0,所以得到質心

          

可見,這個質心的選擇是合乎數(shù)學原理的。

 

K-means方法的缺點是聚類結果依賴于初始選擇的幾個質點位置,看下面這個例子:

如果使用2-means方法,初始選擇d2和d5那么得到的聚類結果就是{d1,d2,d3}{d4,d5,d6},這不是一個合理的聚類結果

解決這種初始種子問題的方案:

1)去處一些游離在外層的文檔后再選擇

2)多選一些種子,取結果好的(RSS?。┑腒個類繼續(xù)算法

3)用層次聚類的方法選擇種子。我認為這不是一個合適的方法,因為對初始N個文檔進行層次聚類代價非常高。

以上的討論都是基于K是已知的,但是我們怎么能從隨機的文檔集合中選擇這個k值呢?

我們可以對k去1~N分別執(zhí)行k-means,得到RSS關于K的函數(shù)下圖:

當RSS由顯著下降到不是那么顯著下降的K值就可以作為最終的K,如圖可以選擇4或9。

四、算法及示例

k 均值算法的計算過程非常直觀:
1、從D 中隨機取k 個元素,作為k 個簇的各自的中心。
2、分別計算剩下的元素到k 個簇中心的相異度,將這些元素分別劃歸到相異度最低的簇。
3、根據(jù)聚類結果,重新計算k 個簇各自的中心,計算方法是取簇中所有元素各自維度的算術平均數(shù)。

4、將D 中全部元素按照新的中心重新聚類。

5、重復第4 步,直到聚類結果不再變化。

6、將結果輸出。

由于算法比較直觀,沒有什么可以過多講解的。下面,我們來看看k-means 算法一個有趣的應用示例:中國男足近幾年到底在亞洲處于幾流水平?

今年中國男足可算是杯具到家了,幾乎到了過街老鼠人人喊打的地步。對于目前中國男足在亞洲的地位,各方也是各執(zhí)一詞,有人說中國男足亞洲二流,有人說三流,還有人說根本不入流,更有人說其實不比日韓差多少,是亞洲一流。既然爭論不能解決問題,我們就讓數(shù)據(jù)告訴我們結果吧。

下圖是采集的亞洲15 只球隊在2005 年-2010 年間大型杯賽的戰(zhàn)績(由于澳大利亞是后來加入亞足聯(lián)的,所以這里沒有收錄)。

其中包括兩次世界杯和一次亞洲杯。我提前對數(shù)據(jù)做了如下預處理:對于世界杯,進入決賽圈則取其最終排名,沒有進入決賽圈的,打入預選賽十強賽賦予40,預選賽小組未出線的賦予50。對于亞洲杯,前四名取其排名,八強賦予5,十六強賦予9,預選賽沒出現(xiàn)的賦予17。這樣做是為了使得所有數(shù)據(jù)變?yōu)闃肆?,便于后續(xù)聚類。

下面先對數(shù)據(jù)進行[0,1]規(guī)格化,下面是規(guī)格化后的數(shù)據(jù):

其中包括兩次世界杯和一次亞洲杯。我提前對數(shù)據(jù)做了如下預處理:對于世界杯,進入決賽圈則取其最終排名,沒有進入決賽圈的,打入預選賽十強賽賦予40,預選賽小組未出線的賦予50。對于亞洲杯,前四名取其排名,八強賦予5,十六強賦予9,預選賽沒出現(xiàn)的賦予17。這樣做是為了使得所有數(shù)據(jù)變?yōu)闃肆?,便于后續(xù)聚類。

下面先對數(shù)據(jù)進行[0,1]規(guī)格化,下面是規(guī)格化后的數(shù)據(jù):

從做到右依次表示各支球隊到當前中心點的歐氏距離,將每支球隊分到最近的簇,可對各支球隊做如下聚類:

中國C,日本A,韓國A,伊朗A,沙特A,伊拉克C,卡塔爾C,阿聯(lián)酋C,烏茲別克斯坦B,泰國C,越南C,阿曼C,巴林B,朝鮮B,印尼C。

第一次聚類結果:

A:日本,韓國,伊朗,沙特;

B:烏茲別克斯坦,巴林,朝鮮;

C:中國,伊拉克,卡塔爾,阿聯(lián)酋,泰國,越南,阿曼,印尼。

下面根據(jù)第一次聚類結果,調整各個簇的中心點。

A 簇的新中心點為: {(0.3+0+0.24+0.3)/4=0.21,(0+0.15+0.76+0.76)/4=0.4175,(0.19+0.13+0.25+0.06)/4=0.1575} = {0.21, 0.4175, 0.1575}

用同樣的方法計算得到B 和C 簇的新中心點分別為{0.7, 0.7333, 0.4167},{1, 0.94,0.40625}。

用調整后的中心點再次進行聚類,得到:

第二次迭代后的結果為:

中國C,日本A,韓國A,伊朗A,沙特A,伊拉克C,卡塔爾C,阿聯(lián)酋C,烏茲別克斯坦B,泰國C,越南C,阿曼C,巴林B,朝鮮B,印尼C。

結果無變化,說明結果已收斂,于是給出最終聚類結果:

亞洲一流:日本,韓國,伊朗,沙特

亞洲二流:烏茲別克斯坦,巴林,朝鮮

亞洲三流:中國,伊拉克,卡塔爾,阿聯(lián)酋,泰國,越南,阿曼,印尼

看來數(shù)據(jù)告訴我們,說國足近幾年處在亞洲三流水平真的是沒有冤枉他們,至少從國際杯賽戰(zhàn)績是這樣的。

其實上面的分析數(shù)據(jù)不僅告訴了我們聚類信息,還提供了一些其它有趣的信息,例如從中可以定量分析出各個球隊之間的差距,例如,在亞洲一流隊伍中,日本與沙特水平最接近,而伊朗則相距他們較遠,這也和近幾年伊朗沒落的實際相符。另外,烏茲別克斯坦和巴林雖然沒有打進近兩屆世界杯,不過憑借預算賽和亞洲杯上的出色表現(xiàn)占據(jù)B 組一席之地,而朝鮮由于打入了2010 世界杯決賽圈而有幸進入B 組,可是同樣奇跡般奪得2007年亞洲杯的伊拉克卻被分在三流,看來亞洲杯冠軍的分量還不如打進世界杯決賽圈重啊。其它有趣的信息,有興趣的朋友可以進一步挖掘。


數(shù)據(jù)分析咨詢請掃描二維碼

若不方便掃碼,搜微信號:CDAshujufenxi

數(shù)據(jù)分析師考試動態(tài)
數(shù)據(jù)分析師資訊
更多

OK
客服在線
立即咨詢
客服在線
立即咨詢
') } function initGt() { var handler = function (captchaObj) { captchaObj.appendTo('#captcha'); captchaObj.onReady(function () { $("#wait").hide(); }).onSuccess(function(){ $('.getcheckcode').removeClass('dis'); $('.getcheckcode').trigger('click'); }); window.captchaObj = captchaObj; }; $('#captcha').show(); $.ajax({ url: "/login/gtstart?t=" + (new Date()).getTime(), // 加隨機數(shù)防止緩存 type: "get", dataType: "json", success: function (data) { $('#text').hide(); $('#wait').show(); // 調用 initGeetest 進行初始化 // 參數(shù)1:配置參數(shù) // 參數(shù)2:回調,回調的第一個參數(shù)驗證碼對象,之后可以使用它調用相應的接口 initGeetest({ // 以下 4 個配置參數(shù)為必須,不能缺少 gt: data.gt, challenge: data.challenge, offline: !data.success, // 表示用戶后臺檢測極驗服務器是否宕機 new_captcha: data.new_captcha, // 用于宕機時表示是新驗證碼的宕機 product: "float", // 產品形式,包括:float,popup width: "280px", https: true // 更多配置參數(shù)說明請參見:http://docs.geetest.com/install/client/web-front/ }, handler); } }); } function codeCutdown() { if(_wait == 0){ //倒計時完成 $(".getcheckcode").removeClass('dis').html("重新獲取"); }else{ $(".getcheckcode").addClass('dis').html("重新獲取("+_wait+"s)"); _wait--; setTimeout(function () { codeCutdown(); },1000); } } function inputValidate(ele,telInput) { var oInput = ele; var inputVal = oInput.val(); var oType = ele.attr('data-type'); var oEtag = $('#etag').val(); var oErr = oInput.closest('.form_box').next('.err_txt'); var empTxt = '請輸入'+oInput.attr('placeholder')+'!'; var errTxt = '請輸入正確的'+oInput.attr('placeholder')+'!'; var pattern; if(inputVal==""){ if(!telInput){ errFun(oErr,empTxt); } return false; }else { switch (oType){ case 'login_mobile': pattern = /^1[3456789]\d{9}$/; if(inputVal.length==11) { $.ajax({ url: '/login/checkmobile', type: "post", dataType: "json", data: { mobile: inputVal, etag: oEtag, page_ur: window.location.href, page_referer: document.referrer }, success: function (data) { } }); } break; case 'login_yzm': pattern = /^\d{6}$/; break; } if(oType=='login_mobile'){ } if(!!validateFun(pattern,inputVal)){ errFun(oErr,'') if(telInput){ $('.getcheckcode').removeClass('dis'); } }else { if(!telInput) { errFun(oErr, errTxt); }else { $('.getcheckcode').addClass('dis'); } return false; } } return true; } function errFun(obj,msg) { obj.html(msg); if(msg==''){ $('.login_submit').removeClass('dis'); }else { $('.login_submit').addClass('dis'); } } function validateFun(pat,val) { return pat.test(val); }