package dist.test;
import dist.MultivariateGaussian;
import shared.DataSet;
import shared.Instance;
import util.linalg.DenseVector;
import util.linalg.RectangularMatrix;
/**
* Testing
* @author Andrew Guillory gtg008g@mail.gatech.edu
* @version 1.0
*/
public class MultivariateGaussianTest {
/**
* The test main
* @param args ignored
*/
public static void main(String[] args) {
Instance[] instances = new Instance[20];
MultivariateGaussian mga = new MultivariateGaussian(new DenseVector(new double[] {100, 100, 100}), RectangularMatrix.eye(3).times(.01));
for (int i = 0; i < instances.length; i++) {
instances[i] = mga.sample();
System.out.println(instances[i]);
}
DataSet set = new DataSet(instances);
MultivariateGaussian mg = new MultivariateGaussian();
mg.estimate(set);
System.out.println(mg);
System.out.println("Most likely " + mg.mode(null));
for (int i = 0; i < 10; i++) {
System.out.println(mg.sample(null));
}
for (int i = 0; i < instances.length; i++) {
System.out.println("Probability of \n" + instances[i]
+ "\n " + mg.p(instances[i]));
}
}
}