function compParams = initRandom(k, X)
  d = size(X,2) ;

  compParams = [];
  for i=1:k
    compParams = [compParams, ...
		  struct('prior',1/k,'mean',normrnd(0,1,1,d),'cov',eye(d))];
  end  