private double createBaseline(double[] time, int[] censor, double[] output)
{
if (regression.BaselineHazard != null)
{
IFittingOptions options = null;
if (regression.BaselineHazard is IFittableDistribution<double, EmpiricalHazardOptions>)
{
// Compute an estimate of the cumulative Hazard
// function using the Nelson-Aalen estimator
options = new EmpiricalHazardOptions()
{
Censor = censor,
Output = output,
Estimator = HazardEstimator.BreslowNelsonAalen
};
}
else if (regression.BaselineHazard is IFittableDistribution<double, SurvivalOptions>)
{
// Compute an estimate of the cumulative Hazard
// function using the Nelson-Aalen estimator
options = new SurvivalOptions()
{
Censor = censor
};
}
regression.BaselineHazard.Fit(time, options);
}
return 0;
}
