Modeling strategies in longitudinal data analysis: Covariate, variance function and correlation structure selection

A modeling paradigm is proposed for covariate, variance and working correlation structure selection for longitudinal data analysis. Appropriate selection of covariates is pertinent to correct variance modeling and selecting the appropriate covariates and variance function is vital to correlation structure selection. This leads to a stepwise model selection procedure that deploys a combination of different model selection criteria. Although these criteria find a common theoretical root based on approximating the Kullback-Leibler distance, they are designed to address different aspects of model selection and have different merits and limitations. For example, the extended quasi-likelihood information criterion (EQIC) with a covariance penalty performs well for covariate selection even when the working variance function is misspecified, but EQIC contains little information on correlation structures. The proposed model selection strategies are outlined and a Monte Carlo assessment of their finite sample properties is reported. Two longitudinal studies are used for illustration.

Discussion of “Generalized estimating equations: Notes on the choice of the working correlation matrix”

Objective To discuss generalized estimating equations as an extension of generalized linear models by commenting on the paper of Ziegler and Vens "Generalized Estimating Equations. Notes on the Choice of the Working Correlation Matrix". Methods Inviting an international group of experts to comment on this paper. Results Several perspectives have been taken by the discussants. Econometricians have established parallels to the generalized method of moments (GMM). Statisticians discussed model assumptions and the aspect of missing data Applied statisticians; commented on practical aspects in data analysis. Conclusions In general, careful modeling correlation is encouraged when considering estimation efficiency and other implications, and a comparison of choosing instruments in GMM and generalized estimating equations, (GEE) would be worthwhile. Some theoretical drawbacks of GEE need to be further addressed and require careful analysis of data This particularly applies to the situation when data are missing at random.

Working-correlation-structure identification in generalized estimating equations

Selecting an appropriate working correlation structure is pertinent to clustered data analysis using generalized estimating equations (GEE) because an inappropriate choice will lead to inefficient parameter estimation. We investigate the well-known criterion of QIC for selecting a working correlation Structure. and have found that performance of the QIC is deteriorated by a term that is theoretically independent of the correlation structures but has to be estimated with an error. This leads LIS to propose a correlation information criterion (CIC) that substantially improves the QIC performance. Extensive simulation studies indicate that the CIC has remarkable improvement in selecting the correct correlation structures. We also illustrate our findings using a data set from the Madras Longitudinal Schizophrenia Study.

Electronic-structures of h2o.bf3 and related n-v addition-compounds - a combined eels-ups study in vapor-phase

The detailed electronic structure of the n-v addition compound H2O·BF3 has been investigated for the first time by a combined use of electron energy loss spectroscopy (EELS) and UV photoelectron spectroscopy (UPS) augmented by MO calculations. The calculated molecular orbital energies of H2O·BF3 agree well with the UPS results and have been used to assign the electronic transitions obtained from EELS and to construct an orbital correlation diagram. The Journal of Chemical Physics is copyrighted by The American Institute of Physics.

Criteria for working-correlation-structure selection in GEE: Assessment via simulation

Efficiency of analysis using generalized estimation equations is enhanced when intracluster correlation structure is accurately modeled. We compare two existing criteria (a quasi-likelihood information criterion, and the Rotnitzky-Jewell criterion) to identify the true correlation structure via simulations with Gaussian or binomial response, covariates varying at cluster or observation level, and exchangeable or AR(l) intracluster correlation structure. Rotnitzky and Jewell's approach performs better when the true intracluster correlation structure is exchangeable, while the quasi-likelihood criteria performs better for an AR(l) structure.

Unbiased estimating equations from working correlation models for irregularly timed repeated measures

The method of generalized estimating equation-, (GEEs) has been criticized recently for a failure to protect against misspecification of working correlation models, which in some cases leads to loss of efficiency or infeasibility of solutions. However, the feasibility and efficiency of GEE methods can be enhanced considerably by using flexible families of working correlation models. We propose two ways of constructing unbiased estimating equations from general correlation models for irregularly timed repeated measures to supplement and enhance GEE. The supplementary estimating equations are obtained by differentiation of the Cholesky decomposition of the working correlation, or as score equations for decoupled Gaussian pseudolikelihood. The estimating equations are solved with computational effort equivalent to that required for a first-order GEE. Full details and analytic expressions are developed for a generalized Markovian model that was evaluated through simulation. Large-sample ".sandwich" standard errors for working correlation parameter estimates are derived and shown to have good performance. The proposed estimating functions are further illustrated in an analysis of repeated measures of pulmonary function in children.

Exocrine pancreatic dysfunction in malnourished Australian Aboriginal children

Pancreatic exocrine dysfunction has been frequently recorded in protein-energy malnutrition in underdeveloped countries. In addition, the pancreas requires optimal nutrition for enzyme synthesis and potentially correctable pancreatic enzyme insufficiency may play a role in the continuation of protein-energy malnutrition. This problem has not been previously evaluated in Australian Aborigines. We have applied a screening test for pancreatic dysfunction (human immunoreactive trypsinogen [IRT] assay) to the study of 398 infants (6-36 months) admitted to the Alice Springs Hospital over a 20-month period. All infants were assessed by anthropometric measures and were assigned to to three nutritional groups (normal, moderate or severely malnourished) and two growth groups (stunted or not stunted). Of the 198 infants who had at least a single serum cationic trypsinogen measurement taken, normal values for serum IRT (with confidence limits) were obtained from 57 children, who were normally nourished. IRT levels were significantly correlated with the degree of underweight but there was no correlation with the degree of stunting or age. Mean IRT levels for the moderate and severely underweight groups were significantly greater than the mean for the normal group (P < 0.01). Seventeen children (8.6%) had trypsinogen levels in excess of the 95th percentile for the normally nourished group, reflecting acinar cell damage or ductal obstruction. We conclude that pancreatic dysfunction may be a common and important overlooked factor contributing to ongoing malnutrition and diseases in malnourished Australian Aboriginal children.

Working correlation structure misspecification, estimation and covariate design: Implications for generalised estimating equations performance

The method of generalised estimating equations for regression modelling of clustered outcomes allows for specification of a working matrix that is intended to approximate the true correlation matrix of the observations. We investigate the asymptotic relative efficiency of the generalised estimating equation for the mean parameters when the correlation parameters are estimated by various methods. The asymptotic relative efficiency depends on three-features of the analysis, namely (i) the discrepancy between the working correlation structure and the unobservable true correlation structure, (ii) the method by which the correlation parameters are estimated and (iii) the 'design', by which we refer to both the structures of the predictor matrices within clusters and distribution of cluster sizes. Analytical and numerical studies of realistic data-analysis scenarios show that choice of working covariance model has a substantial impact on regression estimator efficiency. Protection against avoidable loss of efficiency associated with covariance misspecification is obtained when a 'Gaussian estimation' pseudolikelihood procedure is used with an AR(1) structure.

Electron Mean Energy, Secondary Ionization Coefficients And (E/P)Crit In Binary-Mixtures Of Sf6-N2 And Ccl2f2-N2

Measurements of the ratio of diffusion coefficient to mobility (D/ mu ) of electrons in SF6-N2 and CCl2F2-N2 mixtures over the range 80

Dependence of low-lying energy eigenvalues on the short- and long-range parts of confining potentials

The potential description of a quark-antiquark system seems to work very well in describing a number of hadronic properties. However, the precise form of the potential is unknown. The changes in the low-lying eigenvalues as a result of changes in the long-range part of the potential are investigated in a non-perturbative manner. It is shown by considering a variety of examples that the low-lying eigenvalues are insensitive to the long-range part of the potential.

Performance analysis and control design of a gyro-based wave energy converter

This paper analyses the performance of particular wave-energy converter that uses the gyroscopic effects of a large rotating fly-wheel in combination with a controlled power-take-off device. Controlled gyroscopic forces have been used successfully in the past to reduce the motion of marine structures. With appropriately designed power-take-off elements, gyroscopic forces can be controlled to optimise the extracted energy from the motion of marine structures.

Structure-reactivity correlation in inclusion complexes: deoxycholic acid di-tert-butyl thioketone

Ultraviolet irradiation of crystalline molecular inclusion complexes of deoxycholic acid with di-tert-butyl thioketone results in no reaction. The structure of the above complex has been determined via X-ray diffraction. The absence of expected photoreactions. namely, photoreduction and photooxidation, is rationalized on the basis of the X-ray structure analysis of the complex.

Correlation of plasma IGF-I concentrations and growth rate in aquacultured finfish: a tool for assessing the potential of new diets

A recently developed radioimmunoassay (RIA) for measuring insulin-like growth factor (IGF-I) in a variety of fish species was used to investigate the correlation between growth rate and circulating IGF-I concentrations of barramundi (Lates calcarifer), Atlantic salmon (Salmo salar) and Southern Bluefin tuna (Thunnus maccoyii). Plasma IGF-I concentration significantly increased with increasing ration size in barramundi and IGF-I concentration was positively correlated to growth rates obtained in Atlantic salmon (r2=0.67) and barramundi (r2=0.65) when fed a variety of diet formulations. IGF-I was also positively correlated to protein concentration (r2=0.59). This evidence suggested that measuring IGF-I concentration may provide a useful tool for monitoring fish growth rate and also as a method to rapidly assess different aquaculture diets. However, no such correlation was demonstrated in the tuna study probably due to seasonal cooling of sea surface temperature shortly before blood was sampled. Thus, some recommendations for the design and sampling strategy of nutritional trials where IGF-I concentrations are measured are discussed

Short-term nutritional supplementation during management of pulmonary exacerbations in cystic fibosis: A controlled study, including effects of protein turnover

Effects of nutritional supplements on minimizing weight loss and abnormalities of protein turnover during pulmonary exacerbations in cystic fibrosis (CF) were studied by controlled trial. Patients received pulmonary therapy and either standard diet (n = 10) or adjunctive enteral supplements (n = 12). Initial protein turnover, measured by [15N]glycine kinetics, showed alterations of protein synthesis (P Syn) and catabolism (P Cat), which correlated with the degree of underweight, and negligible net protein deposition (P Dep). With treatment both groups had significant increases in mean body weight and forced expiratory volume in 1 s, expressed as percent predicted value for height (FEV1) by 3 wk, but a significant correlation between initial underweight and subsequent weight gain was observed only in supplemented patients. Mean P Syn and P Dep increased significantly (p < 0.001) only in the supplemented group. Pulmonary exacerbations in CF have important adverse effects on body-protein metabolism, similar to changes in protein-energy malnutrition and infection. These effects are reversed by short-term nutritional support. Strategic nutritional intervention should thus be considered in management, especially in malnourished patients.

Inverse Correlation between Vitamin D and C-Reactive Protein in Newborns

Some studies suggested that adequate vitamin D might reduce inflammation in adults. However, little is known about this association in early life. We aimed to determine the relationship between cord blood 25-hydroxyvitamin D (25(OH)D) and C-reactive protein (CRP) in neonates. Cord blood levels of 25(OH)D and CRP were measured in 1491 neonates in Hefei, China. Potential confounders including maternal sociodemographic characteristics, perinatal health status, lifestyle, and birth outcomes were prospectively collected. The average values of cord blood 25(OH)D and CRP were 39.43 nmol/L (SD = 20.35) and 6.71 mg/L (SD = 3.07), respectively. Stratified by 25(OH)D levels, per 10 nmol/L increase in 25(OH)D, CRP decreased by 1.42 mg/L (95% CI: 0.90, 1.95) among neonates with 25(OH)D <25.0 nmol/L, and decreased by 0.49 mg/L (95% CI: 0.17, 0.80) among neonates with 25(OH)D between 25.0 nmol/L and 49.9 nmol/L, after adjusting for potential confounders. However, no significant association between 25(OH)D and CRP was observed among neonates with 25(OH)D ≥50 nmol/L. Cord blood 25(OH)D and CRP levels showed a significant seasonal trend with lower 25(OH)D and higher CRP during winter-spring than summer-autumn. Stratified by season, a significant linear association of 25(OH)D with CRP was observed in neonates born in winter-spring (adjusted β = −0.11, 95% CI: −0.13, −0.10), but not summer-autumn. Among neonates born in winter-spring, neonates with 25(OH)D <25 nmol/L had higher risk of CRP ≥10 mg/L (adjusted OR = 3.06, 95% CI: 2.00, 4.69), compared to neonates with 25(OH)D ≥25 nmol/L. Neonates with vitamin D deficiency had higher risk of exposure to elevated inflammation at birth.