Updated fiducial distribution of parameters in the associated delta-lognormal population

In this paper we consider a special kind of semicontinous distribution. We try to concern with the situation where the probability of zero observation is associated with the location and scale parameters in lognormal distribution. We first propose a goodness-of-fit test to ensure that the data can be fit by the associated delta-lognormal distribution. Then we define the updated fiducial distributions of the parameters and establish the results that the confidence interval has asymtotically correct level while the significance level of the hypothesis testing is also asymtotically correct. We propose an exact sampling method to sample from the updated fiducial distribution. It can be seen in our simulation study that the inference on the parameters is largely improved. A real data example is also used to illustrate our method.


Introduction
In real applications, such as fisheries research and medical cost analysis, the response variables may be skewed, non-negative and have a non-negligible probability of zero outcomes.These variables are also known as following the semicontinous distribution Fðx; dÞ ¼ d; x ¼ 0; ( where G(x) is the cumulative distribution of a postive variable.In existing researches, δ is usually assumed to be independent of G(x).However, we think that the probability δ is associated with G.For example, consider the precipitation distributions of some certain areas, the areas with larger rainfall per year are more likely to have less dry days.Hence, we can assume that δ is associated with G, say δ = G(a), for some a.In this paper, we try to deal with such assumption by a specified distribution, the delta lognormal distribution.This kind of distributions is first discussed and named by [1].The cumulative distribution function of delta-lognormal distribution is then defined as follow denotes the cumulative distribution function of a lognormal distribution.Since log X follows a normal distribution N(μ, σ), we still refer to μ and σ as the location and scale parameters respectively in the rest of our paper.[2] applied this distribution to deal with the measurement of worker exposure to air contaminants in United States.The use of delta-lognormal distribution to fisheries research was done by [3][4][5].They considered the estimates of the population mean of the delta-lognormal distribution and further studied their robustness.It is easy to calculate the mean of the deltalognormal distribution as Much attention is given to the confidence interval of M by various statisticians.[6] proposed to use the likelihood ratio test to get a better control of the Type I error than the former standard ANOVA F-test and Kruskal-Wallis test.A Bootstrap approach is proposed which is proved to be second-order accurate in [7].[8] considered the case when at least two non-zero observations are observed and modified the profiled loglikelihood function.[9,10] used the generalized pivotal quantities proposed by [11] to construct a generalized pivot for estimating the mean.In their paper, a Beta distribution is used as the generalized pivot for δ.This thought is further developed by [12][13][14].In the papers mentioned above, generalized pivot quantities are proposed for the binomial variable, which is discrete.Meanwhile, the conclusion of [15] on generalized fiducial inference also motivates some new ideas.The recent results are shown in [14], where the authors focus mainly on the improvement of the Beta distribution to approximate the generalized fiducial distribution of δ.
Instead of finding a generalized fiducial distribution, another method is proposed by [16], named "method of variance estimates recovery"(MOVER).This method can be easily applied to many different settings while guarantees the coverage probability of the confidence interval.From the Bayesian perspective, [17] compared the performance of using different prior distributions for both lognormal distribution and delta-lognormal distribution.They further considered the comparison of the means of two lognormal population.
As we can see from the introduction above, the three parameters in delta-lognormal distribution is assumed to be independent.However, in real applications, the probability of zero outcomes may be associated with the location and scale parameters.Consider the case of the spend on children's clothing in [1], a family in a rich community is more likely to be a spender, while the one in a relative poor community may be a nonspender, since it is easy to be influenced by other families in the same community.It is natural to assume that the probabiliy of the nonspender in a community with large μ and σ may be smaller than that of a community with small μ and σ.Similar cases illustrate that in real applications, δ may depend on the other two parameters.We refer to this special kind of distribution as an associated delta-lognormal distribution.Thus, we can learn information about μ and σ from both the nonzero observations and the number of zero observations.Assume that δ is a known function of μ and σ.The unknown parameters of associated delta-lognormal distribution thus become (μ, σ).In this paper we will give the fiducial distributions and infer on the parameters.The idea is that we first obtain the fiducial distributions from the nonzero observations and then update them using the number of nonzero observations which follows a binomial distribution whose success rate is δ(μ, σ).The approach of updating is motivated by the Bayes theorem.The fiducial distributions of (μ, σ) from the nonzero observations is regarded as the "prior distribution", and is combined with the binomial distribution to get the "posterior distribution", which is referred to as the updated fiducial distribution.We further infer on μ, σ and functions of them by this updated fiducial distribution.The updated fiducial distributions of (μ, σ) are not derived from some statistics which are asymptotically normal.The asymptotically results of fiducial distribution given by [15] are no longer applicable here.Coincidentally, the updated fiducial distribution is the posterior distribution under the prior 1/σ.We show that this updated fiducial distribution enjoys the Bernstein-von Mises theorem.Then we show that the marginal fiducial distributions of the parametric functions are asymptotic confidence distributions defined in [18].Therefore, the confidence intervals of the parametric functions have asymptotically correct confidence levels.The significance levels of the hypothesis testings are also asymptotically correct.To deal with the computation, we employ the reject-sampling motivated by the approximate Bayesian computation method, see [19][20][21].Though there are some more superior sampling methods, our method is still promising benifits from its simplicity and exactness.We show in simulation study that our inference can be largely improved, due to the combination of the continous and discrete data.
The rest of the article is organized as follows.In Section 2, we introduce the associated delta-lognormal distribution and propose the updated fiducial distribution of the parameters.We further present approaches of confidence interval estimation and hypothesis testing of the parameters.Their frequentist properties are also given.We conduct simulations in Section 3 and use a real data example to illustrate our method in Section 4. We give our conclusion in the last section.

Methodology: Associated delta-lognormal distribution
In the articles mentioned earlier, three parameters in delta-lognormal distribution are always assumed to be independent.In this section, we consider the case when δ is associated with θ = (μ, σ).We assume that delta is a function of the location and scale parameters, denoted by δ(μ, σ).This means that an observation in the sample generated from the distribution may be 0 with probability δ(μ, σ) and the nonzero observations should follow a lognormal distribution with parameters μ and σ, which is denoted by LN(μ, σ).The cumulative distribution function of the associated delta-lognormal population is Gðx; m; sÞ ¼ dðm; sÞ; x ¼ 0; dðm; sÞ þ ½1 À dðm; sÞ�F LN ðx; m; sÞ; x > 0; ( where F LN (x; μ, σ) is the cumulative distribution function of LN(μ, σ).A sample from this population is denoted by X = (X 1 , X 2 , � � �, X n ).We assume that N 0 observations are zero while the rest N 1 = n − N 0 ones are nonzero.The likelihood function for the number of zero observations N 0 can be given as where n 0 is the observation of N 0 , n 1 = n − n 0 .

Updated fiducial distribution
Without loss of generality, we assume that the first N 1 observations are nonzero, while the rest are 0, that is, X ¼ ðX 1 ; X 2 ; � � � ; X N 1 ; 0; � � � ; 0Þ.Given N 1 = n 1 , the nonzero observations Then the sample mean and variance follow a normal and χ 2 (n 1 − 1) distribution respectively, that is, Let U * N(0, 1) and V * χ 2 (n 1 − 1) be two independent random variables.Then we have Given � Y ¼ � y and S 2 = s 2 , then μ and σ can be regarded as the functions of U and V The joint distribution of (U, V) is Then the joint distribution of (μ, σ) can be calculated as 2 e À ðn 1 À 1Þs 2 2s 2 ; ð2Þ where x obs ¼ ðx 1 ; x 2 ; � � � ; x n 1 ; 0; � � � ; 0Þ.This means that the fiducial distribution of (μ, σ) is If n 1 < 2, we take where � y ¼ 0 when n 1 = 0 and � y ¼ y 1 when n 1 = 1.Then the fiducial density π F (μ, σ|x obs ) in ( 2) is obtained for all n 1 � 0. The fiducial distributions for lognormal distribution is first given by [22].However, there is no common fiducial distribution for binomial variable.A generalized fiducial quantity is proposed by [15], which is a Beta distribution Beta(n 0 , n 1 + 1).Other improvements made on the parameter of the Beta distribution is further proposed by [12,14], which are 0.5[Beta(n 0 , n 1 + 1) + Beta(n 0 + 1, n 1 )] and Beta(n 0 + 0.5, n 1 + 0.5), respectively.Now we consider the problem from the Bayesian perspective, without the need of using generalized fiducial quantities.In Bayesian inference, the prior beliefs about the model parameters θ, say π(θ), are updated by observing data y obs through the likelihood function of the model.We denote the likelihood function by p(y obs |θ) and use the Bayes' theorem to get the posterior distribution The prior distribution is often specified by choosing some tractable distributions that we believe the parameters should obey.For associated delta-lognormal distribution, the prior distributions of (μ, σ) are naturally chosen to be the fiducial distributions (2), and is further updated by the likelihood function (1).We define the updated fiducial distribution of (μ, σ) as where "/" means "proportion to".

Goodness-of-fit test
Let the observation be x 1 , x 2 , � � �, x n .We take δ = G(x 0 ), where x 0 is a preset value and G is the cumulative distribution function of the continuous part.In this paper, we consider the case when G is the lognormal distribution, then In real applications, x 0 maybe known.For example, in Tobit model, see [23], ( then x 0 = y min .When x 0 is unknown, we can obtain x 0 with the following method.
Let n 0 and n 1 be the numbers of zero and nonzero observations, respectively.Without loss of generality, let x 1 ; x 2 ; � � � ; x n 1 be the nonzero ones.Then μ and σ are estimated by ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi Thus, the associated delta can be given by To test the goodness-of-fit, the classical Kolmogorov-Smirnov test is no longer suitable in the zero-inflated model.We consider using the Pearson's chi-square test.The following partition is made on the internal [0, 1), which is 0, (0, a 1 ], (a 1 , a 2 ], � � �, (a k , 1).Let where a 0 = 0. Then p 0 , p 1 , � � �, p k are estimated by Let m i be the number of samples in the interval (a i , a i+1 ), where a k+1 = 1.We can then construct the following test statistic, Given the significance level α, the model of associated delta-lognormal distribution is accepted when T � w 2 1À a ðk À 2Þ:

Inference on functions of parameters
Assume that (μ, σ) follows the updated fiducial distribution π UF (μ, σ|x obs ).Let G = g(μ, σ) which is a random variable.Then we denote the marginal fiducial distribution of g(μ, σ) by Hypothesis testing.For the one-sided hypothesis H 0 : gðm; sÞ � g 0 versus H 1 : gðm; sÞ > g 0 ; The p-value is defined as For the two-sided hypothesis H 0 : gðm; sÞ ¼ g 0 versus H 1 : gðm; sÞ 6 ¼ g 0 ; The p-value is then given by Now we start to investigate the frequenist properties of the confidence interval and the hypothesis testing.First we define the random variable Z i as ( ) are independently identically distributed as f(z, x; μ, σ) given below.The population sample space is then X ¼ f0; 1g � ½0; 1Þ and the dominating measure n ¼ C � LNð0; 1Þ, where C is the counting measure on {0, 1} and LN (0, 1) is the standard log-normal distribution, which has the density as When x = 0, we define the function above as the limit 0. The density f(z, x; μ, σ) with respect to ν is where ðz; xÞ 2 X, θ = (μ, σ) 2 O = (−1, 1) × (0, 1).We first check that f(z, x; μ, σ) is a probability density function.It can be seen that when Z = 1, X = 0, the density is while when Z = 0, X = x, the density becomes �ðlog xÞ : This indicates that f(z, x; μ, σ) is a density function with respect to ν.
Then we show that the family ( 11) is quadratic mean differentiable, which is defined below.

Definition 1 (Quadratic Mean Differentiable)
The family {P θ , θ 2 O} is quadratic mean differentiable at θ 0 if there exists a vector of real-valued functions _ l y 0 such that, as θ !θ 0 , To verify that a family is quadratic mean differentiable, a lemma below is used in this paper.Lemma 1 ([24]).For every θ in an open subset of R k , let p θ be the propbability density.Assume that the map y 7 !s y ðxÞ ¼ ffi ffi ffi ffiffi p y p is continuously differentiable for every x.If the elements of the Fisher information matrix I θ are well defined and continuous in θ, then the density p θ is quadratic mean differentiable.
The proof of this propostion is given in S1 File.

:
Notice that when n 1 � 2, the updated fiducial distribution has the form where y = log x.With simple calculation we can get This means that the updated fiducial distribution can be regarded as a posterior distribution under the prior distribution 1/σ.When n ! 1, Therefore we can apply the famous Bernstein-von Mises Theorem below to the updated fiducial distribution.Lemma 3 (Bernstein-von Mises Theorem, [24]) Let the experiment (P θ : θ 2 O) be differntiable in quadratic mean at θ 0 with nonsigular Fisher information matrix I y 0 , and suppose that for every ε > 0 there exists a sequence of test ψ n such that Furthermore, let the prior measure be absolutely continuous in a neighborhood of θ 0 with a continuous positive density at θ 0 .Then the corresponding posterior distributions satisfy kP ffi ffi n p ðy n À y 0 ÞjX obs À NðD n;y 0 ; I À 1 y 0 Þk !
At the moment we explain notations in (13).The symbol P ffi ffi n p ðy n À y 0 ÞjX obs is the posterior density of ffi ffi ffi n p ðy n À y 0 Þ while NðD n;y 0 ; I À 1 y 0 Þ is a normal distribution with mean which is the L 1 distance between densities f and g.Thus we can obtain the result below.Theorem 4 Under the assumptions of Proposition 2, Bernstein-von Mises theorem holds when the posterior distribution is replaced by the updated fiducial distribution π UF (μ, σ|x obs ).
The proof of Theorem 4 is given in S1 File.
To explore the frequenist properties of the functions of parameters under updated fiducial distribution, we give the definitions of the confidence distribution and asympototic confidence distirbution, which were proposed by [18].
1� is called a confidence distribution for a parameter θ if (i) for each given X n 2 X, H n (�) is a continuous cumulative distribution function; (ii) at the true parameter value θ = θ 0 , H n (�, θ 0 ) = H n (X n , θ 0 ), as a function of the sample X n , has the uniform distribution U(0, 1).The function H n (�) is called asymptotic confidence distribution if requirement (ii) above is replaced by (ii)' : at θ = θ 0 , H n ðy 0 Þ⇝Uð0; 1Þ as n !+ 1, and the continuity requirement on H n (�) is dropped.
The notation "⇝" means convergence in distribution.
Given n 1 � 2, under the fiducial distribution (2), it is well known that the marginal fiducial distributions of μ and σ are confidence distributions.However, under the updated fiducial distribution (6), the fiducial distributions (2) are updated by the discrete variable N 1 .Thus the marginal fiducial distributions are no longer confidence distributions.Except for μ or σ, we consider some functions of them.We have the following theorem.
Theorem 5 Let g(μ, σ) = K(aμ + bσ), where K is a strictly monotone increasing function.Then under the assumptions of Propostion 2, the marginal updated fiducial distribution of g is an asymptotic confidence distribution.
The proof of Theorem 5 is given in S1 File.Apply this theorem to different functions g(μ, σ), we can get the corollary below.

Corollary 6
The marginial updated fiducial distributions of the following functions are all asymptotic confidence distributions: ), the population mean, when The proof of Corollary 6 is given in S1 File.
An example to (v) in Corollary 6 is δ(μ, σ) = F(−μ/σ).We can see that which satisfies (14).The following proposition guanrantees the level of both the confidence interval and the hypothesis testing.

Proposition 7 If the marginal updated fiducial distribution of g(μ, σ) is an asymptotic confidence distribution. Then the level of the confidence interval is asymptotically 1 − α. The significance level of hypothesis testing is asymptotically α.
The proof of Proposition 7 is given in S1 File.From Propostion 7, if g(μ, σ) is taken as in Theorem 5 or Corollary 6, the confidence intervals in (8) and the p-values in ( 9) and ( 10) are asymptotically correct when n ! 1.When the sample size n is moderate, we give simulations in next section.

Sampling from the updated fiducial distribution
To give the confidence intervals of the parameters, we need to compute the γ-quantiles of the updated fiducial distributions.Similarly, to give the p-values of the hypothesis testing, we need to compute the cumumlative distribution functions of the marginal updated fiducial distirbutions at g 0 .However, it is difficult to give the closed forms of them.Fortunately, we can adopt a simple method to produce accurate sample from the updated fiducial distribution, which is known as the reject sampling method.
We can draw parameters from the "prior distribution" and accept the ones that generate the same number of zero as the observed data.This is similar to the reject-ABC method proposed first by [25,26].However, it shall be noticed that there is no approximation error in our sampling method for associated delta-lognormal distribution, since we don't use summary statistics and accept only the parameters which generate the same number of zero.Thus the parameters we accepted are equavilent to sampling from the real posterior distribution (6).
Without loss of generality, assume first that the observation of sample size n is x 1 ; x 2 ; � � � ; x n 1 ; 0; � � � ; 0, where x i > 0 for i = 1, � � �, n 1 and the rest n 0 = n − n 1 ones are zero.A log-transformation is then made to the nonzero observations y 1 ; y 2 ; � � � ; y n 1 .Then the fiducial distributions of μ and σ is given by ( 3), which are ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi where U is the standard normal random variable while V is a χ 2 (n 1 − 1) random variable.
1. Log-transformation is made on the nonzero observation, which is denote by y 1 ; y 2 ; � � � ; y n 1 .
The sample mean and sample variance are calculated and denoted by � y and s 2 .

2.
If n 1 � 2, sample U from the standard normal distribution and V from the w 2 n 1 À 1 distribution, respectively.To sample from the fiducial distribution of the parameters, we simply calculate μ and σ 2 using (15).If n 1 < 2, we draw samples from (4).
4. The process is repeated until we accept a certain number of parameters.
With the sample from the updated fiducial distribution, we then consider the inference on the scalar function g(μ, σ).We first assume that a certain number, say N, parameters are accepted using reject sampling method.We denote these parameters by (μ 1 , σ 1 ), (μ 2 , σ 2 ), � � �, Confidence interval.The confidence interval (8) of g(μ, σ) can be computed as follow.We sort g 0 i s in ascending order g ð1Þ � g ð2Þ � � � � � g ðNÞ : Then we take where [a] is the largest integer not larger than number a. Hypothesis testing.The first hypothesis is testing whether (μ, σ) is in a nondegenerate region.This means that the null hypothesis is (μ, σ) 2 O 0 where O 0 � < × < + .To test this hypothesis, we simply calculate the ratio of (μ i , σ i ) contatining in O 0 as follow and denote this value as the p-value where # A means the number of set A.
We also consider testing the null hypothesis H 0 : θ = θ 0 versus H 1 : θ 6 ¼ θ 0 .The p-value under the null hypothesis is then Thus we reject the null hypothesis when the p-value is not larger than a given level α.

Simulation study
In this section we illustrate the performance of our confidence intervals and hypothesis testing when the sample size is moderate.We take δ(μ, σ) = F((a − μ)/σ).Without loss of generality, we take a = 0. Otherwise, for nonzero observation , which means that we take μ − a as the new location parameter.So we consider δ(μ, σ) = F(−μ/ σ).We can sample from the updated fiducial distribution using the method we proposed.Three simulation studies are conducted in this section.The first simulation study shows the interval estimates of the parameters in associated delta-lognormal distribution, we compare this with that of the fiducial distributions to illustrate the improvements.In the second simulation study, we compare the estimates of population mean of the associated delta-lognormal and the traditional one when δ = 0.5.In the last simulation study, we focus on the estimation and hypothesis testing for δ in associated delta-lognormal distribution, we also compare the result with that of traditional one.

Simulation study I
In this simulation study we consider the estimate of μ and σ.The sample sizes considered are n = 20, 30, 50 and 100.We set the value of σ to 0.5, 1 and 2 while the value of μ is changed to make the corresponding δ = F(−μ/σ) approximately equal to 0.6, 0.5, 0.4, 0.3 and 0.15.Particularly, when σ = 1, the corresponding values of μ are −0.25,0, 0.25, 0.5 and 1.For each parameter setting, we generate 1000 repetitions and for each one we sample n = 4000 pairs of (μ, σ) and calculate the 95 percent confidence intervals of μ and σ using the Eq (16).We compare the results with the estimates obtained from the fiducial distribution and the results are shown in the Figs 1 and 2, where the details are given in the Tables in S1 File.The figures are confidence intervals of μ and σ when σ = 1 and n = 20, 30, 50 and 100.The horizontal coordinates are the values of μ = −0.25,0, 0.25, 0.5 and 1, while the vertical coordinates are the confidence intervals of μ or σ for different μ.The four plots from left to right and from top to bottom denote the cases when n = 20, 30, 50 and 100, respectively.The plot of σ = 1 and 2 are quite similar with that of 0.5, thus we don't put the figures in our context.For this specific δ, we can see that the estimate of μ is largely improved.The lower limits of μ become larger compare to the fiducial distribution while the upper limits are getting smaller.This leads to a significantly smaller confidence interval while retain the coverage probability.However, the impact on σ is not apparently as μ.The average length of the confidence intervals for σ generally get smaller than that of fiducial distribution, with the decreasing of δ and sample size n.The lower limits seems to be always bigger than that of the fiducial distribution while the upper limits gradually become smaller as δ and the sample size increase.We also notice that the distribution of σ is asymmetric, so we suggest to use the 2% and 97% quantile of the sampled σ to construct the 95% confidence interval of σ.

Simulation study II
In this simulation we consider the inference on the log population mean of the associated delta-lognormal distribution which has the form The population mean of the delta lognormal distribution plays a crucial role in statistical analysis and inference.It is a measure of central tendency, providing a summary of the central location of the distribution.For example, in the real data of our paper, we estimate the diagnostic test charges of the patients.The true value of the parameters and the sample sizes are set as we did in the last simulation.We first consider the point estimate of the log population mean, the "posterior mean" and the "posterior median" are considered, the former is approximated by while the latter is approximated by the 0.5 quantile of the N accepted values.We compute the mean bias and the mean squared error of these two estimates and compare with that of [14].
To obtain the estimate of Krishnamoorthy, we compute the mean of "Qtheta" in his paper.The result of the case when σ = 1 is shown in Table 1, the ones for σ = 0.5 and 2 can be found in the S1 File."MB", "MDB" and "GQB" stand for the mean bias of the posterior mean, posterior median and the estimate using the genralized quantity in [14]."MSE" stands for the mean squared error and the subscripts indicate the three estimate.We also use Fig 3 for a better view of the two estimates.It should be noticed that some extreme cases may occur when δ is large and the sample size is small, as is shown in the first plot of Fig 3 where n = 20.In these extreme cases, there are only three or less nonzero observations, making the estimates far from the true value, thus the mean bias and mean squared error become meaningles.So we use the blanks to indicate such problem.However, we can see that the posterior median seems to be a better point estimate of the population mean.The mean bias and the mean squared error are generally smaller, especially when σ is large.

Simulation III
In this simulation we consider the case when δ = 0.5, which happens when μ = 0. We fix μ to 0 while σ = 0.5, 1 and 2. The sample sizes range from 20 to 100.We show in Table 2 the asymptotic 95% confidence intervals of δ.The estimate of δ is compared with that of the generalized fiducial distribution proposed by Hannig, which is a Beta distribution.It can be seen that the   average length is smaller, which means that the estimate becomes more accurate.To illustrate this idea, we also test the hypothesis of δ = 0.1 to 0.9 for the case δ = F(−μ/σ) and calculate the p-value under the null hypothesis.In fact, we can consider any function of μ and σ after drawing pairs of parameters from the posterior distribution.The null hypothesis is set to σ = 0.1, 0.3, 0.5, 0.7 and 0.9.We calculate the p-value for both associated delta-lognormal and compare the result with the traditional one, which use the Beta distribution Beta(n 0 + .5, n 1 + .5)as the generalized fiducial distribution for δ.For each given set we generate 10000 samples and accept N = 10000 pairs of parameters.We calculate p i = F(−μ i /σ i ) for i = 1 to N and calculate the pvalue for δ = p 0 , which is The result is shown in Table 3.A and D in the column named method represent the associated delta-lognormal distribution and the traditional one, respectively.We can see that the pvalue of the same null hypothesis for associated delta-lognormal is more centralized than the traditional delta-lognormal distribution.This means that we are more likely to reject the null hypothesis of the associated delta-lognormal than the traditional ones when the null hypothesis is false.

A real data example
In this section, we use the data set of diagnostic test charges in [27]'s study, see Table 4.This data set is analysed by [7], who showed that the postive part fit a lognormal distribution.The data set is further studied by [9,14].This data set contains 40 patients, but 10 of them had no diagnostic tests during the study period.We assume that the data set comes from an associated delta-lognormal population, where dðm; sÞ ¼ F log x 0 À m s � � : It can be calculated that m ¼ 6:854; ŝ ¼ 1:367; x 0 ¼ 376:657; log x 0 ¼ 5:931.We assume that the data are drawn from the associated delta-lognormal distribution below, To test the goodness-of-fit, we choose k = 4 and create the partition, where a 1 , a 2 , a 3 , a 4 are 250, 500, 900 and 3000, respectively.Given the level α = 0.05, the test statistic ( 7) is 3.916, which is smaller than w 2 0:95 ð2Þ ¼ 5:991.Thus the assumption of the model is accepted.We give the confidence interval of the population mean using the method we proposed in this paper.We accept N = 10000 pairs of (μ, σ) and calculate the 2.5% and 97.5% quantiles.As we have mentioned in last section, the 2% and 97% quantiles are also considered since the distribution of σ is asymmetric.The result is compared with the Fiducial method proposed by [14] and the "MOVER" proposed by [16], see Table 5.It can be seen that the confidence interval is largely improved.