def test_Fit_Exponential_2P():
dist = Exponential_Distribution(Lambda=5, gamma=500)
rawdata = dist.random_samples(20, seed=5)
data = make_right_censored_data(data=rawdata, threshold=dist.mean)
MLE = Fit_Exponential_2P(failures=data.failures,
right_censored=data.right_censored,
method='MLE',
show_probability_plot=False,
print_results=False)
assert_allclose(MLE.Lambda, 7.062867654421206, rtol=rtol, atol=atol)
assert_allclose(MLE.gamma, 500.016737532126, rtol=rtol, atol=atol)
assert_allclose(MLE.AICc, -23.939665128347745, rtol=rtol, atol=atol)
assert_allclose(MLE.BIC, -22.65408293418094, rtol=rtol, atol=atol)
assert_allclose(MLE.loglik, 14.322773740644461, rtol=rtol, atol=atol)
assert_allclose(MLE.AD, 29.413655089419287, rtol=rtol, atol=atol)
LS = Fit_Exponential_2P(failures=data.failures,
right_censored=data.right_censored,
method='LS',
show_probability_plot=False,
print_results=False)
assert_allclose(LS.Lambda, 6.4445633542175, rtol=rtol, atol=atol)
assert_allclose(LS.gamma, 500.01368943066706, rtol=rtol, atol=atol)
assert_allclose(LS.AICc, -23.031777273560103, rtol=rtol, atol=atol)
assert_allclose(LS.BIC, -21.7461950793933, rtol=rtol, atol=atol)
assert_allclose(LS.loglik, 13.86882981325064, rtol=rtol, atol=atol)
assert_allclose(LS.AD, 29.33840933641424, rtol=rtol, atol=atol)
def test_Fit_Exponential_1P():
dist = Exponential_Distribution(Lambda=5)
rawdata = dist.random_samples(20, seed=5)
data = make_right_censored_data(data=rawdata, threshold=dist.mean)
MLE = Fit_Exponential_1P(failures=data.failures,
right_censored=data.right_censored,
method='MLE',
show_probability_plot=False,
print_results=False)
assert_allclose(MLE.Lambda, 6.101198944227536, rtol=rtol, atol=atol)
assert_allclose(MLE.AICc, -22.032339191099148, rtol=rtol, atol=atol)
assert_allclose(MLE.BIC, -21.25882913976738, rtol=rtol, atol=atol)
assert_allclose(MLE.loglik, 12.127280706660684, rtol=rtol, atol=atol)
assert_allclose(MLE.AD, 29.59913306667145, rtol=rtol, atol=atol)
LS = Fit_Exponential_1P(failures=data.failures,
right_censored=data.right_censored,
method='LS',
show_probability_plot=False,
print_results=False)
assert_allclose(LS.Lambda, 5.776959885774546, rtol=rtol, atol=atol)
assert_allclose(LS.AICc, -21.988412212242917, rtol=rtol, atol=atol)
assert_allclose(LS.BIC, -21.214902160911148, rtol=rtol, atol=atol)
assert_allclose(LS.loglik, 12.10531721723257, rtol=rtol, atol=atol)
assert_allclose(LS.AD, 29.52124203457833, rtol=rtol, atol=atol)
def test_Fit_Exponential_2P():
dist = Exponential_Distribution(Lambda=5, gamma=500)
rawdata = dist.random_samples(20, seed=5)
data = make_right_censored_data(data=rawdata, threshold=dist.mean)
fit = Fit_Exponential_2P(failures=data.failures, right_censored=data.right_censored, show_probability_plot=False, print_results=False)
assert_allclose(fit.Lambda, 7.00351280734533,rtol=rtol,atol=atol)
assert_allclose(fit.gamma, 500.015837532126,rtol=rtol,atol=atol)
assert_allclose(fit.AICc, -23.686473231109936,rtol=rtol,atol=atol)
assert_allclose(fit.loglik, 14.196177792025557,rtol=rtol,atol=atol)
def test_Fit_Exponential_1P():
dist = Exponential_Distribution(Lambda=5)
rawdata = dist.random_samples(20, seed=5)
data = make_right_censored_data(data=rawdata, threshold=dist.mean)
fit = Fit_Exponential_1P(failures=data.failures, right_censored=data.right_censored, show_probability_plot=False, print_results=False)
assert_allclose(fit.Lambda, 6.101199434421275,rtol=rtol,atol=atol)
assert_allclose(fit.gamma, 0,rtol=rtol,atol=atol)
assert_allclose(fit.AICc, -22.032339191099254,rtol=rtol,atol=atol)
assert_allclose(fit.loglik, 12.127280706660738,rtol=rtol,atol=atol)