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.

Estimation of growth parameters from multiple-recapture data

This article develops a method for analysis of growth data with multiple recaptures when the initial ages for all individuals are unknown. The existing approaches either impute the initial ages or model them as random effects. Assumptions about the initial age are not verifiable because all the initial ages are unknown. We present an alternative approach that treats all the lengths including the length at first capture as correlated repeated measures for each individual. Optimal estimating equations are developed using the generalized estimating equations approach that only requires the first two moment assumptions. Explicit expressions for estimation of both mean growth parameters and variance components are given to minimize the computational complexity. Simulation studies indicate that the proposed method works well. Two real data sets are analyzed for illustration, one from whelks (Dicathais aegaota) and the other from southern rock lobster (Jasus edwardsii) in South Australia.

Estimating equations for parameters in stochastic growth models from tag-recapture data

James (1991, Biometrics 47, 1519-1530) constructed unbiased estimating functions for estimating the two parameters in the von Bertalanffy growth curve from tag-recapture data. This paper provides unbiased estimating functions for a class of growth models that incorporate stochastic components and explanatory variables. a simulation study using seasonal growth models indicates that the proposed method works well while the least-squares methods that are commonly used in the literature may produce substantially biased estimates. The proposed model and method are also applied to real data from tagged rack lobsters to assess the possible seasonal effect on growth.

Guaiac versus immunochemical tests: Faecal occult blood test screening for colorectal cancer in a rural community

Objective: To describe patient participation and clinical performance in a colorectal cancer (CRC) screening program utilising faecal occult blood test (FOBT). Methods: A community-based intervention was conducted in a small, rural community in north Queensland, 2000/01. One of two FOBT kits – guaiac (Hemoccult-ll) or immunochemical (Inform) – was assigned by general practice and mailed to participants (3,358 patients aged 50–74 years listed with the local practices). Results: Overall participation in FOBT screening was 36.3%. Participation was higher with the immunochemical kit than the guaiac kit (OR=1.9, 95% Cl 1.6-2.2). Women were more likely to comply with testing than men (OR=1.4, 95% Cl 1.2-1.7), and people in their 60s were less likely to participate than those 70–74 years (OR=0.8, 95% Cl 0.6-0.9). The positivity rate was higher for the immunochemical (9.5%) than the guaiac (3.9%) test (χ2=9.2, p=0.002), with positive predictive values for cancer or adenoma of advanced pathology of 37.8% (95% Cl 28.1–48.6) for !nform and 40.0% (95% Cl 16.8–68.7) for Hemoccult-ll. Colonoscopy follow-up was 94.8% with a medical complication rate of 2–3%. Conclusions: An immunochemical FOBT enhanced participation. Higher positivity rates for this kit did not translate into higher false-positive rates, and both test types resulted in a high yield of neoplasia. Implications: In addition to type of FOBT, the ultimate success of a population-based screening program for CRC using FOBT will depend on appropriate education of health professionals and the public as well as significant investment in medical infrastructure for colonoscopy follow-up.

PMU measurement based dynamic load modeling using SVC devices in online environment

Online dynamic load modeling has become possible with the availability of Static Voltage Compensator (SVC) and Phasor Measurement Unit (PMU) devices. The power of the load response to the small random bounded voltage fluctuations caused from SVC can be measured by PMU for modelling purposes. The aim of this paper is to illustrate the capability of identifying an aggregated load model from high voltage substation level in the online environment. The induction motor is used as the main test subject since it contributes the majority of the dynamic loads. A test system representing simple electromechanical generator model serving dynamic loads through the transmission network is used to verify the proposed method. Also, dynamic load with multiple induction motors are modeled to achieve a better realistic load representation.

Replication of the association of common rs9939609 variant of FTO with increased BMI in an Australian adult twin population but no evidence for gene by environment (G x E) interaction

OBJECTIVE: To further investigate a common variant (rs9939609) in the fat mass- and obesity-associated gene (FTO), which recent genome-wide association studies have shown to be associated with body mass index (BMI) and obesity. DESIGN: We examined the effect of this FTO variant on BMI in 3353 Australian adult male and female twins. RESULTS: The minor A allele of rs9939609 was associated with an increased BMI (P=0.0007). Each additional copy of the A allele was associated with a mean BMI increase of approximately 1.04 kg/m(2) (approximately 3.71 kg). Using variance components decomposition, we estimate that this single-nucleotide polymorphism accounts for approximately 3% of the genetic variance in BMI in our sample (approximately 2% of the total variance). By comparing intrapair variances of monozygotic twins of different genotypes we were able to perform a direct test of gene by environment (G x E) interaction in both sexes and gene by parity (G x P) interaction in women, but no evidence was found for either. CONCLUSIONS: In addition to supporting earlier findings that the rs9939609 variant in the FTO gene is associated with an increased BMI, our results indicate that the associated genetic effect does not interact with environment or parity.

Bayesian model comparison in cosmology with Population Monte Carlo

We use Bayesian model selection techniques to test extensions of the standard flat LambdaCDM paradigm. Dark-energy and curvature scenarios, and primordial perturbation models are considered. To that end, we calculate the Bayesian evidence in favour of each model using Population Monte Carlo (PMC), a new adaptive sampling technique which was recently applied in a cosmological context. The Bayesian evidence is immediately available from the PMC sample used for parameter estimation without further computational effort, and it comes with an associated error evaluation. Besides, it provides an unbiased estimator of the evidence after any fixed number of iterations and it is naturally parallelizable, in contrast with MCMC and nested sampling methods. By comparison with analytical predictions for simulated data, we show that our results obtained with PMC are reliable and robust. The variability in the evidence evaluation and the stability for various cases are estimated both from simulations and from data. For the cases we consider, the log-evidence is calculated with a precision of better than 0.08. Using a combined set of recent CMB, SNIa and BAO data, we find inconclusive evidence between flat LambdaCDM and simple dark-energy models. A curved Universe is moderately to strongly disfavoured with respect to a flat cosmology. Using physically well-motivated priors within the slow-roll approximation of inflation, we find a weak preference for a running spectral index. A Harrison-Zel'dovich spectrum is weakly disfavoured. With the current data, tensor modes are not detected; the large prior volume on the tensor-to-scalar ratio r results in moderate evidence in favour of r=0.

In vitro demonstration of the heavy-atom effect for photodynamic therapy

Photodynamic therapy (PDT) is an emerging treatment modality for a range of disease classes, both cancerous and noncancerous. This has brought about an active pursuit of new PDT agents that can be optimized for the unique set of photophysical characteristics that are required for a successful clinical agent. We now describe a totally new class of PDT agent, the BF2-chelated 3,5-diaryl-1H-pyrrol-2-yl-3,5-diarylpyrrol-2-ylideneamines (tetraarylazadipyrromethenes). Optimized synthetic procedures have been developed to facilitate the generation of an array of specifically substituted derivatives to demonstrate how control of key therapeutic parameters such as wavelength of maximum absorbance and singlet-oxygen generation can be achieved. Photosensitizer absorption maxima can be varied within the body's therapeutic window between 650 and 700 nm, with high extinction coefficients ranging from 75,000 to 85,000 M(-1) cm(-1). Photosensitizer singlet-oxygen generation level was modulated by the exploitation of the heavy-atom effect. An array of photosensitizers with and without bromine atom substituents gave rise to a series of compounds with varying singlet-oxygen generation profiles. X-ray structural evidence indicates that the substitution of the bromine atoms has not caused a planarity distortion of the photosensitizer. Comparative singlet-oxygen production levels of each photosensitizer versus two standards demonstrated a modulating effect on singlet-oxygen generation depending upon substituent patterns about the photosensitizer. Confocal laser scanning microscopy imaging of 18a in HeLa cervical carcinoma cells proved that the photosensitizer was exclusively localized to the cellular cytoplasm. In vitro light-induced toxicity assays in HeLa cervical carcinoma and MRC5-SV40 transformed fibroblast cancer cell lines confirmed that the heavy-atom effect is viable in a live cellular system and that it can be exploited to modulate assay efficacy. Direct comparison of the efficacy of the photosensitizers 18b and 19b, which only differ in molecular structure by the presence of two bromine atoms, illustrated an increase in efficacy of more than a 1000-fold in both cell lines. All photosensitizers have very low to nondeterminable dark toxicity in our assay system.