Next Article in Journal
Cosine Similarity Measures of (m, n)-Rung Orthopair Fuzzy Sets and Their Applications in Plant Leaf Disease Classification
Next Article in Special Issue
A New Statistical Technique to Enhance MCGINAR(1) Process Estimates under Symmetric and Asymmetric Data: Fuzzy Time Series Markov Chain and Its Characteristics
Previous Article in Journal
Reaction of Triazolic Aldehydes with Diisopropyl Zinc: Chirality Dissipation versus Amplification
Previous Article in Special Issue
The Modified-Lomax Distribution: Properties, Estimation Methods, and Application
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Bayesian Estimation for the Difference between Coefficients of Quartile Variation of Delta-Lognormal Distributions: An Application to Rainfall in Thailand

by
Noppadon Yosboonruang
1 and
Sa-Aat Niwitpong
2,*
1
Department of Statistics, School of Science, University of Phayao, Phayao 56000, Thailand
2
Department of Applied Statistics, Faculty of Applied Science, King Mongkut’s University of Technology North Bangkok, Bangkok 10800, Thailand
*
Author to whom correspondence should be addressed.
Symmetry 2023, 15(7), 1383; https://doi.org/10.3390/sym15071383
Submission received: 5 June 2023 / Revised: 2 July 2023 / Accepted: 5 July 2023 / Published: 7 July 2023
(This article belongs to the Special Issue Symmetry in Statistics and Data Science, Volume 2)

Abstract

:
The coefficient of quartile variation is a valuable measure used to assess data dispersion when it deviates from a normal distribution or displays skewness. In this study, we focus specifically on the delta-lognormal distribution. The lognormal distribution is characterized by its asymmetrical nature and comprises exclusively positive values. However, when these values undergo a logarithmic transformation, they conform to a symmetrical (or normal) distribution. Consequently, this research aims to establish confidence intervals for the difference between coefficients of quartile variation within lognormal distributions incorporating zero values. We employ the Bayesian, generalized confidence interval, and fiducial generalized confidence interval methods to construct these intervals, involving data simulation using RStudio software. We evaluate the performance of these methods based on coverage probabilities and average lengths. Our findings indicate that the Bayesian method, employing Jeffreys’ prior, performs well in low variability, while the generalized confidence interval method is more suitable for higher variability. Therefore, we recommend using both approaches to construct confidence intervals for the difference between the coefficients of the quartile variation in lognormal distributions that include zero values. Furthermore, we apply these methods to rainfall data in Thailand to illustrate their alignment with actual and simulated data.

1. Introduction

Climate change significantly impacts weather patterns worldwide, including in Thailand, where changes in rainfall dispersion have become particularly noticeable [1]. Thailand has a tropical climate characterized by a rainy season from May to October and a dry season from November to April [2]. However, recent years have witnessed irregular and extreme weather conditions, causing disruptions to traditional rainfall patterns. Several researchers have investigated rainfall dispersion in Thailand, such as Kumphon et al. [3], Szyniszewska and Waylen [4], and Thodsan et al. [5]. Additionally, studies on statistical inference for rainfall distribution in Thailand have been reported by Maneerat et al. [6], Khooriphan et al. [7], Yosboonruang et al. [8,9], and Thangjai et al. [10]. Previous research on statistical inference has primarily focused on measuring rainfall dispersion using variance and coefficient of variation. However, when rainfall dispersion is highly skewed, the coefficient of quartile variation (CQV) becomes a more appropriate tool for analyzing this data type.
The CQV, which stands for the quartile coefficient of dispersion, is a statistical measure that assesses the relative dispersion or variability within a dataset. It quantifies the spread of the data concerning its central tendency, represented by the median. A higher coefficient indicates a more significant variability or dispersion in the dataset, indicating that the values are more spread out from the median. Conversely, a lower coefficient suggests less variability, indicating that the values are closer and more tightly clustered around the median. When the data follow a non-normal distribution, the CQV performs better than the coefficient of variation in measuring relative dispersion. Additionally, when outlier values are present, the CQV is a more appropriate measure for quantifying data dispersion [11]. The CQV has been applied in various subject areas. For instance, Hussein and Morgan [12] used the CQV to measure intravertebral density heterogeneity. Marcoulaki et al. [13] employed the CQV to assess dispersion in computer simulation data for designing central pipeline systems. Chatterjee et al. [14] compared land surface temperature and radiant temperature images at notable coal fire locations using the CQV. Antonetti et al. [15] incorporated water temperature simulations into a fish habitat model and measured thermal heterogeneity using the CQV. Furthermore, researchers have examined statistical inferences regarding the quartile coefficient of variation. Bonett [16] constructed confidence intervals for the CQV, which applied to normal and non-normal distributions. Ambati et al. [17] introduced ratio- and regression-type estimators for estimating the CQV in finite populations. Javed et al. [18] proposed a class of ratio estimators for estimating population variance, utilizing the CQV of an auxiliary variable. Altunkaynak and Gamgam [19] recommended the bootstrap method to establish confidence intervals for the CQV in non-normal distributions. Singh et al. [20] and Ahmed and Shabbir [21] identified an error in the mean squared error of Ambati et al. [17] and rectified it using auxiliary information. Ahmed and Shabbir [21] also presented the Rao regression-type estimator for estimating the CQV with an auxiliary variable. In 2022, Eppen et al. [22] proposed naïve, Rao, and regression estimators for estimating the CQV with an auxiliary variable. Singh and Usman [23] expanded upon the methods introduced by Ambati et al. [17] to estimate the CQV for missing data. Furthermore, Yosboonruang et al. [24] developed a confidence interval for the CQV of a zero-inflated lognormal distribution.
The lognormal distribution plays a significant role in climate change studies due to its ability to model skewed data, analyze extreme events, quantify uncertainty, and facilitate econometric analyses. According to the rainfall data, several researchers have reported that the data follow a lognormal distribution with zero values, also known as a delta-lognormal distribution [8,9,25,26,27,28]. The delta-lognormal distribution consists of positive values following a lognormal distribution and actual zero values following a binomial distribution. The lognormal distribution is asymmetrical in shape. Nevertheless, if the values of the lognormal random variable undergo a logarithmic transformation, they conform to a symmetrical distribution, commonly known as the normal distribution. This distribution has attracted significant interest from researchers studying statistical inference related to it. For example, Li et al. [29] presented generalized and fiducial inference approaches for estimating the mean of a lognormal distribution with excess zeros, with the fiducial approach demonstrating superior performance. Wu and Hsieh [30] constructed a generalized confidence interval (GCI) for the mean of the delta-lognormal distribution using an asymptotic generalized pivotal quantity (GPQ), and their method showed excellent performance. Hasan and Krishnamoorthy [31] introduced fiducial confidence intervals and the method of variance estimate recovery (MOVER) to estimate the mean of a delta-lognormal distribution, receiving recognition for the effectiveness of their proposed methods. In their 2022 study, Zhang et al. [32] proposed the fiducial generalized pivotal quantity (FGPQ) and employed MOVER with FGPQ to construct simultaneous confidence intervals for ratios of means in zero-inflated lognormal distributions, which were also highly recognized for their effectiveness. Furthermore, Yosboonruang et al. [9] introduced various methods to estimate the ratio of coefficients of variation of lognormal distributions with excess zeros, including the fiducial generalized confidence interval (FGCI), Bayesian methods, and the Wald and Fieller log-likelihood methods, with the Bayesian method proving to be the most effective. Recently, in 2023, Thangjai et al. [10] employed the FGCI, Bayesian, and bootstrap methods to establish confidence intervals for the ratio of percentiles of delta-lognormal distributions, with the Bayesian method demonstrating superior performance. Various studies have addressed estimating parameters for lognormal distributions with excess zero values. In this article, we focus on estimating the dispersion of a dataset that follows a delta-lognormal distribution. One effective method for estimating this dispersion is the CQV, building upon the methodology proposed by Yosboonruang and Niwitpong [24] to examine and compare the dispersion between two datasets. Specifically, we aim to construct the highest posterior density (HPD) and confidence intervals for the difference between the CQVs of two delta-lognormal distributions.
The following section presents Bayesian approaches utilizing multiple priors, GCI, and FGCI to construct HPD and confidence intervals. Section 3 provides the simulation results and an empirical study. The final section encompasses the discussion and conclusions of the study.

2. Materials and Methods

Let Y i j , i = 1 , 2 , j = 1 , 2 , , n i be random variables from n observations of delta-lognormal distributions denoted by Y i j Δ ( μ i , σ i 2 , δ i 0 ) , where μ i , σ i 2 , and  δ i 0 represent the mean, variance, and probability of zero values, respectively. For  Y i j > 0 , X i j = ln ( Y i j ) N ( μ i , σ i 2 ) follows a lognormal distribution while Y i j = 0 follows a binomial distribution. Let n i 0 and n i 1 be the numbers of zero and positive values, respectively, such that n i = n i 0 + n i 1 . Aitchison [33] derived the mean and variance of Y i j as μ Y = δ i 1 exp ( μ i + σ i 2 / 2 ) and σ Y 2 = δ i 1 exp ( 2 μ i + σ i 2 ) exp ( σ i 2 ) δ i 1 , where δ i 1 is the probability of positive values.
The CQV is a descriptive statistic used to measure the dispersion between data sets that have different units or to compare within data sets that have different mean values. The CQV is defined by the first and third quartiles as follows:
φ i = Q 3 i Q 1 i Q 3 i + Q 1 i ,
where Q 1 i and Q 3 i denote the first and third quartiles of Y i j , respectively. The quartiles are determined according to Hasan and Krishnamoorthy [31] as
Q r i = exp μ i + Φ 1 r i 4 1 δ i 1 δ i 1 σ i , δ i 1 > 1 r i 4 ,
where Φ is the cumulative standard normal distribution. Since this study focuses on the difference between CQVs, it is defined as
Ψ = φ 1 φ 2 .
This study introduces Bayesian and GCI methods to establish HPD and confidence intervals for the difference between CQVs.

2.1. Bayesian Method

Nowadays, research on statistical inferences and applications is focused on the Bayesian approach because it relies on the population distribution to estimate the parameter of interest [34]. In Bayesian inference, the parameters of interest are directly illustrated by the probability distribution, which is defined as random variables [35].
Regarding the unknown parameters of the delta-lognormal distributions, namely δ i 0 , μ i , and  σ i 2 where δ i 0 = 1 δ i 1 , the joint likelihood function can be defined as
L δ i 0 , μ i , σ i 2 y i j i = 1 2 δ i 0 n i 0 1 δ i 0 n i 1 j = 1 n i 1 1 σ i exp 1 2 σ i 2 ln y i j μ i 2 .
Using Equation (4), the Fisher information matrix of parameters θ = δ 1 0 , μ 1 , σ 1 2 , δ 2 0 , μ 2 , σ 2 2 can be derived by taking the second-order derivative of the log-likelihood function:
I θ = diag n 1 δ 1 0 1 δ 1 0 n 1 1 σ 1 2 n 1 1 δ 1 0 2 σ 1 2 2 n 2 δ 2 0 1 δ 2 0 n 2 1 σ 2 2 n 2 1 δ 2 0 2 σ 2 2 2 .
To estimate the difference between CQVs, HPD intervals are constructed based on the posterior distribution, which is updated using the concept of Bayes’ theorem defined as p θ i | y i 1 , y i 2 , , y i n i = p θ i p y i 1 , y i 2 , , y i n i | θ i p y i 1 , y i 2 , , y i n i .
Since the parameters of interest in this study are δ i 0 , μ i , and  σ i 2 , the posterior of these parameters is computed by integrating the likelihood function in Equation (4) with the prior density function for a delta-lognormal distribution, p δ i 0 , μ i , σ i 2 | y i j . Therefore, the posterior density of δ i 0 , μ i and σ i 2 can be derived as follows:
p δ i 0 | y i j = p δ i 0 , μ i , σ i 2 | y i j d μ i d σ i 2 ,
p μ i | σ i 2 , y i j = p δ i 0 , μ i , σ i 2 | y i j d δ i 0 ,
and
p σ i 2 | y i j = p δ i 0 , μ i , σ i 2 | y i j d δ i 0 d μ i ,
respectively.
Furthermore, Bayesian deep learning can also serve as an alternative approach for generating posterior distributions. By incorporating Bayesian inference techniques, this approach provides several advantages, including robust uncertainty estimation and a principled approach to mitigating overfitting issues. Estimating posterior distributions allows a more comprehensive understanding of the uncertainty associated with the model’s parameters, given the observed data. For a more in-depth exploration of the Bayesian deep learning method, we recommend referring to the research conducted by Zhuang et al. [36].
This article selected three prior distributions, namely the normal gamma prior, Jeffrey’s prior, and the uniform prior for the Bayesian method because these prior distributions can return the closed-form solutions of the posterior distributions. But whenever choosing priors in other forms, the posterior distributions may not have closed-form solutions or follow regular distributions. Therefore, it is necessary to find alternative inference methods, such as Bayesian variable sampling or Markov chain Monte Carlo sampling (MCMC), etc. [37].

2.1.1. The Normal-Gamma Prior

Choosing hyperparameters in the normal-gamma prior involves determining appropriate values for the mean and precision parameters of the normal distribution and the shape and rate parameters of the gamma distribution. These hyperparameters play a crucial role in shaping the prior distribution and subsequently influence the posterior distribution in Bayesian inference.
Maneerat et al. [6] utilized the conjugate families proposed by DeGroot [38] for a normal random sample to derive the posterior distribution of parameters in the normal-gamma prior. The posterior distributions of δ i 0 , μ i , and  σ i 2 are as follows: p δ i 0 | y i j B e t a n i 0 + d i , n i 1 + d i , where d i = 1 6 2 + z α 2 2 ; p μ i | σ i 2 , y i j t n i 1 1 μ ^ i , 1 2 i = 1 n i 1 ln y i j μ ^ i 2 n i 1 1 2 n i 1 ; and p σ i 2 | y i j I n v G a m m a n i 1 1 / 2 , 1 2 i = 1 n i 1 ln y i j μ ^ i 2 , respectively.

2.1.2. Jeffreys’ Prior

According to the delta-lognormal distribution, which is a mixture of the lognormal and binomial distributions, the parameters of interest are δ i 0 , μ i , and  σ i 2 . Following the concept of Jeffreys [39], the prior distributions for these parameters can be obtained by taking the square root of the determinant of the Fisher information matrix (Equation (5)), resulting in p θ i = I θ . Specifically, the Jeffreys prior for δ i 0 is p δ i 0 = δ i 0 δ i 1 1 2 . For the lognormal distribution, the prior distributions for μ i and σ i 2 are p μ i exp n i 1 2 σ i 2 μ i μ ^ i 2 and p σ i 2 = σ i 2 [40]. Consequently, the posterior distributions of these parameters can be computed using Equations (6) - (8) as follows: p δ i 0 | y i j B e t a n i 0 + 0.5 , n i 1 + 0.5 , p μ i | σ i 2 , y i j N μ ^ i , σ i 2 / n i 1 , and  p σ i 2 | y i j I n v G a m m a n i 1 1 / 2 , n i 1 1 σ ^ i 2 / 2 .

2.1.3. The Uniform Prior

According to the uniform prior, it represents a constant function of an a priori probability that assigns equal probabilities to all possible values [41,42]. Therefore, the uniform priors for the parameters of the binomial and lognormal distributions are proportional to 1. By integrating the prior density function for a delta-lognormal distribution, the posterior distributions of δ i 0 , μ i , and  σ i 2 can be determined as follows: p δ i 0 | y i j B e t a n i 0 + 1 , n i 1 + 1 , p μ i | σ i 2 , y i j N μ ^ i , σ i 2 / n i 1 , and
p σ i 2 | y i j I n v G a m m a n i 1 2 / 2 , n i 1 2 σ ^ i 2 / 2 .
Using the posteriors of δ i 0 , μ i , and  σ i 2 obtained for each prior, we substitute these posterior distribution into Equation (2) to calculate the difference between CQVs. Subsequently, we construct the HPD intervals for all methods using the HDInterval package in the R statistical program, following the outlined Algorithm 1 below.
Algorithm 1 Steps to construct HPD intervals for the Bayesian method.
   Step 1. Generate Y i j , where i = 1 , 2 and j = 1 , 2 , , n i , from the delta-lognormal distributions.
   Step 2. Compute δ ^ i 0 , μ ^ i , and  σ ^ i 2 .
   Step 3. Generate the posterior densities of δ i 0 , μ i , and  σ i 2 using each prior.
   Step 4. Compute Q r i using Equation (2).
   Step 5. Compute φ i using Equation (1).
   Step 6. Compute Ψ using Equation (3).
   Step 7. Repeat Steps 3–6 for a total of 2000 times.
   Step 8. Construct HPD intervals for Ψ using each prior.
   Step 9. Repeat Steps 1–8 for a total of 10,000 times.

2.2. Generalized Confidence Interval

The concept of the GCI was introduced by Weerahandi [43]. It is based on the GPQs of the parameters of interest. Furthermore, the construction of confidence intervals for the model parameters using the generalized inference method is discussed in [44]. In this context, the random variables Y i j , where i = 1 , 2 and j = 1 , 2 , , n i , follow delta-lognormal distributions. Referring to Equations (1)–(3), the parameters of interest are δ i 1 , μ i , and  σ i . Let y i j , where i = 1 , 2 and j = 1 , 2 , , n i are the observed values of the random variables Y i j . The GPQs for these parameters possess two important properties: (1) the distribution of GPQs is free from all unknown parameters, and (2) the observed values of GPQs do not depend on the nuisance parameter. Following Wu and Hsieh [30], they computed the variance stabilizing transformation of a binomial distribution using the concept of DasGupta [45]. The GPQ for δ i 1 is given by
R δ i 1 g c i = sin 2 arcsin δ ^ i 1 1 2 n i Z i ,
where Z i N 0 , 1 . Moreover, they used the idea of Krishnamoorthy and Mathew [46] to compute the GPQs for μ i and σ i 2 as follows:
R μ i = μ ^ i Z i n i 1 1 σ ^ i 2 n i 1 χ n i 1 1 2 ,
where Z i N 0 , 1 , and 
R σ i = n i 1 1 σ ^ i 2 χ n i 1 1 2 ,
The pivotal quantities R δ i 1 g c i , R μ i , and  R σ i are consistent with the properties of GPQs. Therefore, we can express R Q r i g c i as follows:
R Q r i g c i = exp R μ i + Φ 1 r i 4 1 R δ i 1 g c i R δ i 1 g c i R σ i , R δ i 1 g c i > 1 R δ i 1 g c i 4 .
By substituting the pivotal quantity from Equation (12) into Equation (1), we obtain
R φ i g c i = R Q 3 i g c i R Q 1 i g c i R Q 3 i g c i + R Q 1 i g c i .
Hence, the pivotal quantity for the difference between CQVs is given by
R Ψ g c i = R φ 1 g c i R φ 2 g c i .
Consequently, the  1 α 100 % confidence interval for Ψ can be expressed as
C I Ψ g c i = R Ψ g c i α / 2 , R Ψ g c i 1 α / 2 ,
where R Ψ g c i α / 2 and R Ψ g c i 1 α / 2 represent the 100 α / 2 -th and 100 1 α / 2 -th percentiles of R Ψ , respectively. The steps for constructing confidence intervals using the GCI method are presented in Algorithm 2.
Algorithm 2 Steps to construct confidence interval for the GCI method.
   Step 1. Generate Y i j , where i = 1 , 2 and j = 1 , 2 , , n i , from the delta-lognormal distributions.
   Step 2. Compute the estimates δ ^ i 1 , μ ^ i , and  σ ^ i 2 .
   Step 3. Generate random variables Z i N 0 , 1 and χ n i 1 1 2 .
   Step 4. Compute the pivotal quantities R δ i 1 g c i , R μ i , and  R σ i .
   Step 5. Repeat Steps 3–4 for a total of 2000 times.
   Step 6. Construct the confidence interval for Ψ .
   Step 7. Repeat Steps 1–6 for a total of 10,000 times.

2.3. Fiducial Generalized Confidence Interval

We extended the FGCI method, as proposed by Yosboonruang et al. [24], to handle the construction of the confidence interval for the difference between CQVs in the delta-lognormal distribution. Consequently, the parameters of interest, based on Equation (2), are δ i 1 , μ i , and  σ i . The fiducial quantities for these parameters can be expressed as R δ i 1 f g c i B e t a n i n i 1 , n i 1 + 0.5 , while R μ i and R σ are represented by Equations (10) and (11), respectively. Furthermore, by utilizing Equation (2), we can determine the fiducial quantity for Q r i as follows:
R Q r i f g c i = exp R μ i + Φ 1 r i 4 1 R δ i 1 f g c i R δ i 1 f g c i R σ i , R δ i 1 f g c i > 1 R δ i 1 f g c i 4 .
Substituting R Q r i f g c i into Equation (1), we can obtain the fiducial quantity for φ i :
R φ i f g c i = R Q 3 i f g c i R Q 1 i f g c i R Q 3 i f g c i + R Q 1 i f g c i .
Accordingly, the fiducial quantity for the difference between CQVs, denoted as R Ψ f g c i , can be represented as
R Ψ f g c i = R φ 1 f g c i R φ 2 f g c i .
Consequently, the  1 α 100 % confidence interval for Ψ can be expressed as
C I Ψ f g c i = R Ψ f g c i α / 2 , R Ψ f g c i 1 α / 2 ,
where R Ψ f g c i α / 2 and R Ψ f g c i 1 α / 2 represent the 100 α / 2 -th and 100 1 α / 2 -th percentiles of R Ψ , respectively. Algorithm 3 outlines the procedure for constructing confidence intervals using the FGCI method.
Algorithm 3 Steps to construct confidence interval for the FGCI method.
   Step 1. Generate Y i j , where i = 1 , 2 and j = 1 , 2 , , n i , from the delta-lognormal distributions.
   Step 2. Compute the estimates δ ^ i 1 , μ ^ i , and  σ ^ i 2 .
   Step 3. Generate random variables Z i N 0 , 1 and χ n i 1 1 2 .
   Step 4. Compute the pivotal quantities R δ i 1 f g c i , R μ i , and  R σ i .
   Step 5. Repeat Steps 3–4 for a total of 2000 times.
   Step 6. Construct the confidence interval for Ψ .
   Step 7. Repeat Steps 1–6 for a total of 10,000 times.

3. Results

3.1. Simulation Study

To compare the performance of various methods, including Bayesian methods based on three different priors (normal-gamma prior—B.NG, Jeffreys’ prior—B.J, and uniform prior—B.U), as well as the GCI and FGCI methods, we conducted simulation studies using Monte Carlo simulation in RStudio version 2022.12.0+353. The goal was to compute the coverage probabilities and average lengths of the HPD and confidence intervals. A favorable method would have coverage probabilities close to or greater than the nominal confidence level of 0.95 and the narrowest intervals. The simulation settings were as follows: sample sizes of 5, 10, 15, 30, 50, 100 and 200; mean values following Hasan and Krishnamoorthy [31] as σ 2 / 2 ; probabilities of positive values as 0.80, 0.85, 0.90, and 0.95 (aligned with Equation (2)); and population variances of 1, 2, 3, 5, and 10. We performed 10,000 replicates for all parameter combinations and 2000 replicates for the Bayesian and GCI methods.
The results in Table 1 and Table 2 provide the HPD and confidence intervals for the difference between CQVs. In cases with small sample sizes and low variance of all probabilities of positive values using the B.NG and B.J methods, the coverage probabilities fall below the nominal confidence level of 0.95. However, the coverage probabilities are either close to or exceed 0.95 for other cases. Figure 1 displays the coverage probabilities of GCI for cases with equal sample sizes. It consistently demonstrates coverage probabilities around 0.95, which are lower than the probabilities obtained from other methods when the sample sizes exceed 5, irrespective of the probabilities of positive values. However, for a small sample size of 5, the B.J method outperforms the other methods.
Moving on to Figure 2, which depicts cases with unequal sample sizes, reveals that the coverage probabilities of GCI are either close to or higher than 0.95 for at least one sample size larger than 5, accompanied by probabilities equal to 0.8 or 0.85. Moreover, as the sample size increases and the probability of positive values reaches or exceeds 0.90, the performance of B.J becomes comparable to or better than that of GCI.
In Figure 3 and Figure 4, the Bayesian method based on Jeffrey’s prior exhibits good performance regarding coverage probabilities for cases with small variances. In contrast, the coverage probabilities of GCI are significantly lower than those of other methods for cases with large variances. Furthermore, both the GCI and FGCI methods consistently show similar coverage probabilities. Additionally, the average lengths of intervals show similar results across all methods. Furthermore, it is worth noting that in Figure 1, Figure 2, Figure 3 and Figure 4, the average lengths tend to narrow for all methods as the variances or sample sizes increase.

3.2. Empirical Study

The rainfall data collected from Bang Klam and Singhanakhon, Songkhla, Thailand, between 2018 and 2022 conform to a delta-lognormal distribution. These data, consisting of positive and zero values, are particularly interesting for the study. The Southern-East Coast Meteorological Center provided the monthly rainfall data presented in Table 3. Both datasets include positive and zero values, with the proportion of positive rainfall data in Bang Klam and Singhanakhon being 53 out of 60 and 48 out of 60, respectively. Descriptive statistics were calculated for the data from Bang Klam and Singhanakhon, resulting in the following values: sample sizes ( n 1 , n 2 ) = (60, 60), proportions of positive values ( δ ^ 1 1 , δ ^ 2 1 ) = (0.88, 0.80), means ( μ ^ 1 , μ ^ 2 ) = (89.3717, 137.5783), and variances ( σ ^ 1 2 , σ ^ 2 2 ) = (18,551.1790, 38,863.9817).
The positive values in both rainfall datasets exhibited right-skewness, as shown in Figure 5. The minimum Akaike information criteria (AIC) and the Bayesian information criteria (BIC) were utilized to confirm the distribution of the data. The results in Table 4 demonstrate that both datasets have lower AIC and BIC values for the lognormal distribution, conforming their adherence to this distribution. Additionally, density and normal Q-Q plots of the log-transformed rainfall data, as depicted in Figure 6, further validate their lognormal distribution.
Table 5 displays the 95% HPD and confidence intervals for the difference between CQVs of the monthly rainfall data from Bang Klam and Singhanakhon. The interval lengths of all methods are similar, consistent with the findings of the simulation study. However, it should be noted that the HPD and confidence intervals generated by all methods exhibit conflicting lower and upper bounds, indicating that these intervals contain zero values. Therefore, there is no significant difference in the monthly rainfall dispersion between both areas.

4. Discussion and Conclusions

This study aimed to construct HPD and confidence intervals for the difference between CQVs of delta-lognormal distributions. We proposed the Bayesian method based on three priors: normal gamma, Jeffreys, and uniform priors, along with the GCI and FGCI methods. To evaluate the performance of these methods, we assessed their coverage probabilities and average lengths under various simulation scenarios.
The findings indicate that the Bayesian approach based on Jeffreys’ prior is suitable for cases with slight variances or small sample sizes. Conversely, as the variance increases, the GCI method outperforms the others. However, it is essential to note that this study focuses on quartiles and specifically applies when the probabilities of positive values are 0.8 or higher. Furthermore, we computed the HPD and confidence intervals for the CQVs of rainfall data from two areas, which followed a delta-lognormal distribution. The empirical results align with the findings from the simulation study, demonstrating that the interval lengths are similar across all methods. In conclusion, we recommend using the Bayesian approach based on Jeffreys’ prior and the GCI method to construct HPD and confidence intervals for the difference between CQVs of delta-lognormal distributions.

Author Contributions

Conceptualization, N.Y. and S.-A.N.; methodology, N.Y. and S.-A.N.; software, N.Y.; validation, N.Y. and S.-A.N.; formal analysis, N.Y. and S.-A.N.; investigation, S.-A.N.; resources, N.Y.; data curation, N.Y.; writing—original draft preparation, N.Y.; writing—review and editing, N.Y. and S.-A.N.; visualization, N.Y.; supervision, S.-A.N.; project administration, S.-A.N.; funding acquisition, N.Y. All authors have read and agreed to the published version of the manuscript.

Funding

This work (Grant No. RGNS 64-131) was financially supported by Office of the Permanent Secretary, Ministry of Higher Education, Science, Research and Innovation, Thailand Science Research and Innovation.

Institutional Review Board Statement

Not applicable.

Informed Consent Statement

Not applicable.

Data Availability Statement

The rainfall data from Bang Klam and Singhanakhon, Songkhla, Thailand, were collected by the Southern-East Coast Meteorological Center.

Acknowledgments

This work (Grant No. RGNS 64-131) was supported by Office of the Permanent Secretary, Ministry of Higher Education, Science, Research and Innovation (OPS MHESI), Thailand Science Research and Innovation (TSRI) and University of Phayao.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. UNDP&FAO. Integrating Agriculture in National Adaptation Plans: Thailand Case Study; UNDP: Bangkok, Thailand, 2019; p. 4. [Google Scholar]
  2. Walalite, T.; Dekker, S.C.; Schot, P.P.; Kardel, I.; Okruszko, T.; Wassen, M.J. Nutrients in tropical and temperate rivers and floodplains—Comparison of the Rivers Songkhram (Thailand) and Narew (Poland). Wetlands 2023, 43, 46. [Google Scholar] [CrossRef]
  3. Kumphon, B.; Kaew-Man, A.; Seenoi, P. A Rainfall Distribution for the Lampao Site in the Chi River Basin, Thailand. J. Water Resour. Prot. 2012, 4, 866–869. [Google Scholar] [CrossRef] [Green Version]
  4. Szyniszewska, A.M.; Waylen, P.R. Determining the daily rainfall characteristics from the monthly rainfall totals in central and northeastern Thailand. Appl. Geogr. 2012, 35, 377–393. [Google Scholar] [CrossRef]
  5. Thodsan, T.; Wu, F.; Torsri, K.; Khampuenson, T.; Yang, G. Impact of the assimilation of multi-platform observations on heavy rainfall forecasts in Kong-Chi Basin, Thailand. Atmosphere 2021, 12, 1497. [Google Scholar] [CrossRef]
  6. Maneerat, P.; Niwitpong, S.-A.; Niwitpong, S. A Bayesian approach to construct confidence intervals for comparing the rainfall dispersion in Thailand. PeerJ 2020, 8, e8502. [Google Scholar] [CrossRef] [Green Version]
  7. Khooriphan, W.; Niwitpong, S.-A.; Niwitpong, S.; Smithpreecha, N. Confidence Intervals for the Ratio of Variances of Delta Gamma Distribution with Applications. Axioms 2022, 11, 689. [Google Scholar] [CrossRef]
  8. Yosboonruang, N.; Niwitpong, S.-A.; Niwitpong, S. Bayesian computation for the common coefficient of variation of delta-lognormal distributions with application to common rainfall dispersion in Thailand. PeerJ 2022, 10, e12858. [Google Scholar] [CrossRef]
  9. Yosboonruang, N.; Niwitpong, S.-A.; Niwitpong, S. Confidence intervals for rainfall dispersions using the ratio of two coefficients of variation of lognormal distributions with excess zeros. PLoS ONE 2022, 17, e0265875. [Google Scholar] [CrossRef]
  10. Thangjai, W.; Niwitpong, S.-A.; Niwitpong, S.; Smithpreecha, N. Confidence Interval Estimation for the Ratio of the Percentiles of Two Delta-Lognormal Distributions with Application to Rainfall Data. Symmetry 2023, 15, 794. [Google Scholar] [CrossRef]
  11. Botta-Dukát, Z. Quartile coefficient of variation is more robust than CV for traits calculated as a ratio. Sci. Rep. 2023, 13, 4671. [Google Scholar] [CrossRef]
  12. Hussein, A.I.; Morgan, E.F. The effect of intravertebral heterogeneity in microstructure on vertebral strength and failure patterns. Osteoporos. Int. 2013, 24, 979–989. [Google Scholar] [CrossRef] [Green Version]
  13. Marcoulaki, E.C.; Papazoglou, I.A.; Pixopoulou, N. Integrated framework for the design of pipeline systems using stochastic optimisation and GIS tools. Chem. Eng. Res. Des. 2012, 24, 2209–2222. [Google Scholar] [CrossRef]
  14. Chatterjee, R.S.; Singh, N.; Thapa, S.; Sharma, D. Retrieval of land surface temperature (LST) from landsat TM6 and TIRS data by single channel radiative transfer algorithm using satellite and ground-based inputs. Int. J. Appl. Earth Obs. Geoinf. 2017, 24, 264–277. [Google Scholar] [CrossRef]
  15. Antonetti, M.; Hoppler, L.; Tonolla, D.; Vanzo, D.; Schmid, M.; Doering, M. Integrating two-dimensional water temperature simulations into a fish habitat model to improve hydro- and thermopeaking impact assessment. River Res. Appl. 2023, 39, 501–521. [Google Scholar] [CrossRef]
  16. Bonett, D.G. Confidence interval for a coefficient of quartile variation. Comput. Stat. Data Anal. 2006, 50, 2953–2957. [Google Scholar] [CrossRef]
  17. Ambati, S.K.; Sedory, S.A.; Singh, S. Estimation of coefficient of dispersion using auxiliary information. Commun. Stat.-Simul. Comput. 2018, 47, 1996–2009. [Google Scholar] [CrossRef]
  18. Javed, K.; Jamal, N.; Hanif, M.; Ali, M.; Shahzad, U.; Luengo, A.V.G. Improved estimator of finite population variance using coefficient of quartile deviation. Asian J. Adv. Res. Rep. 2018, 1, 1–6. [Google Scholar] [CrossRef]
  19. Altunkaynak, B.; Gamgam, H. Bootstrap confidence intervals for the coefficient of quartile variation. Commun. Stat.-Simul. Comput. 2019, 50, 2138–2146. [Google Scholar] [CrossRef]
  20. Singh, R.; Vishwakarma, G.K.; Mishra, P. A note on the estimators for coefficient of dispersion using auxiliary information. Commun. Stat.-Simul. Comput. 2020, 49, 2347–2356. [Google Scholar] [CrossRef]
  21. Ahmed, A.; Shabbir, J. On estimation of coefficient of dispersion using the auxiliary information. Commun. Stat.-Simul. Comput. 2021, 50, 3590–3606. [Google Scholar] [CrossRef]
  22. Eappen, C.V.; Sedory, S.A.; Singh, S. Ratio and regression type estimators of a new measure of coefficient of dispersion. Commun. Stat.-Simul. Comput. 2022, 51, 1899–1920. [Google Scholar] [CrossRef]
  23. Singh, G.N.; Usman, M. Estimation of coefficient of quartile deviation in case of missing data. Commun. Stat.-Theory Methods 2022, 51, 8027–8052. [Google Scholar] [CrossRef]
  24. Yosboonruang, N.; Niwitpong, S.-A. Confidence intervals for the coefficient of quartile variation of a zero-inflated lognormal distribution. Emerg. Sci. J. 2021, 5, 457–470. [Google Scholar] [CrossRef]
  25. Fukuchi, H. Correlation properties of rainfall rates in the United Kingdom. IEE Proc. 1988, 135, 83–88. [Google Scholar] [CrossRef]
  26. Shimizu, K. A bivariate mixed lognormal distribution with an analysis of rainfall data. Am. Meteorol. Soc. 1993, 32, 161–171. [Google Scholar] [CrossRef]
  27. Yue, S. The bivariate lognormal distribution to model a multivariate flood episode. Hydrol. Process. 2000, 14, 2575–2588. [Google Scholar] [CrossRef]
  28. Yosboonruang, N.; Niwitpong, S.-A.; Niwitpong, S. Simultaneous confidence intervals for all pairwise differences between the coefficients of variation of rainfall series in Thailand. PeerJ 2021, 9, e11651. [Google Scholar] [CrossRef]
  29. Li, X.; Zhou, X.; Tian, L. Interval estimation for the mean of lognormal data with excess zeros. Stat. Probab. Lett. 2013, 83, 2447–2453. [Google Scholar] [CrossRef]
  30. Wu, W.-H.; Hsieh, H.-N. Generalized confidence interval estimation for the mean of delta-lognormal distribution: An application to New Zealand trawl survey data. J. Appl. Stat. 2014, 41, 1471–1485. [Google Scholar] [CrossRef]
  31. Hasan, M.S.; Krishnamoorthy, K. Confidence intervals for the mean and a percentile based on zero-inflated lognormal data. J. Stat. Comput. Simul. 2018, 88, 1499–1514. [Google Scholar] [CrossRef]
  32. Zhang, Q.; Xu, J.; Zhao, J.; Liang, H.; Li, X. Simultaneous confidence intervals for ratios of means of zero-inflated log-normal populations. J. Stat. Comput. Simul. 2022, 92, 1113–1132. [Google Scholar] [CrossRef]
  33. Aitchison, J. On the distribution of a positive random variable having a discrete probability mass at the origin. Am. Stat. Assoc. 1955, 50, 901–908. [Google Scholar]
  34. Hespanhol, L.; Vallio, C.S.; Costa, L.M.; Saragiotto, B.T. Understanding and interpreting confidence and credible intervals around effect estimates. Braz. J. Phys. Ther. 2019, 23, 290–301. [Google Scholar] [CrossRef]
  35. O’Neill, P.D. A tutorial introduction to Bayesian inference for stochastic epidemic models using Markov chain Monte Carlo methods. Braz. J. Phys. Ther. 2019, 23, 290–301. [Google Scholar] [CrossRef] [PubMed]
  36. Zhuang, L.; Xu, A.; Wang, X.-L. A prognostic driven predictive maintenance framework based on Bayesian deep learning. Reliab. Eng. Syst. Saf. 2023, 234, 109181. [Google Scholar] [CrossRef]
  37. Zhou, S.; Xu, A.; Tang, Y.; Shen, L. Fast Bayesian inference of reparameterized gamma process with random effects. IEEE Trans. Reliab. 2023. [Google Scholar] [CrossRef]
  38. DeGroot, M.H. Optimal Statistical Decisions; McGraw-Hill, Inc.: New York, NY, USA, 1970. [Google Scholar]
  39. Jeffreys, H. An invariant form for the prior probability in estimation problems. Proc. R. Soc. A Math. Phys. Eng. Sci. 1946, 186, 453–461. [Google Scholar]
  40. Ghosh, J.K.; Delampady, M.; Samanta, T. An Introduction to Bayesian Analysis: Theory and Methods; Springer: New York, NY, USA, 2006. [Google Scholar]
  41. Stone, J.V. Bayes’ Rule: A Tutorial Introduction to Bayesian Analysis; Sebtel Press: Sheffield, UK, 2013. [Google Scholar]
  42. Bolstad, T.A.; Curran, J.M. Introduction to Bayesian Statistics; John Wiley & Sons: Hoboken, NJ, USA, 2017. [Google Scholar]
  43. Weerahandi, S. Generalized confidence intervals. J. Am. Stat. Assoc. 1993, 88, 899–905. [Google Scholar] [CrossRef]
  44. Luo, C.; Shen, L.; Xu, A. Modelling and estimation of system reliability under dynamic operating environments and lifetime ordering constraints. Reliab. Eng. Syst. Saf. 2022, 218, 108136. [Google Scholar] [CrossRef]
  45. DasGupta, A. Asymptotic Theory of Statistics and Probability; Springer: New York, NY, USA, 2008. [Google Scholar]
  46. Krishnamoorthy, K.; Mathew, T. Inferences on the means of lognormal distributions using generalized p-Values generalized confidence intervals. J. Stat. Plan. Inference 2003, 115, 103–121. [Google Scholar] [CrossRef]
Figure 1. Comparison of the coverage probabilities and average lengths of proposed methods for cases with varying sample sizes and probabilities of positive values, assuming equal sample sizes.
Figure 1. Comparison of the coverage probabilities and average lengths of proposed methods for cases with varying sample sizes and probabilities of positive values, assuming equal sample sizes.
Symmetry 15 01383 g001
Figure 2. Comparison the coverage probabilities and average lengths of proposed methods for cases with varying sample sizes and probabilities of positive values, assuming unequal sample sizes.
Figure 2. Comparison the coverage probabilities and average lengths of proposed methods for cases with varying sample sizes and probabilities of positive values, assuming unequal sample sizes.
Symmetry 15 01383 g002
Figure 3. Comparison of the coverage probabilities and average lengths of proposed methods for cases with varying variances and probabilities of positive values.
Figure 3. Comparison of the coverage probabilities and average lengths of proposed methods for cases with varying variances and probabilities of positive values.
Symmetry 15 01383 g003
Figure 4. Comparison of coverage probabilities and average lengths of proposed methods for cases with varying variances and sample sizes.
Figure 4. Comparison of coverage probabilities and average lengths of proposed methods for cases with varying variances and sample sizes.
Symmetry 15 01383 g004aSymmetry 15 01383 g004bSymmetry 15 01383 g004c
Figure 5. Density and normal Q-Q plots of monthly rainfall data in Bang Klam and Singhanakhon, Songkhla, Thailand (2018–2022).
Figure 5. Density and normal Q-Q plots of monthly rainfall data in Bang Klam and Singhanakhon, Songkhla, Thailand (2018–2022).
Symmetry 15 01383 g005
Figure 6. Density and normal Q-Q plots of log-transformed monthly rainfall data in Bang Klam and Singhanakhon, Songkhla, Thailand (2018–2022).
Figure 6. Density and normal Q-Q plots of log-transformed monthly rainfall data in Bang Klam and Singhanakhon, Songkhla, Thailand (2018–2022).
Symmetry 15 01383 g006
Table 1. The coverage probabilities and average lengths of two-sided 95% HPD and confidence intervals for the difference between CQVs in delta-lognormal distributions are presented for cases with equal sample sizes.
Table 1. The coverage probabilities and average lengths of two-sided 95% HPD and confidence intervals for the difference between CQVs in delta-lognormal distributions are presented for cases with equal sample sizes.
n 1 , n 2 δ 1 , δ 2 σ 1 2 , σ 2 2 B.NGB.JB.UGCIFGCI
CPALCPALCPALCPALCPAL
5, 50.80, 0.801, 10.93500.83780.93060.84520.96930.93230.95890.86710.96050.8649
2, 20.94050.71030.93550.72500.97250.66980.95500.74800.95780.7408
3, 30.94660.61360.94280.63210.97680.67970.95530.65460.95700.6449
5, 50.96390.48400.95770.50480.98520.53870.95390.52670.95630.5157
10, 100.98080.30290.97870.32240.99290.34190.95650.34390.95770.3328
0.85, 0.851, 10.93810.81470.93310.81540.96710.89590.95580.81960.96060.8292
2, 20.94370.69650.93810.70820.96770.76030.95490.71700.95940.7192
3, 30.95050.60470.94380.62140.97510.66000.95330.63150.95810.6305
5, 50.96620.47400.95860.49520.98390.51940.95690.50630.96050.5019
10, 100.98350.29630.97910.31760.99280.32860.95580.33060.96050.3247
0.90, 0.901, 10.93960.81610.93630.81660.96640.89780.95720.82080.96080.8306
2, 20.94260.69600.93640.70780.96940.76000.95530.71660.96020.7191
3, 30.95190.60730.94240.62370.97640.66280.95620.63390.95980.6331
5, 50.96210.47660.95520.49770.98230.52230.95450.50880.95830.5048
10, 100.98440.29180.98000.31340.99480.32400.95750.32620.96330.3203
0.95, 0.951, 10.94000.81680.93590.81720.96560.89850.95530.82140.96000.8310
2, 20.94540.69750.93770.70930.96940.76110.95740.71820.96140.7205
3, 30.94990.60790.94440.62460.97510.66330.95350.63460.95800.6338
5, 50.96350.47480.95700.49590.98220.52040.95460.50750.95810.5032
10, 100.98130.29560.97860.31700.99110.32780.95660.32980.96040.3239
10, 100.80, 0.801, 10.95600.65760.94940.65710.96440.68570.96460.66920.96750.6678
2, 20.96110.53980.95340.54740.96960.56170.96610.56140.96980.5559
3, 30.96690.44380.96290.45480.97530.46240.96670.46940.96890.4618
5, 50.97920.31590.97460.32860.98640.33010.96680.34290.96960.3342
10, 100.99250.15990.98970.17050.99390.16840.96440.18340.96690.1762
0.85, 0.851, 10.95910.64300.95330.63550.96600.67040.96280.63620.96740.6438
2, 20.95810.53910.94920.54220.96620.56020.96040.54810.96540.5490
3, 30.96560.44920.96040.45760.97510.46720.96190.46570.96780.4632
5, 50.97950.32390.97270.33590.98440.33770.96370.34570.96920.3405
10, 100.99120.16850.98740.17960.99420.17670.96510.18950.97120.1841
0.90, 0.901, 10.95900.63840.95270.62860.96830.66540.96170.62640.96700.6367
2, 20.96100.53640.95230.53850.96990.55780.96190.54230.96770.5452
3, 30.96780.44970.95810.45750.97370.46770.96210.46390.96810.4629
5, 50.97790.32560.97000.33710.98340.33950.96040.34590.96710.3416
10, 100.99210.16990.98780.18090.99500.17830.96290.19030.96910.1853
0.95, 0.951, 10.96100.61550.95200.59580.96890.64360.95630.57840.96490.6029
2, 20.95930.53160.94930.52740.96610.55290.95380.52140.96380.5335
3, 30.96100.45050.95120.45390.97120.46870.95180.45300.96280.4588
5, 50.97170.33240.96160.34180.97840.34620.95220.34560.96180.3458
10, 100.98820.17680.98030.18780.99200.18520.95120.19480.96130.1915
15, 150.80, 0.801, 10.96610.55990.96160.55170.97150.57410.96900.55370.97240.5588
2, 20.96340.46360.95640.46420.96910.47460.96410.47010.96830.4698
3, 30.97000.37720.96210.38240.97610.38670.96700.38970.97020.3868
5, 50.98020.26140.97420.26940.98380.26830.96700.27740.97030.2730
10, 100.99320.12240.99010.12950.99470.12620.96670.13670.97090.1327
0.85, 0.851, 10.96410.55150.95770.54040.96990.56560.96510.53920.97000.5470
2, 20.96590.46000.95920.45930.96960.47110.96410.46290.96890.4646
3, 30.97310.37900.96560.38310.97720.38830.96840.38860.97290.3873
5, 50.97920.26460.97120.27180.98370.27170.96730.27870.97090.2753
10, 100.98950.12670.98630.13370.99280.13060.96350.14030.96830.1367
15, 150.90, 0.901, 10.96870.53420.95890.51820.97180.54910.96130.51100.96840.5243
2, 20.96600.45610.95760.45220.97030.46770.96220.45150.96830.4572
3, 30.96700.37950.95840.38160.97210.38920.95940.38380.96700.3854
5, 50.97910.26890.97110.27520.98290.27600.96360.28010.96960.2784
10, 100.99230.13030.98590.13720.99390.13430.96450.14290.97000.1400
0.95, 0.951, 10.96280.50330.95480.48120.96980.51930.95460.46580.96300.4866
2, 20.96060.44310.95200.43470.96760.45550.95390.42810.96380.4394
3, 30.96210.37550.95370.37420.97010.38580.95290.37220.96290.3782
5, 50.97260.27210.96320.27650.97760.27970.95320.27860.96210.2796
10, 100.98780.13490.97740.14140.98920.13900.95030.14570.95860.1439
30, 300.80, 0.801, 10.97750.45700.97430.45190.97830.46160.97960.45510.98060.4572
2, 20.97770.36180.97390.36240.97960.36560.97860.36770.98040.3665
3, 30.98040.28310.97640.28580.98280.28580.98050.29160.98100.2890
5, 50.98700.18000.98270.18380.98800.18190.97940.18920.98170.1863
10, 100.99520.07090.99280.07380.99550.07170.98100.07760.98410.0755
0.85, 0.851, 10.97570.43620.97290.42760.97840.44110.97560.42700.97850.4326
2, 20.97160.35650.96760.35470.97390.36030.97220.35730.97670.3586
3, 30.97680.28400.97320.28530.97760.28710.97480.28900.97760.2883
5, 50.98020.18590.97770.18920.98190.18800.97270.19340.97500.1914
10, 100.99130.07660.98800.07960.99260.07760.97260.08290.97520.0811
0.90, 0.901, 10.96920.40770.96420.39700.97120.41330.96490.39290.96820.4015
2, 20.97020.34550.96500.34210.97120.34980.96880.34200.97150.3456
3, 30.96900.28130.96430.28120.97130.28460.96710.28300.97130.2839
5, 50.97570.18990.97060.19230.97730.19220.96590.19540.96990.1944
10, 100.98710.08140.98470.08430.98940.08250.96560.08720.96990.0857
0.95, 0.951, 10.96380.35740.95850.34640.96700.36320.95760.34070.96270.3497
2, 20.95840.32280.95200.31780.96180.32770.95420.31550.95970.3208
3, 30.95860.27270.95330.27140.96140.27660.95450.27100.95950.2737
5, 50.96760.19170.96040.19340.97060.19450.95690.19490.96280.1952
10, 100.97870.08770.97290.09040.98080.08900.95610.09260.96100.0916
50, 500.80, 0.801, 10.98250.39790.98050.39400.98290.39970.98360.39740.98450.3987
2, 20.98340.30350.98040.30370.98290.30520.98490.30820.98660.3072
3, 30.98600.22900.98400.23050.98610.23010.98630.23490.98730.2331
5, 50.98810.13690.98450.13910.98760.13760.98260.14290.98390.1408
10, 100.99500.04830.99350.04980.99550.04860.98270.05200.98440.0508
0.85, 0.851, 10.97760.36280.97700.35600.97900.36520.97610.35520.97890.3603
2, 20.97660.29390.97450.29190.97770.29580.97730.29340.97980.2950
3, 30.97990.22940.97720.22960.98070.23080.97870.23200.98050.2319
5, 50.98300.14450.98150.14610.98500.14540.97940.14880.98120.1477
10, 100.99100.05460.98880.05600.99160.05490.98020.05790.98180.0570
0.90, 0.901, 10.97000.31840.96890.31160.97280.32120.96890.30870.97070.3149
2, 20.96880.27460.96330.27160.96820.27660.96550.27130.96920.2743
3, 30.97230.22320.97020.22260.97330.22490.97120.22360.97450.2248
5, 50.97310.14690.96930.14790.97480.14790.96780.14970.97110.1494
10, 100.98240.05930.97830.06070.98320.05980.96730.06230.97010.0615
0.95, 0.951, 10.96290.27170.95770.26670.96440.27450.95900.26500.96200.2692
2, 20.96180.25000.95650.24750.96330.25220.95960.24710.96210.2498
3, 30.95990.21180.95770.21100.96150.21370.95780.21160.96060.2130
5, 50.96030.14760.95780.14840.96170.14880.95280.14970.95570.1498
10, 100.97120.06500.96950.06630.97420.06550.95830.06760.96000.0671
100, 1000.80, 0.801, 10.98770.32890.98710.32620.98800.32980.98710.32870.98850.3300
2, 20.98560.23960.98460.23940.98580.24010.98740.24220.98810.2419
3, 30.98790.17330.98620.17390.98800.17360.98880.17650.98890.1757
5, 50.99170.09620.99020.09720.99090.09650.98910.09930.99010.0983
10, 100.99450.02930.99330.02990.99470.02940.98850.03090.98890.0304
0.85, 0.851, 10.98090.27360.97920.26950.98130.27440.97850.26900.97950.2726
2, 20.97650.22050.97430.21900.97680.22120.97650.21980.97790.2212
3, 30.98120.16940.97950.16910.98130.16990.98050.17030.98160.1708
5, 50.98230.10190.98200.10250.98260.10220.98080.10380.98160.1035
10, 100.98700.03500.98550.03550.98700.03510.97930.03640.98000.0360
0.90, 0.901, 10.96800.22100.96520.21840.96750.22200.96570.21810.96650.2206
2, 20.96460.19450.96390.19340.96620.19530.96680.19370.96840.1951
3, 30.96970.15870.96840.15830.97060.15920.97080.15900.97230.1597
5, 50.97140.10320.97000.10360.97100.10360.96900.10460.97020.1045
10, 100.97840.03980.97710.04030.97910.03990.97050.04100.97120.0407
0.95, 0.951, 10.95970.19160.95770.18980.96100.19260.95880.19010.95960.1916
2, 20.95540.17760.95280.17670.95720.17850.95770.17750.95850.1784
3, 30.95860.15070.95700.15040.95870.15130.95830.15130.95940.1518
5, 50.95870.10400.95670.10430.96060.10440.95620.10530.95840.1053
10, 100.96340.04370.96160.04420.96530.04390.95520.04490.95570.0447
200, 2000.80, 0.801, 10.99190.26750.99070.26550.99200.26770.99110.26730.99140.2686
2, 20.98720.18800.98740.18750.98860.18820.98950.18940.98930.1896
3, 30.98910.13090.98890.13110.98970.13100.99000.13270.99080.1324
5, 50.99150.06830.99120.06870.99200.06840.99130.06980.99170.0694
10, 100.99390.01850.99320.01870.99290.01850.98980.01920.99040.0189
0.85, 0.851, 10.98040.19290.97900.19130.98100.19310.98010.19160.98070.1932
2, 20.97710.15750.97660.15680.97860.15770.97800.15740.97850.1582
3, 30.97960.12120.97870.12090.98080.12130.98040.12180.98110.1221
5, 50.98250.07200.98090.07220.98180.07210.98160.07290.98230.0728
10, 100.98680.02340.98510.02360.98510.02350.98080.02400.98120.0239
0.90, 0.901, 10.96860.15520.96840.15440.96780.15560.96840.15490.97040.1558
2, 20.96350.13770.96500.13720.96510.13800.96570.13810.96690.1386
3, 30.96790.11230.96690.11220.96830.11250.96840.11290.96930.1132
5, 50.96700.07270.96590.07280.96700.07270.96760.07350.96700.0734
10, 100.97540.02710.97360.02730.97580.02720.96910.02770.97070.0276
0.95, 0.951, 10.95820.13470.95680.13410.95720.13500.95880.13480.95930.1353
2, 20.95400.12580.95160.12550.95570.12610.95490.12630.95560.1266
3, 30.95610.10680.95610.10670.95810.10710.95860.10750.95990.1077
5, 50.95980.07340.95750.07350.95920.07360.95820.07420.95980.0742
10, 100.96080.03010.96110.03020.96190.03010.95740.03060.95850.0305
Note: The values shown in bold represent the shortest expected lengths.
Table 2. The coverage probabilities and average lengths of two-sided 95% HPD and confidence intervals for the difference between CQVs in delta-lognormal distributions are presented for cases with unequal sample sizes.
Table 2. The coverage probabilities and average lengths of two-sided 95% HPD and confidence intervals for the difference between CQVs in delta-lognormal distributions are presented for cases with unequal sample sizes.
n 1 , n 2 δ 1 , δ 2 σ 1 2 , σ 2 2 B.NGB.JB.UGCIFGCI
CPALCPALCPALCPALCPAL
5, 150.80, 0.801, 10.94570.71450.94070.71790.96960.77040.96270.73420.96540.7325
2, 20.95180.59350.94740.60400.97330.63710.96310.62200.96490.6156
3, 30.95980.49550.95560.50940.97900.53250.96550.52810.96740.5194
5, 50.96920.37130.96330.38720.98420.40150.95960.40600.96020.3967
10, 100.98720.21240.98480.22650.99400.23300.96480.24590.96690.2374
0.85, 0.851, 10.94760.69460.94330.69160.96990.74430.95910.69600.96300.7027
2, 20.95290.58570.94730.59340.97050.62350.95970.60160.96370.6020
3, 30.96080.49620.95310.50860.97610.52830.96280.51870.96590.5160
5, 50.97110.36910.96240.38500.98580.39490.96160.39670.96670.3920
10, 100.98570.20890.98190.22420.99190.22610.95950.23770.96430.2324
0.90, 0.901, 10.94540.68380.94120.67790.96390.73330.95690.67830.96170.6884
2, 20.94630.58330.93960.58890.96480.62130.95520.59490.96100.5980
3, 30.95290.49600.94720.50700.97050.52810.95760.51520.96380.5146
5, 50.96670.37400.95850.38900.97930.39970.95750.39960.96260.3960
10, 100.98500.21510.97950.23040.99220.23260.95860.24330.96300.2385
0.95, 0.951, 10.93770.66720.93450.65900.95800.71630.95230.65600.95770.6689
2, 20.93730.57540.93150.57960.95820.61350.94880.58300.95480.5882
3, 30.94390.49560.94100.50520.96260.52790.95200.51140.95690.5129
5, 50.95710.37670.94880.39070.96980.40250.95060.40020.95720.3979
10, 100.97750.21750.96990.23230.98710.23490.95090.24460.95610.2403
10, 500.80, 0.801, 10.96530.54160.95720.54080.97210.55810.96640.54960.97270.5490
2, 20.96660.42920.95960.43460.97220.44140.96670.44610.97050.4416
3, 30.97210.34180.96480.34950.97720.35150.96800.36220.97160.3565
5, 50.98270.22830.97590.23700.98590.23550.96490.25060.97060.2444
10, 100.99460.10340.99220.10990.99600.10730.96850.12240.97150.1177
0.85, 0.851, 10.96790.51770.96330.51100.97500.53390.96940.51150.97510.5178
2, 20.96550.42730.95790.42970.97230.43960.96540.43520.97160.4356
3, 30.96840.34790.96080.35380.97370.35770.96360.36150.97040.3596
5, 50.97940.23950.97310.24800.98340.24680.96520.25780.97250.2541
10, 100.99060.11100.98890.11820.99260.11500.96680.12840.97280.1247
0.90, 0.901, 10.95170.49690.95110.48900.96010.51400.95870.48700.96490.4956
2, 20.95400.41900.95100.42030.96120.43130.96280.42400.96790.4260
3, 30.96110.34580.95810.35110.96630.35580.96490.35740.97090.3568
5, 50.96520.24040.96420.24830.97280.24760.96300.25730.96900.2543
10, 100.98120.11450.98100.12160.98440.11850.96550.13120.96990.1278
0.95, 0.951, 10.95000.46880.94780.45460.95690.48660.95400.44560.96040.4621
2, 20.94660.40790.94780.40580.95480.42090.95730.40320.96430.4109
3, 30.95130.34330.94860.34650.95780.35370.95880.34820.96340.3515
5, 50.95290.24620.95240.25330.95940.25360.95420.25900.96080.2585
10, 100.96810.12250.96850.12980.97130.12660.95450.13780.96160.1355
15, 1000.80, 0.801, 10.95630.45610.94160.45030.96110.46440.94500.45280.95460.4561
2, 20.95950.36240.94310.36360.96530.36880.94500.36880.95710.3684
3, 30.96380.28520.95030.28940.97020.29040.94510.29650.95570.2943
5, 50.97740.18610.96700.19170.98180.18950.94340.20050.95560.1973
10, 100.99330.07640.98650.08090.99510.07800.94220.08920.95420.0866
0.85, 0.851, 10.96470.43070.95140.42290.96830.43950.95640.42210.96560.4282
2, 20.96630.35570.95390.35560.96960.36240.95800.35870.96630.3600
3, 30.96850.28560.95910.28900.97160.29090.95870.29440.96550.2934
5, 50.97760.19070.96940.19600.97930.19430.96120.20370.96860.2012
10, 100.98860.08160.98570.08610.99080.08330.95900.09370.96830.0914
15, 1000.90, 0.901, 10.96340.40410.95670.39350.96720.41360.96270.39000.96910.3992
2, 20.95910.34570.95220.34360.96450.35280.95810.34400.96640.3477
3, 30.96170.28470.95700.28690.96710.29030.96010.28990.96960.2909
5, 50.96970.19330.96510.19820.97310.19700.96100.20410.96830.2028
10, 100.98040.08640.97790.09110.98170.08830.96040.09780.96830.0958
0.95, 0.951, 10.95290.37560.94790.36100.95870.38600.94950.35320.95910.3669
2, 20.94930.33460.94440.32950.95380.34250.95260.32560.96040.3331
3, 30.95040.28060.94630.28060.95700.28680.95140.28030.96050.2840
5, 50.95600.19640.95330.20040.95950.20060.95440.20410.96300.2045
10, 100.96770.09210.96770.09660.96980.09410.95660.10230.96250.1009
50, 2000.80, 0.801, 10.98430.34420.98180.34220.98600.34550.98270.34620.98600.3463
2, 20.98340.25130.98040.25190.98340.25220.98200.25620.98450.2549
3, 30.98480.18230.98240.18370.98600.18290.98180.18790.98410.1864
5, 50.98790.10230.98440.10400.98820.10270.98030.10780.98350.1063
10, 100.99340.03210.99230.03310.99460.03230.98120.03540.98360.0346
0.85, 0.851, 10.97910.29050.97640.28600.97910.29200.97600.28670.98050.2903
2, 20.97620.23610.97230.23480.97690.23740.97450.23620.97720.2373
3, 30.97640.18210.97360.18230.97620.18290.97470.18430.97820.1844
5, 50.97980.11170.97750.11280.98040.11210.97470.11540.97750.1147
10, 100.98670.03950.98510.04050.98700.03970.97310.04250.97690.0419
0.90, 0.901, 10.96780.24880.96570.24420.96930.25060.96580.24360.96900.2478
2, 20.96130.21680.96060.21470.96340.21800.96410.21480.96600.2169
3, 30.96390.17600.96200.17560.96490.17700.96440.17640.96730.1773
5, 50.97210.11450.97080.11530.97160.11520.96960.11700.97220.1169
10, 100.97300.04430.97330.04520.97370.04460.96720.04690.97040.0465
0.95, 0.951, 10.95880.21370.95530.21040.96160.21550.95860.20990.96190.2127
2, 20.95530.19760.95020.19590.95630.19920.95510.19610.95830.1978
3, 30.95450.16720.95220.16660.95500.16840.95400.16740.95620.1682
5, 50.96000.11570.95800.11630.96040.11650.95850.11760.96100.1177
10, 100.96390.04890.96160.04980.96450.04930.95620.05120.95910.0509
100, 2000.80, 0.801, 10.98900.30390.98890.30160.98920.30410.98920.30420.98950.3052
2, 20.98770.21590.98730.21570.98780.21630.98810.21840.98910.2181
3, 30.98450.15290.98430.15340.98590.15310.98460.15580.98650.1551
5, 50.98900.08220.98860.08290.98940.08230.98640.08480.98840.0841
10, 100.99440.02340.99320.02380.99410.02340.98820.02470.98930.0243
0.85, 0.851, 10.98320.23620.98020.23330.98240.23700.98020.23370.98160.2363
2, 20.97880.19170.97800.19060.97930.19220.97980.19120.98060.1925
3, 30.98170.14720.97930.14690.98240.14750.98140.14800.98220.1484
5, 50.98130.08780.97920.08820.98120.08800.97950.08940.98160.0892
10, 100.98660.02930.98590.02970.98620.02940.97950.03050.98020.0302
0.90, 0.901, 10.96740.19060.96620.18890.96740.19140.96660.18900.96800.1907
2, 20.97080.16850.96890.16760.97130.16910.97100.16820.97170.1693
3, 30.96850.13740.96670.13710.96770.13780.96810.13790.97000.1384
5, 50.96890.08900.96680.08930.97060.08930.96750.09020.96850.0902
10, 100.97470.03380.97410.03410.97560.03380.96980.03480.97070.0346
0.95, 0.951, 10.95960.16550.95780.16430.96050.16620.96010.16480.96230.1658
2, 20.95750.15390.95550.15320.95780.15450.95950.15400.95980.1547
3, 30.95590.13060.95280.13040.95570.13110.95660.13120.95680.1316
5, 50.95660.08980.95560.09000.95660.09010.95650.09090.95770.0909
10, 100.96500.03730.96290.03770.96290.03740.95800.03830.96020.0381
Note: The values shown in bold represent the shortest expected lengths.
Table 3. Summary of monthly rainfall data for Bang Klam and Singhanakhon, Songkhla, Thailand (2018–2022).
Table 3. Summary of monthly rainfall data for Bang Klam and Singhanakhon, Songkhla, Thailand (2018–2022).
MonthBang KlamSinghanakhon
2022202120202019201820222021202020192018
Jan.69390113.871.1121.349.40146.6237.5
Feb.21102400193.1050.200
Mar.14.3190010141.247.5009.2
Apr.24.4101730.527.5290.421.20045
May82.5112.49.504725028.2012.539.8
Jun.75.878.323.710.54869289.33022.6
Jul.32.752.593567.575.7175.628890
Aug.428.510932.59660.235.810.728.3
Sept.21.7103744.51761.602028.441.9
Oct.202.23871.1122.418331327.6135.3196.9373.3
Nov.353.6500.6442.8236.9109.5344.5738.2729.7150.2383.7
Dec.732.5122.1233.862309.6933.6232543.954.7279.9
Table 4. AIC and BIC values for monthly rainfall data in Bang Klam and Singhanakorn, Songkhla, Thailand (2018–2022).
Table 4. AIC and BIC values for monthly rainfall data in Bang Klam and Singhanakorn, Songkhla, Thailand (2018–2022).
AreasCriteriaDistributions
NormalCauchyLognormalGammaWeibull
Bang KlamAIC677.8336627.4157586.7384598.2644596.5048
BIC681.7742631.3562590.6790602.2050600.4454
SinghanakornAIC651.0329632.7343590.8083593.3023592.6424
BIC654.7753636.4767594.5507597.0447596.3848
Note: The values shown in bold represent the minimum AIC and BIC.
Table 5. The 95% HPD and confidence intervals for the difference in CQVs of monthly rainfall data between Bang Klam and Singhanakorn, Songkhla, Thailand.
Table 5. The 95% HPD and confidence intervals for the difference in CQVs of monthly rainfall data between Bang Klam and Singhanakorn, Songkhla, Thailand.
MethodsLower BoundUpper BoundLength
B.NG−0.24890.05910.3080
B.J−0.24720.07750.3247
B.U−0.26200.05450.3165
GCI−0.25420.06610.3203
FGCI−0.25600.07430.3303
Disclaimer/Publisher’s Note: The statements, opinions and data contained in all publications are solely those of the individual author(s) and contributor(s) and not of MDPI and/or the editor(s). MDPI and/or the editor(s) disclaim responsibility for any injury to people or property resulting from any ideas, methods, instructions or products referred to in the content.

Share and Cite

MDPI and ACS Style

Yosboonruang, N.; Niwitpong, S.-A. Bayesian Estimation for the Difference between Coefficients of Quartile Variation of Delta-Lognormal Distributions: An Application to Rainfall in Thailand. Symmetry 2023, 15, 1383. https://doi.org/10.3390/sym15071383

AMA Style

Yosboonruang N, Niwitpong S-A. Bayesian Estimation for the Difference between Coefficients of Quartile Variation of Delta-Lognormal Distributions: An Application to Rainfall in Thailand. Symmetry. 2023; 15(7):1383. https://doi.org/10.3390/sym15071383

Chicago/Turabian Style

Yosboonruang, Noppadon, and Sa-Aat Niwitpong. 2023. "Bayesian Estimation for the Difference between Coefficients of Quartile Variation of Delta-Lognormal Distributions: An Application to Rainfall in Thailand" Symmetry 15, no. 7: 1383. https://doi.org/10.3390/sym15071383

APA Style

Yosboonruang, N., & Niwitpong, S. -A. (2023). Bayesian Estimation for the Difference between Coefficients of Quartile Variation of Delta-Lognormal Distributions: An Application to Rainfall in Thailand. Symmetry, 15(7), 1383. https://doi.org/10.3390/sym15071383

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