Next Article in Journal
Drying of a Colloidal Suspension Deposited on a Substrate: Experimental and Numerical Studies
Next Article in Special Issue
Experiments and Mechanical Simulation on Bubble Concrete: Studies on the Effects of Shape and Position of Hollow Bodies Mixed in Concrete
Previous Article in Journal
Tunable Low Frequency Band Gap and Waveguide of Phononic Crystal Plates with Different Filling Ratio
Previous Article in Special Issue
Comparison of Machine Learning Approaches with Traditional Methods for Predicting the Compressive Strength of Rice Husk Ash Concrete
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

On Modeling Concrete Compressive Strength Data Using Laplace Birnbaum-Saunders Distribution Assuming Contaminated Information

by
Farouq Mohammad A. Alam
1,* and
Mazen Nassar
1,2
1
Department of Statistics, Faculty of Science, King Abdulaziz University, P.O. Box 80200, Jeddah 201589, Saudi Arabia
2
Department of Statistics, Faculty of Commerce, Zagazig University, Zagazig 44519, Egypt
*
Author to whom correspondence should be addressed.
Crystals 2021, 11(7), 830; https://doi.org/10.3390/cryst11070830
Submission received: 30 June 2021 / Revised: 14 July 2021 / Accepted: 15 July 2021 / Published: 17 July 2021

Abstract

:
Compressive strength is a well-known measurement to evaluate the endurance of a given concrete mixture to stress factors, such as compressive loads. A suggested approach to assess compressive strength of concrete is to assume that it follows a probability model from which its reliability is calculated. In reliability analysis, a probability distribution’s reliability function is used to calculate the probability of a specimen surviving to a certain threshold without damage. To approximate the reliability of a subject of interest, one must estimate the corresponding parameters of the probability model. Researchers typically formulate an optimization problem, which is often nonlinear, based on the maximum likelihood theory to obtain estimates for the targeted parameters and then estimate the reliability. Nevertheless, there are additional nonlinear optimization problems in practice from which different estimators for the model parameters are obtained once they are solved numerically. Under normal circumstances, these estimators may perform similarly. However, some might become more robust under irregular situations, such as in the case of data contamination. In this paper, nine frequentist estimators are derived for the parameters of the Laplace Birnbaum-Saunders distribution and then applied to a simulated data set and a real data set. Afterwards, they are compared numerically via Monte Carlo comparative simulation study. The resulting estimates for the reliability based on these estimators are also assessed in the latter study.

1. Introduction

In reality, concrete is the most widely used construction material in the world. Concrete compressive strength is a measure used in determining the amount of resistance a structural element can offer to deformation. Compressive strength is a widely used measure to access the performance of a given concrete mixture. Considering this approach of concrete is vital due to it is the main measure deciding how well concrete can withstand loads that influence its measure. It precisely tells us whether a particular mix is suitable to meet the necessities of a particular venture. Concrete can astoundingly stand up to compressive loading. This is often why it is reasonable for constructing arches, columns, dams, foundations, and tunnel linings among other structures.
Researchers from different fields of science may attempt to describe phenomena of interest, such as concrete compressive strength, using either mathematical or probabilistic models. For instance, scientists describe the life of an object using a probabilistic model called a lifetime probability distribution. Such practice is common in the scientific community, especially in many science fields that involve reliability and reliability analyses. For example, in material science, the two-parameter Birnbaum-Saunders (BS) lifetime distribution can be used by analysts to model fatigue of materials due to periodic cyclic loading [1,2]. This non-negative lifetime model has unimodal skewed probability density and hazard rate curves. Unimodal hazard rate curves are common in practice; see, for example, Langlands et al. [3]. A non-negative continuous random variable T is said to follow the BS distribution if the corresponding cumulative distribution function (CDF) is given by:
F ( t ; α , β ) = Φ 1 α t β β t , t > 0 , α , β > 0 ,
where α is the shape parameter, β is the scale parameter, and Φ ( · ) is the CDF of the standard normal distribution. Desmond [4] provided a more general derivation for the BS distribution assuming a biological model which reinforced the physical rationalization for the use of the BS distribution by relaxing the original presumptions of [1,2]. The BS distribution has desirable aspects and a close relationship to the normal distribution. Consequently, at least a couple of hundred papers and a single research monograph have already appeared describing all properties and developments of this distribution; see, for example, the comprehensive review by Balakrishnan and Kundu [5] in this connection. Examples of recent applications of the BS distribution are such as Bourguignon et al. [6], Hassani et al. [7], and Kannan et al. [8], among others. The BS distribution belongs to a generalized family of distributions called the generalized BS distribution [9]. The generalization is obtained by replacing the Gaussian kernel in Equation (1) with kernels of symmetrical distributions such as the Laplace and logistic distributions. Sampling plans from truncated life tests assuming the generalized BS distribution were developed in [10], while the generalized BS distribution was used to analyze air pollutant concentration in [11]. For further details in this connection, see [12]. The Laplace BS (LBS) distribution is a BS distribution based on Laplace kernel. Its properties and some associated estimation methods were studied by Zhu and Balakrishnan [13]. This paper expands upon their work by discussing additional methods to estimate the parameter of the LBS distribution assuming data contamination. A positive random variable is said to follow the LBS distribution if the corresponding CDF, reliability function, and the probability density functions are given by:
F ( t ; α , β ) = 0.5 + 0.5 sgn ( t β ) 1 exp sgn ( t β ) A ( t ; β ) α ,
S ( t ; α , β ) = 0.5 0.5 sgn ( t β ) 1 exp sgn ( t β ) A ( t ; β ) α ,
and
f ( t ; α , β ) = a ( t ; β ) 4 α β t exp sgn ( t β ) A ( t ; β ) α , x , α , β > 0 ,
respectively, such that α is a positive shape parameter, while β is a positive scale parameter, where:
A ( t ; β ) = t β β t
and
a ( t ; β ) = t + β β t .
The LBS distribution and other lifetime models are characterized by essential statistical properties such as hazard (failure) rates and reliability. The latter concept represents the probability of a specimen continuing to exist for a certain amount of time without loss. Approximating the aspects of the life of objects of interest (e.g., reliability) is associated with model parameters estimation. Estimating the probability distribution parameters has been of great interest to scientists and has received much attention in the statistical literature. In practice, one can obtain various estimators for the model parameters. Therefore, several researchers conducted comparative Monte Carlo simulation studies to numerically assess the estimators from different statistical and computation perspectives; see, for example, Gupta and Kundu [14], Alkasasbeh and Raqab [15], Mazucheli et al. [16], do Espirito Santo and Mazucheli [17], and Balakrishnan and Alam [18], among other research contributions.
Due to the advancement of human civilization, researchers currently deal with various types and large amounts of data originated from the targeted phenomena. In practice, there is no guarantee that part of the data may be missing or contaminated with unusual information values called outliers or extremes. In both cases, estimation efficiency is negatively impacted since the quality of data became questionable. Under data contamination, the maximum likelihood method is not robust since the deviations caused by existing outliers can negatively affect the likelihood. Data contamination motivated many researchers to propose alternative estimators to those obtained by the maximum likelihood theory for various distributions; see, for example, Lawson et al. [19], Boudt et al. [20], Agostinelli et al. [21], Wang et al. [22], among other papers.
The aim of this paper is twofold. Firstly, the LBS distribution is fitted to a real concrete compressive strength data using nine estimation methods. Secondly, the performances of the resulting estimators from the considered estimation methods; namely, modified moments estimators (MMEs), maximum likelihood estimators (MLEs), the maximum product of spacings estimators (MPSEs), the least-squares estimators (LSEs), the weighted least-squares estimators (WLSEs), the percentile estimators (PCEs), the Cramér-von Mises estimators (CVMEs), the Anderson-Darling estimators (ADEs), and the right-tailed Anderson-Darling estimators (RADEs) are investigated using numerical applications and Monte Carlo simulations. The remaining parts of this article are organized as follows. Section 2 discusses the nine estimation methods of interest. Section 3 illustrates the practical application of the discussed estimators using a simulated data set and a real data set. Section 4 reports the outcomes of extensive Monte Carlo simulation study to compare the performance of each estimator under different settings. Finally, remarks and future research are used to conclude this paper in Section 5.

2. Reliability and Model Parameters Estimation

In this section, eight nonlinear optimization problems are formulated to obtain eight estimators for the α and β parameters mentioned in the introductory section. The ninth estimators are the MMEs which have closed-form expression [13]. Before establishing the targeted optimization problems, one must discussing some computational considerations to solve them.

2.1. Computational Considerations for the Optimization Process

The first important computational consideration is finding suitable starting values to solve the optimization problems. One of the important and practical aspects of the LBS distribution is that the scale parameter β is actually the median of the population. Hence, a reasonable starting value for this parameter is the sample median, i.e., β ^ ( 0 ) = median ( t 1 , , t 2 ) . Regarding the starting value for the shape parameter α , one can acquire such value by using the relationship between the LBS distribution and the standard Laplace distribution. In fact, if T follows the LBS distribution with model parameters α and β , then:
α Z = T β β T
follows the Laplace distribution with the location parameter equals to 0, and the scale parameter equals to α . Hence, a starting value for α is obtained by determining the sample Z values; say, z 1 , , z n , based on the observed data t 1 , , t n and the corresponding sample median β ^ ( 0 ) , and then finding the mean absolute deviation (MAD) of z 1 , , z n as a starting value for α . That is,
α ^ ( 0 ) = 1 n i = 1 n z i median ( z 1 , , z 2 ) ,
such that
z i = t i β ^ ( 0 ) β ^ ( 0 ) t i .
The second computational consideration is the choice of optimization algorithm from which one can obtain the estimators for the LBS distribution. Zhu and Balakrishnan [13] showed that the probability density function of the LBS distribution Equation (4) is continuous, but it is not differentiable at β . Consequently, estimators like the MLEs require an optimization method free of derivatives to be obtained, such as the Nelder-Mead algorithm [23]. Note that the latter algorithm is utilized to acquire the remaining estimators to avoid any potential bias in the competitive computations.

2.2. Maximum Likelihood Estimation

Zhu and Balakrishnan [13] have obtained MLEs for α and β and prove their existence and uniqueness. One approach to obtain such estimators is by solving the following maximization algorithm:
maximize n log α + i = 1 n log a ( t i , β ) 1 α i = 1 n A ( t i , β ) s . t . α , β R + 2 ,
such that t 1 , , t n are the observed random sample.

2.3. Least-Squares-Based Estimations

Given an observed random sample t 1 , , t n from the LBS distribution with parameters α and β . Suppose that t 1 : n < , t n : n are the corresponding observed sample order statistics, and consider the following minimization problem:
minimize i = 1 n w i 0.5 + 0.5 sgn ( t i : n β ) 1 e sgn ( t i : n β ) A ( t i : n ; β ) α i n + 1 2 s . t . α , β R + 2 ,
Hence, the solutions of the above optimization problem are the LSEs for α and β given that w 1 = = w n = 1 . However, if w i = [ ( n + 1 ) 2 ( n + 1 ) ] / [ i ( n i + 1 ) ] , then the solutions of the above minimization problem are the WLSEs of α and β [24].

2.4. Percentile Estimation

PCEs of α and β are acquired by fitting a linear model to the theoretical percentiles and the sample percentiles [25,26]. This method requires closed-form cumulative distribution and quantile functions. In the case of the LBS distribution, the cumulative distribution function was given by Equation (2), while the quantile function is defined as:
Q ( u ; α , β ) = β 4 α Q L ( u ) + 4 + α Q L ( u ) 2 2 ,
such that 0 < u < 1 , α , β > 0 , and Q L ( u ) is the quantile function of the standard Laplace distribution, i.e., Q L ( u ) = sgn ( u 0.5 ) log ( 1 2 | u 0.5 | ) . The PCEs of α and β are the solutions for the following minimization problem:
minimize 1 n i = 1 n t i : n β 4 α Q L ( p i : n ) + 4 + α Q L ( p i : n ) 2 2 2 s . t . α , β R + 2 ,
where t 1 : n , , t n : n are the observed sample order statistics, and p i : n = i / ( n + 1 ) .

2.5. Maximum Product of Spacing Estimation

Another estimators of α and β that depends on solving a maximization problem are the MPSEs. Recent research indicates that such estimates compete with MLEs in terms of estimation efficiency and asymptotic properties [27,28,29]. Given the observed sample order statistics t 1 : n < < t n : n , then the MPSEs for the model parameters are determined numerically by solving the following maximization problem:
maximize 1 n + 1 i = 1 n + 1 log Δ i s . t . α , β R + 2 ,
such that
Δ i = 0.5 + 0.5 sgn ( t 1 : n β ) 1 e sgn ( t 1 : n β ) A ( t i : n ; β ) α if i = 1 0.5 sgn ( t i : n β ) 1 e sgn ( t i : n β ) A ( t i : n ; β ) α sgn ( t i 1 : n β ) 1 e sgn ( t i 1 : n β ) A ( t i 1 : n ; β ) α if 1 < i n 0.5 0.5 sgn ( t n : n β ) 1 e sgn ( t n : n β ) A ( t n : n ; β ) α if i = n + 1 .

2.6. Goodness-of-Fit Estimations

The remaining three estimators are based on the idea of minimizing goodness-of-fit statistics, i.e., minimizing the difference between the estimated cumulative distribution function and an empirical counterpart. Examples of such statistics are the Cramér-von Mises, the Anderson-Darling, and the right-tailed Anderson-Darling statistics. The CVMEs of α and β are obtained by evaluating the following minimization problem:
minimize 1 12 n + i = 1 n 0.5 + 0.5 sgn ( t i : n β ) 1 e sgn ( t i : n β ) A ( t i : n ; β ) α 2 i 1 2 n 2 s . t . α , β R + 2 ,
On the other hand, ADEs and RADEs of α and β are acquired as solutions for the following minimization problems:
minimize n 1 n i = 1 n ( 2 i 1 ) log 0.5 + 0.5 sgn ( t i : n β ) 1 e sgn ( t i : n β ) A ( t i : n ; β ) α 1 n i = 1 n ( 2 i 1 ) log 0.5 0.5 sgn ( t n i + 1 : n β ) 1 e sgn ( t n i + 1 : n β ) A ( t i : n ; β ) α s . t . α , β R + 2 ,
and
minimize n 2 2 i = 1 n 0.5 + 0.5 sgn ( t i : n β ) 1 e sgn ( t i : n β ) A ( t i : n ; β ) α 1 n i = 1 n ( 2 i 1 ) log 0.5 0.5 sgn ( t n i + 1 : n β ) 1 e sgn ( t n i + 1 : n β ) A ( t i : n ; β ) α s . t . α , β R + 2 ,
respectively.

3. Numerical Applications

In this section, a simulated data is analyzed; afterwards, a real data is analyzed for the sake of illustration.

3.1. Simulated Data Analysis

Suppose 15 random observations are generated from the LBS distribution with α = β = 1 as shown in Table 1.
The simulated data in Table 1 are obtained according to the following algorithm:
  • Generate a random sample U 1 , , U n from the standard uniform distribution (i.e., U i Uniform ( 0 , 1 , i .)
  • For i = 1 , , n , obtain the desired simulated random sample X 1 , , X n from the LBS distribution with model parameters α and β by using Equation (7), i.e.,
    X i = β 4 α Q L ( U i ) + 4 + α Q L ( U i ) 2 2 ,
    such that Q L ( U i ) = sgn ( U i 0.5 ) log ( 1 2 | U i 0.5 | ) .
Before obtaining the estimates for α and β , one must check their existence and uniqueness. Mathematically proving these requirements is beyond the scope of this study; nevertheless, one may prove them using graphical devices. Using extensive Monte Carlo simulations, a three-dimensional (3D) profile plot for each objective function in the preceding section is established, as shown in Figure 1. The 3D charts clearly indicate that there are areas in which global extrema exist and are expected to be unique for each objective function. The MMEs are obtained first for α and β as shown in [13], and then the remaining estimators are obtained successively as shown in Table 2. The outcomes of the latter table are obtained assuming no contamination in the data, and assuming contamination in the upper 20% order statistics. On the other hand, based on the obtained estimates under both assumptions, Table 3 provides the actual reliability probabilities calculated from the complement of Equation (2) vs. their approximated counterparts using different estimators. Both true and approximated reliability probabilities are evaluated at the sample minimum, the sample three quartiles ( Q 1 , Q 2 , Q 3 ), and the sample maximum. From the previous tables, one can easily observe that the PCEs for α and β provided the farthest approximations for the model parameters and the reliability probabilities. The performances of the remaining estimates are further assessed later.

3.2. Real Data Analysis

To illustrate the application of the considered estimation methods in practice, the concrete compressive strength data of [30] is considered for analysis. This data set established from 17 different sources to check the reliability of a suggested strength model. The data gathered concrete comprising cement alongside fly ash, blast furnace slag, and superplasticizer. The data set consisted of a single response variable; namely, the compressive strength of concrete (in MPa), and 8 covariates. Using this data, various estimated models are obtained and compared by the means of Kolmogorov-Smirnov (KS) test. The latter one-sample testing procedure is used to test the null that the distribution function of a given data set is that of the probability distribution of interest. To obtain the KS statistic, one must consider the following steps:
  • Obtain the estimates of the parameters α and β , denoted by α ^ and β ^ .
  • Compute u i : n = F ( x i : n , α ^ , β ^ ) , such that x i : n is the observed ith sample order statistics, where i = 1 , , n . Here, F ( · , · , · ) is given by Equation (2).
  • Calculate the value of KS statistic as follows:
    KS = max 1 i n u i : n i 1 : n n , i n u i : n ,
    and accordingly calculate the p-value to make a decision about the hypotheses.
Furthermore, since ties exist and the model parameters were estimated, the p-values for the KS statistics were obtained using B = 1 , 000 parametric bootstrapping samples. The steps to obtain the bootstrapping p-value for KS test are as follow:
  • For each method, obtain the estimates of the model parameters α and β ; say, α ^ and β ^ .
  • Use the estimates in the previous step and the algorithm in the preceding section to generate a random sample X 1 * , , X n * from the LBS distribution with shape parameter α ^ and scale parameter β ^ .
  • Compute the KS statistics for each bootstrap sample as discussed before, i.e., repeat Steps 2 and 3, B times to obtain KS ( 1 ) , , KS ( B ) .
  • Calculate the p-value as follows:
    p - value * = 1 B j = 1 B I j ,
    where I j is an indicator function, such that I j = 1 if KS ( j ) > KS , and zero otherwise, for j = 1 , , B and KS is the KS statistic obtained from the original data set.
Finally, it is important to mention that since ties were observed in the data, the MPSEs cannot be acquired directly. When ties exist, one may use a generalization of the maximum product of spacings method to obtain the required estimators; see Murage et al. [31] for additional details. Alongside the estimated parameters and the goodness-of-fit statistics, the reliability probability is calculated at 17 MPa, 28 MPa, and 70 MPa by substituting these values in Equation (3) and replacing the model parameters with the corresponding estimates. In practice, concrete compressive strength can fluctuate between 17 MPa and 28 MPa for residential concrete, while in can be higher as 70 MPa in the case of commercial constructions [32].
Table 4, Table 5, Table 6 and Table 7 respectively summarize the analyses of the considered data set assuming no contamination in the data, 20% of upper data contamination (i.e., the upper 20% of order statistics are multiplied by 5), 20% of lower data contamination (i.e., the lower 20% of order statistics are divided by 5), and 40% two-tailed data contamination which is a mixture of the previous data contamination cases. From the latter tables, one can note the following observations based on the values of the KS test statistics:
  • When there is no data contamination, both MLEs and MPSEs performed well in terms of goodness-of-fit.
  • In the case of upper data contamination, ADEs outperformed both MLEs and MPSEs which took second and third place, respectively.
  • On the other hand, both MLEs and MPSEs maintained their performance followed by ADEs in the case of lower data contamination.
  • In contrast, WLSEs have perform better than MPSEs and MLEs when two-tailed data contamination exists.
  • Overall, MLEs and MPSEs provided the best results in terms of goodness-of-fit, and they both have endured data contamination unlike their counterparts. This is most likely due to the fact that the sample size is large (1000+ units). Furthermore, PCEs and MMEs did not perform well among compared to their counterparts in all considered settings.
  • Finally, according to the reliability proportions estimated by MLEs and MPSEs, one can conclude that the sampled specimens of [30] were suitable for residential buildings.

4. Simulation Outcomes

This section presents the outcomes of Monte Carlo simulation experiments based on 1000 random samples from the LBS distribution with different combinations of values for the shape parameter and sample sizes assuming the following scenarios:
  • Model 1: A model with no contamination.
  • Model 2: A model with 10% of severe upper contamination, i.e., the upper 10% of order statistics are multiplied by 5.
  • Model 3: A model with 10% of severe lower contamination, i.e., the lower 10% of order statistics are multiplied by 1/5.
  • Model 4: A model with 20% of severe two-tailed contamination, i.e., the upper 10% of order statistics are multiplied by 5, while the lower 10% of order statistics are multiplied by 1/5.
For each scenario, the simulation study assumes n = 10 ( 10 ) 100 , α = 0.5 ( 0.5 ) 2.0 , and β = 1 , without loss of any generality. To measure estimation efficiency, the simulated bias and simulated root mean-squared-error (RMSE) are calculated as
Bias α = 1 N i = 1 N α ^ i α , Bias β = 1 N i = 1 N β ^ i β ,
RMSE α = 1 N i = 1 N α ^ i α 2 and RMSE β = 1 N i = 1 N β ^ i β 2 ,
such that N = 1000 , while α ^ i ( β ^ i ) is an estimate of the model parameter α ( β ) based on simulation repetition i. Furthermore, to measure the goodness-of-fit of the fitted model parameters based on the nine estimators, the average absolute difference between the true and estimated reliability function ( D abs ), and the maximum absolute difference between the true and estimated reliability function ( D max ) are determined as
D abs = 1 n × N i = 1 N j = 1 n S ( t j ; α , β ) S ( t j ; α ^ i , β ^ i ) ,
and
D max = 1 N i = 1 N max j = 1 , , n S ( t j ; α , β ) S ( t j ; α ^ i , β ^ i ) ,
respectively, such that S ( t ; α , β ) = 1 F ( t ; α , β ) , while α ^ i ( β ^ i ) is an estimate of the model parameter α ( β ) based on simulation repetition i.
From a statistical perspective, an estimator is computationally consistent when its simulated bias tends to 0 as the sample size increases. Furthermore, when the simulated bias neither increases nor decreases when data contamination exists, then one can conclude that the estimator is computationally robust. Here, Figure 2 and Figure 3 clearly indicate that the most consistent and robust estimators for the model parameters α and β are the MPSEs and CVMEs regardless of the sample size and the true value of α .
In computational statistics, an estimator is computationally efficient when the simulated RMSE tends to 0 as the sample size increases regardless of the existence of data contamination. Figure 4 and Figure 5 suggest that the MPSEs, CVMEs, and the LSEs of α and β are the most efficient estimators compared to the other ones regardless of the simulation settings.
In goodness-of-fit analysis, whenever a pair of estimators yields D abs and D max that tends to 0 as the sample size increases, and are not negatively affected by data contamination, then this pair of estimates provides the least difference between the true and estimated reliability function. This is very important in practice since the aim is to find the best approximation for the reliability. Figure 6 and Figure 7 again indicate that MPSEs, CVMEs, and the LSEs of α and β are the estimators that performed well in terms of goodness-of-fit.

5. Conclusions

In this paper, the estimation problem of the parameters and the reliability function of the Laplace Birnbaum-Saunders lifetime distribution is considered. Besides the method of maximum likelihood, eight classical frequentist estimation methods have been discussed for this purpose; namely, modified moments, maximum product of spacings, least-squares, weighted least-squares, percentile, Cramér-von Mises, Anderson-Darling and right-tailed Anderson-Darling estimation methods. Based on the assumption that the invariance property is exist for the different estimation methods, the reliability function is also estimated using the different estimation methods. To compare the performance of the different estimators a Monte Carlo simulation study is conducted. The practical application of the estimators is illustrated by analyzing a simulated data set and one real data set belongs to compressive strength of concrete. Both data analyses and the Monte Carlo simulation study indicated that all methods perform well when there is no contamination in the data. Once there is some contamination in the data, maximum product of spacings, least-squares, and Cramér-von Mises estimates are notably robust compared to the other estimators and the performance of the other method improve as the sample size increases. Data contamination is not the only problem that faces researchers in practice. Data censoring is another practical challenge that needs to be addressed in future research since it negatively impacts estimation efficiency and robustness. Another important research direction is to compare the studied frequentist estimators to Bayesian estimation in terms of performance based on additional real experimental results which are available in literature.

Author Contributions

F.M.A.A. and M.N. contributed equally to this work. Both the authors have read and agreed to the published version of the manuscript.

Funding

This research received no external funding.

Institutional Review Board Statement

Not applicable.

Informed Consent Statement

Not applicable.

Data Availability Statement

No new data were created or analyzed in this study.

Acknowledgments

The authors would like to thank the Editorial Board and two anonymous reviewers for their constructive feedback that improved the final version of the paper.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Birnbaum, Z.; Saunders, S. A new family of life distributions. J. Appl. Probab. 1969, 6, 319–327. [Google Scholar] [CrossRef]
  2. Birnbaum, Z.W.; Saunders, S.C. Estimation for a Family of Life Distributions with Applications to Fatigue. J. Appl. Probab. 1969, 6, 328–347. [Google Scholar] [CrossRef]
  3. Langlands, A.O.; Pocock, S.J.; Kerr, G.R.; Gore, S.M. Long–term survival of patients with breast cancer: A study of the curability of the disease. BMJ 1979, 2, 1247–1251. [Google Scholar] [CrossRef] [Green Version]
  4. Desmond, A. Stochastic models of failure in random environments. Can. J. Stat. 1985, 13, 171–183. [Google Scholar] [CrossRef]
  5. Balakrishnan, N.; Kundu, D. Birnbaum-Saunders distribution: A review of models, analysis, and applications. Appl. Stoch. Model. Bus. Ind. 2019, 35, 4–49. [Google Scholar] [CrossRef] [Green Version]
  6. Bourguignon, M.; Ho, L.L.; Fernandes, F.H. Control charts for monitoring the median parameter of Birnbaum-Saunders distribution. Qual. Reliab. Eng. Int. 2020, 36, 1333–1363. [Google Scholar] [CrossRef]
  7. Hassani, H.; Kalantari, M.; Entezarian, M.R. A new five-parameter Birnbaum–Saunders distribution for modeling bicoid gene expression data. Math. Biosci. 2020, 319, 108275. [Google Scholar] [CrossRef] [PubMed]
  8. Kannan, G.; Jeyadurga, P.; Balamurali, S. Economic design of repetitive group sampling plan based on truncated life test under Birnbaum–Saunders distribution. Commun. Stat. Simul. Comput. 2020, 1–17. [Google Scholar] [CrossRef]
  9. Diáz-Garciá, J.A.; Leiva-Sánchez, V. A new family of life distributions based on the elliptically contoured distributions. J. Stat. Plan. Inference 2005, 128, 445–457. [Google Scholar] [CrossRef]
  10. Balakrishnan, N.; Leiva, V.; López, J. Acceptance Sampling Plans from Truncated Life Tests Based on the Generalized Birnbaum–Saunders Distribution. Commun. Stat. Simul. Comput. 2007, 36, 643–656. [Google Scholar] [CrossRef]
  11. Leiva, V.; Barros, M.; Paula, G.A.; Sanhueza, A. Generalized Birnbaum-Saunders distributions applied to air pollutant concentration. Environmetrics 2008, 19, 235–249. [Google Scholar] [CrossRef]
  12. Sanhueza, A.; Leiva, V.; Balakrishnan, N. The Generalized Birnbaum–Saunders Distribution and Its Theory, Methodology, and Application. Commun. Stat. Theory Methods 2008, 37, 645–670. [Google Scholar] [CrossRef]
  13. Zhu, X.; Balakrishnan, N. Birnbaum-Saunders distribution based on Laplace kernel and some properties and inferential issues. Stat. Probab. Lett. 2015, 101, 1–10. [Google Scholar] [CrossRef]
  14. Gupta, R.D.; Kundu, D. Generalized exponential distribution: Different method of estimations. J. Stat. Comput. Simul. 2001, 69, 315–337. [Google Scholar] [CrossRef]
  15. Alkasasbeh, M.R.; Raqab, M.Z. Estimation of the generalized logistic distribution parameters: Comparative study. Stat. Methodol. 2009, 6, 262–279. [Google Scholar] [CrossRef]
  16. Mazucheli, J.; Louzada, F.; Ghitany, M. Comparison of estimation methods for the parameters of the weighted Lindley distribution. Appl. Math. Comput. 2013, 220, 463–471. [Google Scholar] [CrossRef]
  17. do Espirito Santo, A.; Mazucheli, J. Comparison of estimation methods for the Marshall–Olkin extended Lindley distribution. J. Stat. Comput. Simul. 2014, 85, 3437–3450. [Google Scholar] [CrossRef] [Green Version]
  18. Balakrishnan, N.; Alam, F.M.A. Maximum likelihood estimation of the parameters of Student’s t Birnbaum-Saunders distribution: A comparative study. Commun. Stat. Simul. Comput. 2019, 1–30. [Google Scholar] [CrossRef]
  19. Lawson, C.; Keats, J.; Montgomery, D. Comparison of robust and least-squares regression in computer-generated probability plots. IEEE Trans. Reliab. 1997, 46, 108–115. [Google Scholar] [CrossRef]
  20. Boudt, K.; Caliskan, D.; Croux, C. Robust explicit estimators of Weibull parameters. Metrika 2009, 73, 187–209. [Google Scholar] [CrossRef] [Green Version]
  21. Agostinelli, C.; Marazzi, A.; Yohai, V.J. Robust Estimators of the Generalized Log-Gamma Distribution. Technometrics 2014, 56, 92–101. [Google Scholar] [CrossRef]
  22. Wang, M.; Park, C.; Sun, X. Simple robust parameter estimation for the Birnbaum-Saunders distribution. J. Stat. Distrib. Appl. 2015, 2. [Google Scholar] [CrossRef] [Green Version]
  23. Nelder, J.A.; Mead, R. A Simplex Method for Function Minimization. Comput. J. 1965, 7, 308–313. [Google Scholar] [CrossRef]
  24. Swain, J.J.; Venkatraman, S.; Wilson, J.R. Least-squares estimation of distribution functions in johnson’s translation system. J. Stat. Comput. Simul. 1988, 29, 271–297. [Google Scholar] [CrossRef]
  25. Kao, J.H.K. Computer Methods for Estimating Weibull Parameters in Reliability Studies. IRE Trans. Reliab. Qual. Control. 1958, PGRQC-13, 15–22. [Google Scholar] [CrossRef]
  26. Kao, J.H.K. A Graphical Estimation of Mixed Weibull Parameters in Life-Testing of Electron Tubes. Technometrics 1959, 1, 389–407. [Google Scholar] [CrossRef]
  27. Cheng, R.C.H.; Amin, N.A.K. Maximum Product-of-Spacings Estimation with Applications to the Lognormal Distribution; Technical Report 1; Department of Mathematics, University of Wales IST: Cardiff, UK, 1979. [Google Scholar]
  28. Cheng, R.C.H.; Amin, N.A.K. Estimating Parameters in Continuous Univariate Distributions with a Shifted Origin. J. R. Stat. Soc. Ser. B 1983, 45, 394–403. [Google Scholar] [CrossRef]
  29. Ranneby, B. The Maximum Spacing Method. An Estimation Method Related to the Maximum Likelihood Method. Scand. J. Stat. 1984, 11, 93–112. [Google Scholar]
  30. Yeh, I.C. Modeling of strength of high performance concrete using artificial neural networks. Cem. Concr. Res. 1998, 28, 1797–1808. [Google Scholar] [CrossRef]
  31. Murage, P.; Mungátu, J.; Odero, E. Optimal Threshold Determination for the Maximum Product of Spacing Methodology with Ties for Extreme Events. Open J. Model. Simul. 2019, 7, 149–168. [Google Scholar] [CrossRef] [Green Version]
  32. NRMCA. Concrete in Practice: What, Why and How? CIP 35—Testing Compressive Strength of Concrete; National Ready Mixed Concrete Association: Silver Spring, MD, USA, 2020. [Google Scholar]
Figure 1. 3D plots of the objective functions based on data in Table 1 for (a) MLE, (b) LSE, (c) WLSE, (d) PCE, (e) MPSE, (f) CVME, (g) ADE and (h) RADE. (cont.)
Figure 1. 3D plots of the objective functions based on data in Table 1 for (a) MLE, (b) LSE, (c) WLSE, (d) PCE, (e) MPSE, (f) CVME, (g) ADE and (h) RADE. (cont.)
Crystals 11 00830 g001aCrystals 11 00830 g001b
Figure 2. Simulated biases for the estimators of α .
Figure 2. Simulated biases for the estimators of α .
Crystals 11 00830 g002
Figure 3. Simulated biases for the estimators of β .
Figure 3. Simulated biases for the estimators of β .
Crystals 11 00830 g003
Figure 4. Simulated RMSEs for the estimators of α .
Figure 4. Simulated RMSEs for the estimators of α .
Crystals 11 00830 g004
Figure 5. Simulated RMSEs for the estimators of β .
Figure 5. Simulated RMSEs for the estimators of β .
Crystals 11 00830 g005
Figure 6. Simulated average absolute difference between the true and estimated reliability function.
Figure 6. Simulated average absolute difference between the true and estimated reliability function.
Crystals 11 00830 g006
Figure 7. Simulated maximum absolute difference between the true and estimated reliability function.
Figure 7. Simulated maximum absolute difference between the true and estimated reliability function.
Crystals 11 00830 g007
Table 1. Simulated data from LBS distribution with α = β = 1 .
Table 1. Simulated data from LBS distribution with α = β = 1 .
0.57909022.30479770.81822923.85027556.3717917
0.13144551.05955034.11755311.11460640.9132579
7.87323250.90670381.54539501.16976770.2345327
Table 2. Estimators of α and β based on data in Table 1.
Table 2. Estimators of α and β based on data in Table 1.
No Data ContaminationUpper 20% Data Contamination
Method α ^ β ^ α ^ β ^
MME0.8825371.2363731.4710212.242812
MLE0.9212951.1146061.5337571.169768
LSE1.1484991.2654671.3068021.302364
WLSE1.1566981.2809521.5668341.431192
PCE0.9448541.5222631.9623232.392981
MPSE1.1447151.1093651.8884771.348811
CVME0.9992141.2523081.0471221.262961
ADE1.031281.2703031.7487561.553576
RADE1.0644021.2601562.1202231.493293
Table 3. True reliability vs. approximated reliability based on estimates in Table 2.
Table 3. True reliability vs. approximated reliability based on estimates in Table 2.
No Data Contamination
SourceMin. Q 1 Q 2 Q 3 Max.
α , β 0.9544440.5688250.4485650.1529860.043167
MME0.9776020.6682690.5554440.1715890.044892
MLE0.9692280.621780.50.1582570.042022
LSE0.9555830.6426440.5523520.2247770.080792
WLSE0.9556230.645610.5566970.2286820.083095
PCE0.9813860.7281720.6409780.2337180.071739
MPSE0.9466190.5989420.4979450.1971740.067713
CVME0.9685060.6564950.5550490.1971840.060605
ADE0.9664810.6573230.559580.2060880.066026
RADE0.9630670.65060.5544870.2101020.069557
Upper 10% Data Contamination
SourceMin. Q 1 Q 2 Q 3 Max.
α , β 0.9544440.5688250.4485650.1529860.043167
MME0.9644430.7453420.6921760.4028730.034084
MLE0.9110420.5904190.5155030.2595740.012739
LSE0.9426590.6360590.5562070.2536840.008556
WLSE0.9261490.6393530.5739180.3030640.019867
PCE0.9359470.7094720.6644690.4396850.071761
MPSE0.8918350.6062030.5480990.3190370.031561
CVME0.964750.653410.5562790.2076020.002869
ADE0.9173160.6446170.5868330.3356480.031491
RADE0.8826890.615310.5646430.3528550.048659
Table 4. Statistical analysis assuming no data contamination.
Table 4. Statistical analysis assuming no data contamination.
Source α β KSp-Value Pr ( X > 17 ) Pr ( X > 28 ) Pr ( X > 70 )
MME0.40930.6790.1390.5170.8840.6000.063
MLE0.43733.9430.0640.6020.9000.6780.092
LSE0.46833.7880.0720.5270.8880.6660.102
WLSE0.39433.2610.0740.5880.9120.6770.072
PCE0.26234.0570.1390.3530.9670.7640.030
MPSE0.44433.7190.0650.6580.8960.6710.093
CVME0.46733.7890.0720.5080.8880.6660.101
ADE0.46233.3690.0690.6820.8870.6580.097
RADE0.39034.2110.0750.5170.9200.7010.077
Table 5. Statistical analysis assuming 20% of upper data contamination.
Table 5. Statistical analysis assuming 20% of upper data contamination.
Source α β KSp-Value Pr ( X > 17 ) Pr ( X > 28 ) Pr ( X > 70 )
MME0.85548.9190.2870.4900.8730.6540.116
MLE0.84636.4460.1650.5490.8770.6660.123
LSE0.54134.1100.1940.4780.8770.6590.114
WLSE0.54834.0910.1930.4900.8770.6540.111
PCE0.55570.3140.5280.3640.9250.7450.109
MPSE0.84236.2420.1660.4490.8710.6540.119
CVME0.54034.1100.1940.4970.8780.6590.114
ADE0.95936.8810.1600.5020.8690.6510.121
RADE1.24635.8970.1800.4760.8680.6500.120
Table 6. Statistical analysis assuming 20% of lower data contamination.
Table 6. Statistical analysis assuming 20% of lower data contamination.
Source α β KSp-Value Pr ( X > 17 ) Pr ( X > 28 ) Pr ( X > 70 )
MME0.93217.9760.3830.5160.7830.4660.084
MLE0.90632.2450.1750.4820.8660.6470.122
LSE0.45033.9470.1980.4610.8880.6660.102
WLSE0.36633.8060.1990.4870.9150.6820.071
PCE0.29331.8510.2000.3540.9620.7450.031
MPSE0.92331.8720.1760.5460.8590.6380.124
CVME0.45033.9480.1980.4810.8880.6660.102
ADE1.05029.1620.1810.6040.8640.6360.113
RADE0.37934.4130.1990.4910.9200.7020.076
Table 7. Statistical analysis assuming 40% of two-tailed data contamination.
Table 7. Statistical analysis assuming 40% of two-tailed data contamination.
Source α β KSp-Value Pr ( X > 17 ) Pr ( X > 28 ) Pr ( X > 70 )
MME1.37528.5310.1990.5030.7900.5420.127
MLE1.32533.4000.1650.5980.8410.6320.147
LSE0.54934.1190.1950.4610.8760.6580.115
WLSE1.02433.6830.1560.5250.8730.6510.114
PCE0.56568.7600.5140.3630.9210.7360.108
MPSE1.33133.0620.1630.7470.8370.6270.148
CVME0.54734.1190.1950.4770.8760.6580.115
ADE1.56332.9720.1800.7930.8450.6300.138
RADE1.43633.4890.1750.6090.8670.6490.120
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

Alam, F.M.A.; Nassar, M. On Modeling Concrete Compressive Strength Data Using Laplace Birnbaum-Saunders Distribution Assuming Contaminated Information. Crystals 2021, 11, 830. https://doi.org/10.3390/cryst11070830

AMA Style

Alam FMA, Nassar M. On Modeling Concrete Compressive Strength Data Using Laplace Birnbaum-Saunders Distribution Assuming Contaminated Information. Crystals. 2021; 11(7):830. https://doi.org/10.3390/cryst11070830

Chicago/Turabian Style

Alam, Farouq Mohammad A., and Mazen Nassar. 2021. "On Modeling Concrete Compressive Strength Data Using Laplace Birnbaum-Saunders Distribution Assuming Contaminated Information" Crystals 11, no. 7: 830. https://doi.org/10.3390/cryst11070830

APA Style

Alam, F. M. A., & Nassar, M. (2021). On Modeling Concrete Compressive Strength Data Using Laplace Birnbaum-Saunders Distribution Assuming Contaminated Information. Crystals, 11(7), 830. https://doi.org/10.3390/cryst11070830

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop