The likelihood
module
- class pylife.materialdata.woehler.likelihood.Likelihood(fatigue_data)[source]
Calculate the likelihood a fatigue dataset matches with Wöhler curve parameters.
- likelihood_infinite(SD, TS)[source]
Produces the likelihood functions that are needed to compute the endurance limit and the scatter in load direction. The likelihood functions are represented by a cummalative distribution function. The likelihood function of a runout is 1-Li(fracture).
- Parameters:
SD – Endurnace limit start value to be optimzed, unless the user fixed it.
TS – The scatter in load direction 1/TS to be optimzed, unless the user fixed it.
- Returns:
Sum of the log likelihoods. The negative value is taken since optimizers in statistical packages usually work by minimizing the result of a function. Performing the maximum likelihood estimate of a function is the same as minimizing the negative log likelihood of the function.
- Return type:
neg_sum_lolli
- likelihood_total(SD, TS, k_1, ND, TN)[source]
Produces the likelihood functions that are needed to compute the parameters of the woehler curve. The likelihood functions are represented by probability and cummalative distribution functions. The likelihood function of a runout is 1-Li(fracture). The functions are added together, and the negative value is returned to the optimizer.
- Parameters:
SD – Endurnace limit start value to be optimzed, unless the user fixed it.
TS – The scatter in load direction 1/TS to be optimzed, unless the user fixed it.
k_1 – The slope k_1 to be optimzed, unless the user fixed it.
ND – Load-cycle endurance start value to be optimzed, unless the user fixed it.
TN – The scatter in load-cycle direction 1/TN to be optimzed, unless the user fixed it.
- Returns:
Sum of the log likelihoods. The negative value is taken since optimizers in statistical packages usually work by minimizing the result of a function. Performing the maximum likelihood estimate of a function is the same as minimizing the negative log likelihood of the function.
- Return type:
neg_sum_lolli