Binomial Regression Models with a Flexible Generalized Logit Link Function

: In binomial regression, a link function is used to join the linear predictor variables and the expectation of the response variable. This paper proposes a ﬂexible link function from a new class of generalized logistic distribution, namely a ﬂexible generalized logit (glogit) link. This approach considers both symmetric and asymmetric models, including the cases of lighter and heavier tails, as compared to standard logistic. The glogit is created from the inverse cumulative distribution function of the exponentiated-exponential logistic (EEL) distribution. Using a Bayesian framework, we conduct a simulation study to investigate the model performance compared to the most commonly used link functions, e.g., logit, probit, and complementary log–log. Furthermore, we compared the proposed model with several other asymmetric models using two previously published datasets. The results show that the proposed model outperforms the existing ones and provides ﬂexibility ﬁtting the experimental dataset. Another attractive aspect of the model are analytically tractable and can be easily implemented under a Bayesian approach. logit, probit, and cloglog in the Bayesian framework. The simulation result show that the glogit model always gave the best ﬁt for the datasets generated from the glogit link function and the second-best ﬁt for the datasets generated from another correspond link function. It also shows the second-best ﬁt model when the datasets generated from the glogit was the model that related to the shape of EEL p.d.f. curve determined from α and λ value.


Introduction
Binomial regression is a general term of a regression model in which the dependent variable has binomial distribution. The binomial distribution is formed from a number of successes in some repeated trials of a binomial experiment, where each trial has a probability of success. Binomial regression belongs to the generalized linear models (GLM) class, which uses a link function to connect linear predictor variables to the expectation of response variables. In binomial regression, predictor variables relate the probability of success from a series of independent Bernoulli trials. The link function for the binomial regression is formed from inverse cumulative distribution functions (c.d.f.) of a random distribution, for instance, the logit link function that corresponds to the inverse c.d.f. of the standard logistic distribution, the probit link function uses the inverse c.d.f. of standard normal distribution, and the complementary log-log (cloglog) link function is formed from the inverse c.d.f. of the Gumbel distribution. For more information about GLM and binomial regression, see McCullagh and Nelder [1] or Agresti [2].
The logit, probit, and cloglog links are the three commonly used link functions in a binomial regression. The first link is preferable due to its easy interpretation of the regression coefficients and the odds ratio. The logit model is a linear model for the canonical or natural parameter of the underlying exponential family that was obtained, and it has a closed form [3]. However, the logit model cannot always guarantee a good fit for all binomial regressions, especially for asymmetric models [4]. Czado and Santner [5], Tiku and Vaughan [6], Dobson and Barnett [7] indicated that misspecification of the link function in binomial regression could produce significant bias and increase mean squared error of parameter estimates and the predicted probabilities.
Both logit and probit links are symmetric links. Symmetric links assume that the probability of a given binomial response approaches 0 at the same rate as it approaches 1. This implies that the probability curve has a symmetric form of about 0.5. It is supported by the symmetrical probability density function (p.d.f.), which corresponds to the inverse c.d.f. of the link function. This assumption may not be reasonable in various cases. Many studies found that these symmetric link functions have some limitations. Nagler [8] indicated how sensitive the inferences are if a symmetric link was incorrectly used in the direction of an asymmetric model. Chen et al. [9] concluded that when the probability of a given binomial response approaches 0 at a different rate than it approaches 1 (depending on the covariate) or vice versa, the symmetric link function is inappropriate.
Collet [10] showed that an asymmetric link function could be more fitting than symmetric for binomial regression in some specific situations. One of the asymmetric links that is frequently used in binomial regression is the cloglog link function. However, the cloglog link has a fixed left or negative skewness, which means the probability of a given binomial response rises slowly from 0 but then tapers off more quickly as it approaches 1 and not the other way around [11]. As a result, the cloglog may produce an inaccurate inference for a symmetric model and a positive skewness model, and hence, it is not a flexible link function for binomial regression, although it is an asymmetric link. Binomial regression will be better modeled with flexible link functions that allow for symmetric and asymmetric (both positive and negative skewness) models and that allows the model to determine the required skewness level.
There have been many types of studies proposed that include skewness into the link function. Stukel [12] presented a two-parameter class of generalized logistic models that can approximate many link functions, including probit, logit, and complementary log-log. However, in a Bayesian analysis, it may not be straightforward to be implemented, because of a problem of multiple covariates and noninformative improper priors. Chen et al. [9] proposed skew-probit links using a latent variable approach of Albert and Chib [13], which can lead to proper posterior distributions using standard improper priors. Kim et al. [14] introduced a new link function based on the class of generalized skewed t-link models with the constraint on the shape parameter. Bazan et al. [15] presented the skewed probit link function and some variants with different parameterizations. Naranjo et al. [16] presented asymmetric link function based on the asymmetric exponential power (AEP) distribution. More recently, Caron et al. [17] proposed a skewed Weibull link function for categorical response data arising from the binomial model as well as the multinomial model and showed that logit, probit, and complementary log-log can be obtained as limiting cases.
In this paper, we propose a flexible generalized logit (glogit) link function for the binomial regression model. The link function is created from inverse cumulative distribution function (c.d.f.) of a new class of generalized logistic based on the exponentiated-exponential logistic (EEL) distribution that was introduced by Ghosh and Alzaatreh [18]. The EEL distribution has similar characteristics to Type IV generalized logistic distribution (GLD). The distributions consider both symmetric and asymmetric models, including the case of lighter and heavier tails, simultaneously as compared to standard logistic. The logistic distribution, Type I GLD, and Type II GLD can be obtained as special cases of both distributions. However, the EEL distribution differs from the Type IV GLD in the c.d.f. which is algebraically more tractable, so that it is more easily to be implemented. Some of the researchers use the previous GLD as the link function in binary regression, e.g., Tiku and Vaughan [6] that use Type I GLD and Oral [19] extended the work of Tiku and Vaughan [6] with a stochastic covariate, Valle et al. [20] use the Type III GLD and show that the model has ability to capture heavy and light tails, and Prentice [21] applies Type IV GLD model for Dose-Response Curves.
A Bayesian estimation technique using Markov chain Monte Carlo (MCMC) is provided for implementations of the proposed model. We present a simple syntax for the computation of the model using the JAGS (Just Another Gibbs Sampler) program shown as an attractive aspect. A simulation study to investigate the model performance of the proposed model in comparison with the most commonly used link functions, i.e., logit, probit, and complementary log-log is conducted. We also compare the proposed model with several other asymmetric models using two previously published datasets, i.e., the beetle mortality dataset [22] that was analyzed by Naranjo et al. [16] and potency of three different poisons dataset [23] analyzed by Caron et al. [17].
The rest of the paper is organized as follows: Section 2 describes the binomial regression model and the role of the link function in the model. In Section 3, we introduce the flexible generalized logit (glogit) link based on the c.d.f. corresponding to the exponentiated-exponential logistic (EEL) distribution, and we show the graphs of the distribution in various schemes. Section 4 demonstrates the flexibility of the glogit model compared with logit, probit, and complementary log-log by simulated datasets. Section 5 contains the application and comparison of the glogit model using the beetle mortality and potency of three different poisons datasets. The conclusion and discussion are given in Section 6.

Binomial Regression Model
The general notation of the binomial regression model is given as follows. Suppose y = (y 1 , y 2 , . . . , y n ) is an n × 1 vector of independent binary random variables that took value 1 for success and 0 for failure. Let X denote n × (k + 1) design matrix with rows x i , where x i = (1, x i1 , x i2 , . . . , x ik ) is a (k + 1) × 1 vector of the predictor variable, where value 1 corresponds to an intercept. Furthermore, we consider β = (β 0 , β 1 , . . . , β k ) as a (k + 1) × 1 corresponding the intercept and regression coefficient.
In binomial regression models, the probability p i = P(y i = 1) is predicted by applying link functions to a linear combination of predictor variables that leads to the following model: where F −1 (·) is called as a link function. Link function F −1 can be formed from inverse c.d.f. of a random distribution F. When F is a c.d.f. of a symmetric distribution, the resulting link is also symmetric. The most common symmetric links include logit and probit, whereas an asymmetric link is obtained when F is an asymmetric c.d.f., for example, the cloglog link. The likelihood function L β y, X is derived from the likelihood of binomial distribution as follows: thus, the log-likelihood function β y, X is given by: The β parameters in the model (3) can be estimated by using a maximum likelihood method or using the Bayesian method. The numerical method, such as the Newton-Raphson algorithm or the Nelder-Mead algorithm, can be applied for estimation with maximum likelihood estimation (MLE). While in the Bayesian framework, inferences about model parameters are typically summarized by random draws from a posterior density, which is proportional to the likelihood function for the data multiplied by a prior for the parameters.

Flexible Generalized Logit Link
The flexible generalized logit link function for the binomial regression model proposed in this paper is developed from inverse c.d.f. of a new class of generalized logistic. We start this section with a brief review of the new class of generalized logistic based on the exponentiated-exponential logistic (EEL) distribution given by Ghosh and Alzaatreh [18]. The EEL distribution is obtained using the technique of mixing two distributions of exponentiated-exponential and logistic distribution. For a more detailed discussion about this method, see Alzaatreh et al. [24] as well as Gupta and Kundu [25].

Exponentiated-Exponential Logistic Distribution
A random variable X is an EEL distribution, with parameter θ = (µ, s, α, λ), where µ ∈ and s > 0 are the locations and the scale parameter of logistic distribution, α > 0 and λ > 0 are the shapes and the scale parameter of exponentiated-exponential distribution. The probability density function (p.d.f.) of X is given by: We denote hereafter as X ∼ EEL(µ, s, α, λ). (4) is reduced to the p.d.f. of the logistic distribution. In the special case of λ = 1, the p.d.f. in (4) is the Type I GLD. Moreover, when α = 1, the p.d.f. reduces to the Type II GLD. The illustrations of some typical standard EEL p.d.f.s for different values of α and λ in µ = 0 and s = 1 are provided in Figure 1. Figure 1 shows that the density function of the EEL distribution can take different shapes. It has negative skewness for α < λ, positive skewness for α > λ, and symmetric for α = λ. The greater skewness can be achieved when the difference between α and λ is higher. Furthermore, the standard EEL p.d.f. has lighter tails compared to standard logistic p.d.f. when α > 1 and λ > 1, otherwise if α < 1 and λ < 1 then heavier tails are achieved. more detailed discussion about this method, see Alzaatreh et al. [24] as well as Gupta and Kundu [25].

Exponentiated-Exponential Logistic Distribution
A random variable X is an EEL distribution, with parameter = ( , , , ), where ∈ ℜ and > 0 are the locations and the scale parameter of logistic distribution, > 0 and > 0 are the shapes and the scale parameter of exponentiated-exponential distribution. The probability density function (p.d.f.) of X is given by: We denote hereafter as ~EEL( , , , ). (4) is reduced to the p.d.f. of the logistic distribution. In the special case of = 1, the p.d.f. in (4) is the Type I GLD. Moreover, when = 1, the p.d.f. reduces to the Type II GLD. The illustrations of some typical standard EEL p.d.f.s for different values of and in = 0 and = 1 are provided in Figure 1. Figure 1 shows that the density function of the EEL distribution can take different shapes. It has negative skewness for < , positive skewness for > , and symmetric for = . The greater skewness can be achieved when the difference between and is higher. Furthermore, the standard EEL p.d.f. has lighter tails compared to standard logistic p.d.f. when > 1 and > 1, otherwise if < 1 and < 1 then heavier tails are achieved. Based on the p.d.f. in (4), the c.d.f. for the EEL distribution is given by: Based on the p.d.f. in (4), the c.d.f. for the EEL distribution is given by: Suppose F −1 gl (p) denotes the quantile function for the EEL distribution, where 0 < p < 1. Then, the quantile function, which is an inverse of the EEL c.d.f. can be written as The result is obtained by using F gl F −1 gl (p) = p in (4) and then solving it for F −1 gl (p).

Binomial Regression Model with Glogit Link Function
A binomial regression model, applying the inverse of the EEL c.d.f. as a link function, namely glogit, is presented in this section. The glogit link function is formed from the quantile function of the standard EEL distribution at µ = 0 and s = 1. The glogit link function F −1 gl can be defined as: When α = λ = 1, the glogit link function in (7) will be reduced to the logit link function.

Binomial Regression Model with Glogit Link Function
A binomial regression model, applying the inverse of the EEL c.d.f. as a link function, namely glogit, is presented in this section. The glogit link function is formed from the quantile function of the standard EEL distribution at = 0 and = 1. The glogit link function −1 can be defined as: When = = 1, the glogit link function in (7) will be reduced to the logit link function.  The glogit class for the binomial regression model is obtained by considering −1 in the specification of the binomial regression model given in (1), which is expressed as: We can rewrite the formula in (8) where (•) is the c.d.f. of the EEL distribution given in (5)  The glogit class for the binomial regression model is obtained by considering F −1 gl in the specification of the binomial regression model given in (1), which is expressed as: We can rewrite the formula in (8) as: . . , n.
where F gt (·) is the c.d.f. of the EEL distribution given in (5) in case µ = 0 and s = 1. The likelihood function for the glogit model can be formed by substituting (9) in (2) as: Thus, the log-likelihood function (β, α, λ y, X) is given as: where β, α, and λ in the model (11) are the parameters to be estimated by using a maximum likelihood estimation (MLE) method, or a Bayesian method as previously mentioned in Section 2.
In this paper, we conduct a Bayesian method for the simulation and applications. The Bayesian analysis requires the joint posterior density π(β, α, λ y, X) that is defined by: π(β, α, λ y, X) ∝ L(β, α, λ y, X)π(β)π(α)π(λ), where π(β) is the join prior of regression coefficients, π(α) is prior for the α parameter, and π(λ) is prior for the λ parameter. For the priors of the intercept and regression coefficients, we can adopt the typical proper priors for β, a normal prior, and thus β j ∼ N µ * j , τ * j , where j = 0, 1, . . . , k, µ * j = 0, and τ * j denotes the precision (1/σ 2 * j ) with a small value representing the common situation that is a lack of information available for the priors. Alternatively, informative priors are obtained when µ * j and τ * j are elicited based on historical information, expert opinion or validation from datasets. For α and λ, we can use a proper lognormal prior with the mean 0 and variance 1 ( α ∼ LN(0, 1); λ ∼ LN(0, 1)), for the reasons that their value is positive and around to 1 (get near to the logit model).
With the model setup of prior distributions as shown above, the joint posterior distribution of the model is given by: Using the proportionality follows from Bayes' theorem, the joint posterior distribution (13) can be rewritten as: The analysis of Bayesian models is using MCMC, which sampling from the posterior is done using the Gibbs sampling algorithm. The full conditional distributions for all parameters must be derived to implement the algorithm for the joint posterior distribution (14). The full conditional distributions for β are easy to obtain: where β (−j) = β 0 , β 1 , . . . , β j−1 , β j+1 , . . . , β k and the subscript (−j) denotes the j-th element has been removed. The full conditional distributions for α and λ are given by: Each draw can be done using the adaptive rejection metropolis algorithm implemented in JAGS. JAGS is a program for Bayesian graphical modeling, which has compatibility with classic BUGS (Bayesian inference using Gibbs sampling), developed by Martyn Plummer [26]. We can express the models for k = 2 in JAGS using the following simple syntax:  ])),-lambda)),alpha) } #Prior beta0~dnorm(0, 0.001) beta1~dnorm(0, 0.001) beta2~dnorm(0, 0.001) alpha~dlnorm(0, 1) lambda~dlnorm(0, 1) } We implemented the models using the function autojags in the jagsUI package developed by Kellner [27] in the R software. The autojags applies the convergence diagnostic based on Gelman and Rubin [28] and Brooks and Gelman [29] called Gelman-Rubin diagnostic tests. To compute the Gelman-Rubin statistic, we must run two chains in JAGS using two different sets of initial values. If the MCMC chain has reached convergence, the Gelman-Rubin statistics, namely Rhat, are lower than 1.05, and it concludes that burn-in is adequate, whereas the value in excess of 1.05 indicates a lack of convergence. Another convergence diagnostic method of Raftery and Lewis [30] was considered in autojags. The method estimates and applies the number of iterations needed for a given level of precision in posterior samples when quantiles are the posterior summaries of interest.

Simulation
A simulation study was carried out to evaluate the performance of the glogit link function for the binomial regression models compared to logit, probit, and cloglog. For this purpose, several datasets were generated by the following steps:

•
Step 2: Generate x i1 ∼ Normal(0, 1), then we create x i = (1, x i1 ) ; • Step 3: Compute p i from a given link, so that p i = F(x i β); • Step 4: Generate n i ∼ Poisson(δ) with δ = 100, and Y i ∼ Binomial(n i , p i ), which is the binomial distribution with p i as the probability of success and n i is the sample size for the repeated Bernoulli trials.
We consider seven simulations to accommodate into seven scenarios for symmetric and asymmetric models and lighter and heavier tails (compared to the logit model). The scenarios were considered by applying the different inverse of link functions F(·). The scenarios are as follows: The comparison of graph visualizations for scenarios 1 to 5 represented specifically for the glogit can be seen in Figures 1 and 2. For each scenario, we generated 100 datasets, and for each dataset, we generated m = 100 binomial response variables, thus, y i ∼ Binomial(n i , p i ) for i = 1, . . . , 100. Furthermore, we fit all of the models for each simulated dataset, respectively, to find the best performance in terms of model comparisons. We conducted a Bayesian analysis on the above seven simulated datasets and set with the uninformative prior for each model given in Section 3.2.
To compare the performance of the models under different link functions, we use two model comparisons. The first is the deviance information criterion (DIC) proposed by Spiegelhalter et al. [31], which represents a goodness-of-fit criterion. It is an estimate of expected predictive error. Therefore, the lower value of DIC is better. DIC is given by DIC = D(η) +p D, where D(η) = −2 log L(η) is deviance of model, L(η) is the likelihood, and D(η) = E D(η) data is the deviance posterior mean. The number of parameters is given as p, andp D = D(η) − D(η) is defined as the effective number of parameters, while D(η) is the deviance posterior means of the parameters of interest, with η = E η data . Another criterion for the model comparison is absolute errors (AE). Absolute error is defined as the absolute value of the difference between the prediction value and the true value of a measurement and is commonly given as the maximum possible error given a model's degree of accuracy. AE is defined by AE = n i=1 p i −p i . The smaller values of AE indicate better models. Table 1 shows the result of the simulation for the seven scenarios. Scenarios 1 to 3 are symmetrical models in which the first one is generated from glogit(α = 1; λ = 1), or same with the logit. The second one represents the lighter tails compared to the logit, and the third one represents the heavier tails compared to the logit. The fourth and fifth scenarios are used to represent the asymmetric models in which the fourth one has a positive direction of skewness, while the fifth one has a negative direction. In addition, the sixth and seventh scenarios belong to the cloglog and probit link function. We evaluate each scenario by DIC and AE values to compare the models, i.e., glogit, logit, cloglog, and probit.
First, we focus on the results of the first scenario. When the datasets generated from the logit link that equaled the glogit link for α = λ = 1, then the logit model is the best, followed by the glogit. However, the DIC and AE values of both models are very close compared to the other. For the second to fifth scenario that used datasets generated from the glogit model for the several values of α and λ, the glogit model performs best. In the second scenario, which has lighter tails compared to the standard logistic, the second-best fit model is the probit model. For the third scenario, the logit model gives almost identical values of DIC and AE, with the result that it is the second-best fit model. In the fourth scenario, the second-best fit is the probit model, whereas, in the fifth one, the second-best fit is the cloglog model. This result is consistent with the fact that the cloglog model also has negative skewness. When the datasets are generated from cloglog and probit models (Scenario 6 and 7), the second-best fit is the glogit model after them. However, the differences in the model comparison values are quite similar. In summary, the simulation results show that the glogit model is quite flexible and performs better than the standard link function.

Application
This section discusses the result of applying the proposed method to the empirical cases. It aims to study the performance of the flexible generalized logit link function, compared with other recent asymmetric link functions. Two datasets are considered, i.e., beetle mortality dataset of Bliss [22] and potency of three different poison datasets of Finney [23].

Beetle Mortality Dataset
The beetle mortality dataset in Table 2 is taken from Bliss [22]. The data consist of the toxicological experiment results based on the number of flour beetles dead after several hours of carbon disulfide (CS 2 ) exposure at various concentrations (transformed to ten base logarithm of dose). More information about the data can be found in Bliss [22]. The dataset in Table 2 was previously analyzed by Naranjo et al. [16], who proposed the AEP-based link model to analyze a binomial regression by using the inverse of the Asymmetric Exponential Power (AEP) c.d.f. as the link function. Furthermore, they compared the model with some symmetric and asymmetric link models, such as the probit and t-link (T(8)) models that are proposed by Albert and Chib [13], the skewed probit (SP) model that are proposed by Chen et al. [9], and the skewed generalized t-link (SGT) model proposed by Kim et al. [14]. The model is given by where the response variable is binary with (y i = 1) representing the dead and (y i = 0) representing the beetles that survived, p i is the proportion of dead beetles (ndead i ) to the number of beetles exposed (n i ) under the i th dose. A comparison of the five models was considered by using different link functions F −1 . Five models were fitted to the data, and four different criteria of model comparison were used, e.g., the deviance information criterion (DIC), the Bayesian information criterion (BIC), the Akaike information criterion (AIC), and the absolute errors (AE). However, in this paper, only DIC and AE are used to compare with the flexible generalized logit (glogit) model because the DIC, BIC, and AIC give the same conclusions. The formulas of DIC and AE are given in Section 4. All of the estimation and comparison criteria values for the models, except the glogit model, are adopted from Naranjo et al. [16] that applied Bayesian analysis to estimate the parameters. For the glogit model, we also used a Bayesian framework, whose join posterior is given in (12), including the settings for the prior.
The posterior estimates (mean and standard deviation), the DIC value, and the AE value for the fitted models are reported in Table 3. The glogit model is fitted to the dataset, and the result is compared with those for the probit, T(8), SP, SGT, and AEP-based link models processed by Naranjo et al. [16]. These results presented in Table 3 show that the glogit model provides the best fit, followed by the AEP-based link model indicated by the smallest value of DIC. For criterion AE, we note that the AEP-based link model has the smallest value so that the AEP-based link model is superior, followed by the glogit model. Furthermore, the glogit and AEP-based link models have similar values of model comparison compared with the other in this dataset.

The Potency of Three Different Poisons Dataset
In the second application, we take the potency of three different poisons dataset from Finney [23] and the analysis results from Caron et al. [17]. The data consist of the test results of the different doses and the different poisons to know the potency of the poisons to kill insects. The poisons are rotenone, deguelin, and a mixture of rotenone and deguelin. The data are provided in Table 4, where n is the number of observations. More information about this data can be found in Finney [23].
Caron et al. [17] proposed a skewed Weibull (SW) link function and analyzed the dataset in Table 4 in MLE and Bayesian framework. In the Bayesian method, they compared the proposed model with another skewed link function, such as the AEP-based link (AEP) model presented by Naranjo et al. [16] and the skewed probit (SP) model proposed by Chen et al. [9]. The model that was used is: where the p i is the proportion of dead insects (dead i ) to the number of observation (n i ) under the i th log(dose), rotenone dummy variable, and deguelin dummy variable. The mixture variable is considered as the reference poison (rotenone = 0 and deguelin = 0). A comparison to the models was considered by applying different F −1 (link functions). They are the SP, AEP, SW, and glogit models. In this section, we applied the DIC and the mean absolute error (MAE) for model comparison. The MAE is defined as MAE = 1 ⁄n n i=1 p i −p i . The estimation and the model comparison values for the models, except the glogit model, were taken from Caron et al. [17]. More detail about their methods can be found in Caron et al.'s study [17]. For the glogit model, we used the model presented in (12), including the settings for the prior.
The model comparison results of DIC and MAE for the fitted models are presented in Table 5. The glogit model is fitted to the dataset, and the result is compared with SP, AEP, and SW models analyzed by Caron et al. [17]. The results reported in Table 5 show that the glogit model provides the best fit, followed by the AEP-based link model, indicated by the smallest value of DIC and MAE. Overall, we see the superiority of the glogit against other link functions for the datasets. The mean and standard deviation of posterior estimates values for all of the models are also given in Table 5.

Conclusions
A flexible generalized logit (glogit) link function for analyzing binomial regression is proposed in this paper. The link function is based on a new class of generalized logistic distribution, exponentiated-exponential logistic (EEL) distribution. This distribution considers both symmetric and asymmetric models to accommodate positive or negative skewness. Moreover, this allows for the case of a lighter or heavier tails behavior compared to the standard logistic.
The paper provides a simulation of binomial regression using seven scenarios to indicate the flexibility and superiority of the glogit compared with common link functions, i.e., logit, probit, and cloglog in the Bayesian framework. The simulation result show that the glogit model always gave the best fit for the datasets generated from the glogit link function and the second-best fit for the datasets generated from another correspond link function. It also shows the second-best fit model when the datasets generated from the glogit was the model that related to the shape of EEL p.d.f. curve determined from α and λ value.
In addition, applying the model for the datasets from previous studies indicates that the glogit fits the data well and outperforms the other asymmetric link functions, such as the probit and t-link (T(8)) models proposed by Albert and Chib [13], the skewed probit (SP) model proposed by Chen et al. [9], the skewed generalized t-link (SGT) model proposed by Kim et al. [14], AEP-based link (AEP) model presented by Naranjo et al. [16], and finally, a skewed Weibull (SW) model presented by Caron et al. [17]. The proposed model is evidently flexible and capable of handling the datasets. The performance of the glogit is very good, indicated by the smallest value of DIC, followed by the AEP-based link model for both of the datasets. We are convinced that the glogit model is a good option and could be used in practice due to its flexibility. Other attractive aspects of the model are that it is analytically tractable and can be easily implemented under a Bayesian approach using the simple syntax.
This study adopted the typical priors for the Bayesian analysis with non-informative priors. Priors with other distributions can also be considered. The prior selection is an important activity that was not analyzed in this paper and hence, can be considered in future research.
Author Contributions: R.B.P., H.K., N.I., and B.S.S.U. understood and designed the research; R.B.P. collected and analyzed the data, and drafted the paper. All authors critically read and revised the draft and approved the final paper. All authors have read and agreed to the published version of the manuscript.