public void RunTest4()
{
// Data from: http://www.sph.emory.edu/~cdckms/CoxPH/prophaz2.html
// with added tied times
double[,] data =
{
{ 50, 1, 1 },
{ 60, 1, 1 },
{ 40, 1, 1 },
{ 51, 1, 1 },
{ 70, 2, 1 },
{ 45, 3, 0 },
{ 35, 5, 0 },
{ 62, 7, 1 },
{ 50, 11, 0 },
{ 45, 4, 0 },
{ 57, 6, 0 },
{ 32, 8, 0 },
{ 57, 9, 1 },
{ 60, 10, 1 },
};
var regression = new ProportionalHazards(1);
double[][] inputs = data.GetColumn(0).ToJagged();
double[] time = data.GetColumn(1);
int[] output = data.GetColumn(2).ToInt32();
var target = new ProportionalHazardsNewtonRaphson(regression);
double error = target.Run(inputs, time, output);
double log = -2 * regression.GetPartialLogLikelihood(inputs, time, output);
Assert.AreEqual(0.04863, regression.Coefficients[0], 1e-4);
Assert.IsFalse(Double.IsNaN(regression.Coefficients[0]));
Assert.AreEqual(0.04186, regression.StandardErrors[0], 1e-4);
Assert.IsFalse(Double.IsNaN(regression.StandardErrors[0]));
}