Efficient parameter estimation in longitudinal data analysis using a hybrid GEE method

The method of generalized estimating equations (GEEs) provides consistent estimates of the regression parameters in a marginal regression model for longitudinal data, even when the working correlation model is misspecified (Liang and Zeger, 1986). However, the efficiency of a GEE estimate can be seriously affected by the choice of the working correlation model. This study addresses this problem by proposing a hybrid method that combines multiple GEEs based on different working correlation models, using the empirical likelihood method (Qin and Lawless, 1994). Analyses show that this hybrid method is more efficient than a GEE using a misspecified working correlation model. Furthermore, if one of the working correlation structures correctly models the within-subject correlations, then this hybrid method provides the most efficient parameter estimates. In simulations, the hybrid method's finite-sample performance is superior to a GEE under any of the commonly used working correlation models and is almost fully efficient in all scenarios studied. The hybrid method is illustrated using data from a longitudinal study of the respiratory infection rates in 275 Indonesian children.

Weighted rank regression for clustered data analysis

We consider ranked-based regression models for clustered data analysis. A weighted Wilcoxon rank method is proposed to take account of within-cluster correlations and varying cluster sizes. The asymptotic normality of the resulting estimators is established. A method to estimate covariance of the estimators is also given, which can bypass estimation of the density function. Simulation studies are carried out to compare different estimators for a number of scenarios on the correlation structure, presence/absence of outliers and different correlation values. The proposed methods appear to perform well, in particular, the one incorporating the correlation in the weighting achieves the highest efficiency and robustness against misspecification of correlation structure and outliers. A real example is provided for illustration.

Process Analytical Technology Approach on Fluid Bed Granulation and Drying : Identifying Critical Relationships and Constructing the Design Space

Fluid bed granulation is a key pharmaceutical process which improves many of the powder properties for tablet compression. Dry mixing, wetting and drying phases are included in the fluid bed granulation process. Granules of high quality can be obtained by understanding and controlling the critical process parameters by timely measurements. Physical process measurements and particle size data of a fluid bed granulator that are analysed in an integrated manner are included in process analytical technologies (PAT). Recent regulatory guidelines strongly encourage the pharmaceutical industry to apply scientific and risk management approaches to the development of a product and its manufacturing process. The aim of this study was to utilise PAT tools to increase the process understanding of fluid bed granulation and drying. Inlet air humidity levels and granulation liquid feed affect powder moisture during fluid bed granulation. Moisture influences on many process, granule and tablet qualities. The approach in this thesis was to identify sources of variation that are mainly related to moisture. The aim was to determine correlations and relationships, and utilise the PAT and design space concepts for the fluid bed granulation and drying. Monitoring the material behaviour in a fluidised bed has traditionally relied on the observational ability and experience of an operator. There has been a lack of good criteria for characterising material behaviour during spraying and drying phases, even though the entire performance of a process and end product quality are dependent on it. The granules were produced in an instrumented bench-scale Glatt WSG5 fluid bed granulator. The effect of inlet air humidity and granulation liquid feed on the temperature measurements at different locations of a fluid bed granulator system were determined. This revealed dynamic changes in the measurements and enabled finding the most optimal sites for process control. The moisture originating from the granulation liquid and inlet air affected the temperature of the mass and pressure difference over granules. Moreover, the effects of inlet air humidity and granulation liquid feed rate on granule size were evaluated and compensatory techniques used to optimize particle size. Various end-point indication techniques of drying were compared. The ∆T method, which is based on thermodynamic principles, eliminated the effects of humidity variations and resulted in the most precise estimation of the drying end-point. The influence of fluidisation behaviour on drying end-point detection was determined. The feasibility of the ∆T method and thus the similarities of end-point moisture contents were found to be dependent on the variation in fluidisation between manufacturing batches. A novel parameter that describes behaviour of material in a fluid bed was developed. Flow rate of the process air and turbine fan speed were used to calculate this parameter and it was compared to the fluidisation behaviour and the particle size results. The design space process trajectories for smooth fluidisation based on the fluidisation parameters were determined. With this design space it is possible to avoid excessive fluidisation and improper fluidisation and bed collapse. Furthermore, various process phenomena and failure modes were observed with the in-line particle size analyser. Both rapid increase and a decrease in granule size could be monitored in a timely manner. The fluidisation parameter and the pressure difference over filters were also discovered to express particle size when the granules had been formed. The various physical parameters evaluated in this thesis give valuable information of fluid bed process performance and increase the process understanding.

Rank-based regression for analysis of repeated measures

We consider rank-based regression models for repeated measures. To account for possible withinsubject correlations, we decompose the total ranks into between- and within-subject ranks and obtain two different estimators based on between- and within-subject ranks. A simple perturbation method is then introduced to generate bootstrap replicates of the estimating functions and the parameter estimates. This provides a convenient way for combining the corresponding two types of estimating function for more efficient estimation.

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.

A generalized estimating equations approach for analysis of the impact of new technology on a trawl fishery

The article describes a generalized estimating equations approach that was used to investigate the impact of technology on vessel performance in a trawl fishery during 1988-96, while accounting for spatial and temporal correlations in the catch-effort data. Robust estimation of parameters in the presence of several levels of clustering depended more on the choice of cluster definition than on the choice of correlation structure within the cluster. Models with smaller cluster sizes produced stable results, while models with larger cluster sizes, that may have had complex within-cluster correlation structures and that had within-cluster covariates, produced estimates sensitive to the correlation structure. The preferred model arising from this dataset assumed that catches from a vessel were correlated in the same years and the same areas, but independent in different years and areas. The model that assumed catches from a vessel were correlated in all years and areas, equivalent to a random effects term for vessel, produced spurious results. This was an unexpected finding that highlighted the need to adopt a systematic strategy for modelling. The article proposes a modelling strategy of selecting the best cluster definition first, and the working correlation structure (within clusters) second. The article discusses the selection and interpretation of the model in the light of background knowledge of the data and utility of the model, and the potential for this modelling approach to apply in similar statistical situations.

Influence of environmental and instrumental variables on the non-invasive prediction of Brix in pineapple using near infrared spectroscopy

The Brix content of pineapple fruit can be non-invasively predicted from the second derivative of near infrared reflectance spectra. Correlations obtained using a NIRSystems 6500 spectrophotometer through multiple linear regression and modified partial least squares analyses using a post-dispersive configuration were comparable with that from a pre-dispersive configuration in terms of accuracy (e.g. coefficient of determination, R2, 0.73; standard error of cross validation, SECV, 1.01°Brix). The effective depth of sample assessed was slightly greater using the post-dispersive technique (about 20 mm for pineapple fruit), as expected in relation to the higher incident light intensity, relative to the pre-dispersive configuration. The effect of such environmental variables as temperature, humidity and external light, and instrumental variables such as the number of scans averaged to form a spectrum, were considered with respect to the accuracy and precision of the measurement of absorbance at 876 nm, as a key term in the calibration for Brix, and predicted Brix. The application of post-dispersive near infrared technology to in-line assessment of intact fruit in a packing shed environment is discussed.

Non-invasive assessment of pineapple and mango fruit quality using near infra-red spectroscopy

The potential of near infra-red (NIR) spectroscopy for non-invasive measurement of fruit quality of pineapple (Ananas comosus var. Smooth Cayenne) and mango (Magnifera indica var. Kensington) fruit was assessed. A remote reflectance fibre optic probe, placed in contact with the fruit skin surface in a light-proof box, was used to deliver monochromatic light to the fruit, and to collect NIR reflectance spectra (760–2500 nm). The probe illuminated and collected reflected radiation from an area of about 16 cm2. The NIR spectral attributes were correlated with pineapple juice Brix and with mango flesh dry matter (DM) measured from fruit flesh directly underlying the scanned area. The highest correlations for both fruit were found using the second derivative of the spectra (d2 log 1/R) and an additive calibration equation. Multiple linear regression (MLR) on pineapple fruit spectra (n = 85) gave a calibration equation using d2 log 1/R at wavelengths of 866, 760, 1232 and 832 nm with a multiple coefficient of determination (R2) of 0.75, and a standard error of calibration (SEC) of 1.21 °Brix. Modified partial least squares (MPLS) regression analysis yielded a calibration equation with R2 = 0.91, SEC = 0.69, and a standard error of cross validation (SECV) of 1.09 oBrix. For mango, MLR gave a calibration equation using d2 log 1/R at 904, 872, 1660 and 1516 nm with R2 = 0.90, and SEC = 0.85% DM and a bias of 0.39. Using MPLS analysis, a calibration equation with R2 = 0.98, SEC = 0.54 and SECV = 1.19 was obtained. We conclude that NIR technology offers the potential to assess fruit sweetness in intact whole pineapple and DM in mango fruit, respectively, to within 1° Brix and 1% DM, and could be used for the grading of fruit in fruit packing sheds.

A Genetic Analysis of Lysosomal Enzyme Activities in Brahman Cattle

Analyses of variance and co variance were carried out on the activities of three lysosomal enzymes in mononuclear blood cells from Brahman cattle. These were hexosaminidase (HEX), beta-D-galacto-sidase (GAL) and acid alpha-glucosidase (GLU) which had been measured in blood mononuclear cells from 1752 cattle from 6 herds in a Pompe's disease control programme. Herd of origin and date of bleeding significantly affected the level of activity of all enzymes. In addition, HEX and GAL were affected by age and HEX by the sex of the animal bled. Estimates of heritability from sire variances were 0.29:t 0.09 for HEX, 0.31 :t 0.09 for GAL and 0.44:t 0.09 for GLU. Genetic correlations between all enzymes were positive. The data indicate the existence of a major gene causing Pompe's disease and responsible for 16% of the genetic variation in GLU. One standard deviation of selection differential for high GLU should almost eliminate Pompe's disease from the population. The effi-ciency of selection would be aided by estimating the breeding value for GLU using measurements of HEX and GLU and taking account of an animal's sex, age, date of bleeding and herd of origin.

Plasma Biochemical and Haematological Changes in Steers of Different Temperament during Feedlotting

Correlations between liveweight gains and temperament during feedlotting (Voisinet et al., 1997) may be an indication that cattle of different temperament suffer different degrees of stress. This paper describes the effect of grouping into feedlot pens, cattle of good, poor and mixed (some good and some poor) temperament on blood parameters that are indicators of stress (Bennett et al., 1989). Animal production for a consuming world : proceedings of 9th Congress of the Asian-Australasian Association of Animal Production Societies [AAAP] and 23rd Biennial Conference of the Australian Society of Animal Production [ASAP] and 17th Annual Symposium of the University of Sydney, Dairy Research Foundation, [DRF]. 2-7 July 2000, Sydney, Australia.

Relationships between Patterns of Behaviour and Liveweight Gains of Feedlot Steers

There is considerable individual variation in the performance of feedlot cattle even when they are of the same genotype and sourced from the same property (Hasker et al., 1996). It is possible that behavioural differences between individuals may account for some of this variation. This paper reports the correlations between some behavioural patterns and average daily gain (ADG) of Brahman steers that were feedlot-fed for 100 days.

Genetic Parameters for Staple Length and Staple Strength of Merino Wool Produced in Central and North West Queensland

Estimates of genetic parameters are presented for staple length and staple strength for 15 month old, medium Peppin sheep at Longreach and Julia Creek Queensland. The effects of birth type, sex and year of birth are shown. There were significant interactions for sex by site and for sex by year of birth. Heritability of staple length and strength were respectively 0.75 and 0.37 for the Longreach flock and 0.70 and 0.23 for the Julia Creek flock. The heritability of staple strength agrees with other published data however the estimate for staple length is very high. Phenotypic and genetic correlations with greasy fleece weight, yield, clean fleece weight, average fibre diameter and liveweight are in general agreement with other published estimates. Animal production for a consuming world : proceedings of 9th Congress of the Asian-Australasian Association of Animal Production Societies [AAAP] and 23rd Biennial Conference of the Australian Society of Animal Production [ASAP] and 17th Annual Symposium of the University of Sydney, Dairy Research Foundation, [DRF]. 2-7 July 2000 University of New South Wales, Sydney, Australia.

Genetic Parameters for Physiological Characters in Merino Rams in Central and North West Queensland

Genetic and phenotypic parameters for respiration rate (RR) and rectal temperature (RT) are presented for weaner and hogget Merino rams, at Longreach and Julia Creek, Queensland. Heritability estimates for RT and RR at both sites and at both ages ranged from moderate to very high. Phenotypic and genetic correlations between these characters are also reported. AAABG 14th Conference; Proceedings of the Association for the Advancement of Animal Breeding and Genetics. AAABG

Distribution and abundance of white grubs (Coleoptera: Scarabaeidae) on groundnut in southern India

Scarab species associated with groundnuts were surveyed in Andhra Pradesh, Karnataka and Tamil Nadu, southern India, between 1995 and 2001. Scarab adults were collected from trees on which they were feeding and/or mating, and larvae (white grubs) from groundnut fields. Holotrichia species, especially H. reynaudi and H. serrata were the major species associated with groundnut. H. reynaudi predominated in the central Deccan area, while H. serrata was most abundant in areas to the south and west. A new, undescribed, Holotrichia species near H. consanguinea was collected south and south-west of Hyderabad in mixed populations with H. reynaudi. However, the full extent of this new species’ distribution remains uncertain. H. rufoflava was rarely associated with groundnut, but was common as an adult at some locations. Other genera encountered during surveys were Anomala, Adoretus, Schizonycha, Autoserica. In survey data, densities of Holotrichia larvae and ‘all other white grubs’ were both very highly correlated with % of damaged groundnut plants. These correlations in combination with concurrent observations of plant damage establish a causal link between white grubs and plant damage and death in southern Indian groundnut. Ranking of preferred host trees for adults were developed from field observations for four Holotrichia species and Schizonycha spp. and will assist grower-initiated surveys of pest occurrence. In combination with insecticide efficacy data published elsewhere, the survey provides the basis for an environmentally friendly and economically viable pest-management system for white grubs on groundnut in southern India.

Many-body effects in electrochemisorption : Part III. Beyond mean field and hartree-fock

Several channels provided by many-body couplings — both fermion-fermion and fermion-boson — for the evolution of the chemisorption system are discussed. This provides an opportunity of a systematic study of the effects of correlations reflected through the intricate pole structure of the absorbate Green functions. The results of Newns, Anda and others in the context of chemisorption are generalized.