Univariate and bivariate distribution of growth traits in beef buffaloes from Brazil

The aim of this study was to analyze the weight at birth (BW) and adjusted at 205 (W205), 365 (W365) and 550 (W55O) days in beef buffaloes from Brazil, using two approaches: parametric, by normal distribution, and non-parametric, by kernel function, and thus estimating the genetic, environmental and phenotypic correlation among traits. Information of 5,169 animals at birth (BW), 3,792 at 205 days (W205), 3.883 at 365 days (W365) and 1,524 at 550 days of age (W550) were used. The birth weight distribution presented an evident discrepancy in relation to the normal distribution. However, W205, W365 and W550 presented normal distributions. The birth weight presented weak genetic, environmental, and phenotypic associations with the other weight measurements. on the other hand, the weight traits at 205, 365, 550 days of age showed a high genetic correlation.

Robust Bayesian analysis of heavy-tailed stochastic volatility models using scale mixtures of normal distributions

Fundação de Amparo à Pesquisa do Estado de São Paulo (FAPESP)

Power law distribution in education: University entrance examination

We studied the statistical distribution of candidate's performance which is measured through their marks in university entrance examination (Vestibular) of UNESP (Universidade Estadual Paulista) for years 1998, 1999, and 2000. All students are divided in three groups: Physical, Biological and Humanities. We paid special attention to the examination of Portuguese language which is common for all and examinations for the particular area. We observed long ubiquitous power law tails in Physical and Biological sciences. This indicate the presence of strong positive feedback in sciences. We are able to explain completely these statistical distributions through Gradually Truncated Power law distributions which we developed recently to explain statistical behavior of financial market. The statistical distribution in case of Portuguese language and humanities is close to normal distribution. We discuss the possible reason for this peculiar behavior.

A história da origem da curva normal

Pós-graduação em Educação Matemática - IGCE

Ecological distribution of the shrimp Pleoticus muelleri (Bate, 1888) and Artemesia longinaris Bate, 1888 (Decapoda, Penaeoidea) in the southeastern Brazilian littoral

Fundação de Amparo à Pesquisa do Estado de São Paulo (FAPESP)

Parameter recovery for a skew-normal IRT model under a Bayesian approach: hierarchical framework, prior and kernel sensitivity and sample size

Item response theory (IRT) comprises a set of statistical models which are useful in many fields, especially when there is an interest in studying latent variables (or latent traits). Usually such latent traits are assumed to be random variables and a convenient distribution is assigned to them. A very common choice for such a distribution has been the standard normal. Recently, Azevedo et al. [Bayesian inference for a skew-normal IRT model under the centred parameterization, Comput. Stat. Data Anal. 55 (2011), pp. 353-365] proposed a skew-normal distribution under the centred parameterization (SNCP) as had been studied in [R. B. Arellano-Valle and A. Azzalini, The centred parametrization for the multivariate skew-normal distribution, J. Multivariate Anal. 99(7) (2008), pp. 1362-1382], to model the latent trait distribution. This approach allows one to represent any asymmetric behaviour concerning the latent trait distribution. Also, they developed a Metropolis-Hastings within the Gibbs sampling (MHWGS) algorithm based on the density of the SNCP. They showed that the algorithm recovers all parameters properly. Their results indicated that, in the presence of asymmetry, the proposed model and the estimation algorithm perform better than the usual model and estimation methods. Our main goal in this paper is to propose another type of MHWGS algorithm based on a stochastic representation (hierarchical structure) of the SNCP studied in [N. Henze, A probabilistic representation of the skew-normal distribution, Scand. J. Statist. 13 (1986), pp. 271-275]. Our algorithm has only one Metropolis-Hastings step, in opposition to the algorithm developed by Azevedo et al., which has two such steps. This not only makes the implementation easier but also reduces the number of proposal densities to be used, which can be a problem in the implementation of MHWGS algorithms, as can be seen in [R.J. Patz and B.W. Junker, A straightforward approach to Markov Chain Monte Carlo methods for item response models, J. Educ. Behav. Stat. 24(2) (1999), pp. 146-178; R. J. Patz and B. W. Junker, The applications and extensions of MCMC in IRT: Multiple item types, missing data, and rated responses, J. Educ. Behav. Stat. 24(4) (1999), pp. 342-366; A. Gelman, G.O. Roberts, and W.R. Gilks, Efficient Metropolis jumping rules, Bayesian Stat. 5 (1996), pp. 599-607]. Moreover, we consider a modified beta prior (which generalizes the one considered in [3]) and a Jeffreys prior for the asymmetry parameter. Furthermore, we study the sensitivity of such priors as well as the use of different kernel densities for this parameter. Finally, we assess the impact of the number of examinees, number of items and the asymmetry level on the parameter recovery. Results of the simulation study indicated that our approach performed equally as well as that in [3], in terms of parameter recovery, mainly using the Jeffreys prior. Also, they indicated that the asymmetry level has the highest impact on parameter recovery, even though it is relatively small. A real data analysis is considered jointly with the development of model fitting assessment tools. The results are compared with the ones obtained by Azevedo et al. The results indicate that using the hierarchical approach allows us to implement MCMC algorithms more easily, it facilitates diagnosis of the convergence and also it can be very useful to fit more complex skew IRT models.

Distribution pattern of 68Ga-DOTATATE in disease-free patients

68Ga-DOTA-DPhe1,Tyr3-octreotate (68Ga-DOTATATE) is a somatostatin analogue that shows high affinity for somatostatin receptor subtype 2 (sst2) and has been used for imaging neuroendocrine tumours. However, normal uptake patterns and potential pitfalls have not been described with this high-sensitivity radiotracer. The aim of this study was therefore to outline the normal distribution pattern of 68Ga-DOTATATE in disease-free patients, to provide standardized uptake values (SUVs) of various organs and to compare our results with the current knowledge on sst2 receptor expression in vitro.

Model Evaluation Based on the Distribution of Estimated Absolute Prediction Error

The construction of a reliable, practically useful prediction rule for future response is heavily dependent on the "adequacy" of the fitted regression model. In this article, we consider the absolute prediction error, the expected value of the absolute difference between the future and predicted responses, as the model evaluation criterion. This prediction error is easier to interpret than the average squared error and is equivalent to the mis-classification error for the binary outcome. We show that the distributions of the apparent error and its cross-validation counterparts are approximately normal even under a misspecified fitted model. When the prediction rule is "unsmooth", the variance of the above normal distribution can be estimated well via a perturbation-resampling method. We also show how to approximate the distribution of the difference of the estimated prediction errors from two competing models. With two real examples, we demonstrate that the resulting interval estimates for prediction errors provide much more information about model adequacy than the point estimates alone.

Modeling absolute differences in life expectancy with a censored skew-normal regression approach.

Parameter estimates from commonly used multivariable parametric survival regression models do not directly quantify differences in years of life expectancy. Gaussian linear regression models give results in terms of absolute mean differences, but are not appropriate in modeling life expectancy, because in many situations time to death has a negative skewed distribution. A regression approach using a skew-normal distribution would be an alternative to parametric survival models in the modeling of life expectancy, because parameter estimates can be interpreted in terms of survival time differences while allowing for skewness of the distribution. In this paper we show how to use the skew-normal regression so that censored and left-truncated observations are accounted for. With this we model differences in life expectancy using data from the Swiss National Cohort Study and from official life expectancy estimates and compare the results with those derived from commonly used survival regression models. We conclude that a censored skew-normal survival regression approach for left-truncated observations can be used to model differences in life expectancy across covariates of interest.

Size frequency distribution of the β-amyloid (aβ) deposits in dementia with Lewy bodies with associated Alzheimer’s disease pathology

The objective is to study beta-amyloid (Abeta) deposition in dementia with Lewy bodies (DLB) with Alzheimer's disease (AD) pathology (DLB/AD). The size frequency distributions of the Abeta deposits were studied and fitted by log-normal and power-law models. Patients were ten clinically and pathologically diagnosed DLB/AD cases. Size distributions had a single peak and were positively skewed and similar to those described in AD and Down's syndrome. Size distributions had smaller means in DLB/AD than in AD. Log-normal and power-law models were fitted to the size distributions of the classic and diffuse deposits, respectively. Size distributions of Abeta deposits were similar in DLB/AD and AD. Size distributions of the diffuse deposits were fitted by a power-law model suggesting that aggregation/disaggregation of Abeta was the predominant factor, whereas the classic deposits were fitted by a log-normal distribution suggesting that surface diffusion was important in the pathogenesis of the classic deposits.

Factors determining the size frequency distribution of β-amyloid (Aβ) deposits in Alzheimer's disease

The size frequency distributions of discrete β-amyloid (Aβ) deposits were studied in single sections of the temporal lobe from patients with Alzheimer's disease. The size distributions were unimodal and positively skewed. In 18/25 (72%) tissues examined, a log normal distribution was a good fit to the data. This suggests that the abundances of deposit sizes are distributed randomly on a log scale about a mean value. Three hypotheses were proposed to account for the data: (1) sectioning in a single plane, (2) growth and disappearance of Aβ deposits, and (3) the origin of Aβ deposits from clusters of neuronal cell bodies. Size distributions obtained by serial reconstruction through the tissue were similar to those observed in single sections, which would not support the first hypothesis. The log normal distribution of Aβ deposit size suggests a model in which the rate of growth of a deposit is proportional to its volume. However, mean deposit size and the ratio of large to small deposits were not positively correlated with patient age or disease duration. The frequency distribution of Aβ deposits which were closely associated with 0, 1, 2, 3, or more neuronal cell bodies deviated significantly from a log normal distribution, which would not support the neuronal origin hypothesis. On the basis of the present data, growth and resolution of Aβ deposits would appear to be the most likely explanation for the log normal size distributions.

A quantitative study of abnormally enlarged neurons in cognitively normal brain and neurodegenerative disease

Abnormally enlarged neurons (AEN) occur in many neurodegenerative diseases. To define AEN more objectively, the frequency distribution of the ratio of greatest cell diameter(CD) to greatest nuclear diameter (ND) was studied in populations of cortical neurons in tissue sections of seven cognitively normal brains. The frequency distribution of CD/ND deviated from a normal distribution in 15 out of 18 populations of neurons studied and hence, the 95th percentile (95P) was used to define a limit of the CD/ND ratio excluding the5% most extreme observations. The 95P of the CD/ ND ratio varied from 2.0 to 3.0 in different cases and regions and a value of 95P = 3.0 was chosen to define the limit for normalneurons under non-pathological conditions. Based on the 95P = 3.0 criterion, the proportion of AEN with a CD/ND ≥ 3 varied from 2.6% in Alzheimer's disease (AD) to 20.3% in Pick's disease (PiD). The data suggest: (1) that a CL/ND ≥ 3.0 may be a useful morphological criterion for defining AEN, and (2) AEN were most numerous in PiD and corticobasal degeneration (CBD) and least abundant in AD and in dementia with Lewy bodies (DLB). © 2013 Dustri-Verlag Dr. K. Feistle.

Statnote 35: are the data log-normal?

In many of the Statnotes described in this series, the statistical tests assume the data are a random sample from a normal distribution These Statnotes include most of the familiar statistical tests such as the ‘t’ test, analysis of variance (ANOVA), and Pearson’s correlation coefficient (‘r’). Nevertheless, many variables exhibit a more or less ‘skewed’ distribution. A skewed distribution is asymmetrical and the mean is displaced either to the left (positive skew) or to the right (negative skew). If the mean of the distribution is low, the degree of variation large, and when values can only be positive, a positively skewed distribution is usually the result. Many distributions have potentially a low mean and high variance including that of the abundance of bacterial species on plants, the latent period of an infectious disease, and the sensitivity of certain fungi to fungicides. These positively skewed distributions are often fitted successfully by a variant of the normal distribution called the log-normal distribution. This statnote describes fitting the log-normal distribution with reference to two scenarios: (1) the frequency distribution of bacterial numbers isolated from cloths in a domestic environment and (2), the sizes of lichenised ‘areolae’ growing on the hypothalus of Rhizocarpon geographicum (L.) DC.

Statnote 36: Do the data fit the Poisson distribution

In previous Statnotes, many of the statistical tests described rely on the assumption that the data are a random sample from a normal or Gaussian distribution. These include most of the tests in common usage such as the ‘t’ test ), the various types of analysis of variance (ANOVA), and Pearson’s correlation coefficient (‘r’) . In microbiology research, however, not all variables can be assumed to follow a normal distribution. Yeast populations, for example, are a notable feature of freshwater habitats, representatives of over 100 genera having been recorded . Most common are the ‘red yeasts’ such as Rhodotorula, Rhodosporidium, and Sporobolomyces and ‘black yeasts’ such as Aurobasidium pelculans, together with species of Candida. Despite the abundance of genera and species, the overall density of an individual species in freshwater is likely to be low and hence, samples taken from such a population will contain very low numbers of cells. A rare organism living in an aquatic environment may be distributed more or less at random in a volume of water and therefore, samples taken from such an environment may result in counts which are more likely to be distributed according to the Poisson than the normal distribution. The Poisson distribution was named after the French mathematician Siméon Poisson (1781-1840) and has many applications in biology, especially in describing rare or randomly distributed events, e.g., the number of mutations in a given sequence of DNA after exposure to a fixed amount of radiation or the number of cells infected by a virus given a fixed level of exposure. This Statnote describes how to fit the Poisson distribution to counts of yeast cells in samples taken from a freshwater lake.

Income distribution dependence of poverty measure:a theoretical analysis

Using a modified deprivation (or poverty) function, in this paper, we theoretically study the changes in poverty with respect to the 'global' mean and variance of the income distribution using Indian survey data. We show that when the income obeys a log-normal distribution, a rising mean income generally indicates a reduction in poverty while an increase in the variance of the income distribution increases poverty. This altruistic view for a developing economy, however, is not tenable anymore once the poverty index is found to follow a pareto distribution. Here although a rising mean income indicates a reduction in poverty, due to the presence of an inflexion point in the poverty function, there is a critical value of the variance below which poverty decreases with increasing variance while beyond this value, poverty undergoes a steep increase followed by a decrease with respect to higher variance. Identifying this inflexion point as the poverty line, we show that the pareto poverty function satisfies all three standard axioms of a poverty index [N.C. Kakwani, Econometrica 43 (1980) 437; A.K. Sen, Econometrica 44 (1976) 219] whereas the log-normal distribution falls short of this requisite. Following these results, we make quantitative predictions to correlate a developing with a developed economy. © 2006 Elsevier B.V. All rights reserved.