Quantile regression for longitudinal data with a working correlation model

This paper proposes a linear quantile regression analysis method for longitudinal data that combines the between- and within-subject estimating functions, which incorporates the correlations between repeated measurements. Therefore, the proposed method results in more efficient parameter estimation relative to the estimating functions based on an independence working model. To reduce computational burdens, the induced smoothing method is introduced to obtain parameter estimates and their variances. Under some regularity conditions, the estimators derived by the induced smoothing method are consistent and have asymptotically normal distributions. A number of simulation studies are carried out to evaluate the performance of the proposed method. The results indicate that the efficiency gain for the proposed method is substantial especially when strong within correlations exist. Finally, a dataset from the audiology growth research is used to illustrate the proposed methodology.

Study of Pinna nobilis growth from inner record: How biased are posterior adductor muscle scars estimates?

Previous studies have shown that the external growth records of the posterior adductor muscle scar (PAMS) of the bivalve Pinna nobilis are incomplete and do not produce accurate age estimations. We have developed a new methodology to study age and growth using the inner record of the PAMS, which avoids the necessity of costly in situ shell measurements or isotopic studies. Using the inner record we identified the positions of PAMS previously obscured by nacre and estimated the number of missing records in adult specimens with strong abrasion of the calcite layer in the anterior portion of the shell. The study of the PAMS and inner record of two shells that were 6 years old when collected showed that only 2 and 3 PAMS were observed, while 6 inner records could be counted, thus confirming our working methodology. Growth parameters of a P. nobilis population located in Moraira, Spain (western Mediterranean) were estimated with the new methodology and compared to those obtained using PAMS data and in situ measurements. For the comparisons, we applied different models considering the data alternatively as length-at-age (LA) and tag-recapture (TR). Among every method we tested to fit the Von Bertalanffy growth model, we observed that LA data from inner record fitted to the model using non-linear mixed effects and the estimation of missing records using the calcite width was the most appropriate. The equation obtained with this method, L = 573*(1 - e(-0.16(t-0.02))), is very similar to that calculated previously from in situ measurements for the same population.

Fatigue crack growth rates and fracture toughness in electroslag refined En 52 steel

En 52 steel has been electroslag refined and the resultant effects of refining on its mechanical properties have been assessed. It was found that refining caused a decrease in fatigue crack growth rates and increases in fatigue strength, fracture toughness, Charpy fracture energy and tensile ductility. Fatigue crack growth rates in region I and in region III were found to be considerably lower in the electroslag refined steel: they were unaffected in region II. The fracture toughness values for the electroslag refined steel are nearly twice those estimated for the unrefined steel. Measurements on heat-treated samples have shown that the electroslag refined steel has a better response to heat-treatment. The improvement in the mechanical properties is explained in terms of the removal of nonmetallic inclusions and a reduction in the sulphur content of the steel.

Deep inelastic electroproduction structure functions for polarized and unpolarized proton targets

Following Ioffe's method of QCD sum rules the structure functions F2(x) for deep inelastic ep and en scattering are calculated. Valence u-quark and d-quark distributions are obtained in the range 0.1 less, approximate x <0.4 and compared with data. In the case of polarized targets the structure function g1(x) and the asymmetry Image Full-size image are calculated. The latter is in satisfactory agreement in sign and magnitude with experiments for x in the range 0.1< x < 0.4.

A maximum likelihood approach for estimating growth from tag–recapture data

The Fabens method is commonly used to estimate growth parameters k and l infinity in the von Bertalanffy model from tag-recapture data. However, the Fabens method of estimation has an inherent bias when individual growth is variable. This paper presents an asymptotically unbiassed method using a maximum likelihood approach that takes account of individual variability in both maximum length and age-at-tagging. It is assumed that each individual's growth follows a von Bertalanffy curve with its own maximum length and age-at-tagging. The parameter k is assumed to be a constant to ensure that the mean growth follows a von Bertalanffy curve and to avoid overparameterization. Our method also makes more efficient use nf thp measurements at tno and recapture and includes diagnostic techniques for checking distributional assumptions. The method is reasonably robust and performs better than the Fabens method when individual growth differs from the von Bertalanffy relationship. When measurement error is negligible, the estimation involves maximizing the profile likelihood of one parameter only. The method is applied to tag-recapture data for the grooved tiger prawn (Penaeus semisulcatus) from the Gulf of Carpentaria, Australia.

Accounting for individual variability in the von Bertalanffy growth model

Estimation of von Bertalanffy growth parameters has received considerable attention in fisheries research. Since Sainsbury (1980, Can. J. Fish. Aquat. Sci. 37: 241-247) much of this research effort has centered on accounting for individual variability in the growth parameters. In this paper we demonstrate that, in analysis of tagging data, Sainsbury's method and its derivatives do not, in general, satisfactorily account for individual variability in growth, leading to inconsistent parameter estimates (the bias does not tend to zero as sample size increases to infinity). The bias arises because these methods do not use appropriate conditional expectations as a basis for estimation. This bias is found to be similar to that of the Fabens method. Such methods would be appropriate only under the assumption that the individual growth parameters that generate the growth increment were independent of the growth parameters that generated the initial length. However, such an assumption would be unrealistic. The results are derived analytically, and illustrated with a simulation study. Until techniques that take full account of the appropriate conditioning have been developed, the effect of individual variability on growth has yet to be fully understood.

A simple method for estimating growth parameters from multiple length-frequency data in presence of continuous recruitment

The extended recruitment season for short-lived species such as prawns biases the estimation of growth parameters from length-frequency data when conventional methods are used. We propose a simple method for overcoming this bias given a time series of length-frequency data. The difficulties arising from extended recruitment are eliminated by predicting the growth of the succeeding samples and the length increments of the recruits in previous samples. This method requires that some maximum size at recruitment can be specified. The advantages of this multiple length-frequency method are: it is simple to use; it requires only three parameters; no specific distributions need to be assumed; and the actual seasonal recruitment pattern does not have to be specified. We illustrate the new method with length-frequency data on the tiger prawn Penaeus esculentus from the north-western Gulf of Carpentaria, Australia.

Robust estimating functions and bias correction for longitudinal data analysis

Robust methods are useful in making reliable statistical inferences when there are small deviations from the model assumptions. The widely used method of the generalized estimating equations can be "robustified" by replacing the standardized residuals with the M-residuals. If the Pearson residuals are assumed to be unbiased from zero, parameter estimators from the robust approach are asymptotically biased when error distributions are not symmetric. We propose a distribution-free method for correcting this bias. Our extensive numerical studies show that the proposed method can reduce the bias substantially. Examples are given for illustration.

Effects of variance-function misspecification in analysis of longitudinal data

The approach of generalized estimating equations (GEE) is based on the framework of generalized linear models but allows for specification of a working matrix for modeling within-subject correlations. The variance is often assumed to be a known function of the mean. This article investigates the impacts of misspecifying the variance function on estimators of the mean parameters for quantitative responses. Our numerical studies indicate that (1) correct specification of the variance function can improve the estimation efficiency even if the correlation structure is misspecified; (2) misspecification of the variance function impacts much more on estimators for within-cluster covariates than for cluster-level covariates; and (3) if the variance function is misspecified, correct choice of the correlation structure may not necessarily improve estimation efficiency. We illustrate impacts of different variance functions using a real data set from cow growth.

Growth curves with time-dependent explanatory variables

A class of growth models incorporating time-dependent factors and stochastic perturbations are introduced. The proposed model includes the existing growth models used in fisheries as special cases. Particular attention is given to growth of a population (in average weight or length) from which observations are taken randomly each time and the analysis of tag-recapture data. Two real data sets are used for illustration: (a) to estimate the seasonal effect and population density effect on growth of farmed prawn (Penaeus monodon) from weight data and (b) to assess the effect of tagging on growth of barramundi (Lates calcarifer)

Maximum likelihood estimation of mortality and growth with individual variability from multiple length-frequency data

We consider estimation of mortality rates and growth parameters from length-frequency data of a fish stock and derive the underlying length distribution of the population and the catch when there is individual variability in the von Bertalanffy growth parameter L-infinity. The model is flexible enough to accommodate 1) any recruitment pattern as a function of both time and length, 2) length-specific selectivity, and 3) varying fishing effort over time. The maximum likelihood method gives consistent estimates, provided the underlying distribution for individual variation in growth is correctly specified. Simulation results indicate that our method is reasonably robust to violations in the assumptions. The method is applied to tiger prawn data (Penaeus semisulcatus) to obtain estimates of natural and fishing mortality.

Population structure, mortality and growth of Pinna nobilis Linnaeus, 1758 (Mollusca, Bivalvia) at different depths in Moraira bay (Alicante, Western Mediterranean)

An investigation to characterize the causes of Pinna nobilis population structure in Moraira bay (Western Mediterranean) was developed. Individuals of two areas of the same Posidonia meadow, located at different depths (A1, -13 and A2, -6 m), were inventoried, tagged, their positions accurately recorded and monitored from July 1997 to July 2002. On each area, different aspects of population demography were studied (i.e. spatial distribution, size structure, displacement evidences, mortality, growth and shell orientation). A comparison between both groups of individuals was carried out, finding important differences between them. In A1, the individuals were more aggregated and mean and maximum size were higher (A1, 10.3 and A2, 6 individuals/100 m(2); A1, x = 47.2 +/- 9.9; A2, x = 29.8 +/- 7.4 cm, P < 0.001, respectively). In A2, growth rate and mortality were higher, the latter concentrated on the largest individuals, in contrast to A1, where the smallest individuals had the higher mortality rate [A1, L = 56.03(1 - e(-0.17t)); A2, L = 37.59(1 - e(-0.40t)), P < 0.001; mean annual mortality A1: 32 dead individuals out of 135, 23.7% and A2: 16 dead individuals out of 36, 44.4%, and total mortality coefficients (z), z(A1(-30)) = 0.28, z(A1(31-45)) = 0.05, z(A1(46-)) = 0.08; z(A2(-30)) = 0.15, z(A2(31-45)) = 0.25]. A common shell orientation N-S, coincident with the maximum shore exposure, was observed in A2. Spatial distribution in both areas showed not enough evidence to discard a random distribution of the individuals, despite the greater aggregation on the deeper area (A1) (A1, chi(2) = 0.41, df = 3, P > 0.5, A2, chi(2)= 0.98, df = 2 and 0.3 < P < 0.5). The obtained results have demonstrated that the depth-related size segregation usually shown by P. nobilis is mainly caused by differences in mortality and growth among individuals located at different depths, rather than by the active displacement of individuals previously reported in the literature. Furthermore, dwarf individuals are observed in shallower levels and as a consequence, the relationship between size and age are not comparable even among groups of individuals inhabiting the same meadow at different depths. The final causes of the differences on mortality and growth are also discussed.

Kallikrein-related peptidases in prostate cancer: From molecular function to clinical application

Prostate cancer is a leading contributor to male cancer-related deaths worldwide. Kallikrein-related peptidases (KLKs) are serine proteases that exhibit deregulated expression in prostate cancer, with KLK3, or prostate specific antigen (PSA), being the widely-employed clinical biomarker for prostate cancer. Other KLKs, such as KLK2, show promise as prostate cancer biomarkers and, additionally, their altered expression has been utilised for the design of KLK-targeted therapies. There is also a large body of in vitro and in vivo evidence supporting their role in cancer-related processes. Here, we review the literature on studies to date investigating the potential of other KLKs, in addition to PSA, as biomarkers and in therapeutic options, as well as their current known functional roles in cancer progression. Increased knowledge of these KLK-mediated functions, including degradation of the extracellular matrix, local invasion, cancer cell proliferation, interactions with fibroblasts, angiogenesis, migration, bone metastasis and tumour growth in vivo, may help define new roles as prognostic biomarkers and novel therapeutic targets for this cancer.

Evaluation of growth and changes in body composition following neonatal diagnosis of cystic fibrosis

Early deficits in nutritional status that might require specific treatment and early response to nutritional therapy were studied longitudinally in 25 infants with cystic fibrosis (CF) diagnosed by neonatal screening, using anthropometric and research body composition methodology, and evaluation of pancreatic function. At the time of confirmed diagnosis (mean 5.4 weeks), body mass, length, total body fat (TBF), and total body potassium (TBK) were all significantly reduced. Following diagnosis and commencement of therapy there was a normalization of weight, length, and TBK by 6-12 months of age, indicating catch-up growth. But in some individuals the response was incomplete, and as a group, mean total body fat remained significantly lower than normal at 1 year of age. Seven of 25 (28%) were pancreatic sufficient at diagnosis, and all but one had evidence of declining pancreatic function requiring the institution of pancreatic enzyme therapy during the next 1-9 months. The median age of commencement of enzyme therapy was 10 weeks (range 5 weeks to 11 months). These longitudinal assessments emphasize the dynamic changes occurring in absorptive function, body composition, and nutritional status following neonatal diagnosis of cystic fibrosis and may reflect previously described abnormalities of energy metabolism in this age group. Abnormal body composition is evident in most CF infants following diagnosis by neonatal screening but pancreatic damage may still be evolving. We suggest that early active nutritional therapy and surveillance for changes in pancreatic function are warranted in CF infants diagnosed by neonatal screening.