computeSystemSurvivalSignature {ReliabilityTheory}

Compute the survival signature of a system

Description

The system survival signature (Coolen and Coolen-Maturi, 2012) is a generalisation of the signature to systems with multiple component types. This function automatically computes the survival signature of the specified system. Here, system implies components (as opposed to links) are unreliable.

Usage

computeSystemSurvivalSignature(sys, cutsets=NULL, frac=FALSE)

Arguments

sys

a system object representing the system whose types are to be set. This should have been created by a call to createSystem.

cutsets

if the cut-sets of the system or network are already known they may be passed in as a list of numeric vectors. This can save time because cut-set computation is the slowest part of the algorithm. Leaving as NULL causes the function to find the cut sets itself.

frac

if TRUE then the function prints out survival signature probabilities as fractions rather than decimals.

Details

The survival signature of a system with K types of component is the functional \Phi(l_1, \dots, l_K) giving the probability that the system works given exactly l_k of the components of type k are working. See Coolen and Coolen-Maturi (2012) for details. Thus, the survival signature can be represented by a table with K+1 columns, the first K being the number of each type of component which is working and the final column being the probability the system works.

The system or network is specified by means of a system object, whereby each end of the system is denoted by nodes named s and t which are taken to be perfectly reliable. It is easy to construct the appropriate reliability block diagram representation using the function createSystem. Note that each physically distinct component should be separately numbered when constructing this object.

Once the topology of the system has been defined (or at definition time), one must indicate the type of each component (if not done when initially calling createSystem it can later be modified using setCompTypes). The Examples section below features the full computation of the survival signature for Figure 1 in Coolen and Coolen-Maturi (2012) and Figure 2 in Coolen et al (2013) to make this clear.

Value

computeSystemSurvivalSignature returns a data frame with K+1 columns. The first K columns represent the function inputs, l_1, \dots, l_K and the final column is the probability that the system works given the corresponding numbers of each component which are working.

Note

Please feel free to email [e-mail] with any queries or if you encounter errors when running this function.

Author(s)

Louis J. M. Aslett [e-mail] (https://www.louisaslett.com/)

References

Coolen, F. P. A. and Coolen-Maturi, T. (2012), Generalizing the signature to systems with multiple types of components, in 'Complex Systems and Dependability', Springer, pp. 115-130.

Coolen, F. P. A., Coolen-Maturi, T., Al-nefaiee, A. H. and Aboalkhair, A. M. (2013), ‘Recent advances in system reliability using the survival signature’, Proceedings of Advances in Risk and Reliability Technology Symposium, Loughborough.

See Also

computeSystemSignature

Examples

## EXAMPLE 1
## Figure 1 in Coolen and Coolen-Maturi (2012)

# First, define the structure, ensuring that each physically separate component
# is separately numbered
fig1 <- createSystem(s -- 1 -- 2:3 -- 4 -- 5:6 -- t, 2 -- 5, 3 -- 6)

# Second, specify the type of each of those numbered components
# (leaving s,t with no type)
fig1 <- setCompTypes(fig1,
                     list("Type 1" = c("1","2","5"),
                          "Type 2" = c("3","4","6")))

# Third, compute the survival signature (getting fractions rather than decimals)
computeSystemSurvivalSignature(fig1, frac = TRUE)



## EXAMPLE 2
## Figure 3 in Coolen et al (2013)

# First, define the structure, ensuring that each physically separate component
# is separately numbered.
# For this example, we demonstrate how to define the component types at system
# creation time
fig3 <- createSystem(s -- 1:4 -- 2:5 -- 3:6 -- t, s -- 7:8, 8 -- 9, 7:9 -- t,
                     types = list("Type 1" = "1",
                                  "Type 2" = c("2","3","4","7"),
                                  "Type 3" = c("5","6","8","9")))

# Third, compute the survival signature (getting fractions rather than decimals)
computeSystemSurvivalSignature(fig3, frac=TRUE)