Mathematica 教學 FindClusters and scatterplot

R-Project首頁圖示的散佈圖，覺得蠻有用的！

(* 產生隨機變數 *)

GaussianRandomData[n_Integer,p_,sigma_]:=
Table[p+{Re[#],Im[#]}&[RandomReal[NormalDistribution[0,sigma]]
E^(I*RandomReal[{0,2\[Pi]}])],{n}];

datapairs=BlockRandom[SeedRandom[1234];
Join[GaussianRandomData[100,{2,1},.3],
GaussianRandomData[100,{1,1.5},.2],
GaussianRandomData[100,{1,1.1},.4],
GaussianRandomData[100,{1.75,1.75},0.1]]];

(* 資料分群 *)

cl=FindClusters[datapairs];


(* 計算各樣本點與總體平均的距離，距離越大，點也越大 *)

pointdata=Transpose[{{Red,Blue,Green,Brown},
Flatten[{PointSize[EuclideanDistance[Mean@datapairs,#]/50],
Point[#]}&/@cl[[#]]]&/@Range[4]}];

Plot[0,{x,Min@datapairs[[All,1]],1.2Max@datapairs[[All,2]]},
AxesOrigin->Mean@datapairs,
PlotRange->{{0,1.2Max@datapairs[[All,1]]},
{0,1.2Max@datapairs[[All,2]]}},
Epilog->pointdata]

Tags: Mathematica 教學, Mathematica 教學 畫圖, Teaching

• 編輯文章
• 此文章網址
• 轉寄此文章
• 目前有0 意見