Stable monomers in the ancestral sequence reconstruction of the last opisthokont common ancestor of dimeric triosephosphate isomerase.

Function and structure are strongly coupled in obligated oligomers such as Triosephosphate isomerase (TIM). In animals and fungi, TIM monomers are inactive and unstable. Previously, we used ancestral sequence reconstruction to study TIM evolution and found that before these lineages diverged, the last opisthokonta common ancestor of TIM (LOCATIM) was an obligated oligomer that resembles those of extant TIMs. Notably, calorimetric evidence indicated that ancestral TIM monomers are more structured than extant ones. To further increase confidence about the function, structure, and stability of the LOCATIM, in this work, we applied two different inference methodologies and the worst plausible case scenario for both of them, to infer four sequences of this ancestor and test the robustness of their physicochemical properties. The extensive biophysical characterization of the four reconstructed sequences of LOCATIM showed very similar hydrodynamic and spectroscopic properties, as well as ligand-binding energetics and catalytic parameters. Their 3D structures were also conserved. Although differences were observed in melting temperature, all LOCATIMs showed reversible urea-induced unfolding transitions, and for those that reached equilibrium, high conformational stability was estimated (ΔGTot = 40.6-46.2 kcal/mol). The stability of the inactive monomeric intermediates was also high (ΔGunf = 12.6-18.4 kcal/mol), resembling some protozoan TIMs rather than the unstable monomer observed in extant opisthokonts. A comparative analysis of the 3D structure of ancestral and extant TIMs shows a correlation between the higher stability of the ancestral monomers with the presence of several hydrogen bonds located in the "bottom" part of the barrel.

Hueso nasal como marcador de aneuploidías y resultado perinatal / Nasal bone as a marker of aneuploidies and perinatal outcome

Objetivo: Describir el resultado perinatal de los embarazos en función de la evaluación del hueso nasal como marcador de aneuploidía. Métodos: De 1006 embarazadas, 607 cumplieron con los criterios de inclusión para este estudio prospectivo, descriptivo, correlacional no causal donde se correlacionó la ausencia/presencia de hueso nasal con la presencia de síndrome de Down a través de cariotipo fetal prenatal y/o posnatal, así como examen clínico neonatal. Los datos fueron analizados mediantes frecuencias absolutas, porcentajes, capacidad diagnóstica del hueso nasal (índice de Youden), sensibilidad, especificidad, valor predictivo positivo, valor predictivo negativo y cocientes de probabilidad, positivo y negativo. Resultados: La prevalencia de síndrome de Down fue de 1,48 %, la ausencia del hueso nasal como marcador aislado, obtuvo un índice de Youden de 0,55 (0,23 - 0,88), sensibilidad de 55,56 %, especificidad de 99,50 %, valor predictivo positivo de 62,5 %, valor predictivo negativo de 99,33 %, cocientes de probabilidad positivo (hueso nasal ausente) 111 (IC 95 % 31 - 394) y cocientes de probabilidad negativo (hueso nasal presente) de 0,45 (IC 95 % 0,22 -0,93). Conclusión: La ausencia de hueso nasal en primer trimestre aumenta el riesgo de síndrome de Down en 111 veces y la presencia del mismo lo disminuye, sin valor como prueba diagnóstica sino de pesquisa debe considerarse como un marcador secundario(AU)

Objective: To know the perinatal outcome based on nasal bone evaluation as an aneuploidy marker. Methods: From 1006 pregnant women, 607 met the inclusion criteria for this prospective, descriptive, correlational not causal research correlating nasal bone absence / presence with Down syndrome through prenatal / postnatal fetal karyotype and neonatal clinical examination. Absolute frequencies and percentages, nasal bone performance as a diagnostic test (Youden índex), sensitivity, specificity, positive predictive value, negative predictive value, likelihood ratios positive and negative, were calculated. Results: 1.48 % was the Down syndrome prevalence on the sample. The nasal bone absence as an isolated marker obtained an 0,55 Youden index (0.23 to 0.88 ), sensitivity 55,56%, specificity 99,50%, positive predictive value 62,5%, negative predictive value 99,33%, likelihood ratios positive (absent nasal bone) 111, (95% CI 31-394) and likelihood ratios negative (nasal bone present ) 0,45 (95% CI 0 22 -0.93 ). Conclusion: The nasal bone absence in first trimester increases Down syndrome risk 111 times and nasal bone presence decreases it with poor performance as a diagnostic test, so it should be considered a screening test and a secondary marker. Recommendations correlate these results with other markers to improve detection rates and quantify nasal bone measurements in order to make nasal bone nomograms in first trimester pregnancies(AU)

Conformal normal curvature and detection of masked observations in multivariate null intercept measurement error models.

Measurement errors occur very commonly in practice. After fitting the model, influence diagnostics is an important step in statistical data analysis. The most frequently used diagnostic method for measurement error models is the local influence. However, this methodology may fail to detect masked influential observations. To overcome this limitation, we propose the use of the conformal normal curvature with the forward search algorithm. The results are presented through easy to interpret plots considering different perturbation schemes. The proposed methodology is illustrated with three real data sets and one simulated data set, two of which have been previously analyzed in the literature. The third data set deals with the stability of the hygroscopic solid dosage in pharmaceutical processes to ensure the maintenance of product safety quality. In this application, the analytical mass balance is subject to measurement errors, which require attention in the modeling process and diagnostic analysis.

A targeted likelihood estimation comparing cefepime and piperacillin/tazobactam in critically ill patients with community-acquired pneumonia (CAP).

Cefepime and piperacillin/tazobactam are antimicrobials recommended by IDSA/ATS guidelines for the empirical management of patients admitted to the intensive care unit (ICU) with community-acquired pneumonia (CAP). Concerns have been raised about which should be used in clinical practice. This study aims to compare the effect of cefepime and piperacillin/tazobactam in critically ill CAP patients through a targeted maximum likelihood estimation (TMLE). A total of 2026 ICU-admitted patients with CAP were included. Among them, (47%) presented respiratory failure, and (27%) developed septic shock. A total of (68%) received cefepime and (32%) piperacillin/tazobactam-based treatment. After running the TMLE, we found that cefepime and piperacillin/tazobactam-based treatments have comparable 28-day, hospital, and ICU mortality. Additionally, age, PTT, serum potassium and temperature were associated with preferring cefepime over piperacillin/tazobactam (OR 1.14 95% CI [1.01-1.27], p = 0.03), (OR 1.14 95% CI [1.03-1.26], p = 0.009), (OR 1.1 95% CI [1.01-1.22], p = 0.039) and (OR 1.13 95% CI [1.03-1.24], p = 0.014)]. Our study found a similar mortality rate among ICU-admitted CAP patients treated with cefepime and piperacillin/tazobactam. Clinicians may consider factors such as availability and safety profiles when making treatment decisions.

Clear cell likelihood score may improve diagnosis and management of renal masses.

The detection of solid renal masses has increased over time due to incidental findings during imaging studies conducted for unrelated medical conditions. Approximately 20% of lesions measuring less than 4 cm are benign and 80% are malignant. Clear cell renal cell carcinoma (ccRCC) is the most frequent among renal carcinomas, responsible for 65-80% of cases. The increased detection of renal masses facilitates early diagnosis and treatment. However, it also leads to more invasive interventions, which result in higher morbidity and costs. Currently, only histological analysis can offer an accurate diagnosis. Surgical nephron loss significantly elevates morbidity and mortality rates. Active surveillance represents a conservative management approach for patients diagnosed with a solid renal mass that is endorsed by both American Urological Association and the European Society for Medical Oncology. However, active surveillance is used in a minority of patients and varies across institutions. The lack of clinical studies using a standardized approach to incidentally detected small renal masses precludes the widespread use of active surveillance. Hence, there is an urgent need for better patient selection, distinguishing those who require surgery from those suitable for active surveillance. The clear cell likelihood score (ccLS) represents a novel MRI tool for assessing the probability of a renal mass being a ccRCC. In this study, we present a comprehensive review of renal masses and their evaluation using the ccLS to facilitate shared decision between urologists and patients.

Family of bivariate distributions on the unit square: theoretical properties and applications.

We introduce the bivariate unit-log-symmetric model based on the bivariate log-symmetric distribution (BLS) defined in Vila et al. [25] as a flexible family of bivariate distributions over the unit square. We then study its mathematical properties such as stochastic representations, quantiles, conditional distributions, independence of the marginal distributions and marginal moments. Maximum likelihood estimation method is discussed and examined through Monte Carlo simulation. Finally, the proposed model is used to analyze some soccer data sets.

Continuous-Time Model Identification of the Subglottal System.

Mathematical models that accurately simulate the physiological systems of the human body serve as cornerstone instruments for advancing medical science and facilitating innovative clinical interventions. One application is the modeling of the subglottal tract and neck skin properties for its use in the ambulatory assessment of vocal function, by enabling non-invasive monitoring of glottal airflow via a neck surface accelerometer. For the technique to be effective, the development of an accurate building block model for the subglottal tract is required. Such a model is expected to utilize glottal volume velocity as the input parameter and yield neck skin acceleration as the corresponding output. In contrast to preceding efforts that employed frequency-domain methods, the present paper leverages system identification techniques to derive a parsimonious continuous-time model of the subglottal tract using time-domain data samples. Additionally, an examination of the model order is conducted through the application of various information criteria. Once a low-order model is successfully fitted, an inverse filter based on a Kalman smoother is utilized for the estimation of glottal volume velocity and related aerodynamic metrics, thereby constituting the most efficient execution of these estimates thus far. Anticipated reductions in computational time and complexity due to the lower order of the subglottal model hold particular relevance for real-time monitoring. Simultaneously, the methodology proves efficient in generating a spectrum of aerodynamic features essential for ambulatory vocal function assessment.

A novel flexible exponent power-X family of distributions with applications to COVID-19 mortality rate in Mexico and Canada.

This paper aims to introduce a novel family of probability distributions by the well-known method of the T-X family of distributions. The proposed family is called a "Novel Generalized Exponent Power X Family" of distributions. A three-parameters special sub-model of the proposed method is derived and named a "Novel Generalized Exponent Power Weibull" distribution (NGEP-Wei for short). For the proposed family, some statistical properties are derived including the hazard rate function, moments, moment generating function, order statistics, residual life, and reverse residual life. The well-known method of estimation, the maximum likelihood estimation method is used for estimating the model parameters. Besides, a comprehensive Monte Carlo simulation study is conducted to assess the efficacy of this estimation method. Finally, the model selection criterion such as Akaike information criterion (AINC), the correct information criterion (CINC), the Bayesian information criterion (BINC), the Hannan-Quinn information criterion (HQINC), the Cramer-von-Misses (CRMI), and the ANDA (Anderson-Darling) are used for comparison purpose. The comparison of the NGEP-Wei with other rival distributions is made by Two COVID-19 data sets. In terms of performance, we show that the proposed method outperforms the other competing methods included in this study.

A model for bimodal rates and proportions.

The beta model is the most important distribution for fitting data with the unit interval. However, the beta distribution is not suitable to model bimodal unit interval data. In this paper, we propose a bimodal beta distribution constructed by using an approach based on the alpha-skew-normal model. We discuss several properties of this distribution, such as bimodality, real moments, entropies and identifiability. Furthermore, we propose a new regression model based on the proposed model and discuss residuals. Estimation is performed by maximum likelihood. A Monte Carlo experiment is conducted to evaluate the performances of these estimators in finite samples with a discussion of the results. An application is provided to show the modelling competence of the proposed distribution when the data sets show bimodality.

On the Nuisance Parameter Elimination Principle in Hypothesis Testing.

The Non-Informative Nuisance Parameter Principle concerns the problem of how inferences about a parameter of interest should be made in the presence of nuisance parameters. The principle is examined in the context of the hypothesis testing problem. We prove that the mixed test obeys the principle for discrete sample spaces. We also show how adherence of the mixed test to the principle can make performance of the test much easier. These findings are illustrated with new solutions to well-known problems of testing hypotheses for count data.

Performance of tree-building methods using a morphological dataset and a well-supported Hexapoda phylogeny.

Recently, many studies have addressed the performance of phylogenetic tree-building methods (maximum parsimony, maximum likelihood, and Bayesian inference), focusing primarily on simulated data. However, for discrete morphological data, there is no consensus yet on which methods recover the phylogeny with better performance. To address this lack of consensus, we investigate the performance of different methods using an empirical dataset for hexapods as a model. As an empirical test of performance, we applied normalized indices to effectively measure accuracy (normalized Robinson-Foulds metric, nRF) and precision, which are measured via resolution, one minus Colless' consensus fork index (1-CFI). Additionally, to further explore phylogenetic accuracy and support measures, we calculated other statistics, such as the true positive rate (statistical power) and the false positive rate (type I error), and constructed receiver operating characteristic plots to visualize the relationship between these statistics. We applied the normalized indices to the reconstructed trees from the reanalyses of an empirical discrete morphological dataset from extant Hexapoda using a well-supported phylogenomic tree as a reference. Maximum likelihood and Bayesian inference applying the k-state Markov (Mk) model (without or with a discrete gamma distribution) performed better, showing higher precision (resolution). Additionally, our results suggest that most available tree topology tests are reliable estimators of the performance measures applied in this study. Thus, we suggest that likelihood-based methods and tree topology tests should be used more often in phylogenetic tree studies based on discrete morphological characters. Our study provides a fair indication that morphological datasets have robust phylogenetic signal.

An Optimal Integral Controller for Adaptive Optics Systems.

Integral controllers are commonly employed in astronomical adaptive optics. This work presents a novel tuning procedure for integral controllers in adaptive optics systems which relies on information about the measured disturbances. This tuning procedure consists of two main steps. First, it models and identifies measured disturbances as continuous-time-damped oscillators using Whittles´s likelihood and the wavefront sensor output signal. Second, it determines the integral controller gain of the adaptive optics system by minimizing the output variance. The effectiveness of this proposed method is evaluated through theoretical examples and numerical simulations conducted using the Object-Oriented Matlab Adaptive Optics toolbox. The simulation results demonstrate that this approach accurately estimates the disturbance model and can reduce the output variance. Our proposal results in improved performance and better astronomical images even in challenging atmospheric conditions. These findings significantly contribute to adaptive optics system operations in astronomical observatories and establish our procedure as a promising tool for fine-tuning integral controllers in astronomical adaptive optics systems.

Negative log-binomial model with optimal robust variance to estimate the prevalence ratio, in cross-sectional population studies.

BACKGROUND: Cross-sectional studies are useful for the estimation of prevalence of a particular event with concerns in specific populations, as in the case of diseases or other public health interests. Most of these studies have been carried out with binary binomial logistic regression model which estimates OR values that could be overestimated due to the adjustment of the model. Thus, the selection of the best multivariate model for cross-sectional studies is a priority to control the overestimation of the associations. METHODS: We compared the precision of the estimates of the prevalence ratio (PR) of the negative Log-binomial model (NLB) with Mantel-Haenszel (MH) and the regression models Cox, Log-Poisson, Log-binomial, and the OR of the binary logistic regression in population-based cross-sectional studies. The prevalence from a previous cross-sectional study carried out in Colombia about the association of mental health disorders with the consumption of psychoactive substances (e.g., cocaine, marijuana, cigarette, alcohol and risk of consumption of psychoactive substances) were used. The precision of the point estimates of the PR was evaluated for the NLB model with robust variance estimates, controlled with confounding variables, and confidence interval of 95%. RESULTS: The NLB model adjusted with robust variance showed accuracy in the measurements of crude PRs, standard errors of estimate and its corresponding confidence intervals (95%CI) as well as a high precision of the PR estimate and standard errors of estimate after the adjustment of the model by grouped age compared with the MH PR estimate. Obtained PRs and 95%CI entre NLB y MH were: cocaine consumption (2.931,IC95%: 0.723-11.889 vs. 2.913, IC95%: 0.786-12.845), marijuana consumption (3.444, IC95%: 1.856-6.391 vs. 3.407, IC95%: 1.848, 6.281), cigarette smoking (2.175,IC95%: 1.493, 3.167 vs. 2.209, IC95%: 1.518-3.214), alcohol consumption (1.243,IC95%: 1.158-1.334 vs. 1.241, IC95%: 1.157-1.332), and risk of consumption of psychoactive substances (1.086, IC95%: 1.047-1.127 vs. 1.086, IC95%: 1.047, 1.126). The NLB model adjusted with robust variance showed mayor precision when increasing the prevalence, then the other models with robust variance with respect to MH. CONCLUSIONS: The NLB model with robust variance was shown as a powerful strategy for the estimation of PRs for cross-sectional population-based studies, as high precision levels were identified for point estimators, standard errors of estimate and its corresponding confidence intervals, after the adjustment of confounding variables. In addition, it does not represent convergence issues for high prevalence cases (as it occur with the Log-binomial model) and could be considered in cases of overdispersion and with greater precision and goodness of fit than the other models with robust variance, as it was shown with the data set of the cross-sectional study used in here.

CBIR-SAR System Using Stochastic Distance.

This article proposes a system for Content-Based Image Retrieval (CBIR) using stochastic distance for Synthetic-Aperture Radar (SAR) images. The methodology consists of three essential steps for image retrieval. First, it estimates the roughness (α^) and scale (γ^) parameters of the GI0 distribution that models SAR data in intensity. The parameters of the model were estimated using the Maximum Likelihood Estimation and the fast approach of the Log-Cumulants method. Second, using the triangular distance, CBIR-SAR evaluates the similarity between a query image and images in the database. The stochastic distance can identify the most similar regions according to the image features, which are the estimated parameters of the data model. Third, the performance of our proposal was evaluated by applying the Mean Average Precision (MAP) measure and considering clippings from three radar sensors, i.e., UAVSAR, OrbiSaR-2, and ALOS PALSAR. The CBIR-SAR results for synthetic images achieved the highest MAP value, retrieving extremely heterogeneous regions. Regarding the real SAR images, CBIR-SAR achieved MAP values above 0.833 for all polarization channels for image samples of forest (UAVSAR) and urban areas (ORBISAR). Our results confirmed that the proposed method is sensitive to the degree of texture, and hence, it relies on good estimates. They are inputs to the stochastic distance for effective image retrieval.

Likelihood-Based Tests of Species Tree Hypotheses.

Likelihood-based tests of phylogenetic trees are a foundation of modern systematics. Over the past decade, an enormous wealth and diversity of model-based approaches have been developed for phylogenetic inference of both gene trees and species trees. However, while many techniques exist for conducting formal likelihood-based tests of gene trees, such frameworks are comparatively underdeveloped and underutilized for testing species tree hypotheses. To date, widely used tests of tree topology are designed to assess the fit of classical models of molecular sequence data and individual gene trees and thus are not readily applicable to the problem of species tree inference. To address this issue, we derive several analogous likelihood-based approaches for testing topologies using modern species tree models and heuristic algorithms that use gene tree topologies as input for maximum likelihood estimation under the multispecies coalescent. For the purpose of comparing support for species trees, these tests leverage the statistical procedures of their original gene tree-based counterparts that have an extended history for testing phylogenetic hypotheses at a single locus. We discuss and demonstrate a number of applications, limitations, and important considerations of these tests using simulated and empirical phylogenomic data sets that include both bifurcating topologies and reticulate network models of species relationships. Finally, we introduce the open-source R package SpeciesTopoTestR (SpeciesTopology Tests in R) that includes a suite of functions for conducting formal likelihood-based tests of species topologies given a set of input gene tree topologies.

Fast inference for robust nonlinear mixed-effects models.

The interest for nonlinear mixed-effects models comes from application areas as pharmacokinetics, growth curves and HIV viral dynamics. However, the modeling procedure usually leads to many difficulties, as the inclusion of random effects, the estimation process and the model sensitivity to atypical or nonnormal data. The scale mixture of normal distributions include heavy-tailed models, as the Student-t, slash and contaminated normal distributions, and provide competitive alternatives to the usual models, enabling the obtention of robust estimates against outlying observations. Our proposal is to compare two estimation methods in nonlinear mixed-effects models where the random components follow a multivariate scale mixture of normal distributions. For this purpose, a Monte Carlo expectation-maximization algorithm (MCEM) and an efficient likelihood-based approximate method are developed. Results show that the approximate method is much faster and enables a fairly efficient likelihood maximization, although a slightly larger bias may be produced, especially for the fixed-effects parameters. A discussion on the robustness aspects of the proposed models are also provided. Two real nonlinear applications are discussed and a brief simulation study is presented.

Genetic association between feed efficiency, growth, scrotal circumference, and carcass traits in Guzerat cattle.

The objective of this study was to estimate the genetic parameters for feed efficiency-related traits and their genetic correlations with growth, male fertility, and carcass traits using multi-trait analysis in Guzerat cattle. Further, it aimed to predict the direct and correlated responses for feed efficiency traits when selection was applied for growth, male fertility, and carcass traits. The evaluated traits were adjusted weight at 120 (W120), 210 (W210), 365 (W365), and 450 days of age (W450), adjusted scrotal circumference at 365 days of age (SC365) and at 450 days of age (SC450), scrotal circumference, ribeye area (REA), backfat thickness (BFT), rump fat thickness (RFT), residual feed intake (RFI), and dry matter intake (DMI). The genetic parameters were obtained by the restricted maximum likelihood method (REML), using an animal model in multi-trait analyses. The heritability estimates for W120, W210, W365, W450, SC365, and SC450 varied from low to high (0.17 to 0.39). The carcass traits, REA, BFT, and RFT, displayed low to moderate heritability estimates, 0.27, 0.10, and 0.31, respectively. The heritability estimates for RFI (0.15) and DMI (0.23) were low and moderate, respectively. The RFI showed low genetic correlations with growth traits, ranging from - 0.07 to 0.22, from 0.03 to 0.05 for scrotal circumference, and from - 0.35 to 0.16 for carcass, except for DMI, which ranged from 0.42 to 0.46. The RFI and DMI presented enough additive genetic variability to be used as selection criteria in Guzerat breed genetic improvement program. Additionally, the response to selection for RFI would be higher when selection is performed directly for this trait. The selection for residual feed intake would not promote unfavorable correlated responses for scrotal circumference, carcass (yield and finish), and growth traits. Therefore, the selection for more efficient animals would not compromise the productive, reproductive, and carcass performance, contributing to reduce the production costs, increasing the profitability and sustainability of beef cattle production in tropical areas.

The unit generalized half-normal quantile regression model: formulation, estimation, diagnostics, and numerical applications.

In this paper, we propose and derive a new regression model for response variables defined on the open unit interval. By reparameterizing the unit generalized half-normal distribution, we get the interpretation of its location parameter as being a quantile of the distribution. In addition, we can evaluate effects of the explanatory variables in the conditional quantiles of the response variable as an alternative to the Kumaraswamy quantile regression model. The suitability of our proposal is demonstrated with two simulated examples and two real applications. For such data sets, the obtained fits of the proposed regression model are compared with that provided by a Kumaraswamy regression model.

Comparison of the Effective Reproduction Number (Rt) Estimation Methods of COVID-19 Using Simulation Data Based on Available Data from Iran, USA, UK, India, and Brazil.

BACKGROUND: Accurate determination of the effective reproduction number (Rt) is a very important strategy in the epidemiology of contagious diseases, including coronavirus disease 2019 (COVID-19). This study compares different methods of estimating the Rt of susceptible population to identify the most accurate method for estimating Rt. STUDY DESIGN: A secondary study. METHODS: The value of Rt was estimated using attack rate (AR), exponential growth (EG), maximum likelihood (ML), time-dependent (TD), and sequential Bayesian (SB) methods, for Iran, the United States, the United Kingdom, India, and Brazil from June to October 2021. In order to accurately compare these methods, a simulation study was designed using forty scenarios. RESULTS: The lowest mean square error (MSE) was observed for TD and ML methods, with 15 and 12 cases, respectively. Therefore, considering the estimated values of Rt based on the TD method, it was found that Rt values in the United Kingdom (1.33; 95% CI: 1.14-1.52) and the United States (1.25; 95% CI: 1.12-1.38) substantially have been more than those in other countries, such as Iran (1.07; 95% CI: 0.95-1.19), India (0.99; 95% CI: 0.89-1.08), and Brazil (0.98; 95% CI: 0.84-1.14) from June to October 2021. CONCLUSION: The important result of this study is that TD and ML methods lead to a more accurate estimation of Rt of population than other methods. Therefore, in order to monitor and determine the epidemic situation and have a more accurate prediction of the incidence rate, as well as control COVID-19 and similar diseases, the use of these two methods is suggested to more accurately estimate Rt.

The modified slash Lindley-Weibull distribution with applications to nutrition data.

This work presents an extension of the slash Lindley-Weibull distribution, of which it can be considered a modification. The new family is obtained by using the quotient of two independent random variables: a two-parameter Lindley-Weibull distribution divided by a power of the exponential distribution with parameter equal to 2. We present the pdf and cdf of the new distribution, analyzing their risk functions. Some statistical properties are studied and the moments and coefficients of asymmetry and kurtosis are shown. The parameter estimation problem is carried out by the maximum likelihood method. The method is assessed by a Monte Carlo simulation study. We use nutrition data, which are characterized by high kurtosis, to illustrate the usefulness of the proposed model.