Next Article in Journal
TFD-IIS-CRMCB: Telecom Fraud Detection for Incomplete Information Systems Based on Correlated Relation and Maximal Consistent Block
Next Article in Special Issue
Detecting Structural Change Point in ARMA Models via Neural Network Regression and LSCUSUM Methods
Previous Article in Journal
Distance in Information and Statistical Physics III
Previous Article in Special Issue
Asymptotics of Subsampling for Generalized Linear Regression Models under Unbounded Design
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Comparing Several Gamma Means: An Improved Log-Likelihood Ratio Test

Department of Mathematics and Statistics, York University, Toronto, ON M3J 1P3, Canada
Submission received: 20 November 2022 / Revised: 28 December 2022 / Accepted: 3 January 2023 / Published: 5 January 2023
(This article belongs to the Special Issue Recent Advances in Statistical Theory and Applications)

Abstract

:
The two-parameter gamma distribution is one of the most commonly used distributions in analyzing environmental, meteorological, medical, and survival data. It has a two-dimensional minimal sufficient statistic, and the two parameters can be taken to be the mean and shape parameters. This makes it closely comparable to the normal model, but it differs substantially in that the exact distribution for the minimal sufficient statistic is not available. A Bartlett-type correction of the log-likelihood ratio statistic is proposed for the one-sample gamma mean problem and extended to testing for homogeneity of k 2 independent gamma means. The exact correction factor, in general, does not exist in closed form. In this paper, a simulation algorithm is proposed to obtain the correction factor numerically. Real-life examples and simulation studies are used to illustrate the application and the accuracy of the proposed method.

1. Introduction

Consider a sample of ( x 1 , , x n ) from the two-parameter gamma model with mean μ and shape λ . The joint density is
f ( x 1 , , x n ; μ , λ ) = Γ n ( λ ) λ μ n λ exp λ t λ μ s i = 1 n 1 x i
where ( s , t ) = i = 1 n x i , i = 1 n log x i is a minimal sufficient statistic. The two-parameter gamma distribution is often used to model the non-negative data with right-skewed distribution. Moreover, depending on the values of the parameters, it can have a decreasing failure rate, a constant failure rate, or an increasing failure rate. This makes it a valuable model for analyzing data arising from engineering, environmental, meteorological, and medical studies.
Similar to the normal distribution, the two-parameter gamma distribution has a two-dimensional minimal sufficient statistic ( s , t ) . Another version of the minimal sufficient statistic is ( r , s ) , where r is the log offset of the arithmetic mean from the geometric mean
r = log i = 1 n x i n log i = 1 n x i 1 / n = log s n t n .
Notice that the density of log x i has location model form for a fixed λ . It follows from [1] that the conditional density for s given r and the marginal density for r take the form
f ( r ; λ ) = Γ ( n λ ) Γ n ( λ ) n n λ + 1 / 2 exp { n λ r } h n ( r ) f ( s | r ; μ , λ ) = Γ 1 ( n λ ) λ μ n λ exp n λ log s λ s μ 1 s ,
where h n ( r ) appears in the transformed measure
i = 1 n 1 x i d x i = 1 s d s × n h n ( r ) d r
which requires ( n 2 ) -dimensional integration. Hence, the joint density for ( s , r ) is
f ( s , r ; μ , λ ) = f ( s | r ; μ , λ ) f ( r ) .
By change of variable, the joint density for (s, t) takes the form
f ( s , t ; μ , λ ) = Γ n ( λ ) λ μ n λ exp λ t λ μ s 1 s 1 n h n log s n t n .
The same result can be obtained by using the properties of the exponential transformation model by [2] or the conditional argument in [3]. Note that h n ( · ) requires ( n 2 ) -dimensional integration, and it is available exactly only for small value of n (see [3,4]). Unlike the normal model, where inference for the normal mean can be obtained explicitly, inference for the gamma mean is a complicated and challenging problem. Many asymptotics inferential methods for the gamma mean exist in statistical literature. Moreover, most of the existing asymptotics methods are likelihood-based methods. Some are very simple but not very accurate, and others are very accurate but mathematically complicated and also computationally intensive. Furthermore, only limited methods can be applied to the problem of comparing the means of k > 2 independent gamma distributions.
Ref. [5] considered the log-likelihood function obtained from (1), which takes the form
( μ , λ ) = n log Γ ( λ ) + n λ log λ n λ log μ + λ t λ μ s .
The maximum likelihood estimate (MLE) is ( μ ^ , λ ^ ) , where μ ^ = s n , and λ ^ satisfies
ψ ( λ ^ ) + log λ ^ log s n + t n = 0
with ψ ( · ) being the digamma function. In addition, the observed information matrix evaluated at MLE is
j ^ = n λ ^ μ ^ 2 + 2 λ ^ μ ^ 3 s n μ ^ s μ ^ 2 n μ ^ s μ ^ 2 n λ ^ + n ψ ( λ ^ ) .
It is well-known that the variance-covariance matrix of the maximum likelihood estimators can be approximated by j ^ 1 . Hence, the standardized maximumm likelihood estimator is
Q ( μ ) = μ ^ μ v a r ^ ( μ ^ ) ,
where v a r ^ ( μ ^ ) can be approximated by the ( 1 , 1 ) entry of j ^ 1 . With the regularity conditions stated in [6] and also in [7], for large n,
Q ( μ ) d N ( 0 , 1 )
with first order accuracy, O ( n 1 / 2 ) . Thus, inference for μ can be obtained based on the limiting distribution. This method is generally known as the Wald method or the asymptotic MLE method.
Another commonly used method to obtain inference for μ is the log-likelihood ratio method. Let λ ˜ μ be the constrained MLE, which maximizes ( μ , λ ) for a fixed μ . In this case, λ ˜ μ must satisfy
ψ ( λ ˜ μ ) + log λ ˜ μ log μ + 1 s n μ + t n = 0 .
Then the log-likelihood ratio statistic is
W ( μ ) = 2 [ ( μ ^ , λ ^ ) ( μ , λ ˜ μ ) ] .
Again, with the regularity conditions stated in [6] and also in [7], for large n,
W ( μ ) d χ 1 2 .
with first order accuracy, O ( n 1 / 2 ) . Hence, inference for μ can be approximated based on the limiting χ 1 2 distribution. This method is also known as the Wilks method.
To improve the accuracy of the log-likelihood ratio method, ref. [8] applied the Bartlett correction to the log-likelihood ratio statistic (see [9]). The resulting Bartlett corrected log-likelihood ratio statistic takes the form
W ( μ ) = W ( μ ) 1 + B ( μ ) / n
where
B ( μ ) = 1 6 λ ˜ μ + 1 2 D 3 ( λ ˜ μ ) [ D 2 ( λ ˜ μ ) ] 2 + 1 4 D 2 ( λ ˜ μ ) D 2 ( λ ) = λ 2 [ λ 1 ψ ( λ ) ] D 3 ( λ ) = λ 3 [ λ 2 ψ ( λ ) ] ,
and B ( · ) is known as the Bartlett correction factor. Ref. [9] showed that the Bartlett corrected log-likelihood ratio statistic converges to χ 1 2 distribution with fourth order accuracy. Hence, inference for μ can be approximated based on the limiting χ 1 2 distribution.
Refs. [3,4] showed that the exact form of h n ( · ) in (2) is only available when n is small. Jensen used the fact that the model is an exponential-transformation model and applied the saddlepoint method to approximate h n ( · ) and derived an inference procedure for μ with third order accuracy. However, due to the complexity of the method, ref. [3] only provided tables for 1, 2.5, 97.5, and 99 percentiles of μ for sample sizes 10, 20, 40, and , which were obtained by extensive iterative calculations. On the other hand, ref. [4] proposed another third order method to obtain inference for μ . This method is asymptotically equivalent to Jensen’s method with the exception that it involves direct implementation of the method derived in [10].
Note that Gross and Clark’s method is very simple but not very accurate. The log-likelihood ratio method is slightly more complicated because of the calculation of the constrained MLE, and it is still not very accurate. The Bartlett corrected log-likelihood ratio method by [8] gives very accurate results. It is also relatively straightforward because [8] derived all the necessary equations. The method presented in [4] is also very accurate, but it is computational intensive. Gross and Clark’s method, Jensen’s method and Fraser, Reid and Wong’s method are not applicable to the problem of testing homogeneity of k > 2 independent gamma means. The log-likelihood ratio method can be extended to this problem but it has only first order accuracy. Ref. [8] also derived the explicit Bartlett correction factor for testing equality of two independent gamma means, but, due to the complexity of the method, did not derive the explicit Bartlett correction factor for testing homogeneity of k > 2 independent gamma means.
In this paper, a Bartlett-type correction of the log-likelihood ratio statistic is proposed in Section 2. The proposed Bartlett-type correction factor is numerically obtained by simulations. The proposed method is then applied testing homogeneity of k 2 independent gamma means in Section 3 and Section 4, respectively. Some concluding remarks are given in Section 5. Real-life examples and simulation studies results are presented to compare the accuracy of the proposed method with the existing methods.

2. Main Results

Let ( θ ) be the log-likelihood function with a p-dimensional parameter θ . With the regularity conditions stated in [6], the log-likelihood ratio statistic,
W ( θ ) = 2 [ ( θ ^ ) ( θ ) ] ,
is asymptotically distributed as χ p 2 with first order accuracy, where θ ^ is the overall MLE, which maximizes ( θ ) . Ref. [9] showed that the mean of W ( θ ) can be expressed as
E [ W ( θ ) ] = p 1 + B ( θ ) n + O ( n 2 ) ,
where n is the size of the observed sample and B ( · ) is the Bartlett correction factor. Hence, the Bartlett corrected log-likelihood ratio statistic is
W ( θ ) = W ( θ ) 1 + B ( θ ) / n d χ p 2
and has mean p with fourth order accuracy.
The above method can be generalized to the case when ψ = ψ ( θ ) is the parameter of interest and dimension of ψ is m < p . With the regularity conditions stated in [6], the log-likelihood ratio statistic
W ( ψ ) = 2 [ ( θ ^ ) ( θ ˜ ψ ) ]
is asymptotically distributed as χ m 2 with first order accuracy. Note that θ ^ is the overall MLE, which maximizes ( θ ) , and θ ˜ ψ is the constrained MLE, which maximizes ( θ ) for a given value of ψ . The Bartlett corrected log-likelihood ratio statistic is then
W ( ψ ) = W ( ψ ) 1 + B ( ψ ) / n d χ m 2
and W ( ψ ) has mean m with fourth order accuracy.
Theoretically, the Bartlett correction method gives extremely accurate results, even for small sample sizes. However, obtaining the explicit closed form expression of the Bartlett correction factor is a very difficult problem. There exists only limited problems in statistical literature that the explicit closed form, or even the asymptotic form of the Bartlett correction is available. For example, Ref. [8] obtained the Bartlett correction factor for the one-sample gamma mean problem as well as for the equality of the two independent gamma means problem only. However, they did not discuss the case for testing homogeneity of k > 2 independent gamma means. The aim of this paper is to propose a systematic way of approximating a Bartlett-type correction factor.
Since the log-likelihood ratio statistic W ( ψ ) has the limiting distribution χ m 2 , similar to the Bartlett correction method, we want to find a scale transformation of W ( ψ ) such that the transformed statistic has the exact mean m. An obvious transformation is
W ( ψ ) = W ( ψ ) E [ W ( ψ ) ] / m .
However, calculating E [ W ( ψ ) ] is extremely complicated, if not impossible. If we have an observed sample of W ( ψ ) , we can then apply the method of moments to estimate E [ W ( ψ ) ] . The primary task is to obtain such a sample. We propose to employ simulations to create such a sample. The main idea to to generate samples from the original model but with the parameters being the constrained MLE obtained from the original observed sample. We summarized the idea into the following algorithm.
Assume: ( x 1 , , x n ) is a sample from a model with density f ( · ; θ ) , and ψ is the parameter of interest.
Have:The log-likelihood function ( θ ) is given in (3).
From the log-likelihood function, we can obtain θ ^ , θ ˜ ψ , and W ( ψ ) .
Step 1:Generate a sample of size n , ( x 1 s , , x n s ) , from the density f ( · ; θ ˜ ψ ) .
Step 2:From the simulated sample, obtain the log-likelihood ratio statistic and denote it as W s ( ψ ) .
Step 3:Repeat Steps 1 to 2 N times, where N is large.
As a result, we have ( W 1 s ( ψ ) , , W N s ( ψ ) ) .
Step 4:Obtain
W ¯ s ( ψ ) = i = 1 N W i s ( ψ ) N .
Step 5:By the method of moments, W ¯ s ( ψ ) is an consistent estimate of E [ W ( ψ ) ] . Hence,
W ( ψ ) = W ( ψ ) W ¯ s ( ψ ) / m
has mean m, and its limiting distribution is χ m 2 .

3. One-Sample Gamma Mean Problem

Consider the one-sample gamma mean problem with the log-likelihood given in (3). Ref. [5] proposed to use Wald statistic given in (4) to obtain inference for μ , whereas [8] recommended to use the Bartlett corrected log-likelihood ratio statistic given in (6) to obtain inference for μ . In this paper, a Bartlett-typed corrected log-likelihood ratio statistic based on the algorithm given in Section 2 is proposed as an alternative approach to obtain inference for μ . To compare the results obtained by these methods, we consider the data set given in [5], which is the survival time of 20 mice exposed to 240 rad of gamma radiation:
152 152 115 109 137 88 94 77 160 165
125 40 128 123 136 101 62 153 83 69
Table 1 recorded the 95% confidence intervals for μ and the p-values for testing H 0 : μ = 133 vs. H a : μ 133 obtained by the methods discussed in this paper. From Table 1, we observed that results obtained by Jensen and Kristensen’s method and by the proposed method are almost identical. However, they are very different from the results obtained by Gross and Clark’s method and the standard log-likelihood ratio statistic method.
Simulation studies are performed to compare the accuracy of the methods discussed in this paper. In particular, 5000 simulated samples are obtained for each combination of μ , λ , and n. Moreover, we use N = 100 for each of the simulated sample to approximate the mean of the log-likelihood ratio statistic. The proportion of samples that was rejected at 5% significance level are recorded in Table 2. Theoretically, the true percentage of samples that will be rejected is 5% with a standard deviation of 0.31%. Extensive simulation studies were performed. All results are very similar and, therefore, only a subset of them are presented in Table 2. Results from other combinations of the parameters are available from the authors upon request.
We observed that results by Gross and Clark’s method are significantly different from the nominal 5% level, but the accuracy is improving slowly as the sample size n increases. Results by the log-likelihood ratio method are slightly better, and the accuracy improves much faster as the sample size increases. Results by both Jensen and Kristensen’s method and the proposed method are very accurate and always within 3 standard deviation of the nominal 5% level, even when the sample size is as small as 5.

4. Testing Homogeneity of k Independent Gamma Means

In this section, the proposed method is extended to testing homogeneity of k 2 independent gamma means problem. Let ( x i 1 , x i n i ) be a sample from the two-parameter gamma model with mean μ i and shape λ i , where i = 1 , , k and k 2 . Moreover, assume the k two-parameter gamma models are independent. Let i ( μ i , λ i ) be the log-likelihood function from the ith model. Then the joint log-likelihood function is
( μ 1 , , μ k , λ 1 , , λ k ) = i = 1 k i ( μ i , λ i ) = i = 1 k n i λ i log λ i n i λ i log μ i n i log Γ ( λ i ) + λ i t i λ i μ i s i
where
s i = j = 1 n i x i j and t i = j = 1 n i log x i j .
Since the k models are independent, the overall MLE is ( μ ^ 1 , , μ ^ k , λ ^ 1 , , λ ^ k ) where μ ^ i = s i n i and λ ^ i must satisfy
ψ ( λ ^ i ) + log ( λ ^ i ) log s i n i + t i n i = 0 .
Hence, the log-likelihood function is evaluated at MLE ( μ ^ 1 , , μ ^ k , λ ^ 1 , , λ ^ k ) .
To test homogeneity of k gamma means, the null and alternative hypotheses are
H 0 : μ 1 = = μ k = μ vs . H a : not all equal .
The constrained log-likelihood function is
( μ , , μ , λ 1 , , λ k ) = i = 1 k n i λ i log λ i n i λ i log μ n i log Γ ( λ i ) + λ i t i λ i μ s i .
The constrained MLE can must satisfy
μ ˜ = i = 1 k s i λ ˜ i i = 1 k n i λ ˜ i and ψ ( λ ˜ i ) + log ( λ ˜ i ) + 1 log μ ˜ log s i n i μ ˜ + t i n i = 0 .
Thus, the log-likelihood function evaluated at the constrained MLE is ( μ ˜ , , μ ˜ , λ ˜ 1 , , λ ˜ k ) . Finally, the log-likelihood ratio statistic is
W ( μ ˜ ) = 2 ( μ ^ 1 , , μ ^ k , λ ^ 1 , , λ ^ k ) ( μ ˜ , , μ ˜ , λ ˜ 1 , , λ ˜ k ) ,
and it is asymptotically distributed as χ k 1 2 with first order accuracy. Note that for this problem, the Wald method is not applicable. Moreover, ref. [8] did not derive the exact Bartlett correction factor for k > 2 . However, the proposed method is still applicable.
To illustrate the application of the log-likelihood ratio method and the proposed method for this problem, we consider the intervals in service-hours between failures of the air-conditioning equipment in 10 Boeing 720 jet aircrafts reported in “Example T” from [11]. It is assumed that the reported times for each aircraft is distributed as a two-parameter gamma distribution. The question of interest is whether the ten aircrafts have the same mean intervals in service hours between failure. In other words, we are testing
H 0 : μ 1 = = μ 10 = μ vs . H a : not all equal .
The log-likelihood ratio method gives a p-value of 0.0871, whereas the proposed method gives a p-value of 0.1295 (using N = 100,000). At 10% level of significance, the two methods give contradictory results with the log-likelihood ratio method rejecting H 0 , and the proposed method failing to reject H 0 .
As in the one sample case, to compare the accuracy of the two methods, extensive simulation studies were performed. For each combination of k , n 1 , , n k , ( μ 1 , λ 1 ) , , ( μ k , λ k ) , 5000 simulated samples are obtained. Moreover, for each simulated sample, N = 100 is used to estimate the mean of the log-likelihood ratio statistic. The proportion of samples that reject the null hypothesis of homogeneity of the k means at 5% significance levels are recorded. Since all results are very similar, only results from the 8 cases listed in Table 3 are reported.
Results in Table 4 and Table 5 demonstrate that the log-likelihood ratio method gives unsatisfactory results, especially when the sample sizes are small. However, the accuracy of the results improve as the sample sizes increase. In comparison, the results from the proposed method are very accurate, and they are always within 3 standard deviation of the nominal 5% level, regardless of the sample sizes.

5. Conclusions

In this paper, a Bartlett-type corrected log-likelihood ratio method for comparing the means of several independent gamma distributions is proposed. The method can easily be applied in statistics software, such as R. Simulation results demonstrate the log-likelihood ratio method does not give satisfactory results, especially when the sample sizes are small. However, the proposed method is extremely accurate even when the sample sizes are small. One advantage of the proposed method is that it is not restricted to the gamma means problem, as it is applicable to any parametric models.

Funding

This research was funded by Natural Sciences and Engineering Researh Council of Canada Grant RGPIN 2017-05179.

Institutional Review Board Statement

Not applicable.

Data Availability Statement

Data is contained within the article.

Acknowledgments

The authors would like to thank the editor and the referees whose comments have improved the quality of this paper.

Conflicts of Interest

The author declares no conflict of interest.

References

  1. Fraser, D.A.S. The Structure of Inference; Wiley: New York, NY, USA, 1968. [Google Scholar]
  2. Barndorff-Nielsen, O.E.; Blasild, P.; Jensen, J.L.; Jorgensen, B. Exponential tranformation models. Proc. R. Soc. Lond. Ser. A 1982, 39, 41–65. [Google Scholar]
  3. Jensen, J.L. Inference for the mean of a gamma distribution with unknown shape parameter. Scand. J. Stat. 1986, 13, 135–151. [Google Scholar]
  4. Fraser, D.A.S.; Reid, N.; Wong, A. Simple and accurate inference for the mean of the gamma model. Can. J. Stat. 1997, 25, 91–99. [Google Scholar] [CrossRef] [Green Version]
  5. Gross, A.J.; Clark, V.A. Survival Distributions: Reliability Applications in the Biomedical Sciences; Wiley: New York, NY, USA, 1975. [Google Scholar]
  6. Cox, D.R.; Hinkley, D.V. Theoretical Statistics; CRC Press: London, UK, 1979. [Google Scholar]
  7. Barndorff-Nielsen, O.E.; Cox, D.R. Inference and Asymptotics; Chapman and Hall: London, UK, 1994. [Google Scholar]
  8. Jensen, J.L.; Kristensen, I.B. Saddlepoint approximation to exact test and improved likelihood ratio tests for the gamma distribution. Commun. Stat. Theory Methods 1991, 20, 1515–1532. [Google Scholar] [CrossRef]
  9. Bartlett, M.S. Properties of sufficiency and statistical tests. Proc. R. Soc. Lond. Ser. A 1937, 160, 268–282. [Google Scholar]
  10. Fraser, D.A.S.; Reid, N. Ancillaries and third order significance. Util. Math. 1995, 7, 33–53. [Google Scholar]
  11. Cox, D.R.; Snell, E.J. Applied Statistics: Principles and Examples; Chapman and Hall: London, UK, 1981. [Google Scholar]
Table 1. 95% confidence interval for μ and p-value for testing H 0 : μ = 133 vs. H a : μ 133 .
Table 1. 95% confidence interval for μ and p-value for testing H 0 : μ = 133 vs. H a : μ 133 .
MethodConfidence Intervalp-Value
Gross and Clark(96.7, 130.2)0.0223
Log-likelihood ratio(97.6, 133.0)0.0499
Jensen and Kristensen(97.0, 133.9)0.0586
Proposed(97.0, 133.9)0.0596
Table 2. Observed proportion of samples (in percentage) that was rejected μ = μ 0 at 5% significance levels of 5000 simulated samples.
Table 2. Observed proportion of samples (in percentage) that was rejected μ = μ 0 at 5% significance levels of 5000 simulated samples.
μ 0 λ 0 nGross and ClarkLog-Likelihood RatioJensen and KristensenProposed
1.50.5521.829.705.705.20
1014.566.825.125.04
1511.845.944.684.98
2010.905.985.105.30
2.02.0517.429.365.765.02
1010.547.425.245.64
159.086.645.625.78
207.385.544.625.04
4.01.0518.769.265.485.02
1012.666.965.045.10
1510.806.205.165.32
208.765.604.664.76
Table 3. Various combinations of k , ( μ 1 , λ 1 ) , , ( μ k , λ k ) .
Table 3. Various combinations of k , ( μ 1 , λ 1 ) , , ( μ k , λ k ) .
Casek ( μ 1 , λ 1 ) , , ( μ k , λ k )
(1)3 ( 2 , 2 ) , ( 2 , 2 ) , ( 2 , 2 )
(2)3 ( 1 , 2 ) , ( 1 , 4 ) , ( 1 , 5 )
(3)3 ( 3 , 1 3 ) , ( 3 , 2 3 ) , ( 3 , 4 3 )
(4)3 ( 1 , 1 ) , ( 1 , 1 2 ) , ( 1 , 1 3 )
(5)5 ( 2 , 2 ) , ( 2 , 2 ) , ( 2 , 2 ) , ( 2 , 2 ) , ( 2 , 2 )
(6)5 ( 1 , 2 ) , ( 1 , 3 ) , ( 1 , 5 ) , ( 1 , 7 ) , ( 1 , 8 )
(7)5 ( 1 , 1 3 ) , ( 1 , 1 2 ) , ( 1 , 1 4 ) , ( 1 , 3 4 ) , ( 1 , 2 3 )
(8)5 ( 2 , 1 ) , ( 2 , 3 ) , ( 2 , 1 4 ) , ( 2 , 2 3 ) , ( 2 , 2 )
Table 4. Observed proportion of samples (in percentage) that rejected the null hypothesis of homogeneity of means at 5% significance levels.
Table 4. Observed proportion of samples (in percentage) that rejected the null hypothesis of homogeneity of means at 5% significance levels.
Case (1)
n 1 n 2 n 3 Log-likelihood ratioProposed
5550.10820.0490
58150.09200.0544
812160.07320.0496
1015200.06920.0542
3030300.05920.0552
Case (2)
n 1 n 2 n 3 Log-likelihood ratioProposed
5550.10900.0532
58150.09980.0546
812160.07380.0546
1015200.07060.0552
3030300.06180.0554
Case (3)
n 1 n 2 n 3 Log-likelihood ratioProposed
5550.11180.0480
58150.09860.0564
812160.07580.0486
1015200.06560.0460
3030300.06180.0550
Case (4)
n 1 n 2 n 3 Log-likelihood ratioProposed
5550.10880.0522
58150.09620.0536
812160.07140.0474
1015200.06800.0522
3030300.05620.0514
Table 5. Observed proportion of samples (in percentage) that was rejected homogeneity of means at 5% significance levels of 5000 simulated samples.
Table 5. Observed proportion of samples (in percentage) that was rejected homogeneity of means at 5% significance levels of 5000 simulated samples.
Case (5)
n 1 n 2 n 3 n 4 n 5 Log-likelihood ratioProposed
555550.14500.0530
57810150.09440.0502
8121622300.07060.0478
10131621250.07440.0508
30303030300.05620.0488
Case (6)
n 1 n 2 n 3 n 4 n 5 Log-likelihood ratioProposed
555550.13200.0464
57810150.09800.0504
8121622300.08260.0560
10131621250.07680.0550
30303030300.05840.0512
Case (7)
n 1 n 2 n 3 n 4 n 5 Log-likelihood ratioProposed
555550.13840.0462
57810150.10340.0480
8121622300.07740.0514
10131621250.08040.0544
30303030300.06320.0544
Case (8)
n 1 n 2 n 3 n 4 n 5 Log-likelihood ratioProposed
555550.13620.0482
57810150.10140.0478
8121622300.07500.0534
10131621250.07540.0564
30303030300.05760.0492
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

Wong, A. Comparing Several Gamma Means: An Improved Log-Likelihood Ratio Test. Entropy 2023, 25, 111. https://0-doi-org.brum.beds.ac.uk/10.3390/e25010111

AMA Style

Wong A. Comparing Several Gamma Means: An Improved Log-Likelihood Ratio Test. Entropy. 2023; 25(1):111. https://0-doi-org.brum.beds.ac.uk/10.3390/e25010111

Chicago/Turabian Style

Wong, Augustine. 2023. "Comparing Several Gamma Means: An Improved Log-Likelihood Ratio Test" Entropy 25, no. 1: 111. https://0-doi-org.brum.beds.ac.uk/10.3390/e25010111

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