Fitting mixtures of Kent distributions to aid in joint set identification

When examining a rock mass, joint sets and their orientations can play a significant role with regard to how the rock mass will behave. To identify joint sets present in the rock mass, the orientation of individual fracture planer can be measured on exposed rock faces and the resulting data can be examined for heterogeneity. In this article, the expectation-maximization algorithm is used to lit mixtures of Kent component distributions to the fracture data to aid in the identification of joint sets. An additional uniform component is also included in the model to accommodate the noise present in the data.

X-ray scattering studies of mixed langmuir monolayers and Langmuir-Blodgett films of a noncentrosymmetric porphyrin with cadmium arachidate

The structures of mixed Langmuir (floating) monolayers and Langmuir-Blodgett (LB) films of a phenanthroline-porphyrin with cadmium arachidate (PhenPor + CdAr) have been investigated by synchrotron X-ray grazing incidence diffraction (GIXD) and specular X-ray reflectivity (SXR). GIXD measurements of the floating monolayers showed only one peak, arising from the CdAr domains in the films, at a scattering angle of 21.5 degrees. This is consistent with a hexagonal structure (alpha = 4.77 Angstrom). The correlation length in these domains is 250 Angstrom. GMD measurements of the LB films, however, show two sets of diffraction features: one arises from CdAr domains with a rectangular in-plane structure (alpha = 7.44 Angstrom and b = 4.90 Angstrom) and a correlation length of 85 Angstrom; the other is from porphyrin domains with an oblique in-plane structure (alpha (p) 15.2 Angstrom, b(p) = 8.86 Angstrom, and gamma (p) = 80 degrees) and a correlation length of 105 Angstrom. These dimensions are consistent with the surface pressure-area isotherm measurements and indicate that the two components are immiscible. The thickness of the bilayer is 57 Angstrom, and there is no correlation between the bilayers. Introduction of a trigger compound does not alter the structure of the films but slightly increases the bilayer thickness. The SXR measurements of the floating monolayers also support the suggested immiscibility of the two components in the films.

A methodology for analysis of sugarcane productivity trends - I. Analysis across districts

Historically, few articles have addressed the use of district level mill production data for analysing the effect of varietal change on sugarcane productivity trends. This appears to be due to lack of compiled district data sets and appropriate methods by which to analyse these data. Recently, varietal data on tonnes of sugarcane per hectare (TCH), sugar content (CCS), and their product, tonnes of sugar content per hectare (TSH) on a district basis, have been compiled. This study was conducted to develop a methodology for regular analysis of such data from mill districts to assess productivity trends over time, accounting for variety and variety x environment interaction effects for 3 mill districts (Mulgrave, Babinda, and Tully) from 1958 to 1995. Restricted maximum likelihood methodology was used to analyse the district level data and best linear unbiased predictors for random effects, and best linear unbiased estimates for fixed effects were computed in a mixed model analysis. In the combined analysis over districts, Q124 was the top ranking variety for TCH, and Q120 was top ranking for both CCS and TSH. Overall production for TCH increased over the 38-year period investigated. Some of this increase can be attributed to varietal improvement, although the predictors for TCH have shown little progress since the introduction of Q99 in 1976. Although smaller gains have been made in varietal improvement for CCS, overall production for CCS decreased over the 38 years due to non-varietal factors. Varietal improvement in TSH appears to have peaked in the mid-1980s. Overall production for TSH remained stable over time due to the varietal increase in TCH and the non-varietal decrease in CCS.

Adaptation to chronic eccentric exercise in humans: the influence of contraction velocity

We compared changes in muscle fibre composition and muscle strength indices following a 10 week isokinetic resistance training programme consisting of fast (3.14 rad(.)s(-1)) or slow (0.52 rad(.)s(-1)) velocity eccentric muscle contractions. A group of 20 non-resistance trained subjects were assigned to a FAST (n = 7), SLOW (n = 6) or non-training CONTROL (n = 7) group. A unilateral training protocol targeted the elbow flexor muscle group and consisted of 24 maximal eccentric isokinetic contractions (four sets of six repetitions) performed three times a week for 10 weeks. Muscle biopsy samples were obtained from the belly of the biceps brachii. Isometric torque and concentric and eccentric torque at 0.52 and 3.14 rad(.)s(-1) were examined at 0, 5 and 10 weeks. After 10 weeks, the FAST group demonstrated significant [mean (SEM)] increases in eccentric [29.6 (6.4)%] and concentric torque [27.4 (7.3) %] at 3.14 rad(.)s(-1), isometric torque [21.3 (4.3)%] and eccentric torque [25.2 (7.2) %] at 0.52 rad(.)s(-1). The percentage of type I fibres in the FAST group decreased from [53.8 (6.6)% to 39.1 (4.4)%] while type lib fibre percentage increased from [5.8 (1.9)% to 12.9 (3.3)%; P < 0.05]. In contrast. the SLOW group did not experience significant changes in muscle fibre type or muscle torque. We conclude that neuromuscular adaptations to eccentric training stimuli may be influenced by differences in the ability to cope with chronic exposure to relatively fast and slow eccentric contraction velocities. Possible mechanisms include greater cumulative damage to contractile tissues or stress induced by slow eccentric muscle contractions.

A total-evidence phylogeny of ticks provides insights into the evolution of life cycles and biogeography

We inferred the phylogeny of 33 species of ticks from the subfamilies Rhipicephalinae and Hyalomminae from analyses of nuclear and mitochondrial DNA and morphology. We used nucleotide sequences from 12S rRNA, cytochrome c oxidase I, internal transcribed spacer 2 of the nuclear rRNA, and 18S rRNA. Nucleotide sequences and morphology were analyzed separately and together in a total-evidence analysis. Analyses of the five partitions together (3303 characters) gave the best-resolved and the best-supported hypothesis so far for the phylogeny of ticks in the Rhipicephalinae and Hyalomminae, despite the fact that some partitions did not have data for some taxa. However, most of the hidden conflict (lower support in the total-evidence analyses compared to that in the individual analyses) was found in those partitions that had taxa without data. The partitions with complete taxonomic sampling had more hidden support (higher support in the total-evidence analyses compared to that in the separate-partition analyses) than hidden conflict. Mapping of geographic origins of ticks onto our phylogeny indicates an African origin for the Rhipicephalinae sensu lato (i.e., including Hyalomma spp.), the Rhipicephalus-Boophilus lineage, the Dermacentor-Anocentor lineage, and the Rhipicephalus-Booophilus-Nosomma-Hyalomma-Rhipicentor lineage. The Nosomma-Hyalomma lineage appears to have evolved in Asia. Our total-evidence phylogeny indicates that (i) the genus Rhipicephalus is paraphyletic with respect to the genus Boophilus, (ii) the genus Dermacentor is paraphyletic with respect to the genus Anocentor, and (iii) some subgenera of the genera Hyalomma and Rhipicephalus are paraphyletic with respect to other subgenera in these genera. Study of the Rhipicephalinae and Hyalomminae over the last 7 years has shown that analyses of individual datasets (e.g., one gene or morphology) seldom resolve many phylogenetic relationships, but analyses of more than one dataset can generate well-resolved phylogenies for these ticks. (C) 2001 Academic Press.

Associations involving delays (particularly long delays) between certain weather parameters and geomagnetic activity

Four sunspot-minimum periods (1963-1966, 1971-1977, 1983-1987 and 1992-1997) have been examined for the results which are presented. Using several different weather parameters, tropospheric gravity waves, enhanced cold fronts and two rainfall data sets in Eastern Australia, associations at reasonably high levels of significance have been found with enhanced geomagnetic activity (EGA). Statistically this EGA involved either short delays of several days or long delays of about 20 days. The geomagnetic parameters used were (a) the AE index (b) the hourly H component for a number of stations and (c) the daily K-P-sum value. The K-P-sum analyses have shown that the EGA associated with the delays form part of four or five cycles of recurrent geomagnetic activity for 27-day periodicities. Furthermore statistically two recurrent cycles are found to exist concurrently, one apparently related to the short delays and the other to the long delays. Periodicities of 13.5 days are created because the two sets are displaced from each other by approximately this interval. A brief reference is made to the 13.5 periodicity known to exist for geomagnetic activity and the evidence in the literature for active regions on the sun to be displaced by 180 degrees of solar longitude.

Comments on alternative calculations of the broadening of spectral lines of neutral sodium by H-atom collisions

With the exception of the sodium D-lines, recent calculations of line broadening cross sections for several multiplets of sodium by Leininger et al (Leininger T, Gadea F X and Dickinson A 2000 J. Phys. B: At. Mol. Opt. Phys. 33 1805) are in substantial disagreement with cross sections interpolated from the tables of Anstee and O'Mara (Anstee and O'Mara 1995 Mon. Not. R. Astron. Soc. 276 859) and Barklem and O'Mara (Barklem P S and O'Mara B J 1997 Mon. Not. R. Astron. Soc. 290 102). The discrepancy is as large as a factor of 3 for the 3p-4d multiplet. The two theories are tested by using the results of each to synthesize lines in the solar spectrum. It is found that generally the data from the theory of Anstee, Barklem and O'Mara produce the best match to the observed solar spectrum. It is found, using a simple model for reflection of the optical electron by the potential barrier between the two atoms, that the reflection coefficient is too large for avoided crossings with the upper states of subordinate lines to contribute to line broadening, supporting the neglect of avoided ionic crossings by Anstee, Barklem and O'Mara for these lines. The large discrepancies between the two sets of calculations is a result of an approximate treatment of avoided ionic crossings for these lines by Leininger et al (Leininger T, Gadea F X and Dickinson A 2000 J. Phys. B: At. Mol. Opt. Phys. 33 1805).

Investigating technical and scale efficiencies of Australian universities through data envelopment analysis

Performance indicators in the public sector have often been criticised for being inadequate and not conducive to analysing efficiency. The main objective of this study is to use data envelopment analysis (DEA) to examine the relative efficiency of Australian universities. Three performance models are developed, namely, overall performance, performance on delivery of educational services, and performance on fee-paying enrolments. The findings based on 1995 data show that the university sector was performing well on technical and scale efficiency but there was room for improving performance on fee-paying enrolments. There were also small slacks in input utilisation. More universities were operating at decreasing returns to scale, indicating a potential to downsize. DEA helps in identifying the reference sets for inefficient institutions and objectively determines productivity improvements. As such, it can be a valuable benchmarking tool for educational administrators and assist in more efficient allocation of scarce resources. In the absence of market mechanisms to price educational outputs, which renders traditional production or cost functions inappropriate, universities are particularly obliged to seek alternative efficiency analysis methods such as DEA.

Carapace dentition patterns, morphometrics and allozyme differentiation amongst two toothed freshwater crab species (Potamonautes warreni and P-unispinus) (Decapoda : Brachyura : Potamonautidae) from river systems in South Africa

The taxonomic relationship between two toothed South African river crabs, Potamonautes warreni and P. unispinus, is unclear. The problem stems from the widespread variation in carapace dentition patterns amongst P. warreni individuals over its biogeographic range, where single toothed individuals may appear similar in carapace morphology to P. unispinus. Ten populations of P. warreni and 18 populations of P. unispinus were collected and the morphometric and genetic differentiation between the two taxa quantified. Patterns of morphometric and genetic variation were examined using multivariate statistics and protein gel electrophoresis, respectively. Principal component analyses of carapace characters showed that the two species are morphologically indistinguishable. However, discriminate functions analyses and additional statistical results corroborate the morphological distinction between the two taxa. Allozyme electrophoresis of 17 protein coding loci, indicated a close genetic similarity between the two species (I = 0.92). A fixed allelic difference at one locus (LT-2) and extensive genetic variability at another locus (PGM-1) indicate that two gene pools are present and that the two taxa are genetically isolated. Intraspecific genetic I values for both species were > 0.97 and indicated no apparent genetic structuring on a micro or macro-geographic scale. The variation in carapace dentition among P. warreni populations possesses no genetic basis and may possibly toe the product of ecogenesis. The value of dentition patterns in the systematics of river crabs is discussed. Dentition patterns among river crab species appear to be conserved and reliable as species specific diagnostic markers, but should ideally be used in combination with other morphological data sets and genetic evidence.

The GRIP domain is a specific targeting sequence for a population of trans-Golgi network derived tubulo-vesicular carriers

Vesicular carriers for intracellular transport associate with unique sets of accessory molecules that dictate budding and docking on specific membrane domains. Although many of these accessory molecules are peripheral membrane proteins, in most cases the targeting sequences responsible for their membrane recruitment have yet to be identified. We have previously defined a novel Golgi targeting domain (GRIP) shared by a family of coiled-coil peripheral membrane Golgi proteins implicated in membrane trafficking. We show here that the docking site for the GRIP motif of p230 is a specific domain of Golgi. membranes. By immunoelectron microscopy of HeLa cells stably expressing a green fluorescent protein (GFP)-p230(GRIP) fusion protein, we show binding specifically to a subset of membranes of the trans-Golgi network (TGN). Real-time imaging of live HeLa cells revealed that the GFP-p230(GRIP) was associated with highly dynamic tubular extensions of the TGN, which have the appearance and behaviour of transport carriers. To further define the nature of the GRIP membrane binding site, in vitro budding assays were performed using purified rat liver Golgi membranes and cytosol from GFP-p230(GRIP) transfected cells. Analysis of Golgi-derived vesicles by sucrose gradient fractionation demonstrated that GFP-p230(GRIP) binds to a specific population of vesicles distinct from those labelled for beta -COP or gamma -adaptin. The GFP-p230(GRIP) fusion protein is recruited to the same vesicle population as full-length p230, demonstrating that the GRIP domain is solely proficient as a targeting signal for membrane binding of the native molecule. Therefore, p230 GRIP is a targeting signal for recruitment to a highly selective membrane attachment site on a specific population of trans-Golgi network tubulovesicular carriers.

Precipitation hardening of Cu-Fe-Cr alloys - Part I - Mechanical and electrical properties

This research is part of a project whose scope was to investigate the engineering properties of new non-commercial alloy formulations based on the Cu rich corner of the Cu-Fe-Cr ternary system with the primary aim of exploring the development of a new cost-effective high-strength, high-conductivity copper alloy. The literature indicated that Cu rich Cu-Cr and Cu-Fe alloys have been thoroughly investigated. A number of commercial alloys have been developed and these are used for a variety of applications requiring combinations of high-strength, high-conductivity and resistance to softening. Little evidence was found in the literature that the Cu rich corner of the Cu-Fe-Cr system had previously been investigated for the purpose of developing high-strength, high-conductivity copper alloys resistant to softening. The aim of these present investigations was to explore the possibility that new alloys could be developed that combined the properties of both sets of alloys, ie large precipitation hardening response combined with the ability to stabilise cold worked microstructures to high temperatures while at the same maintain high electrical conductivity. To assess the feasibility of this goal the following alloys were chosen for investigation: Cu-0.7wt%Cr-0.3wt%Fe, Cu-0.7wt%Cr-0.8wt%Fe, Cu-0.7wt%Cr-2.0wt%Fe. This paper reports on the mechanical property investigation which indicated that the Cu-0.7wt%Cr-0.3wt%Fe, and Cu-0.7wt%Cr-2.0wt%Fe alloys were worthy of further investigation. (C) 2001 Kluwer Academic Publishers.

Brief note on the variation of constants formula for fuzzy differential equations

This note gives a theory of state transition matrices for linear systems of fuzzy differential equations. This is used to give a fuzzy version of the classical variation of constants formula. A simple example of a time-independent control system is used to illustrate the methods. While similar problems to the crisp case arise for time-dependent systems, in time-independent cases the calculations are elementary solutions of eigenvalue-eigenvector problems. In particular, for nonnegative or nonpositive matrices, the problems at each level set, can easily be solved in MATLAB to give the level sets of the fuzzy solution. (C) 2002 Elsevier Science B.V. All rights reserved.

Riddling and invariance for discontinuous maps preserving Lebesgue measure

In this paper we use the mixture of topological and measure-theoretic dynamical approaches to consider riddling of invariant sets for some discontinuous maps of compact regions of the plane that preserve two-dimensional Lebesgue measure. We consider maps that are piecewise continuous and with invertible except on a closed zero measure set. We show that riddling is an invariant property that can be used to characterize invariant sets, and prove results that give a non-trivial decomposion of what we call partially riddled invariant sets into smaller invariant sets. For a particular example, a piecewise isometry that arises in signal processing (the overflow oscillation map), we present evidence that the closure of the set of trajectories that accumulate on the discontinuity is fully riddled. This supports a conjecture that there are typically an infinite number of periodic orbits for this system.

Theory and applications of fuzzy Volterra integral equations

Formulations of fuzzy integral equations in terms of the Aumann integral do not reflect the behavior of corresponding crisp models. Consequently, they are ill-adapted to describe physical phenomena, even when vagueness and uncertainty are present. A similar situation for fuzzy ODEs has been obviated by interpretation in terms of families of differential inclusions. The paper extends this formalism to fuzzy integral equations and shows that the resulting solution sets and attainability sets are fuzzy and far better descriptions of uncertain models involving integral equations. The investigation is restricted to Volterra type equations with mildly restrictive conditions, but the methods are capable of extensive generalization to other types and more general assumptions. The results are illustrated by integral equations relating to control models with fuzzy uncertainties.

A mixture model-based approach to the clustering of microarray expression data

Motivation: This paper introduces the software EMMIX-GENE that has been developed for the specific purpose of a model-based approach to the clustering of microarray expression data, in particular, of tissue samples on a very large number of genes. The latter is a nonstandard problem in parametric cluster analysis because the dimension of the feature space (the number of genes) is typically much greater than the number of tissues. A feasible approach is provided by first selecting a subset of the genes relevant for the clustering of the tissue samples by fitting mixtures of t distributions to rank the genes in order of increasing size of the likelihood ratio statistic for the test of one versus two components in the mixture model. The imposition of a threshold on the likelihood ratio statistic used in conjunction with a threshold on the size of a cluster allows the selection of a relevant set of genes. However, even this reduced set of genes will usually be too large for a normal mixture model to be fitted directly to the tissues, and so the use of mixtures of factor analyzers is exploited to reduce effectively the dimension of the feature space of genes. Results: The usefulness of the EMMIX-GENE approach for the clustering of tissue samples is demonstrated on two well-known data sets on colon and leukaemia tissues. For both data sets, relevant subsets of the genes are able to be selected that reveal interesting clusterings of the tissues that are either consistent with the external classification of the tissues or with background and biological knowledge of these sets.