Asymptotic efficiency and small sample power of a locally most powerful linear rank test for the log-logistic distribution

Assume that two independent random samples are distributed according to a log-logistic distribution (LLD). In this study, the score functions for the locally most powerful rank test were derived for the location and scale parameters. The Wilcoxon rank-sum test was shown to be locally most powerful rank test for the LLD. The asymptotic efficiency of the Wilcoxon rank-sum test was derived and compared with that of the modified Wilcoxon rank-sum test for the LLD.


Introduction
Testing hypotheses is one of the most important problems in performing nonparametric statistics. Various nonparametric statistics have been shown to be a locally most powerful rank test (LMPRT) for a specific distribution over the course of many years. For example, [9] obtained the LMPRTs for comparing two possibly censored samples for a given alternative, by deriving scores for censored and uncensored observations. In finance literature, the Lévy distribution arises as a special case for describing security price returns by mixtures of distributions. Runde [13] derived the score functions of LMPRTs for the location and scale parameters with the Lévy distribution. Hájek et al. [5] discussed the LMPRTs for various distributions. The most powerful test for correlation is well-known, and the LMPRT uses Fisher-Yates expected normal scores. However, the bivariate normal distribution does not fit some types of data. Conover [2] considered the LMPRTs for correlation with four examples. Pandit [11] investigated the LMPRT for testing independence against a weighted contamination alternative. Few studies have discussed the LMPRT when only a fraction of treated subjects respond to treatment; however, [12] examined the LMPRT for several of these case types.
Let {X ij ; i ¼ 1; 2; j ¼ 1; . . .; n i } be two random samples of size n 1 and n 2 independent observations, each of which have a continuous distribution described as F 1 and F 2 , respectively. Assume that X ij is distributed according to a log-logistic distribution (LLD); see e.g., [14], with the probability density function (pdf) f ðxÞ and the cumulative distribution function (cdf) FðxÞ as follows: respectively. The parameters l; a and b denote the location, scale and shape parameters, respectively. Note that the median of the LLD is a þ l. This distribution is a special case of Burr's type-XII distribution; see e.g., [1]. In addition, the LLD is a special case of the kappa distributions introduced by [6]. The LLD is often applied to economics as a simple model for wealth or income distribution [3]. Note that the LLD is known as the Fisk distribution in economics. Additionally, the LLD is used in survival analysis as a parametric model for events whose hazard rate increases initially and decreases later. On the basis of these data types, we proposed the following hypothesis: For a recent comparison study of many nonparametric tests for scale, see [8]. To test these hypotheses, we developed a linear rank statistic. Let V j ¼ 0 if the jth smallest of the N ¼ n 1 þ n 2 observations is from X 1j ; otherwise, V j ¼ 1. This provides a general two-sample linear rank statistic, as follows: aðR j Þ: Herein, R j denotes the rank of sample X 2j . In Sect. 2, we derive a LMPRT for H L and H S . Since finite sample sizes are used in practice, we investigate a small sample power of the linear rank tests for the LLD in Sect. 3. Finally, conclusions are presented in Sect. 4.

Locally most powerful linear rank test
In this section, we derive the LMPRT for H L and H S . Assume that f is absolutely continuous and that Z f 0 ðxÞ j jdx \ 1: The LMPRT for H L is given by [5], as follows: where U ð1:NÞ \ Á Á Á \ U ðN:NÞ are ordered statistics from the uniform distribution in the interval ½0; 1. In addition, the LMPRT for H S is given by [5] as LMPRT for the location parameter In this section, the score function of the LMPRT is derived for the LLD for the location parameter. The standard LLD is given as follows: From (2), the first derivative of pdf and the inverse of cdf for the standard LLD are, respectively, given by The score function w L ðU ðj:NÞ ; f Þ for H L is obtained from the following: using (1). (4) denotes the inverse function of R j . Thus, the Wilcoxon rank-sum test is the LMPRT for H L with LLDs. Note that the LMPRT is valid in the neighborhood close to the null hypothesis. In many cases, the asymptotic efficiency of the LMPRT is the highest for the adjusted distribution. However, moderate to large sample sizes are required to assume a specific distribution for deriving the LMPRT. Herein, we consider the asymptotic efficiency of both the original and modified Wilcoxon rank-sum tests for the LLD. Tamura [15] proposed the modified Wilcoxon rank-sum test to raise the asymptotic relative efficiency as follows: Note that the test statistic TðpÞ is the original Wilcoxon rank-sum test when p ¼ 1. By applying the idea of [15] to the score function (3), we can obtain the score function of the inverse rank test: The asymptotic efficiency of the score function (6) for the location parameter is given as follows: The asymptotic efficiency AELðp; f Þ and asymptotic relative efficiency AREðp; 1; f Þ to the original Wilcoxon rank-sum test for the LLD are listed in Table 1.
The results shown in Table 1 indicate that the asymptotic efficiency of the TðpÞ test is higher than that of the original Wilcoxon rank-sum test for LLDs when p [ 1.
LMPRT for the scale parameter Consider the score function of the LMPRT for H S . Using a similar procedure as that applied to the location case, we obtain the score function, given below: using (1). Thus, (7) reveals that the Wilcoxon rank-sum test is the LMPRT for H S with a LLD. Note that the score function of the modified Wilcoxon rank-sum test is given as [4]: w Ã S ðU ði:NÞ ; f Þ ¼ ðp þ 1ÞU p ði:NÞ À 1: From this, the asymptotic efficiency of the modified Wilcoxon rank-sum test for the scale parameter is given by xf ðxÞ 2 FðxÞ pÀ1 dx 2

:
The asymptotic efficiency AESðp; f Þ and the asymptotic relative efficiency AREðp; 1; f Þ to the original Wilcoxon rank-sum test for the LLD are listed in Table 2.
Therefore, the results shown in Table 2 reveal that the original Wilcoxon rank-sum test is suitable for testing LLDs.
Herein, we consider another distribution, as follows: x qþ1 f1 þ ð2 q À 1Þx Àq g ð1 þ qÞ=q ; For the case q ¼ 1, this distribution is equivalent to that of the standard LLD. Thus, f Ã and F Ã is another extension of the standard LLD. Using a procedure similar to the one presented earlier in this section, the score function w S ðU ði:NÞ ; f Þ is given by w S ðU ði:NÞ ; f Ã Þ ¼ ðq þ 1ÞU q ði:NÞ À 1: Therefore, the modified Wilcoxon rank-sum test is the LMPRT for H S with f Ã . Additionally, we compared the asymptotic efficiency of the modified Wilcoxon rank-sum test with that of the original Wilcoxon rank-sum test for f Ã in Table 3. Table 3 shows that the modified Wilcoxon rank-sum test is more efficient than is the original Wilcoxon rank-sum test for f Ã .

Powers of small sample sizes
In this section, we present our investigation of the behavior of the original Wilcoxon rank-sum test and the modified Wilcoxon rank-sum test for the LLD f ðxÞ. In a previous section, we obtained asymptotic results. However, because finite sample sizes are used in practice, we investigated a small sample power of the original and modified Wilcoxon rank-sum tests. Generally, the location, scale, locationscale, and shape parameters of the X 1j and X 2j samples are unequal. Here, we considered a two-sample problem in which the hypothesis, H 0 : F 1 ðxÞ ¼ F 2 ðxÞ, was tested against H 1 : F 1 ðxÞ 6 ¼ F 2 ðxÞ. The following assumptions were made for LLDs, as follows:   LLDðl; a; bÞ with location parameter l, scale parameter a and shape parameter b. Table 4 and Table 5 show the analysis results of the power of the TðpÞ for ðn 1 ; n 2 Þ ¼ ð10; 10Þ and (5,10). The significance level was 5 %, and the simulations were repeated 1,000,000 times. In this study, we assumed p ¼ 0:5; 0:75; 1 and 2. The exact probability of the original Wilcoxon rank-sum test Tð1Þ is 0.044 for n 1 ¼ n 2 ¼ 10 under the null hypothesis. Thus, we adopted an exact probability of 0.044 for the TðpÞ test.
The results indicated that a smaller p in the TðpÞ test was more powerful than a larger p in the TðpÞ test for equal and unequal sample sizes, with b 1 for the shifted location, scale, and location-scale parameters. For the case of b [ 1, a larger p in the TðpÞ test was more efficient than a small p in the TðpÞ test for the changed scale parameter, but not the location and location-scale parameters. The results indicated Tð0:5Þ and Tð0:75Þ tests were suitable for a wide range of LLD. Therefore, a smaller p in the TðpÞ test is suitable for parameters associated with LLDs.

Conclusion and discussions
In this paper, we considered the LMPRT for the location and scale parameters with the LLD. We showed that the Wilcoxon rank-sum test is suitable for hypothesis testing with the LLD. In addition, the asymptotic efficiency of the modified Wilcoxon rank-sum test was higher than that of the original Wilcoxon rank-sum test for the location parameter with a LLD. For the scale parameter, the asymptotic efficiency of the original Wilcoxon rank-sum test was highest for the LLD. Additionally, the score function of another extension of the LLD was derived as the modified Wilcoxon rank-sum test. However, because finite sample sizes are used in practice, we investigated a small sample power of the original and modified Wilcoxon rank-sum tests. The results indicated that a smaller p in the TðpÞ test was suitable for parameters associated with LLDs. Future research directions should consider the twosample location-scale problem-i.e., where two null hypotheses H 0 are simultaneously tested, see; e.g., [7,10]-for the LLD.