Intra-cluster correlation structure in longitudinal data analysis: Selection criteria and misspecification tests

Selection criteria and misspecification tests for the intra-cluster correlation structure (ICS) in longitudinal data analysis are considered. In particular, the asymptotical distribution of the correlation information criterion (CIC) is derived and a new method for selecting a working ICS is proposed by standardizing the selection criterion as the p-value. The CIC test is found to be powerful in detecting misspecification of the working ICS structures, while with respect to the working ICS selection, the standardized CIC test is also shown to have satisfactory performance. Some simulation studies and applications to two real longitudinal datasets are made to illustrate how these criteria and tests might be useful.

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.

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.

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.

Induced smoothing for rank regression with censored survival times

Adaptions of weighted rank regression to the accelerated failure time model for censored survival data have been successful in yielding asymptotically normal estimates and flexible weighting schemes to increase statistical efficiencies. However, for only one simple weighting scheme, Gehan or Wilcoxon weights, are estimating equations guaranteed to be monotone in parameter components, and even in this case are step functions, requiring the equivalent of linear programming for computation. The lack of smoothness makes standard error or covariance matrix estimation even more difficult. An induced smoothing technique overcame these difficulties in various problems involving monotone but pure jump estimating equations, including conventional rank regression. The present paper applies induced smoothing to the Gehan-Wilcoxon weighted rank regression for the accelerated failure time model, for the more difficult case of survival time data subject to censoring, where the inapplicability of permutation arguments necessitates a new method of estimating null variance of estimating functions. Smooth monotone parameter estimation and rapid, reliable standard error or covariance matrix estimation is obtained.

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.

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.

Mandatory nap times and group napping patterns in child care: An observational study

Policy provision for naps is typical in child care settings, but there is variability in the practices employed. One practice that might modify children’s early sleep patterns is the allocation of a mandatory nap time in which all children are required to lie on their beds without alternate activity permitted. There is currently limited evidence of the effects of such practices on children’s napping patterns. This study examined the association between duration of mandatory nap times and group-level napping patterns in child care settings. Observations were undertaken in a community sample of 113 preschool rooms with a scheduled nap time (N = 2,114 children). Results showed that 83.5% of child care settings implemented a mandatory nap time (range = 15–145 min) while 14.2% provided alternate activities for children throughout the nap time period. Overall, 31% of children napped during nap times. Compared to rooms with ≤ 30 min of mandatory nap time, rooms with 31–60 min and > 60 min of mandatory nap time had a two-and-a-half and fourfold increase, respectively, in the proportion of children napping. Nap onset latency did not significantly differ across groups. Among preschool children, exposure to longer mandatory nap times in child care may increase incidence of napping.

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.

Spectral analysis of pair-correlation bandwidth: Application to cell biology images

Images from cell biology experiments often indicate the presence of cell clustering, which can provide insight into the mechanisms driving the collective cell behaviour. Pair-correlation functions provide quantitative information about the presence, or absence, of clustering in a spatial distribution of cells. This is because the pair-correlation function describes the ratio of the abundance of pairs of cells, separated by a particular distance, relative to a randomly distributed reference population. Pair-correlation functions are often presented as a kernel density estimate where the frequency of pairs of objects are grouped using a particular bandwidth (or bin width), Δ>0. The choice of bandwidth has a dramatic impact: choosing Δ too large produces a pair-correlation function that contains insufficient information, whereas choosing Δ too small produces a pair-correlation signal dominated by fluctuations. Presently, there is little guidance available regarding how to make an objective choice of Δ. We present a new technique to choose Δ by analysing the power spectrum of the discrete Fourier transform of the pair-correlation function. Using synthetic simulation data, we confirm that our approach allows us to objectively choose Δ such that the appropriately binned pair-correlation function captures known features in uniform and clustered synthetic images. We also apply our technique to images from two different cell biology assays. The first assay corresponds to an approximately uniform distribution of cells, while the second assay involves a time series of images of a cell population which forms aggregates over time. The appropriately binned pair-correlation function allows us to make quantitative inferences about the average aggregate size, as well as quantifying how the average aggregate size changes with time.

Past-day recall of sedentary time: Validity of a self-reported measure of sedentary time in a university population

Background Sedentary behaviour is associated with several deleterious health consequences. Although device-based measures of sedentary time are available, they are costly and do not provide a measure of domain specific sedentary time. High quality self-report measures are necessary to accurately capture domain specific sedentary time, and to provide an alternative to devices when cost is an issue. In this study, the Past-day Adults’ Sedentary Time (PAST) questionnaire, previously shown to have acceptable validity and reliability in a sample of breast cancer survivors, was modified for a university sample and validity of the modified questionnaire was examined compared with activPAL. Methods Participants (n = 58, age = 18–55 years, 48% female, 66% students) were recruited from the University of Queensland (students and staff). They answered the PAST questionnaire, which asked about time spent sitting or lying down for work, study, travel, television viewing, leisure-time computer use, reading, eating, socialising and other purposes, during the previous day. Time reported for these questions was summed to provide a measure of total sedentary time. Participants also wore an activPAL device for the full day prior to completing the questionnaire and recorded their wake and sleep times in an activity log. Total waking sedentary time derived from the activPAL was used as the criterion measure. Correlation (Pearson's r) and agreement (Bland–Altman plots) between PAST and activPAL sedentary time were examined. Results Participants were sedentary (activPAL-determined) for approximately 66% of waking hours. The correlation between PAST and activPAL sedentary time for the whole sample was r = 0.50 [95% confidence interval (CI) = 0.28–0.67]; and higher for non-students (r = 0.63, 95% CI = 0.26–0.84) than students (r = 0.46, 95% CI = 0.16–0.68). Bland–Altman plots revealed that the mean difference between the two measures was 19 min although limits of agreement were wide (95% limits of agreement −4.1 to 4.7 h). Discussion The PAST questionnaire provides an acceptable measure of sedentary time in this population, which included students and adults with high workplace sitting. These findings support earlier research that questionnaires employing past-day recall of sedentary time provide a viable alternative to existing sedentary behaviour questionnaires.

Supporting sleep in early care and education: An assessment of observed sleep-times using a sleep practices optimality index

Aim: The aim was to investigate whether the sleep practices in early childhood education (ECE) settings align with current evidence on optimal practice to support sleep. Background: Internationally, scheduled sleep times are a common feature of daily schedules in ECE settings, yet little is known about the degree to which care practices in these settings align with the evidence regarding appropriate support of sleep. Methods: Observations were conducted in 130 Australian ECE rooms attended by preschool children (Mean = 4.9 years). Of these rooms, 118 had daily scheduled sleep times. Observed practices were scored against an optimality index, the Sleep Environment and Practices Optimality Score, developed with reference to current evidence regarding sleep scheduling, routines, environmental stimuli, and emotional climate. Cluster analysis was applied to identify patterns and prevalence of care practices in the sleep time. Results: Three sleep practices types were identified. Supportive rooms (36%) engaged in practices that maintained regular schedules, promoted routine, reduced environmental stimulation, and maintained positive emotional climate. The majority of ECE rooms (64%), although offering opportunity for sleep, did not engage in supportive practices: Ambivalent rooms (45%) were emotionally positive but did not support sleep; Unsupportive rooms (19%) were both emotionally negative and unsupportive in their practices. Conclusions: Although ECE rooms schedule sleep time, many do not adopt practices that are supportive of sleep. Our results underscore the need for education about sleep supporting practice and research to ascertain the impact of sleep practices in ECE settings on children’s sleep health and broader well-being.

Crude oil and agricultural futures: An analysis of correlation dynamics

Correlations between oil and agricultural commodities have varied over previous decades, impacted by renewable fuels policy and turbulent economic conditions. We estimate smooth transition conditional correlation models for 12 agricultural commodities and WTI crude oil. While a structural change in correlations occurred concurrently with the introduction of biofuel policy, oil and food price levels are also key influences. High correlation between biofuel feedstocks and oil is more likely to occur when food and oil price levels are high. Correlation with oil returns is strong for biofuel feedstocks, unlike with other agricultural futures, suggesting limited contagion from energy to food markets.