Varying the number of Gaussians (anti-annular region)

Generate data from above distributions

Number of points for training and test data

n1 = 1000 ; n2 = 1000 ;

Generate data

data1 = AntiAnnularRegion[c1, rr1, rr2, n1] ; data2 = AntiAnnularRegion[c1, rr1, rr2, n2] ;

Visualize training data

gd1 = PlotData2D[data1, pr, style1, nS] ; gd = Show[gd1, FS] ;

[Graphics:../HTMLFiles/index_46.gif]

Generate sample EM results with varying number of component densities

results = Table[EMAll1[data1, k], {k, 2, 15}] ;

Plot log-likelihood of training data as a function of k

PlotLogLikelihood2[Last /@ results, data1] ;

[Graphics:../HTMLFiles/index_49.gif]

Plot log-likelihood of test data as a function of k

PlotLogLikelihood2[Last /@ results, data2] ;

[Graphics:../HTMLFiles/index_51.gif]


Created by Mathematica  (September 8, 2003)