KrigingPOD

class KrigingPOD(*args)

Kriging based POD.

Available constructor:

KrigingPOD(inputSample, outputSample, detection, noiseThres, saturationThres, boxCox)

Parameters:

inputSample : 2-d sequence of float

Vector of the input values. The first column must correspond with the defect sizes.

outputSample : 2-d sequence of float

Vector of the signals, of dimension 1.

detection : float

Detection value of the signal.

noiseThres : float

Value for low censored data. Default is None.

saturationThres : float

Value for high censored data. Default is None

boxCox : bool or float

Enable or not the Box Cox transformation. If boxCox is a float, the Box Cox transformation is enabled with the given value. Default is False.

Warning

The first column of the input sample must corresponds with the defects sample.

Notes

This class aims at building the POD based on a kriging model. No assumptions are required for the residuals with this method. The POD are computed by simulating conditional prediction. For each, a Monte Carlo simulation is performed. The accuracy of the Monte Carlo simulation is taken into account using the TCL.

The return POD model corresponds with an interpolate function built with the POD values computed for the given defect sizes. The default values are 20 defect sizes between the minimum and maximum value of the defect sample. The defect sizes can be changed using the method setDefectSizes.

The default kriging model is built with a linear basis only for the defect size and constant otherwise. The covariance model is an anisotropic squared exponential model. Parameters are estimated using the TNC algorithm, the initial starting point of the TNC is found thanks to a quasi random search of the best loglikelihood value among 1000 computations.

For advanced use, all parameters can be defined thanks to dedicated set methods. Moreover, if the user has already built a kriging result, it can be given as parameter using the method setKrigingResult, then the POD is computed based on this kriging result.

A progress bar is shown if the verbosity is enabled. It can be disabled using the method setVerbose.

Methods

computeDetectionSize(**kwargs) Compute the detection size for a given probability level.
drawBoxCoxLikelihood([name]) Draw the loglikelihood versus the Box Cox parameter.
drawPOD(**kwargs) Draw the POD curve.
drawValidationGraph(**kwargs) Draw the validation graph of the metamodel.
getBasis() Accessor to the kriging basis.
getBoxCoxParameter() Accessor to the Box Cox parameter.
getCovarianceModel() Accessor to the kriging covariance model.
getDefectSizes() Accessor to the defect size where POD is computed.
getDistribution() Accessor to the parameters distribution.
getInitialStartSize() Accessor to the initial random search size.
getKrigingResult() Accessor to the kriging result.
getPODCLModel([confidenceLevel]) Accessor to the POD model at a given confidence level.
getPODModel() Accessor to the POD model.
getQ2() Accessor to the Q2 value.
getSamplingSize() Accessor to the Monte Carlo sampling size.
getSimulationSize() Accessor to the simulation size.
getVerbose() Accessor to the verbosity.
run() Build the POD models.
setBasis(basis) Accessor to the kriging basis.
setCovarianceModel(covarianceModel) Accessor to the kriging covariance model.
setDefectSizes(size) Accessor to the defect size where POD is computed.
setDistribution(distribution) Accessor to the parameters distribution.
setInitialStartSize(size) Accessor to the initial random search size.
setKrigingResult(result) Accessor to the kriging result.
setSamplingSize(size) Accessor to the Monte Carlo sampling size.
setSimulationSize(size) Accessor to the simulation size.
setVerbose(verbose) Accessor to the verbosity.
computeDetectionSize(**kwargs)

Compute the detection size for a given probability level.

Parameters:

probabilityLevel : float

The probability level for which the defect size is computed.

confidenceLevel : float

The confidence level associated to the given probability level the defect size is computed. Default is None.

Returns:

result : collection of openturns.NumericalPointWithDescription

A NumericalPointWithDescription containing the detection size computed at the given probability level and confidence level if provided.

drawBoxCoxLikelihood(name=None)

Draw the loglikelihood versus the Box Cox parameter.

Parameters:

name : string

name of the figure to be saved with transparent option sets to True and bbox_inches=’tight’. It can be only the file name or the full path name. Default is None.

Returns:

fig : matplotlib.figure

Matplotlib figure object.

ax : matplotlib.axes

Matplotlib axes object.

Notes

This method is available only when the parameter boxCox is set to True.

drawPOD(**kwargs)

Draw the POD curve.

Parameters:

probabilityLevel : float

The probability level for which the defect size is computed. Default is None.

confidenceLevel : float

The confidence level associated to the given probability level the defect size is computed. Default is None.

defectMin, defectMax : float

Define the interval where the curve is plotted. Default : min and max values of the input sample.

nbPt : int

The number of points to draw the curves. Default is 100.

name : string

name of the figure to be saved with transparent option sets to True and bbox_inches=’tight’. It can be only the file name or the full path name. Default is None.

Returns:

fig : matplotlib.figure

Matplotlib figure object.

ax : matplotlib.axes

Matplotlib axes object.

drawValidationGraph(**kwargs)

Draw the validation graph of the metamodel.

Parameters:

name : string

name of the figure to be saved with transparent option sets to True and bbox_inches=’tight’. It can be only the file name or the full path name. Default is None.

Returns:

fig : matplotlib.figure

Matplotlib figure object.

ax : matplotlib.axes

Matplotlib axes object.

getBasis()

Accessor to the kriging basis.

Returns:

basis : openturns.Basis

The basis used as trend in the kriging model. Default is a linear basis for the defect and constant for the other parameters.

getBoxCoxParameter()

Accessor to the Box Cox parameter.

Returns:

lambdaBoxCox : float

The Box Cox parameter used to transform the data. If the transformation is not enabled None is returned.

getCovarianceModel()

Accessor to the kriging covariance model.

Returns:

covarianceModel : openturns.CovarianceModel

The covariance model in the kriging model. Default is an anisotropic squared exponential covariance model.

getDefectSizes()

Accessor to the defect size where POD is computed.

Returns:

defectSize : sequence of float

The defect sizes where the Monte Carlo simulation is performed to compute the POD.

getDistribution()

Accessor to the parameters distribution.

Returns:

distribution : openturns.ComposedDistribution

The input parameters distribution used for the Monte Carlo simulation. Default is a Uniform distribution for all parameters.

getInitialStartSize()

Accessor to the initial random search size.

Returns:

size : int

The size of the initial random search to find the best loglikelihood value to start the TNC algorithm to optimize the covariance model parameters. Default is 1000.

getKrigingResult()

Accessor to the kriging result.

Returns:

result : openturns.KrigingResult

The kriging result.

getPODCLModel(confidenceLevel=0.95)

Accessor to the POD model at a given confidence level.

Parameters:

confidenceLevel : float

The confidence level the POD must be computed. Default is 0.95

Returns:

PODModelCl : openturns.NumericalMathFunction

The function which computes the probability of detection for a given defect value at the confidence level given as parameter.

getPODModel()

Accessor to the POD model.

Returns:

PODModel : openturns.NumericalMathFunction

The function which computes the probability of detection for a given defect value.

getQ2()

Accessor to the Q2 value.

Returns:

Q2 : float

The Q2 value computed analytically using Dubrule (1983) technique.

getSamplingSize()

Accessor to the Monte Carlo sampling size.

Returns:

size : int

The size of the Monte Carlo simulation used to compute the POD for each defect size.

getSimulationSize()

Accessor to the simulation size.

Returns:

size : int

The size of the simulation used to compute the confidence interval.

getVerbose()

Accessor to the verbosity.

Returns:

verbose : bool

Enable or disable the verbosity. Default is True.

run()

Build the POD models.

Notes

This method build the kriging model. First the censored data are filtered if needed. The Box Cox transformation is performed if it is enabled. Then it builds the POD models : conditional samples are simulated for each defect size, then the distributions of the probability estimator (for MC simulation) are built. Eventually, a sample of this distribution is used to compute the mean POD and the POD at the confidence level.

setBasis(basis)

Accessor to the kriging basis.

Parameters:

basis : openturns.Basis

The basis used as trend in the kriging model.

setCovarianceModel(covarianceModel)

Accessor to the kriging covariance model.

Parameters:

covarianceModel : openturns.CovarianceModel

The covariance model in the kriging model.

setDefectSizes(size)

Accessor to the defect size where POD is computed.

Parameters:

defectSize : sequence of float

The defect sizes where the Monte Carlo simulation is performed to compute the POD.

setDistribution(distribution)

Accessor to the parameters distribution.

Parameters:

distribution : openturns.ComposedDistribution

The input parameters distribution used for the Monte Carlo simulation.

setInitialStartSize(size)

Accessor to the initial random search size.

Parameters:

size : int

The size of the initial random search to find the best loglikelihood value to start the TNC algorithm to optimize the covariance model parameters.

setKrigingResult(result)

Accessor to the kriging result.

Parameters:

result : openturns.KrigingResult

The kriging result.

setSamplingSize(size)

Accessor to the Monte Carlo sampling size.

Parameters:

size : int

The size of the Monte Carlo simulation used to compute the POD for each defect size.

setSimulationSize(size)

Accessor to the simulation size.

Parameters:

size : int

The size of the simulation used to compute the confidence interval.

setVerbose(verbose)

Accessor to the verbosity.

Parameters:

verbose : bool

Enable or disable the verbosity.