{"id":463,"date":"2011-06-20T20:08:00","date_gmt":"2011-06-20T12:08:00","guid":{"rendered":"http:\/\/note.systw.net\/note\/?p=463"},"modified":"2023-11-02T20:09:57","modified_gmt":"2023-11-02T12:09:57","slug":"k-means","status":"publish","type":"post","link":"https:\/\/systw.net\/note\/archives\/463","title":{"rendered":"K-means"},"content":{"rendered":"\n

k-means<\/strong>
\u512a\u9ede:\u7c21\u55ae\u6709\u6548\u7387,\u4efb\u4e00\u6f14\u7b97\u6cd5\u6240\u5f97\u4e4b\u7d50\u679c,\u90fd\u53ef\u900f\u904e\u6b64\u65b9\u6cd5\u9032\u884c\u6539\u5584
\u7f3a\u9ede:1\u5404\u7fa4\u7684\u8cc7\u6599\u5206\u4f48\u9808\u5448\u5713\u5f62\u5206\u4f48\u4e14\u6709\u985e\u4f3c\u7684\u5927\u5c0f\u6642,\u624d\u80fd\u5f97\u5230\u8f03\u4f73\u7684\u89e3,2\u521d\u59cb\u503c\u6c92\u8a2d\u597d\u5bb9\u6613\u9677\u5165\u5340\u57df\u6700\u4f73\u89e3<\/p>\n\n\n\n

\u8aaa\u660e<\/strong>
\u5728\u5206\u5272\u5f0f\u65b9\u6cd5\u4e2d\uff0c\u6700\u5e38\u88ab\u4f7f\u7528\u4e5f\u6700\u5177\u4ee3\u8868\u6027\u7684\u65b9\u6cd5\u70baK-means \u5206\u7fa4\u6f14\u7b97\u6cd5\uff0c\u8a72\u65b9\u6cd5\u7531MacQueen\u65bc1967 \u6240\u63d0\u51fa\u4f86\u3002k-means \u5206\u7fa4\u6f14\u7b97\u6cd5\u662f\u4e00\u7a2e\u7c21\u55ae\u6613\u61c2\u4e14\u5ee3\u70ba\u4f7f\u7528\u7684\u975e\u76e3\u7763\u5f0f\u5b78\u7fd2\u7684\u6f14\u7b97\u6cd5\u3002K-means \u5206\u7fa4\u6f14\u7b97\u6cd5\u8981\u5148\u6c7a\u5b9a\u5c07\u8cc7\u6599\u5206\u70bax\u500b\u7fa4\u7d44\uff0c\u5176\u4e2d\u6bcf\u4e00\u500b\u7fa4\u7d44\u5fc5\u9808\u81f3\u5c11\u5305\u542b\u4e00\u9805\u8cc7\u6599\uff0c\u800c\u4efb\u4f55\u4e00\u500b\u8cc7\u6599\u90fd\u61c9\u8a72\u5c6c\u65bc\u67d0\u4e00\u500b\u7fa4\u7d44\uff0c\u4f9d\u64da\u8cc7\u6599\u4e4b\u9593\u7684\u8ddd\u96e2\uff0c\u53d6\u5f97\u7fa4\u7d44\u7684\u4e2d\u5fc3\u4f4d\u7f6e\uff0c\u5207\u5272\u51fa\u7fa4\u96c6\u8207\u7fa4\u96c6\u4e4b\u9593\u7684\u754c\u7dda\u3002
 <\/p>\n\n\n\n

step<\/strong>
1\u6c7a\u5b9a\u7fa4\u6578\u91cf
2\u521d\u59cb\u5316,\u6709\u4ee5\u4e0b\u5169\u7a2e\u505a\u6cd5
\u505a\u6cd51,\u5c07\u5404\u8cc7\u6599\u7269\u4ef6\u6307\u6d3e(ex:\u96a8\u6a5f\u6307\u6d3e)\u5230\u5176\u4e2d\u4e00\u7fa4,\u5728\u8a08\u7b97\u5404\u7fa4\u4e2d\u5fc3
\u505a\u6cd52,\u6307\u6d3e(ex:\u96a8\u6a5f\u6307\u6d3e)\u5176\u4e2d\u4e00\u500b\u8cc7\u6599\u7269\u4ef6\u7576\u5176\u4e2d\u4e00\u500b\u7fa4\u4e2d\u5fc3
3\u5206\u914d\u6bcf\u500b\u8cc7\u6599\u7269\u4ef6\u5230\u64c1\u6709\u6700\u77ed\u77e9\u96e2\u7684\u7fa4(\u6bcf\u500b\u7269\u4ef6\u8981\u7b97\u81ea\u5df1\u8207\u6bcf\u500b\u7fa4\u4e2d\u5fc3\u7684\u8ddd\u96e2,\u82e5\u7b97\u51fa\u4f86\u7684\u8ddd\u96e2\u5047\u8a2d\u662f\u548ca\u7fa4\u6700\u8fd1,\u5247\u8a72\u7269\u4ef6\u5c6c\u65bca\u7fa4)
4\u91cd\u7b97\u6240\u6709\u7fa4\u4e2d\u5fc3
5\u91cd\u89863,4\u76f4\u5230\u6240\u6709\u7fa4\u4e2d\u5fc3\u4e0d\u5728\u6709\u4efb\u4f55\u6539\u8b8a
ps:
2\u7dad\u8ddd\u96e2\u516c\u5f0f sqrt( (1x-2x)^2 + (1y-2y)^2 )
2\u7dadN\u9ede\u4e2d\u5fc3\u516c\u5f0f (x1+x2+x3+…+xN)\/N,(y1+y2+…+yN)\/N
ps:
\u8a2d\u9ede\u67091,2,3,…,N
\u7dad\u5ea6\u6709x,y,z
3\u7dad\u8ddd\u96e2\u516c\u5f0f sqrt( (1x-2x)^2 + (1y-2y)^2 + (1z-2z)^2 )
3\u7dadN\u9ede\u4e2d\u5fc3\u516c\u5f0f (x1+x2+x3+…+xN)\/N,(y1+y2+y3+…+yN)\/N,(z1+z2+z3+…+zN)\/N<\/p>\n\n\n\n

ex:<\/p>\n\n\n\n

 <\/td>x1 <\/td>x2 <\/td><\/tr>
a <\/td>5 <\/td>3 <\/td><\/tr>
b <\/td>-1 <\/td>1 <\/td><\/tr>
c <\/td>1 <\/td>-2 <\/td><\/tr>
d <\/td>-3 <\/td>-2 <\/td><\/tr><\/tbody><\/table><\/figure>\n\n\n\n

1
\u8a2d\u5b9a\u7fa4\u6578\u91cf\u70ba2
\u4e26\u96a8\u6a5f\u5206\u5169\u7fa4\u70baab,cd,
2
\u6c7a\u5b9aab\u7684\u7fa4\u4e2d\u5fc3\u548ccd\u7684\u7fa4\u4e2d\u5fc3,\u5982\u4e0b<\/p>\n\n\n\n

 <\/td>x1 <\/td>x2 <\/td><\/tr>
ab <\/td>(5+(-1))\/2=2 <\/td>(3+1)\/2=2 <\/td><\/tr>
cd <\/td>(1+(-3))\/2=-1 <\/td>(-2+(-2))\/2=-2 <\/td><\/tr><\/tbody><\/table><\/figure>\n\n\n\n

3
\u627e\u51fa\u9ede\u8207ab\u6216cd\u90a3\u500b\u6700\u8fd1,\u4ee5\u5206\u914d\u6bcf\u500b\u8cc7\u6599\u7269\u4ef6\u5230\u6700\u9069\u5408\u7684\u7fa4,\u5982\u4e0b
min{(a,ab),(a,cd)}={10,61}=(a,ab)=sqrt(10)
min{(b,ab),(b,cd)}={10,9}=(b,cd)=sqrt(9) \u56e0b,cd\u5c0f\u65bcb,ab,\u6240\u4ee5b\u5f9eab\u7fa4\u6539\u5230cd\u7fa4
min{(c,ab),(c,cd)}={17,4}=(c,cd)=sqrt(4)
min{(d,ab),(d,cd)}={41,4}=(d,cd)=sqrt(4)
ps:(x1,y1),(x2,y2)\u7684\u8ddd\u96e2=sqrt((x1-x2)(x1-x2)+(y1-y2)(y1-y2))
\u6839\u64da\u4ee5\u4e0a\u7d50\u679c,\u5169\u7fa4\u5df2\u8b8a\u6210a,bcd
4
\u91cd\u7b97\u6240\u6709\u7fa4\u4e2d\u5fc3,\u6c7a\u5b9aa\u7684\u7fa4\u4e2d\u5fc3\u548cb(cd)\u7684\u7fa4\u4e2d\u5fc3,\u5982\u4e0b<\/p>\n\n\n\n

 <\/td>x1 <\/td>x2 <\/td><\/tr>
a <\/td>5 <\/td>3 <\/td><\/tr>
b(cd) <\/td>(-1+(-1))\/2=-1 <\/td>(1+(-2))\/2=-0.5  <\/td><\/tr><\/tbody><\/table><\/figure>\n\n\n\n

5
\u627e\u51fa\u9ede\u8207a\u6216b(cd)\u90a3\u500b\u6700\u8fd1,\u4ee5\u5206\u914d\u6bcf\u500b\u8cc7\u6599\u7269\u4ef6\u5230\u6700\u9069\u5408\u7684\u7fa4,\u5982\u4e0b
min{(a,a),(a,b(cd))}={0,52}=(a,a)=sqrt(0)
min{(b,a),(b,b(cd))}={40,4}=(b,b(cd))=sqrt(4)
min{(c,a),(c,b(cd))}={41,5}=(c,b(cd))=sqrt(5)
min{(d,a),(d,b(cd))}={89,5}=(d,b(cd))=sqrt(5)
\u6839\u64da\u4ee5\u4e0a\u7d50\u679c,\u5169\u7fa4\u4ecd\u7136\u662fa,bcd
6
\u56e0\u7d50\u679c\u4e0d\u8b8a,\u6240\u4ee5\u70ba\u6700\u7d42\u5206\u7fa4\u7d50\u679c<\/p>\n\n\n\n


……………………………….
\u4f7f\u7528\u5de5\u5177
minitab<\/p>\n\n\n\n

k-mean\u6b65\u9a5f<\/strong>
1\u5148\u5c07\u8cc7\u6599\u8cbc\u5230worksheet
2\u9ede\u64castat > multivariate > cluster k-mean ,\u6703\u8df3\u51fa\u5c0d\u8a71\u8996\u7a97
2.1\u5c07\u5de6\u6b04\u8981\u7528\u5230\u7684\u6b04\u4f4d\u9078\u53d6\u597d\u5f8c\u6309select,\u5247\u8981\u7528\u5230\u7684\u6b04\u4f4d\u6703\u51fa\u73fe\u5728\u53f3\u6b04variables\u4e2d
2.2\u9ede\u64canumber of cluster,\u4e26\u8f38\u5165\u7fa4\u7684\u6578\u91cf
2.3\u9ede\u64ca\u4e0b\u65b9storage,\u6703\u8df3\u51fa\u5c0d\u8a71\u8996\u7a97,
2.3.1\u5728cluster membership column\u5167\u8f38\u5165\u4e00\u6b04\u4f4d ex:d
2.3.2\u9ede\u64caok\u5f8cstorage\u5c0d\u8a71\u8996\u7a97\u95dc\u9589
2.4\u9ede\u64caok,cluster k-means\u5c0d\u8a71\u8996\u7a97\u95dc\u9589
4\u7d50\u679c\u986f\u793a
ps:
k-mean\u7684\u8f38\u5165\u8cc7\u6599\u5f62\u5f0f\u901a\u5e38\u70ba\u4e00\u7b46\u7b46\u7684\u8cc7\u6599,\u5982\u4e0b(\u53d6\u81eauci iris data)
5.7,4.4,1.5,0.4,Iris-setosa
5.4,3.9,1.3,0.4,Iris-setosa
5.1,3.5,1.4,0.3,Iris-setosa
ex:
\u4ee5\u4e0a\u8ff0\u7bc4\u4f8b\u70ba\u4f8b,\u5247\u7d50\u679c\u70ba
a (5,3) \u7fa41
b (-1,1) \u7fa42
c (1,-2) \u7fa42
d (-3,-2)\u7fa42
\u5c07\u4e0a\u8ff0\u7bc4\u4f8b\u6539\u70ba3\u500b\u7fa4\u70ba\u4f8b,\u5247\u7d50\u679c\u70ba
a (5,3) \u7fa41
b (-1,1) \u7fa42
c (1,-2) \u7fa43
d (-3,-2)\u7fa42<\/p>\n\n\n\n

…………………………..<\/p>\n\n\n\n


kmean for php \u7a0b\u5f0f\u78bc<\/strong>(beta)
ps:\u597d\u50cf\u9084\u6709\u4e00\u4e9b\u554f\u984c<\/p>\n\n\n\n

\/\/kmean for php
\/\/time record
$startime=caclutime();
if($_GET['debug']==1)$kmean['para']['debug']=1; \/\/\u662f\u5426\u555f\u7528debug mode
$kmean['cluster']['count']=3; \/\/\u8981\u5206\u7684\u7fa4\u6578\u91cf
$datapath='data'; \/\/\u8cc7\u6599\u4f86\u6e90
\/\/\u8b80\u53d6\u8cc7\u6599
input_fromfile($kmean['dot'],$datapath);
\/\/\u57f7\u884ckmean
kmean($kmean);
\/\/$endtime
$endtime=caclutime();
printf('query took:%3.5fsec',($endtime-$startime));
\/\/\/\/\/\/\/\/\/\/\/\/\/\/\/\/\/\/\/\/\/\/\/\/\/\/\/\/\/\/\/\/\/\/\/\/\/\/\/\/\/\/\/\/\/\/\/\/\/\/\/\/\/\/\/\/\/\/\/\/\/\/\/\/\/\/\/\/\/\/\/\/\/\/\/\/\/\/\/\/\/\/\/\/\/\/\/\/\/\/\/\/\/\/\/\/\/\/
\/\/\/\/\/\/\/\/\/\/\/\/\/\/\/\/\/\/\/\/\/\/\/\/\/\/\/\/\/\/\/\/\/\/\/\/\/\/\/\/\/\/\/\/\/\/\/\/\/\/\/\/\/\/\/\/\/\/\/\/\/\/\/\/\/\/\/\/\/\/\/\/\/\/\/\/\/\/\/\/\/\/\/\/\/\/\/\/\/\/\/
\/\/kmean main function start
function kmean($kmean){
\/*
paramenter list
$_GET['init']=1,2,NULL
$_GET['debug']=1,NULL
*\/
\/\/\u8a18\u9304\u6240\u6709\u7fa4\u4e2d\u5fc3\u7684\u4f4d\u7f6e
\/\/$kmean[cluster][center][int:clusterid:1~n][int:dimensions:1~m]
\/\/ex:
\/\/$kmean[cluster][center][1][1]=cluster1 center value1
\/\/$kmean[cluster][center][1][2]=cluster1 center value2
\/\/$kmean[cluster][center][1][n]=cluster1 center valuem
$kmean['cluster']['center'][1][1]=0;
$kmean['cluster']['center'][1][2]=0;
\/\/\u8a18\u9304\u8a72\u7fa4\u76ee\u524d\u6709\u7684\u9ede
\/\/$kmean['cluster']['dot'][int:cluster_id:1~m][]
\/\/$kmean[cluster][dot][1][0]=dot_id \u7b2c\u4e00\u500b\u7fa4\u7684\u7b2c\u4e00\u500b\u9ede\u7684id
\/\/$kmean[cluster][dot][1][1]=dot_id \u7b2c\u4e00\u500b\u7fa4\u7684\u7b2c\u4e8c\u500b\u9ede\u7684id
\/\/\u8a18\u9304\u6240\u6709\u9ede\u7684\u4f4d\u7f6e\u53ca\u540d\u7a31
\/\/$kmean[dot][int:dot_id:1~n][str:name:0]
\/\/$kmean[dot][int:dot_id:1~n][int:dimensions:1~m]
\/\/ex:
\/\/$kmean[dot][1][0]=name
\/\/$kmean[dot][1][1]=value1
\/\/$kmean[dot][1][2]=value2
\/\/$kmean[dot][1][3]=value3
\/\/\u8a18\u9304\u8a72\u9ede\u76ee\u524d\u7684\u7fa4
\/\/$kmean[dotcluster][int:round:1~n][int:dot_id:1~n]
\/\/ex:
\/\/$kmean[dotcluster][0][1]= \u5b58\u653e\u7b2c0\u56de\u5408\u9ede1\u7684\u7fa4
\/\/$kmean[dotcluster[0][2]= \u5b58\u653e\u7b2c0\u56de\u5408\u9ede2\u7684\u7fa4
\/\/$kmean[dotcluster][1][1]= \u5b58\u653e\u7b2c1\u56de\u5408\u9ede1\u7684\u7fa4
\/\/$kmean[dotcluster][1][2]= \u5b58\u653e\u7b2c1\u56de\u5408\u9ede2\u7684\u7fa4
if($_GET['debug']==1)$debug=1; \/\/\u662f\u5426\u555f\u7528debug mode
\/\/\u53d6\u5f97\u6240\u6709\u9ede\u7684\u7e3d\u6578
$kmean['dotsum']=count($kmean['dot']);
\/\/\u521d\u59cb\u5316,\u5148\u5206\u914ddot\u7d66\u6307\u5b9a\u7684\u7fa4\u6578
if($_GET['init']==null)init_dotcentr_order($kmean);
if($_GET['init']==1)init_dotcentr_remainder($kmean);
if($_GET['init']==2)init_dotcentr_rand($kmean);

\/\/\/\/\/\/\/\/\/\/\/\/\/kmean start
while($same<$kmean['dotsum']){
  $same=0; \/\/\u8a18\u9304\u6240\u6709\u7684\u9ede\u5728\u672c\u56de\u5408\u8207\u4e0a\u4e00\u56de\u5408\u76f8\u540c\u7684\u6578\u91cf
  $round++; \/\/\u8a18\u9304\u56de\u5408\u6578
  if($debug==1)echo 'step'.$round.'\/n';
\/\/\u7b97\u7fa4\u4e2d\u5fc3
  for($i=1;$i<=$kmean['cluster']['count'];$i++){
\/\/$kmean[dot]\u662f\u6240\u6709\u9ede\u7684\u8cc7\u6599
\/\/$kmean[cluster][dot][$i]\u662f$i\u7fa4\u7684\u6240\u64c1\u6709\u7684\u9ede
\/\/\u7fa4\u4e2d\u5fc3\u7684\u7d50\u679c\u6703\u5132\u5b58\u5728$kmean[cluster][center][$i]
  computer_clustercenter($kmean['dot'], $kmean['cluster']['dot'][$i] , $kmean['cluster']['center'][$i]);
  if($debug==1){echo 'cluster'.$i.'center=';print_r($kmean['cluster']['center'][$i]);echo '< br>';}
}
unset($kmean['cluster']['dot']);\/\/\u6e05\u7a7a\u5404\u7fa4\u6709\u90a3\u4e9b\u9ede\u7684\u8cc7\u6599
\/\/\u91cd\u65b0\u53d6\u5f97\u5404\u9ede\u8207\u7fa4\u4e2d\u5fc3\u7684\u8ddd\u96e2,$i\u662f\u9ede
for($i=1;$i<=$kmean['dotsum'];$i++){
\/\/\u53d6\u5f97\u8a72\u9ede\u8207\u6240\u6709\u7fa4\u7684\u8ddd\u96e2,$j\u662f\u7fa4
  for($j=1;$j<=$kmean['cluster']['count'];$j++){
\/\/\u53d6\u5f97\u8a72\u9ede\u8207\u8a72\u7fa4\u7684\u8ddd\u96e2
    $dot2centr[$j]=computer_dot2center($kmean['dot'][$i],$kmean['cluster']['center'][$j]);
\/\/\u627e\u51fa\u548c\u9ede\u6700\u8fd1\u7684\u7fa4\u4e2d\u5fc3
    if($dot2centrmin==null){ 
      $dot2centrmin=$dot2centr[$j];
      $kmean['dotcluster'][$round][$i]=$j;
    }else{
      if($dot2centr[$j]<$dot2centrmin){
        $dot2centrmin=$dot2centr[$j];
        $kmean['dotcluster'][$round][$i]=$j; \/\/\u8a2d\u5b9a\u70ba\u9ede\u7684\u65b0\u7fa4
      }
    }
  }
if($debug==1){echo 'dot'.$i.'-clustercenter distance=';print_r($dot2centr);echo '< br>';}
$dot2centrmin=null;
$kmean['cluster']['dot'][$kmean['dotcluster'][$round][$i]][]=$i; \/\/\u8a18\u9304\u5404\u7fa4\u6709\u90a3\u4e9b\u9ede\u7684\u8cc7\u6599
if( $kmean['dotcluster'][$round][$i]== $kmean['dotcluster'][$round-1][$i])$same++;
\/\/\u82e5\u8a72\u56de\u5408\u9ede\u7684\u7fa4\u8207\u4e0a\u4e00\u56de\u5408\u76f8\u540c\u5247\u8a18\u9304\u6b21\u6578,\u5247\u8a18\u9304\u7684\u6b21\u6578\u7b49\u540c\u65bc\u65bc\u9ede\u7684\u6578\u91cf\u6642,\u5c31\u8868\u793a\u672c\u56de\u5408\u8207\u4e0a\u4e00\u56de\u5408\u9ede\u7684\u7fa4\u90fd\u6c92\u8b8a\u904e,\u5c07\u89f8\u767cwhile\u7684\u96e2\u958b\u689d\u4ef6
if($debug==1){if( $kmean['dotcluster'][$round][$i]!=$kmean['dotcluster'][$round-1][$i]){
  $debuginfo.=$i.'='.$kmean['dotcluster'][$round-1][$i].'--->'.$kmean['dotcluster'][$round][$i].'< br>';}}
}
if($debug==1){echo 'cluster dot=';print_r($kmean['cluster']['dot']);echo '< br>'.$debuginfo;}
}
\/\/\/\/\/\/kmean end

\/\/\/\/\/\/\/\/\/\/\/\/\/\/\/\/\/\/
\/\/report
for($i=1;$i<=$kmean['dotsum'];$i++){
while(list($key,$data)=each($kmean['dot'][$i])){
$result.=$data.',';
}
$result.='='.$kmean['dotcluster'][$round][$i];
$result.='
';
}
echo $result;
\/\/other infomation
echo '< br>round='.$round.'< br>';
}
\/\/kmean main function end
\/\/\/\/\/\/\/\/\/\/\/\/\/\/\/\/\/\/\/\/\/\/\/\/\/\/\/\/\/\/\/\/\/\/\/\/\/\/\/\/\/\/\/\/\/\/\/\/\/\/\/\/\/
\/\/\/\/\/\/\/\/\/\/\/\/\/\/\/\/\/\/\/\/\/\/\/\/\/\/\/\/\/\/\/\/\/\/\/\/\/\/\/\/\/\/\/\/\/\/\/\/\/\/\/\/\/\/\/\/\/\/\/\/\/\/\/\/\/\/\/\/\/\/\/\/\/\/\/\/\/\/\/\/\/\/\/\/\/\/\/\/\/\/\/\/\/\/\/
\/\/init method1 use %
function init_dotcentr_remainder(&$kmean){
for($i=1;$i<=$kmean['dotsum'];$i++){
$clusterid=$i%$kmean['cluster']['count']+1; \/\/\u8a08\u7b97\u8a72\u9ede\u5206\u914d\u7684\u7fa4
$kmean['dotcluster'][0][$i]=$clusterid; \/\/\u8a18\u9304\u7b2c0\u56de\u5408\u6642\u8a72\u9ede($i)\u70ba\u90a3\u4e00\u7fa4
$kmean['cluster']['dot'][$clusterid][]=$i; \/\/\u8a18\u9304\u8a72\u7fa4($clusterid)\u6709\u90a3\u4e9b\u9ede($i),\u4e26\u505a\u7d22\u5f15($j)
}
\/\/print_r($kmean['dotcluster'][0]);
\/\/print_r($kmean['cluster']['dot']);
}
\/\/\/\/\/\/\/\/\/\/\/\/\/\/\/\/\/\/\/\/\/\/\/\/\/\/\/\/\/\/\/\/\/\/\/\/\/\/\/\/\/\/\/\/\/\/\/\/\/\/\/\/\/\/\/\/\/\/\/\/\/\/\/\/\/\/\/\/\/\/\/\/\/\/\/\/\/\/\/\/\/\/\/\/\/\/\/\/\/\/\/\/\/\/\/\/
\/\/init method2 use rand
function init_dotcentr_rand(&$kmean){
for($i=1;$i<=$kmean['dotsum'];$i++){
$clusterid=rand(1,$kmean['cluster']['count']); \/\/\u8a08\u7b97\u8a72\u9ede\u5206\u914d\u7684\u7fa4
$kmean['dotcluster'][0][$i]=$clusterid; \/\/\u8a18\u9304\u7b2c0\u56de\u5408\u6642\u8a72\u9ede($i)\u70ba\u90a3\u4e00\u7fa4
$kmean['cluster']['dot'][$clusterid][]=$i; \/\/\u8a18\u9304\u8a72\u7fa4($clusterid)\u6709\u90a3\u4e9b\u9ede($i),\u4e26\u505a\u7d22\u5f15($j)
}
\/\/print_r($kmean['dotcluster'][0]);
\/\/print_r($kmean['cluster']['dot']);
}
\/\/\/\/\/\/\/\/\/\/\/\/\/\/\/\/\/\/\/\/\/\/\/\/\/\/\/\/\/\/\/\/\/\/\/\/\/\/\/\/\/\/\/\/\/\/\/\/\/\/\/\/\/\/\/\/\/\/\/\/\/\/\/\/\/\/\/\/\/\/\/\/\/\/\/\/\/\/\/\/\/\/\/\/\/\/\/\/\/\/\/\/\/\/\/
\/\/init method3 user by order
function init_dotcentr_order(&$kmean){
$groupcount=ceil($kmean['dotsum']\/$kmean['cluster']['count']);
for($i=1;$i<=$kmean['dotsum'];$i++){
if($i%$groupcount==1)$clusterid++;
$kmean['dotcluster'][0][$i]=$clusterid; \/\/\u8a18\u9304\u7b2c0\u56de\u5408\u6642\u8a72\u9ede($i)\u70ba\u90a3\u4e00\u7fa4
$kmean['cluster']['dot'][$clusterid][]=$i; \/\/\u8a18\u9304\u8a72\u7fa4($clusterid)\u6709\u90a3\u4e9b\u9ede($i),\u4e26\u505a\u7d22\u5f15($j)
}
\/\/print_r($kmean['dotcluster'][0]);
\/\/print_r($kmean['cluster']['dot']);
}
\/\/\/\/\/\/\/\/\/\/\/\/\/\/\/\/\/\/\/\/\/\/\/\/\/\/\/\/\/\/\/\/\/\/\/\/\/\/\/\/\/\/\/\/\/\/\/\/\/\/\/\/\/\/\/\/\/\/\/\/\/\/\/\/\/\/\/\/\/\/\/\/\/\/\/\/\/\/\/\/\/\/\/\/\/\/\/\/\/\/\/\/\/\/\/\/
function computer_clustercenter($dot_all,$clusterdot_index,&$center){
\/\/$dot_all[1~n][1~m]
\/\/$clusterdot_index[1~x]
\/\/$center[1~m]
$debug=0;
$clusterdot_total=count($clusterdot_index);
$dimensions_total=count($dot_all[1])-1;\/\/\u8981\u6263\u9664name,\u6240\u4ee5\u8981\u6e1b1
for($i=1;$i<=$dimensions_total;$i++){
$center[$i]=0;
for($j=0;$j<$clusterdot_total;$j++){
$center[$i]+=$dot_all[$clusterdot_index[$j]][$i];
if($debug==1){$debuginfo.=$dot_all[$clusterdot_index[$j]][$i].'+';}
}
$center[$i]=$center[$i]\/$clusterdot_total;
if($debug==1){echo $center[$i].'='.$debuginfo.'\/'.$clusterdot_total.';< br>'; $debuginfo=null;}
}
\/\/ (x1+x2+x3+...+xN)\/N,(y1+y2+y3+...+yN)\/N,(z1+z2+z3+...+zN)\/N
}
\/\/\/\/\/\/\/\/\/\/\/\/\/\/\/\/\/\/\/\/\/\/\/\/\/\/\/\/\/\/\/\/\/\/\/\/\/\/\/\/\/\/\/\/\/\/\/\/\/\/\/\/\/\/\/\/\/\/\/\/\/\/\/\/\/\/\/\/\/\/\/\/\/\/\/\/\/\/\/\/\/\/\/\/\/\/\/\/\/\/\/\/\/\/\/
function computer_dot2center($dot,$center){
$dimensions_total=count($center);
for($i=1;$i<=$dimensions_total;$i++){
$distance+=($dot[$i]-$center[$i])*($dot[$i]-$center[$i]);
}
\/\/sqrt( (1x-2x)^2 + (1y-2y)^2 + (1z-2z)^2 )
return $distance;
}
\/\/kmean end
\/\/\/\/\/\/\/\/\/\/\/\/\/\/\/\/\/\/\/\/\/\/\/\/\/\/\/\/\/\/\/\/\/\/\/\/\/\/\/\/\/\/\/\/\/\/\/\/\/\/\/\/\/\/\/\/\/\/\/\/\/\/\/\/\/\/\/\/\/\/\/\/\/\/\/\/\/\/\/\/\/\/\/\/\/\/\/\/\/
\/\/\/\/\/\/\/\/\/\/\/\/\/\/\/\/\/\/\/\/\/\/\/\/\/\/\/\/\/\/\/\/\/\/\/\/\/\/\/\/\/\/\/\/\/\/\/\/\/\/\/\/\/\/\/\/\/\/\/\/\/\/\/\/\/\/\/\/\/\/\/\/\/\/\/\/\/\/\/\/\/\/\/\/\/\/\/\/\/\/\/\/\/\/\/
\/\/\/\/\/\/\/\/\/\/optional function
\/\/\/\/\/\/\/\/\/\/\/\/\/\/\/\/\/\/\/\/\/\/\/\/\/\/\/
\/\/\/\/\/\/\/\/\/\/\/\/\/\/\/\/\/\/\/\/\/\/\/\/\/\/\/\/\/\/\/\/\/\/\/\/\/\/\/\/\/\/\/\/\/\/\/\/\/\/\/\/\/\/\/\/\/\/\/\/\/\/\/\/\/\/\/\/\/\/\/\/\/\/\/\/\/\/\/\/\/\/\/\/\/\/\/\/\/\/\/\/\/\/
function input_fromfile(&$kmeandot,$sourcedata){
\/*
data fromat description
name1,attribute1_value,attribute2_value,...attributeN_value
name2,attribute1_value,attribute2_value,...attributeN_value
.......
nameN,attribute1_value,attribute2_value,...attributeN_value
ex:
Iris-setosa,5.1,3.5,1.4,0.2
Iris-setosa,4.9,3,1.4,0.2
Iris-setosa,4.7,3.2,1.3,0.2
Iris-setosa,4.6,3.1,1.5,0.2
Iris-setosa,5,3.6,1.4,0.2
ps:attribute must is continue,not discrete
*\/
$inputdata=file($sourcedata);
$i=1;
while(list($key,$dot_id)=each($inputdata)){
$dot_id = trim($dot_id, \" n.\");
$kmeandot[$i++]=explode(',',$dot_id);
}
}
\/\/\/\/\/\/\/\/\/\/\/\/\/\/\/\/\/\/\/\/\/\/\/\/\/\/\/\/\/\/\/\/\/\/\/\/\/\/\/\/\/\/\/\/\/\/\/\/\/\/\/\/\/\/\/\/\/\/\/\/\/\/\/\/\/\/\/\/\/\/\/\/\/\/\/\/\/\/\/\/\/\/\/\/\/\/\/\/\/\/\/\/\/\/
function caclutime(){
$time = explode( \" \", microtime());
$usec = (double)$time[0];
$sec = (double)$time[1];
return $sec + $usec;
}<\/code><\/pre>\n","protected":false},"excerpt":{"rendered":"
k-means\u512a\u9ede:\u7c21\u55ae\u6709\u6548\u7387,\u4efb\u4e00\u6f14\u7b97\u6cd5\u6240\u5f97\u4e4b\u7d50\u679c,\u90fd\u53ef\u900f …<\/p>\n","protected":false},"author":1,"featured_media":0,"comment_status":"closed","ping_status":"closed","sticky":false,"template":"","format":"standard","meta":{"fifu_image_url":"","fifu_image_alt":"","_jetpack_memberships_contains_paid_content":false,"footnotes":"","jetpack_publicize_message":"","jetpack_publicize_feature_enabled":true,"jetpack_social_post_already_shared":false,"jetpack_social_options":{"image_generator_settings":{"template":"highway","default_image_id":0,"font":"","enabled":false},"version":2}},"categories":[13],"tags":[],"class_list":["post-463","post","type-post","status-publish","format-standard","hentry","category-dataanalysis"],"jetpack_publicize_connections":[],"jetpack_featured_media_url":"","jetpack_sharing_enabled":true,"_links":{"self":[{"href":"https:\/\/systw.net\/note\/wp-json\/wp\/v2\/posts\/463","targetHints":{"allow":["GET"]}}],"collection":[{"href":"https:\/\/systw.net\/note\/wp-json\/wp\/v2\/posts"}],"about":[{"href":"https:\/\/systw.net\/note\/wp-json\/wp\/v2\/types\/post"}],"author":[{"embeddable":true,"href":"https:\/\/systw.net\/note\/wp-json\/wp\/v2\/users\/1"}],"replies":[{"embeddable":true,"href":"https:\/\/systw.net\/note\/wp-json\/wp\/v2\/comments?post=463"}],"version-history":[{"count":0,"href":"https:\/\/systw.net\/note\/wp-json\/wp\/v2\/posts\/463\/revisions"}],"wp:attachment":[{"href":"https:\/\/systw.net\/note\/wp-json\/wp\/v2\/media?parent=463"}],"wp:term":[{"taxonomy":"category","embeddable":true,"href":"https:\/\/systw.net\/note\/wp-json\/wp\/v2\/categories?post=463"},{"taxonomy":"post_tag","embeddable":true,"href":"https:\/\/systw.net\/note\/wp-json\/wp\/v2\/tags?post=463"}],"curies":[{"name":"wp","href":"https:\/\/api.w.org\/{rel}","templated":true}]}}