public void ComputeTest2()
{
MultinomialLogisticRegression mlr = createExample1();
double[][] inputs = example1.Submatrix(null, 1, 2).ToJagged();
double[] outputs = example1.Submatrix(null, 0, 0).Reshape(0);
double[] responses;
// Tested against values extracted from predicted probabilities
// table from: http://www.ats.ucla.edu/stat/r/dae/mlogit.htm
responses = mlr.Compute(inputs[0]);
Assert.AreEqual(0.9479577862063925, responses[0], 1e-5);
Assert.AreEqual(0.0502297144022469, responses[1], 1e-5);
Assert.AreEqual(0.0018124993913602, responses[2], 1e-5);
responses = mlr.Compute(inputs[5]);
Assert.AreEqual(0.772875639435192, responses[0], 1e-5);
Assert.AreEqual(0.208690558456066, responses[1], 1e-5);
Assert.AreEqual(0.018433802108742, responses[2], 1e-5);
responses = mlr.Compute(inputs[11]);
Assert.AreEqual(0.772875639435192, responses[0], 1e-5);
Assert.AreEqual(0.208690558456066, responses[1], 1e-5);
Assert.AreEqual(0.018433802108742, responses[2], 1e-5);
responses = mlr.Compute(inputs[12]);
Assert.AreEqual(0.695617266629850, responses[0], 1e-5);
Assert.AreEqual(0.271439833912059, responses[1], 1e-5);
Assert.AreEqual(0.032942899458091, responses[2], 1e-5);
}