Dear all,

 

I am attempting to model some one-dimensional data using Gaussian
mixture model with mclust.  Generally, the data that I have have 3
overlapping populations (with one of them being the majority, and the
other two combining to less than 15%) and for some reason, mclust
consistently ignores the smaller peaks, giving me strange values for the
means (completely anti-intuitive in terms of visual inspection by using
plot(density())).  However, when I use rnorm to generate some fake data,
mclust works consistently well.  As I am quite new to R, I was wondering
whether anyone could help me out.  I am using the em function in mclust,
with parameters that are arbitrarily set 

 

Pro <- c(0.33,0.33,0.33)

Mean<- c(-0.5,0,0.5)

Sigmasq <- c(0.05,0.05,0.05)

 

And a sample of my data set is

 

c(-0.016133,-0.07668,-0.000625,0.031329,-0.011094,0.014199,-0.014141,0.0
36836,-0.007695,-0.011738,-0.021953,0.046565,0.010859,-0.050313,-0.16781
3,0.01543,-0.057598,-0.034336,0.182275,0.032959,0.01918,0.009248,-0.1952
73,-0.00918,-0.017813,0.003828,-0.113867,0.004014,0.031504,0.004277,0.05
2188,-0.030859,-0.214023,-0.329102,-0.07832,-0.008379,0.05833,-0.007285,
-0.036992,0.035768,0.055006,-0.000781,0.005067,-0.025811,0.021016,-0.002
598,-0.036799,-0.03119,-0.004482,-0.024473,-0.108115,-0.116318,0.158008,
0.04252,-0.032129,0.00707,-0.073398,-0.115605,-0.033945,-0.022793,0.0418
55,0.006289,0.250273,0.042607,-0.000449,0.030098,0.041238,-0.028926,-0.1
11895,0.003867,0.015625,-0.018906,-0.00291,-0.027188,0.00957,-0.133369,0
.018652,0.138652,0.038789,-0.050107,0.135908,-0.052725,0.005977,0.030977
,0.005371,-0.179902,-0.008691,0.033711,0.164033,-0.063457,0.022734,-0.04
7227,0.025918,-0.005557,-0.104453,0.021348,-0.054902,-0.069277,-0.115273
,0.038906,0.171211,0.000645,-0.064873,0.014062,-0.00252,-0.017715,-0.000
586,0.174609,-0.056396,0.000937,-0.217148,-0.203105,0.006533,0.015371,-0
.024629,0.015244,0.002949,0.024805,-0.104023,0.007964,-0.198633,-0.00783
3,-0.000518,0.018232,0.000195,0.028575,-0.028145,-0.030098,-0.002148,-0.
035723,-0.005996,0.023027,0.034512,0.009189,0.049252,-0.016641,-0.023262
,-0.013379,-0.013633,0.150996,0.040391,0.153809,0.001182,0.040371,-0.016
191,-0.05097,-0.12041,0.042617,0.018828,-0.002617,-0.043887,-0.025764,-0
.016836,0.023535,0.040625,0.158789,-0.026934,0.02791,-0.108027,0.037979,
-0.011865,0.06127,0.04416,-0.07833,0.019922,0.000685,-0.071885,0.000479,
0.006211,-0.030879,-0.009188,-0.061895,-0.069102,0.032051,-0.082637,-0.2
46484,-0.015586,-0.008555,0.265664,0.050781,0.008242,0.169785,-0.025977,
0.017871,-0.239492,0.005234,0.006865,0.007344,0.237861,-0.110742,0.20800
8,0.189336,0.205469,-0.111729,-0.023438,0.49,-0.028281,0.177988,0.00998,
-0.002402,0.161924,0.220859,0.026455,0.019629,-0.015098,-0.110771,-0.014
414,0.121211,0.028439,-0.026143,0.024989,-0.060801,-0.023124,-0.012734,0
.168398,0.039955,-0.038984,-0.028543,0.157412,-0.015547,-0.012617,0.0316
07,-0.053437,0.027246,0.003906,-0.218613,0.024902,0.020273,0.011914,0.16
2051,0.00582,-0.019189,-0.009029,0.001875,0.015273,0.175303,-0.092441,-0
.086738,-0.022871,0.027852,-0.108809,0.005938,-0.016543,-0.019288,0.2105
66,-0.022813,-0.001748,-0.108574,0.164971,-0.075186)

 

If anyone could help me out, I would be extremely grateful. 

 

Best,

 

Ken Lo

 

 

**********************************************

Ken Lo, PhD
Post-Doctorate Research Affliliate

Department of Cancer Genetics

Buffalo Life Science Complex, L2-104
Roswell Park Cancer Institute
Elm and Carlton Streets
Buffalo, New York 14263
Telephone:  716-845-3941
Fax:         716-845-3940
E-mail:      Ken.Lo@RoswellPark.org <mailto:Ken.Lo@RoswellPark.org> 
Web:        www.RoswellPark.org <http://www.RoswellPark.org>  

Located in the
Buffalo Life Science Complex
on the Buffalo Niagara Medical Campus

**********************************************









This email message may contain legally privileged and/or confidential information.  If you are not the intended recipient(s), or the employee or agent responsible for the delivery of this message to the intended recipient(s), you are hereby notified that any disclosure, copying, distribution, or use of this email message is prohibited.  If you have received this message in error, please notify the sender immediately by e-mail and delete this email message from your computer. Thank you.
	[[alternative HTML version deleted]]

