Designs for generalized linear models with several variables and model uncertainty

Standard factorial designs sometimes may be inadequate for experiments that aim to estimate a generalized linear model, for example, for describing a binary response in terms of several variables. A method is proposed for finding exact designs for such experiments that uses a criterion allowing for uncertainty in the link function, the linear predictor, or the model parameters, together with a design search. Designs are assessed and compared by simulation of the distribution of efficiencies relative to locally optimal designs over a space of possible models. Exact designs are investigated for two applications, and their advantages over factorial and central composite designs are demonstrated.

A flux ratio theorem for multicomponent linear diffusion equations

Ussing [1] considered the steady flux of a single chemical component diffusing through a membrane under the influence of chemical potentials and derived from his linear model, an expression for the ratio of this flux and that of the complementary experiment in which the boundary conditions were interchanged. Here, an extension of Ussing's flux ratio theorem is obtained for n chemically interacting components governed by a linear system of diffusion-migration equations that may also incorporate linear temporary trapping reactions. The determinants of the output flux matrices for complementary experiments are shown to satisfy an Ussing flux ratio formula for steady state conditions of the same form as for the well-known one-component case. (C) 2000 Elsevier Science Ltd. All rights reserved.

Use of confocal microscopy to follow the development of penetrative hyphae during growth of Rhizopus oligosporus in an artificial solid-state fermentation system

Two methods were compared for determining the concentration of penetrative biomass during growth of Rhizopus oligosporus on an artificial solid substrate consisting of an inert gel and starch as the sole source of carbon and energy. The first method was based on the use of a hand microtome to make sections of approximately 0.2- to 0.4-mm thickness parallel to the substrate surface and the determination of the glucosamine content in each slice. Use of glucosamine measurements to estimate biomass concentrations was shown to be problematic due to the large variations in glucosamine content with mycelial age. The second method was a novel method based on the use of confocal scanning laser microscopy to estimate the fractional volume occupied by the biomass. Although it is not simple to translate fractional volumes into dry weights of hyphae due to the lack of experimentally determined conversion factors, measurement of the fractional volumes in themselves is useful for characterizing fungal penetration into the substrate. Growth of penetrative biomass in the artificial model substrate showed two forms of growth with an indistinct mass in the region close to the substrate surface and a few hyphae penetrating perpendicularly to the surface in regions further away from the substrate surface. The biomass profiles against depth obtained from the confocal microscopy showed two linear regions on log-linear plots, which are possibly related to different oxygen availability at different depths within the substrate. Confocal microscopy has the potential to be a powerful tool in the investigation of fungal growth mechanisms in solid-state fermentation. (C) 2003 Wiley Periodicals, Inc.

Finite deformation model of simple shear of fault with microrotations: apparent strain localisation and en-echelon fracture pattern

Strain localisation is a widespread phenomenon often observed in shear and compressive loading of geomaterials, for example, the fault gouge. It is believed that the main mechanisms of strain localisation are strain softening and mismatch between dilatancy and pressure sensitivity. Observations show that gouge deformation is accompanied by considerable rotations of grains. In our previous work as a model for gouge material, we proposed a continuum description for an assembly of particles of equal radius in which the particle rotation is treated as an independent degree of freedom. We showed that there exist critical values of the model parameters for which the displacement gradient exhibits a pronounced localisation at the mid-surface layers of the fault, even in the absence of inelasticity. Here, we generalise the model to the case of finite deformations characteristic for the gouge deformation. We derive objective constitutive relationships relating the Jaumann rates of stress and moment stress to the relative strain and curvature rates, respectively. The model suggests that the pattern of localisation remains the same as in the linear case. However, the presence of the Jaumann terms leads to the emergence of non-zero normal stresses acting along and perpendicular to the shear layer (with zero hydrostatic pressure), and localised along the mid-line of the gouge; these stress components are absent in the linear model of simple shear. These additional normal stresses, albeit small, cause a change in the direction in which the maximal normal stresses act and in which en-echelon fracturing is formed.

Norwegian long-distance reflexives: Lectal variation and speaker attitudes

10 lectal variables were examined with respect to Norwegian speakers' acceptance of long-distance reflexives (LDR), using a questionnaire to elicit grammaticality judgements on 50 potential LDR sentences. A sample of 180 speakers completed the questionnaire. The data was analysed using a general linear model univariate model, and Spearman's correlation. In this sample the results showed that dialect and level of education had significant effects on speakers' acceptance of long-distance reflexives, while sex, age, being a native speaker, having both native-speaker parents, living in the city or the country, and the speaker's attitudes to the two Norwegian writing languages had no influence on speakers' acceptance of long-distance reflexives. It is suggested that the influence of Danish on Norwegian writing and on the southern dialects may be the cause of the observed variation with respect to LDR in Norwegian.

Population subdivision in marine environments: the contributions of biogeography, geographic distance, and discontinuous habitat to genetic differentiation in a blennioid

The relative importance of factors that may promote genetic differentiation in marine organisms is largely unknown. Here, contributions to population structure from biogeography, habitat distribution, and isolation by distance were investigated in Axoclinus nigricaudus, a small subtidal rock reef fish, throughout its range in the Gulf of California. A 408 basepair fragment of the mitochondrial control region was sequenced from 105 individuals. Variation was significantly partitioned between many pairs of populations. Phylogenetic analyses, hierarchical analyses of variance, and general linear models substantiated a major break between two putative biogeographic regions. This genetic discontinuity coincides with an abrupt change in ecological characteristics (including temperature and salinity) but does not coincide with known oceanographic circulation patterns. Geographic distance and the nature of habitat separating populations (continuous habitat along a shoreline, discontinuous habitat along a shoreline, and open water) also contributed to population structure in general linear model analyses. To verify that local populations are genetically stable over time, one population was resampled on four occasions over eighteen months; it showed no evidence of a temporal component to diversity. These results indicate that having a planktonic life stage does not preclude geographically partitioned genetic variation over relatively small geographic distances in marine environments. Moreover, levels of genetic differentiation among populations of Axoclinus nigricaudus cannot be explained by a single factor, but are due to the combined influences of a biogeographic boundary, habitat, and geographic distance.

Short-term spatial diffusion in sigma o+ o- optical molasses

We have measured the spatial diffusion of atoms in a three-dimensional sigma(+)-sigma(-) optical molasses over twenty milliseconds timescale, starting from the initial interaction of the atoms with the molasses. We find that the diffusion constants agree well with a linear model for these short time scales and also compare favourably to other studies of diffusion made over longer time scales. These measurements enable us to quantify the detection method known as freezing molasses. We discuss this method, for detecting and measuring the momentum distribution of cold atoms, which relies on the slow diffusion of atoms in optical molasses to produce a freeze-frame of the spatial distribution of the atoms. This method enables a longer interrogation interval, providing a greatly increased signal-to-noise ratio. (C) 1998 Elsevier Science B.V.

A comparison of lipid variables as predictors of cardiovascular disease in the Asia Pacific Region

PURPOSE: Many guidelines advocate measurement of total or low density lipoprotein cholesterol (LDL), high density lipoprotein cholesterol (HDL), and triglycerides (TG) to determine treatment recommendations for preventing coronary heart disease (CHD) and cardiovascular disease (CVD). This analysis is a comparison of lipid variables as predictors of cardiovascular disease. METHODS: Hazard ratios for coronary and cardiovascular deaths by fourths of total cholesterol (TC), LDL, HDL, TG, non-HDL, TC/HDL, and TG/HDL values, and for a one standard deviation change in these variables, were derived in an individual participant data meta-analysis of 32 cohort studies conducted in the Asia-Pacific region. The predictive value of each lipid variable was assessed using the likelihood ratio statistic. RESULTS: Adjusting for confounders and regression dilution, each lipid variable had a positive (negative for HDL) log-linear association with fatal CHD and CVD. Individuals in the highest fourth of each lipid variable had approximately twice the risk of CHD compared with those with lowest levels. TG and HDL were each better predictors of CHD and CVD risk compared with TC alone, with test statistics similar to TC/HDL and TG/HDL ratios. Calculated LDL was a relatively poor predictor. CONCLUSIONS: While LDL reduction remains the main target of intervention for lipid-lowering, these data support the potential use of TG or lipid ratios for CHD risk prediction. (c) 2005 Elsevier Inc. All rights reserved.

Anthropometric measures in relation to Basal Cell Carcinoma: a longitudinal study

Background: The relationship between anthropometric indices and risk of basal cell carcinoma ( BCC) is largely unknown. We aimed to examine the association between anthropometric measures and development of BCC and to demonstrate whether adherence to World Health Organisation guidelines for body mass index, waist circumference, and waist/ hip ratio was associated with risk of BCC, independent of sun exposure. Methods: Study participants were participants in a community- based skin cancer prevention trial in Nambour, a town in southeast Queensland ( latitude 26 degrees S). In 1992, height, weight, and waist and hip circumferences were measured for all 1621 participants and weight was remeasured at the end of the trial in 1996. Prevalence proportion ratios were calculated using a log- binomial model to estimate the risk of BCC prior to or prevalent in 1992, while Poisson regression with robust error variances was used to estimate the relative risk of BCC during the follow- up period. Results: At baseline, 94 participants had a current BCC, and 202 had a history of BCC. During the 5- year follow- up period, 179 participants developed one or more new BCCs. We found no significant association between any of the anthropometric measures or indices and risk of BCC after controlling for potential confounding factors including sun exposure. There was a suggestion that short- term weight gain may increase the risk of developing BCC for women only. Conclusion: Adherence to World Health Organisation guidelines for body mass index, waist circumference and waist/ hip ratio is not significantly associated with occurrence of basal cell carcinomas of the skin.

Disposition of naproxen, naproxen acyl glucuronide and its rearrangement isomers in the isolated perfused rat liver

1. An isolated perfused rat liver (IPRL) preparation was used to investigate separately the disposition of the non-steroidal anti-inflammatory drug (NSAID) naproxen (NAP), its reactive acyl glucuronide metabolite (NAG) and a mixture of NAG rearrangement isomers (isoNAG), each at 30 mug NAP equivalents ml(-1) perfusate (n = 4 each group). 2. Following administration to the IPRL, NAP was eliminated slowly in a log-linear manner with an apparent elimination half-life (t(1/2)) of 13.4 +/-4.4 h. No metabolites were detected in perfusate, while NAG was the only metabolite present in bile in measurable amounts (3.9 +/-0.8%, of the dose). Following their administration to the IPRL, both NAG and isoNAG were rapidly hydrolysed (t(1/2) in perfusate=57 +/-3 and 75 +/- 14min respectively). NAG also rearranged to isoNAG in the perfusate. Both NAG and isoNAG were excreted intact in bile (24.6 and 14.8% of the NAG and isoNAG doses, respectively). 3. Covalent NAP-protein adducts in the liver increased as the dose changed from NAP to NAG to isoNAG (0.20 to 0.34 to 0.48% of the doses, respectively). Similarly, formation of covalent NAP-protein adducts in perfusate were greater in isoNAG-dosed perfusions. The comparative results Suggest that isoNAG is a better substrate for adduct formation with liver proteins than NAG.

Levene tests of homogeneity of variance for general block and treatment designs

This article develops a weighted least squares version of Levene's test of homogeneity of variance for a general design, available both for univariate and multivariate situations. When the design is balanced, the univariate and two common multivariate test statistics turn out to be proportional to the corresponding ordinary least squares test statistics obtained from an analysis of variance of the absolute values of the standardized mean-based residuals from the original analysis of the data. The constant of proportionality is simply a design-dependent multiplier (which does not necessarily tend to unity). Explicit results are presented for randomized block and Latin square designs and are illustrated for factorial treatment designs and split-plot experiments. The distribution of the univariate test statistic is close to a standard F-distribution, although it can be slightly underdispersed. For a complex design, the test assesses homogeneity of variance across blocks, treatments, or treatment factors and offers an objective interpretation of residual plot.

The effect of a triclosan-containing dentifrice on the progression of periodontal disease in an adult population

Objectives: The aim of the present study was to determine the effect of unsupervised, long-term use of a 0.3% triclosan/2% copolymer dentifrice on the progression of periodontal disease in a general adult population. Methods: Five hundred and four volunteers were enrolled in a double-blind, controlled clinical trial. Participants were matched for disease status, plaque index, age and gender. At the baseline examination, probing pocket depths and relative attachment levels were recorded and participants were assigned to either the test or control group. Re-examinations took place after 6, 12, 24, 36, 48 and 60 months. Subgingival plaque samples were collected at each examination and assayed for Porphyromonas gingivalis , Actinobacillus actinomycetemcomitans and Prevotella intermedia . A generalised linear model was used to analyse the data, with a number of covariates thought to influence the responses included as the possible confounding effects. Results: The triclosan/copolymer dentifrice had a significant effect in subjects with interproximal probing depths greater than or equal to3.5 mm, where it significantly reduced the number of sites with probing depths greater than or equal to3.5 mm at the following examination, when compared with the control group (p

Independent contribution of plaque complexity to myocardial ischemia during Dobutamine stress echocardiography

The influence of complex plaque morphology on the extent of demand-induced ischemia in unselected patients is not well defined. We sought to investigate the functional significance of lesion morphology in patients who underwent coronary angiography and dobutamine stress echocardiography (DSE).,Angiography and DSE were performed within a 6-month period (mean 1 +/- 1 month) in 196 patients. Angiographic assessments involved quantification of stenosis severity, assessment of the extent of jeopardized myocardium, and categorization of plaque morphology according to the Ambrose classification. DSE was interpreted by separate investigators with respect to wall motion score index (WMSI) and number of coronary territories involved. A general linear model was constructed to assess,the independent contribution of patient characteristics and angiographic and DSE results with respect to extent of ischemic myocardium. Complex lesion morphology was seen in 62 patients (32%). Patients with complex lesions were more likely to have had prior myocardial infarction (p < 0.001) and be current smokers (p = 0.03). During angiography, they exhibited a trend toward a greater number of diseased vessels, had a greater coronary jeopardy score (p < 0.001) and more frequent collateral flow (p = 0.03). During echocardiography, patients had a higher stress WMSI (p < 0.001) and were more likely to show ischemia in all 3 arterial territories (p < 0.01). On multivariate regression, the coronary artery jeopardy score and the presence of complex plaque morphology were independent predictors of the extent of ischemic myocardium (R 2 = 34%, p < 0.001). Thus, patients with complex plaque morphology are older, more likely to smoke, and more likely to have had prior myocardial. infarction. They exhibit more extensive disease with higher coronary jeopardy scores and a higher resting and peak stress WMSI. Despite these differences, complex plaque morphology remains an independent predictor of the extent of ischemia during stress. (C) 2003 by Excerpta Medica, Inc.

Prospective evaluation of a D-optimal designed population pharmacokinetic study

Recently, methods for computing D-optimal designs for population pharmacokinetic studies have become available. However there are few publications that have prospectively evaluated the benefits of D-optimality in population or single-subject settings. This study compared a population optimal design with an empirical design for estimating the base pharmacokinetic model for enoxaparin in a stratified randomized setting. The population pharmacokinetic D-optimal design for enoxaparin was estimated using the PFIM function (MATLAB version 6.0.0.88). The optimal design was based on a one-compartment model with lognormal between subject variability and proportional residual variability and consisted of a single design with three sampling windows (0-30 min, 1.5-5 hr and 11 - 12 hr post-dose) for all patients. The empirical design consisted of three sample time windows per patient from a total of nine windows that collectively represented the entire dose interval. Each patient was assigned to have one blood sample taken from three different windows. Windows for blood sampling times were also provided for the optimal design. Ninety six patients were recruited into the study who were currently receiving enoxaparin therapy. Patients were randomly assigned to either the optimal or empirical sampling design, stratified for body mass index. The exact times of blood samples and doses were recorded. Analysis was undertaken using NONMEM (version 5). The empirical design supported a one compartment linear model with additive residual error, while the optimal design supported a two compartment linear model with additive residual error as did the model derived from the full data set. A posterior predictive check was performed where the models arising from the empirical and optimal designs were used to predict into the full data set. This revealed the optimal'' design derived model was superior to the empirical design model in terms of precision and was similar to the model developed from the full dataset. This study suggests optimal design techniques may be useful, even when the optimized design was based on a model that was misspecified in terms of the structural and statistical models and when the implementation of the optimal designed study deviated from the nominal design.

Improved slope adjustment functions for soil erosion prediction

Numerous studies in the last 60 years have investigated the relationship between land slope and soil erosion rates. However, relatively few of these have investigated slope gradient responses: ( a) for steep slopes, (b) for specific erosion processes, and ( c) as a function of soil properties. Simulated rainfall was applied in the laboratory on 16 soils and 16 overburdens at 100 mm/h to 3 replicates of unconsolidated flume plots 3 m long by 0.8 m wide and 0.15 m deep at slopes of 20, 5, 10, 15, and 30% slope in that order. Sediment delivery at each slope was measured to determine the relationship between slope steepness and erosion rate. Data from this study were evaluated alongside data and existing slope adjustment functions from more than 55 other studies from the literature. Data and the literature strongly support a logistic slope adjustment function of the form S = A + B/[1 + exp (C - D sin theta)] where S is the slope adjustment factor and A, B, C, and D are coefficients that depend on the dominant detachment and transport processes. Average coefficient values when interill-only processes are active are A - 1.50, B 6.51, C 0.94, and D 5.30 (r(2) = 0.99). When rill erosion is also potentially active, the average slope response is greater and coefficient values are A - 1.12, B 16.05, C 2.61, and D 8.32 (r(2) = 0.93). The interill-only function predicts increases in sediment delivery rates from 5 to 30% slope that are approximately double the predictions based on existing published interill functions. The rill + interill function is similar to a previously reported value. The above relationships represent a mean slope response for all soils, yet the response of individual soils varied substantially from a 2.5-fold to a 50-fold increase over the range of slopes studied. The magnitude of the slope response was found to be inversely related ( log - log linear) to the dispersed silt and clay content of the soil, and 3 slope adjustment equations are proposed that provide a better estimate of slope response when this soil property is known. Evaluation of the slope adjustment equations proposed in this paper using independent datasets showed that the new equations can improve soil erosion predictions.