Estimation of Heritability in Heteroscedastic One‐Way Unbalanced Random Model

A modified estimator of heritability is proposed under heteroscedastic one way unbalanced random model. The distribution, moments and probability of permissible values (PPV) for conventional and modified estimators are derived. The behaviour of two estimators has been investigated, numerically, to devise a suitable estimator of heritability under variance heterogeneity. The numerical results reveal that under balanced case the heteroscedasticity affects the bias, MSE and PPV of conventional estimator, marginally. In case of unbalanced situations, the conventional estimator underestimates the parameter when more variable group has more observations and overestimates when more variable group has less observations, MSE of the conventional estimator decreases when more variable group has more observations and increases when more variable group has less observations and PPV is marginally decreased. The MSE and PPV are comparable for two estimators while the bias of modified estimator is less than the conventional estimator particularly for small and medium values of the parameter. These results suggest the use of modified estimator with equal or more observations for more variable group in presence of variance heterogeneity.     

The Utilization of Kurtosis in the Estimation of the Parameters of the One‐Way Random Effects Model

This paper deals with the problem of estimating the components of the variance of one‐way random effects model under non‐normality situation using a prior knowledge of coefficient of kurtosis. We have suggested two classes of estimators and for the within and between variances respectively. Optimum estimators in the classes of and are identified with their mean squared errors formulae and compared with that of usual ANOVA unbiased and Shoukri , Tracy and Mian 's (1990) estimators. It is found that the proposed estimators are more efficient than the ANOVA unbiased estimators and Shoukri , Tracy and Mian (1990) estimators.     

The Effect of Design Imbalance on the Power of the F‐test of a Variance Component in the One‐Way Random Model

The ANOVA‐based F‐test used for testing the significance of the random effect variance component is a valid test for an unbalanced one‐way random model. However, it does not have an uniform optimum property. For example, this test is not uniformly most powerful invariant (UMPI). In fact, there is no UMPI test in the unbalanced case (see Khuri , Mathew , and Sinha , 1998). The power of the F‐test depends not only on the design used, but also on the true values of the variance components. As Khuri (1996) noted, we can gain a better insight into the effect of data imbalance on the power of the F‐test using a method for modelling the power in terms of the design parameters and the variance components. In this study, generalized linear modelling (GLM) techniques are used for this purpose. It is shown that GLM, in combination with a method of generating designs with a specified degree of imbalance, is an effective way of studying the behavior of the power of the F‐test in a one‐way random model.     

A Maximum Likelihood Estimator for Shape Parameters of S‐Distributions

The S‐distribution is a four‐parameter distribution that is defined in terms of a differential equation, in which the cumulative is represented as the dependent variable: The article proposes a maximum likelihood estimator for the shape parameters of this distribution.     

Nonparametric Tests in the Unbalanced Multivariate One‐Way Design

A nonparametric model for the multivariate one‐way design is discussed which entails continuous as well as discontinuous distributions and, therefore, allows for ordinal data. Nonparametric hypotheses are formulated by the normalized version of the marginal distribution functions as well as the common distribution functions. The differences between the distribution functions are described by means of the so‐called relative treatment effects, for which unbiased and consistent estimators are derived. The asymptotic distribution of the vector of the effect estimators is derived and under the marignal hypothesis a consistent estimator for the asymptotic covariance matrix is given. Nonparametric versions of the Wald‐type statistic, the ANOVA‐type statistic and the Lawley‐Hotelling statistic are considered and compared by means of a simulation study. Finally, these tests are applied to a psychiatric clinical trial.     

Analysis of the Binary Littermate Data in the One‐Way Layout

In certain experiments with laboratory animals like rodents, it is common to observe the binary littermate data that the fetuses in the same litter provide a set of independent Bernoulli responses conditional on the litter specific abnormal probability. For testing homogeneity of abnormal proportions of several groups, five methods are discussed in the frame of the one‐way layout with litter effects. Simulation results show a superiority of the Rao‐Scott procedure.     

A Method of Moments Estimator for Random Effect Multivariate Meta‐Analysis

Summary Meta‐analysis is a powerful approach to combine evidence from multiple studies to make inference about one or more parameters of interest, such as regression coefficients. The validity of the fixed effect model meta‐analysis depends on the underlying assumption that all studies in the meta‐analysis share the same effect size. In the presence of heterogeneity, the fixed effect model incorrectly ignores the between‐study variance and may yield false positive results. The random effect model takes into account both within‐study and between‐study variances. It is more conservative than the fixed effect model and should be favored in the presence of heterogeneity. In this paper, we develop a noniterative method of moments estimator for the between‐study covariance matrix in the random effect model multivariate meta‐analysis. To our knowledge, it is the first such method of moments estimator in the matrix form. We show that our estimator is a multivariate extension of DerSimonian and Laird’s univariate method of moments estimator, and it is invariant to linear transformations. In the simulation study, our method performs well when compared to existing random effect model multivariate meta‐analysis approaches. We also apply our method in the analysis of a real data example.     

On the Probability of a Negative Estimator of Heteroscedasticity in the Random Coefficient Model

In the estimation of a linear regression model with random coefficients, sometimes negative estimates of variances of random coefficients are obtained–an undesirable feature. In this paper we have obtained the asymptotic bounds of the probability of obtaining negative estimators. A simple illustration is also provided for the purpose.     

A Robust Alternative to the Schemper–Henderson Estimator of Prediction Error

Summary In clinical applications, the prediction error of survival models has to be taken into consideration to assess the practical suitability of conclusions drawn from these models. Different approaches to evaluate the predictive performance of survival models have been suggested in the literature. In this article, we analyze the properties of the estimator of prediction error developed by Schemper and Henderson (2000 , Biometrics 56, 249–255), which quantifies the absolute distance between predicted and observed survival functions. We provide a formal proof that the estimator proposed by Schemper and Henderson is not robust against misspecification of the survival model, that is, the estimator will only be meaningful if the model family used for deriving predictions has been specified correctly. To remedy this problem, we construct a new estimator of the absolute distance between predicted and observed survival functions. We show that this modified Schemper–Henderson estimator is robust against model misspecification, allowing its practical application to a wide class of survival models. The properties of the Schemper–Henderson estimator and its new modification are illustrated by means of a simulation study and the analysis of two clinical data sets.     

Analytical Quantile Solution for the S‐distribution,Random Number Generation and Statistical Data Modeling

The selection of a specific statistical distribution as a model for describing the population behavior of a given variable is seldom a simple problem. One strategy consists in testing different distributions (normal, lognormal, Weibull, etc.), and selecting the one providing the best fit to the observed data and being the most parsimonious. Alternatively, one can make a choice based on theoretical arguments and simply fit the corresponding parameters to the observed data. In either case, different distributions can give similar results and provide almost equivalent models for a given data set. Model selection can be more complicated when the goal is to describe a trend in the distribution of a given variable. In those cases, changes in shape and skewness are difficult to represent by a single distributional form. As an alternative to the use of complicated families of distributions as models for data, the S‐distribution [Voit, E. O. (1992) Biom. J. 7 , 855–878] provides a highly flexible mathematical form in which the density is defined as a function of the cumulative. S‐distributions can accurately approximate many known continuous and unimodal distributions, preserving the well known limit relationships between them. Besides representing well‐known distributions, S‐distributions provide an infinity of new possibilities that do not correspond with known classical distributions. Although the utility and performance of this general form has been clearly proved in different applications, its definition as a differential equation is a potential drawback for some problems. In this paper we obtain an analytical solution for the quantile equation that highly simplifies the use of S‐distributions. We show the utility of this solution in different applications. After classifying the different qualitative behaviors of the S‐distribution in parameter space, we show how to obtain different S‐distributions that accomplish specific constraints. One of the most interesting cases is the possibility of obtaining distributions that acomplish P(X ≤ X_c) = 0. Then, we demonstrate that the quantile solution facilitates the use of S‐distributions in Monte‐Carlo experiments through the generation of random samples. Finally, we show how to fit an S‐distribution to actual data, so that the resulting distribution can be used as a statistical model for them.     

A Bagging‐Based Correction for the Mixture Model Estimator of Population Size

Estimation of a population size by means of capture‐recapture techniques is an important problem occurring in many areas of life and social sciences. We consider the frequencies of frequencies situation, where a count variable is used to summarize how often a unit has been identified in the target population of interest. The distribution of this count variable is zero‐truncated since zero identifications do not occur in the sample. As an application we consider the surveillance of scrapie in Great Britain. In this case study holdings with scrapie that are not identified (zero counts) do not enter the surveillance database. The count variable of interest is the number of scrapie cases per holding. For count distributions a common model is the Poisson distribution and, to adjust for potential heterogeneity, a discrete mixture of Poisson distributions is used. Mixtures of Poissons usually provide an excellent fit as will be demonstrated in the application of interest. However, as it has been recently demonstrated, mixtures also suffer under the so‐called boundary problem, resulting in overestimation of population size. It is suggested here to select the mixture model on the basis of the Bayesian Information Criterion. This strategy is further refined by employing a bagging procedure leading to a series of estimates of population size. Using the median of this series, highly influential size estimates are avoided. In limited simulation studies it is shown that the procedure leads to estimates with remarkable small bias. (© 2008 WILEY‐VCH Verlag GmbH & Co. KGaA, Weinheim)     

Nonparametric One‐way Analysis of Variance of Replicated Bivariate Spatial Point Patterns

A common problem in neuropathological studies is to assess the spatial patterning of cells on tissue sections and to compare spatial patterning between disorder groups. For a single cell type, the cell positions constitute a univariate point process and interest focuses on the degree of spatial aggregation. For two different cell types, the cell positions constitute a bivariate point process and the degree of spatial interaction between the cell types is of interest. We discuss the problem of analysing univariate and bivariate spatial point patterns in the one‐way design where cell patterns have been obtained for groups of subjects. A bootstrapping procedure to perform a nonparametric one‐way analysis of variance of the spatial aggregation of a univariate point process has been suggested by Diggle, Lange and Bene? (1991). We extend their replication‐based approach to allow the comparison of the spatial interaction of two cell types between groups, to include planned comparisons (contrasts) and to assess whole groups against complete spatial randomness and spatial independence. We also accommodate several replicate tissue sections per subject. An advantage of our approach is that it can be applied when processes are not stationary, a common problem in brain tissue sections since neurons are arranged in cortical layers. We illustrate our methods by applying them to a neuropathological study to investigate abnormalities in the functional relationship between neurons and astrocytes in HIV associated dementia. (© 2004 WILEY‐VCH Verlag GmbH & Co. KGaA, Weinheim)     

On the Effect of Unbalancedness and Heteroscedasticity on the ANOVA Estimator of Group Variance Component in One-Way Random Model

The expression for rth cumulant of ANOVA estimator of group variance component is derived in the One-way unbalanced random model under heteroscedasticity. The expression is used to study the effect of unbalancedness and heteroscedasticity on the mean and variance of the estimator, numerically. The computed results reveal that the unbalancedness and heteroscedasticity have a combined effect on the mean and variance of the estimator. For certain situations of unequal group sizes and error variances, the mean and variance of the estimator are increased and for certain other situations the values are decreased.     

Crystal structure of an S‐formylglutathione hydrolase from Pseudoalteromonas haloplanktis TAC125

S‐formylglutathione hydrolases (FGHs) constitute a family of ubiquitous enzymes which play a key role in formaldehyde detoxification both in prokaryotes and eukaryotes, catalyzing the hydrolysis of S‐formylglutathione to formic acid and glutathione. While a large number of functional studies have been reported on these enzymes, few structural studies have so far been carried out. In this article we report on the functional and structural characterization of PhEst, a FGH isolated from the psychrophilic bacterium Pseudoalteromonas haloplanktis. According to our functional studies, this enzyme is able to efficiently hydrolyze several thioester substrates with very small acyl moieties. By contrast, the enzyme shows no activity toward substrates with bulky acyl groups. These data are in line with structural studies which highlight for this enzyme a very narrow acyl‐binding pocket in a typical α/β‐hydrolase fold. PhEst represents the first cold‐adapted FGH structurally characterized to date; comparison with its mesophilic counterparts of known three‐dimensional structure allowed to obtain useful insights into molecular determinants responsible for the ability of this psychrophilic enzyme to work at low temperature. © 2010 Wiley Periodicals, Inc. Biopolymers 93: 669–677, 2010. This article was originally published online as an accepted preprint. The “Published Online” date corresponds to the preprint version. You can request a copy of the preprint by emailing the Biopolymers editorial office at biopolymers@wiley.com     

The role of S‐acylation in protein trafficking

Protein S‐acylation, also known as palmitoylation, consists of the addition of a lipid molecule to one or more cysteine residues through a thioester bond. This modification, which is widespread in eukaryotes, is thought to affect over 12% of the human proteome. S‐acylation allows the reversible association of peripheral proteins with membranes or, in the case of integral membrane proteins, modulates their behavior within the plane of the membrane. This review focuses on the consequences of protein S‐acylation on intracellular trafficking and membrane association. We summarize relevant information that illustrates how lipid modification of proteins plays an important role in dictating precise intracellular movements within cells by regulating membrane‐cytosol exchange, through membrane microdomain segregation, or by modifying the flux of the proteins by means of vesicular or diffusional transport systems. Finally, we highlight some of the key open questions and major challenges in the field.      

Identification and characterization of the zebrafish glutathione S‐transferase Pi‐1

Zebrafish has in recent years emerged as a popular vertebrate model for use in pharmacological and toxicological studies. While there have been sporadic studies on the zebrafish glutathione S‐transferases (GSTs), the zebrafish GST gene superfamily still awaits to be fully elucidated. We report here the identification of 15 zebrafish cytosolic GST genes in NCBI GenBank database and the expression, purification, and enzymatic characterization of the zebrafish cytosolic GST Pi‐1 (GSTP1). The cDNA encoding the zebrafish GSTP1 was cloned from a 3‐month‐old female zebrafish, expressed in Eschelichia coli host cells, and purified. Purified GSTP1 displayed glutathione‐conjugating activity toward 1‐chloro‐2,4‐dinitrobenzene as a representative substrate. The enzymatic characteristics of the zebrafish GSTP1, including pH‐dependency, effects of metal cations, and kinetic parameters, were studied. Moreover, the expression of zebrafish GSTP1 at different developmental stages during embryogenesis, throughout larval development, onto maturity was examined.