A Gaussian pseudolikelihood approach for quantile regression with repeated measurements

To enhance the efficiency of regression parameter estimation by modeling the correlation structure of correlated binary error terms in quantile regression with repeated measurements, we propose a Gaussian pseudolikelihood approach for estimating correlation parameters and selecting the most appropriate working correlation matrix simultaneously. The induced smoothing method is applied to estimate the covariance of the regression parameter estimates, which can bypass density estimation of the errors. Extensive numerical studies indicate that the proposed method performs well in selecting an accurate correlation structure and improving regression parameter estimation efficiency. The proposed method is further illustrated by analyzing a dental dataset.

Repeated measurements of learned irrelevance by a novel within-subject paradigm in humans

Learned irrelevance (LIrr) refers to the retardation of classical conditioning following preexposure of the to-be-associated stimuli. Healthy volunteers have been tested on three occasions with a new LIrr paradigm avoiding methodological problems which afflict traditional paradigms. A significant LIrr effect was demonstrated on each occasion. Thus, the new paradigm enables repeated measurements of LIrr and might be useful in evaluating long-term effects of medication in psychiatric disorders exhibiting aberrant LIrr.

Repeated measurements of cerebral blood flow in the left superior temporal gyrus reveal tonic hyperactivity in patients with auditory verbal hallucinations: a possible trait marker

Background: The left superior temporal gyrus (STG) has been suggested to play a key role in auditory verbal hallucinations (AVH) in patients with schizophrenia. Methods: Eleven medicated subjects with schizophrenia and medication-resistant AVH and 19 healthy controls underwent perfusion magnetic resonance (MR) imaging with arterial spin labeling (ASL). Three additional repeated measurements were conducted in the patients. Patients underwent a treatment with transcranial magnetic stimulation (TMS) between the first 2 measurements. The main outcome measure was the pooled cerebral blood flow (CBF), which consisted of the regional CBF measurement in the left STG and the global CBF measurement in the whole brain. Results: Regional CBF in the left STG in patients was significantly higher compared to controls (p < 0.0001) and to the global CBF in patients (p < 0.004) at baseline. Regional CBF in the left STG remained significantly increased compared to the global CBF in patients across time (p < 0.0007), and it remained increased in patients after TMS compared to the baseline CBF in controls (p < 0.0001). After TMS, PANSS (p = 0.003) and PSYRATS (p = 0.01) scores decreased significantly in patients. Conclusions: This study demonstrated tonically increased regional CBF in the left STG in patients with schizophrenia and auditory hallucinations despite a decrease in symptoms after TMS. These findings were consistent with what has previously been termed a trait marker of AVH in schizophrenia.

Increasing the power of genome wide association studies in natural populations using repeated measures : evaluation and implementation

1. Genomewide association studies (GWAS) enable detailed dissections of the genetic basis for organisms' ability to adapt to a changing environment. In long-term studies of natural populations, individuals are often marked at one point in their life and then repeatedly recaptured. It is therefore essential that a method for GWAS includes the process of repeated sampling. In a GWAS, the effects of thousands of single-nucleotide polymorphisms (SNPs) need to be fitted and any model development is constrained by the computational requirements. A method is therefore required that can fit a highly hierarchical model and at the same time is computationally fast enough to be useful. 2. Our method fits fixed SNP effects in a linear mixed model that can include both random polygenic effects and permanent environmental effects. In this way, the model can correct for population structure and model repeated measures. The covariance structure of the linear mixed model is first estimated and subsequently used in a generalized least squares setting to fit the SNP effects. The method was evaluated in a simulation study based on observed genotypes from a long-term study of collared flycatchers in Sweden. 3. The method we present here was successful in estimating permanent environmental effects from simulated repeated measures data. Additionally, we found that especially for variable phenotypes having large variation between years, the repeated measurements model has a substantial increase in power compared to a model using average phenotypes as a response. 4. The method is available in the R package RepeatABEL. It increases the power in GWAS having repeated measures, especially for long-term studies of natural populations, and the R implementation is expected to facilitate modelling of longitudinal data for studies of both animal and human populations.

Corneal topography with Scheimpflug imaging and videokeratography : comparative study of normal eyes

PURPOSE: To compare the repeatability within anterior corneal topography measurements and agreement between measurements with the Pentacam HR rotating Scheimpflug camera and with a previously validated Placido disk–based videokeratoscope (Medmont E300). ------ SETTING: Contact Lens and Visual Optics Laboratory, School of Optometry, Queensland University of Technology, Brisbane, Queensland, Australia. ----- METHODS: Normal eyes in 101 young adult subjects had corneal topography measured using the Scheimpflug camera (6 repeated measurements) and videokeratoscope (4 repeated measurements). The best-fitting axial power corneal spherocylinder was calculated and converted into power vectors. Corneal higher-order aberrations (HOAs) (up to the 8th Zernike order) were calculated using the corneal elevation data from each instrument. ----- RESULTS: Both instruments showed excellent repeatability for axial power spherocylinder measurements (repeatability coefficients <0.25 diopter; intraclass correlation coefficients >0.9) and good agreement for all power vectors. Agreement between the 2 instruments was closest when the mean of multiple measurements was used in analysis. For corneal HOAs, both instruments showed reasonable repeatability for most aberration terms and good correlation and agreement for many aberrations (eg, spherical aberration, coma, higher-order root mean square). For other aberrations (eg, trefoil and tetrafoil), the 2 instruments showed relatively poor agreement. ----- CONCLUSIONS: For normal corneas, the Scheimpflug system showed excellent repeatability and reasonable agreement with a previously validated videokeratoscope for the anterior corneal axial curvature best-fitting spherocylinder and several corneal HOAs. However, for certain aberrations with higher azimuthal frequencies, the Scheimpflug system had poor agreement with the videokeratoscope; thus, caution should be used when interpreting these corneal aberrations with the Scheimpflug system.

An economic evaluation of the Safe Motherhood programme in Guangxi, China

Maternal and infant mortality is a global health issue with a significant social and economic impact. Each year, over half a million women worldwide die due to complications related to pregnancy or childbirth, four million infants die in the first 28 days of life, and eight million infants die in the first year. Ninety-nine percent of maternal and infant deaths are in developing countries. Reducing maternal and infant mortality is among the key international development goals. In China, the national maternal mortality ratio and infant mortality rate were reduced greatly in the past two decades, yet a large discrepancy remains between urban and rural areas. To address this problem, a large-scale Safe Motherhood Programme was initiated in 2000. The programme was implemented in Guangxi in 2003. Interventions in the programme included both demand-side and supply side-interventions focusing on increasing health service use and improving birth outcomes. Little is known about the effects and economic outcomes of the Safe Motherhood Programme in Guangxi, although it has been implemented for seven years. The aim of this research is to estimate the effectiveness and cost-effectiveness of the interventions in the Safe Motherhood Programme in Guangxi, China. The objectives of this research include: 1. To evaluate whether the changes of health service use and birth outcomes are associated with the interventions in the Safe Motherhood Programme. 2. To estimate the cost-effectiveness of the interventions in the Safe Motherhood Programme and quantify the uncertainty surrounding the decision. 3. To assess the expected value of perfect information associated with both the whole decision and individual parameters, and interpret the findings to inform priority setting in further research and policy making in this area. A quasi-experimental study design was used in this research to assess the effectiveness of the programme in increasing health service use and improving birth outcomes. The study subjects were 51 intervention counties and 30 control counties. Data on the health service use, birth outcomes and socio-economic factors from 2001 to 2007 were collected from the programme database and statistical yearbooks. Based on the profile plots of the data, general linear mixed models were used to evaluate the effectiveness of the programme while controlling for the effects of baseline levels of the response variables, change of socio-economic factors over time and correlations among repeated measurements from the same county. Redundant multicollinear variables were deleted from the mixed model using the results of the multicollinearity diagnoses. For each response variable, the best covariance structure was selected from 15 alternatives according to the fit statistics including Akaike information criterion, Finite-population corrected Akaike information criterion, and Schwarz.s Bayesian information criterion. Residual diagnostics were used to validate the model assumptions. Statistical inferences were made to show the effect of the programme on health service use and birth outcomes. A decision analytic model was developed to evaluate the cost-effectiveness of the programme, quantify the decision uncertainty, and estimate the expected value of perfect information associated with the decision. The model was used to describe the transitions between health states for women and infants and reflect the change of both costs and health benefits associated with implementing the programme. Result gained from the mixed models and other relevant evidence identified were synthesised appropriately to inform the input parameters of the model. Incremental cost-effectiveness ratios of the programme were calculated for the two groups of intervention counties over time. Uncertainty surrounding the parameters was dealt with using probabilistic sensitivity analysis, and uncertainty relating to model assumptions was handled using scenario analysis. Finally the expected value of perfect information for both the whole model and individual parameters in the model were estimated to inform priority setting in further research in this area.The annual change rates of the antenatal care rate and the institutionalised delivery rate were improved significantly in the intervention counties after the programme was implemented. Significant improvements were also found in the annual change rates of the maternal mortality ratio, the infant mortality rate, the incidence rate of neonatal tetanus and the mortality rate of neonatal tetanus in the intervention counties after the implementation of the programme. The annual change rate of the neonatal mortality rate was also improved, although the improvement was only close to statistical significance. The influences of the socio-economic factors on the health service use indicators and birth outcomes were identified. The rural income per capita had a significant positive impact on the health service use indicators, and a significant negative impact on the birth outcomes. The number of beds in healthcare institutions per 1,000 population and the number of rural telephone subscribers per 1,000 were found to be positively significantly related to the institutionalised delivery rate. The length of highway per square kilometre negatively influenced the maternal mortality ratio. The percentage of employed persons in the primary industry had a significant negative impact on the institutionalised delivery rate, and a significant positive impact on the infant mortality rate and neonatal mortality rate. The incremental costs of implementing the programme over the existing practice were US $11.1 million from the societal perspective, and US $13.8 million from the perspective of the Ministry of Health. Overall, 28,711 life years were generated by the programme, producing an overall incremental cost-effectiveness ratio of US $386 from the societal perspective, and US $480 from the perspective of the Ministry of Health, both of which were below the threshold willingness-to-pay ratio of US $675. The expected net monetary benefit generated by the programme was US $8.3 million from the societal perspective, and US $5.5 million from the perspective of the Ministry of Health. The overall probability that the programme was cost-effective was 0.93 and 0.89 from the two perspectives, respectively. The incremental cost-effectiveness ratio of the programme was insensitive to the different estimates of the three parameters relating to the model assumptions. Further research could be conducted to reduce the uncertainty surrounding the decision, in which the upper limit of investment was US $0.6 million from the societal perspective, and US $1.3 million from the perspective of the Ministry of Health. It is also worthwhile to get a more precise estimate of the improvement of infant mortality rate. The population expected value of perfect information for individual parameters associated with this parameter was US $0.99 million from the societal perspective, and US $1.14 million from the perspective of the Ministry of Health. The findings from this study have shown that the interventions in the Safe Motherhood Programme were both effective and cost-effective in increasing health service use and improving birth outcomes in rural areas of Guangxi, China. Therefore, the programme represents a good public health investment and should be adopted and further expanded to an even broader area if possible. This research provides economic evidence to inform efficient decision making in improving maternal and infant health in developing countries.

Evaluation of the iPhone with an acrylic sleeve versus the Scoliometer for rib hump measurement in scoliosis

Background. Vertebral rotation found in structural scoliosis contributes to trunkal asymmetry which is commonly measured with a simple Scoliometer device on a patient's thorax in the forward flexed position. The new generation of mobile 'smartphones' have an integrated accelerometer, making accurate angle measurement possible, which provides a potentially useful clinical tool for assessing rib hump deformity. This study aimed to compare rib hump angle measurements performed using a Smartphone and traditional Scoliometer on a set of plaster torsos representing the range of torsional deformities seen in clinical practice. Methods. Nine observers measured the rib hump found on eight plaster torsos moulded from scoliosis patients with both a Scoliometer and an Apple iPhone on separate occasions. Each observer repeated the measurements at least a week after the original measurements, and were blinded to previous results. Intra-observer reliability and inter-observer reliability were analysed using the method of Bland and Altman and 95% confidence intervals were calculated. The Intra-Class Correlation Coefficients (ICC) were calculated for repeated measurements of each of the eight plaster torso moulds by the nine observers. Results. Mean absolute difference between pairs of iPhone/Scoliometer measurements was 2.1 degrees, with a small (1 degrees) bias toward higher rib hump angles with the iPhone. 95% confidence intervals for intra-observer variability were +/- 1.8 degrees (Scoliometer) and +/- 3.2 degrees (iPhone). 95% confidence intervals for inter-observer variability were +/- 4.9 degrees (iPhone) and +/- 3.8 degrees (Scoliometer). The measurement errors and confidence intervals found were similar to or better than the range of previously published thoracic rib hump measurement studies. Conclusions. The iPhone is a clinically equivalent rib hump measurement tool to the Scoliometer in spinal deformity patients. The novel use of plaster torsos as rib hump models avoids the variables of patient fatigue and discomfort, inconsistent positioning and deformity progression using human subjects in a single or multiple measurement sessions.

Seasonal variation in cardiovascular disease risk factors in a subarctic population : the Tromso Study 1979-2008

Background Seasonal changes in cardiovascular disease (CVD) risk factors may be due to exposure to seasonal environmental variables like temperature and acute infections or seasonal behavioural patterns in physical activity and diet. Investigating the seasonal pattern of risk factors should help determine the causes of the seasonal pattern in CVD. Few studies have investigated the seasonal variation in risk factors using repeated measurements from the same individual, which is important as individual and population seasonal patterns may differ. Methods The authors investigated the seasonal pattern in systolic and diastolic blood pressure, heart rate, body weight, total cholesterol, triglycerides, high-density lipoprotein cholesterol, C reactive protein and fibrinogen. Measurements came from 38 037 participants in the population-based cohort, the Tromsø Study, examined up to eight times from 1979 to 2008. Individual and population seasonal patterns were estimated using a cosinor in a mixed model. Results All risk factors had a highly statistically significant seasonal pattern with a peak time in winter, except for triglycerides (peak in autumn), C reactive protein and fibrinogen (peak in spring). The sizes of the seasonal variations were clinically modest. Conclusions Although the authors found highly statistically significant individual seasonal patterns for all risk factors, the sizes of the changes were modest, probably because this subarctic population is well adapted to a harsh climate. Better protection against seasonal risk factors like cold weather could help reduce the winter excess in CVD observed in milder climates.

A homogeneous surface-enhanced Raman scattering platform for ultra-trace detection of trinitrotoluene in the environment

A facile and sensitive surface-enhanced Raman scattering substrate was prepared by controlled potentiostatic deposition of a closely packed single layer of gold nanostructures (AuNS) over a flat gold (pAu) platform. The nanometer scale inter-particle distance between the particles resulted in high population of ‘hot spots’ which enormously enhanced the scattered Raman photons. A renewed methodology was followed to precisely quantify the SERS substrate enhancement factor (SSEF) and it was estimated to be (2.2 ± 0.17) × 105. The reproducibility of the SERS signal acquired by the developed substrate was tested by establishing the relative standard deviation (RSD) of 150 repeated measurements from various locations on the substrate surface. A low RSD of 4.37 confirmed the homogeneity of the developed substrate. The sensitivity of pAu/AuNS was proven by determining 100 fM 2,4,6-trinitrotoluene (TNT) comfortably. As a proof of concept on the potential of the new pAu/AuNS substrate in field analysis, TNT in soil and water matrices was selectively detected after forming a Meisenheimer complex with cysteamine.

Ability of commonly used prediction equations to predict resting energy expenditure in children with inflammatory bowel disease

Background: Paediatric onset inflammatory bowel disease (IBD) may cause alterations in energy requirements and invalidate the use of standard prediction equations. Our aim was to evaluate four commonly used prediction equations for resting energy expenditure (REE) in children with IBD. Methods: Sixty-three children had repeated measurements of REE as part of a longitudinal research study yielding a total of 243 measurements. These were compared with predicted REE from Schofield, Oxford, FAO/WHO/UNU, and Harris-Benedict equations using the Bland-Altman method. Results: Mean (±SD) age of the patients was 14.2 (2.4) years. Mean measured REE was 1566 (336) kcal per day compared with 1491 (236), 1441 (255), 1481 (232), and 1435 (212) kcal per day calculated from Schofield, Oxford, FAO/WHO/UNU, and Harris-Benedict, respectively. While the Schofield equation demonstrated the least difference between measured and predicted REE, it, along with the other equations tested, did not perform uniformly across all subjects, indicating greater errors at either end of the spectrum of energy expenditure. Smaller differences were found for all prediction equations for Crohn's disease compared with ulcerative colitis. Conclusions: Of the commonly used equations, the equation of Schofield should be used in pediatric patients with IBD when measured values are not able to be obtained. (Inflamm Bowel Dis 2010;) Copyright © 2010 Crohn's & Colitis Foundation of America, Inc.

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.

Serial Changes in Markers Measuring Coagulation, Fibrinolysis, and Vasoactivity in Patients with Ischemic Stroke

Ischemic stroke (IS) is a heterogeneous disease in which outcome is influenced by many factors. The hemostatic system is activated in association with cerebral ischemia, and thus, markers measuring coagulation, fibrinolysis, and vasoactivity could be useful tools in clinical practice. We investigated whether repeated measurements of these markers reveal patterns that might help in evaluating IS patients, including the early diagnosis of stroke subtypes, in estimating prognosis and risk of recurrence, and in selecting a treatment for secondary prevention of stroke. Vasoconstrictor peptide endothelin-1 (ET-1), homocysteine (Hcy), indicators of thrombin formation and activation (prothrombin fragment 1+2/F1+2, thrombin-antithrombin complex/TAT), indicators of plasmin formation and fibrinolysis (tissue plasminogen activator/t-PA, plasminogen activator inhibitor-1/PAI-1, and D-dimer), and natural anticoagulants (antithrombin/AT, protein C/PC, and protein S/PS) were measured in 102 consecutive mild to moderate IS patients on four occasions: on admission and at 1 week, 1 month, and 3 months after stroke, and once in controls. All patients underwent neurological examination and blood sampling in the same session. Furthermore, 42 IS patients with heterozygous factor V Leiden mutation (FVLm) were selected from 740 IS patients without an obvious etiology, and evaluated in detail for specific clinical, laboratory, and radiological features. Measurements of ET-1 and Hcy levels did not disclose information that could aid in the diagnostic evaluation of IS patients. F1+2 level at 3 months after IS had a positive correlation with recurrence of thromboembolic events, and thus, may be used as a predictive marker of subsequent cerebral events. The D-dimer and AT levels on admission and 1 week after IS were strongly associated with stroke severity, outcome, and disability. The specific analysis of IS patients with FVLm more often revealed a positive family history of thrombosis, a higher prevalence of peripheral vascular disease, and multiple infarctions in brain images, most of which were `silent infarcts´. Results of this study support the view that IS patients with sustained activation of both the fibrinolytic and the coagulation systems and increased thrombin generation may have an unfavorable prognosis. The level of activation may reflect the ongoing thrombotic process and the extent of thrombosis. Changes in these markers could be useful in predicting prognosis of IS patients. A clear need exists for a randomized prospective study to determine whether a subgroup of IS patients with markers indicating activation of fibrinolytic and coagulation systems might benefit from more aggressive secondary prevention of IS.

Context dependency of rewards and services in an Indian ant-plant interaction: southern sites favour the mutualism between plants and ants

Protection-based ant-plant mutualisms may vary in strength due to differences in ant rewards, abundance of protective ants and herbivory pressure. We investigated geographical and temporal variation in host plant traits and herbivory pressure at five sites spanning the distribution range of the myrmecophyte Humboldtia brunonis (Fabaceae) in the Indian Western Ghats. Southern siteshad, onaverage, 2.4 times greater abundance of domatia-bearing individuals, 1.6 times greater extrafloral nectary numbers per leaf, 1.2 times larger extrafloral nectary sizes, 2.2 times greater extrafloral nectar (EFN) volumes and a two-fold increase in total amino acid and total sugar concentrations in EFN compared with northern sites. Astrong protection-based mutualismwith ants occurred at only one southern site where herbivory was highest, suggesting that investments in attracting ants correlate with anti-herbivore benefits gained from the presence of protective ants. Our results confirm a temporally stable north-south gradient in myrmecophytic traits in this ant-plant as several of these traits were re-sampled after a 5-y interval. However, the chemical composition of EFN varied at both spatial and short-term temporal scales suggesting that only repeated measurements of rewards such as EFN can reveal the real spectrum of trait variation in an ant-plant mutualistic system.

Avaliação do consumo alimentar no período reprodutivo e variação de peso no pós-parto

O presente estudo tem como tema central o ganho de peso e o consumo alimentar no período reprodutivo. A tese está dividida em três partes: i) revisão da literatura científica sobre os fatores de risco para o ganho de peso excessivo na gestação e para a sua retenção no pós-parto; ii) análise quantitativa e qualitativa do consumo alimentar (energia, nutrientes e alimentos) do período gestacional para o pós-parto; e iii) avaliação prospectiva da associação entre a dieta hiperproteica (≥ 1.2 g/kg) e a variação de peso pós-gestacional. A revisão da literatura apontou como principais determinantes da variação de peso no pós-parto os seguintes fatores: ganho de peso gestacional, peso pré-gestacional, dieta, atividade física, lactação, idade, escolaridade, renda, paridade e raça. As partes ii e iii correspondem a dois artigos, sendo que a população do estudo empregada nas duas análises foi baseada em uma coorte de mulheres no pós-parto atendidas no Centro Municipal de Saúde (CMS) Marcolino Candau, localizado no Município do Rio de Janeiro. Dentre as 709 mulheres convidadas para participar da pesquisa, 479 ingressaram na coorte. As mulheres foram recrutadas através dos serviços de pré-natal e na rotina pediátrica de imunização do BCG, oferecidos pelo CMS, e imediatamente após o parto, na maternidade central de referência. A coleta de dados ocorreu entre maio de 1999 e abril de 2001, sendo 15 meses de recrutamento e 9 meses de seguimento. As mulheres foram entrevistadas aos 15 dias, 2, 6 e 9 meses do pós-parto. Foram consideradas elegíveis para as duas analises as mulheres com idade entre 18 e 45 anos de idade, entrevistadas até 30 dias do pós-parto, sem doenças crônicas pré-existentes, com gravidez de 35 ou mais semanas gestacionais e sem gestação gemelar. O estudo sobre a composição da dieta do período gestacional para o pós-parto contou com 276 mulheres que responderam aos dois questionários de freqüência de consumo alimentar (QFCA) aplicados aos 15 dias e aos 6 meses do seguimento. Os resultados mostraram que as mulheres consumiram dietas mais energéticas na gravidez e aquelas que mais restringiram o consumo de energia no pós-parto, apresentaram um incremento na densidade de proteínas da dieta. As análises das medidas repetidos do peso foram realizadas para avaliar o efeito da dieta hiperproteica (≥ 1.2 g/kg) sobre a variação de peso no pós-parto, para as quais haviam dados de uma amostra de 430 mulheres. Os resultados mostraram que as mulheres com dieta hiperproteica perderam mais peso do que as mulheres com dieta normo ou hipoproteica (< 1.2 g/kg) (226 g/mês versus 123 g/mês). As duas analises indicam que um pequeno incremento no consumo de proteínas promove maior perda de peso. Estudos longitudinais futuros devem incluir avaliação relativa a segurança das dietas hiperproteicas no pós-parto.