WOLFRAM

Wolfram言語 ＆ システムドキュメントセンター

二変数のデータ集合をクラスタ化する

オールドフェイスフル間欠泉の噴射時間の長さについてのデータを得て，それに基づく分布を構築する：
Wolfram Language code: data = ExampleData[{"Statistics", "OldFaithful"}]; 𝒟 = KernelMixtureDistribution[data, "SheatherJones"];
この分布に基づく噴射のシミュレーションを行う：
Wolfram Language code: rng = RandomVariate[𝒟, 500];
FindClustersを使ってデータを分ける：
Wolfram Language code: Show[SmoothDensityHistogram[rng, ColorFunction -> "TemperatureMap", Mesh -> 10], ListPlot[FindClusters[rng, 2], PlotRange -> {{1, 6}, {30, 110}}, AspectRatio -> 1, Frame -> True, ImageSize -> 220, PlotStyle -> {Directive[PointSize[0.0185], Red], Directive[PointSize[0.0185], Black]}, PlotLabel -> Text[Style["Simulated Eruptions", FontFamily -> "Verdana", Bold]], Background -> Directive[Opacity[.5], White], Axes -> None]]
3Dの可視化を作成する：
Wolfram Language code: p2 = ListPlot[FindClusters[rng, 2], PlotRange -> {{1, 6}, {30, 110}}, AspectRatio -> 1, Frame -> True, ImageSize -> 220, PlotStyle -> {Directive[PointSize[0.0185], Red], Directive[PointSize[0.0185], Black]}, PlotLabel -> Text[Style["Simulated Eruptions", FontFamily -> "Verdana", Bold]], Background -> Directive[Opacity[.5], White], Axes -> None];
Wolfram Language code: p1 = Plot3D[Evaluate[PDF[𝒟, {x, y}]], {x, 1, 6}, {y, 30, 105}, Mesh -> 25, PlotRange -> All, PlotLabel -> Text[Style["Duration vs. Waiting Time", FontFamily -> "Verdana", Large]], BoxRatios -> 1, PlotPoints -> 100, MeshStyle -> Opacity[.3], ColorFunction -> "TemperatureMap", Boxed -> False, ViewPoint -> {2, -2, 1}, Epilog -> Inset[p2, {Right, Bottom}, {Right, Bottom}], ImageSize -> 570]