Robust Mixture Modelling Using the t Distribution

Normal mixture models are being increasingly used to model the distributions of a wide variety of random phenomena and to cluster sets of continuous multivariate data. However, for a set of data containing a group or groups of observations with longer than normal tails or atypical observations, the use of normal components may unduly affect the fit of the mixture model. In this paper, we consider a more robust approach by modelling the data by a mixture of t distributions. The use of the ECM algorithm to fit this t mixture model is described and examples of its use are given in the context of clustering multivariate data in the presence of atypical observations in the form of background noise.

Heterogeneous distribution of isoactins in cultured vascular smooth muscle cells does not reflect segregation of contractile and cytoskeletal domains

We have previously demonstrated that or-smooth muscle (alpha -SM) actin is predominantly distributed in the central region and beta -non-muscle (beta -NM) actin in the periphery of cultured rabbit aortic smooth muscle cells (SMCs). To determine whether this reflects a special form of segregation of contractile and cytoskeletal components in SMCs, this study systematically investigated the distribution relationship of structural proteins using high-resolution confocal laser scanning fluorescent microscopy. Not only isoactins but also smooth muscle myosin heavy chain, alpha -actinin, vinculin, and vimentin were heterogeneously distributed in the cultured SMCs. The predominant distribution of beta -NM actin in the cell periphery was associated with densely distributed vinculin plaques and disrupted or striated myosin and ol-actinin aggregates, which may reflect a process of stress fiber assembly during cell spreading and focal adhesion formation. The high-level labeling of alpha -SM actin in the central portion of stress fibers was related to continuous myosin and punctate alpha -actinin distribution, which may represent the maturation of the fibrillar structures. The findings also suggest that the stress fibers, in which actin and myosin filaments organize into sar-comere-like units with alpha -actinin-rich dense bodies analogous to Z-lines, are the contractile vimentin structures of cultured SMCs that link to the network of vimentin-containing intermediate alpha -actinin filaments through the dense bodies and dense plaques.

Morphological and electrophysiological properties of principal neurons in the rat lateral amygdala in vitro

In this study, we characterize the electrophysiological and morphological properties of spiny principal neurons in the rat lateral amygdala using whole cell recordings in acute brain slices. These neurons exhibited a range of firing properties in response to prolonged current injection. Responses varied from cells that showed full spike frequency adaptation, spiking three to five times, to those that showed no adaptation. The differences in firing patterns were largely explained by the amplitude of the afterhyperpolarization (AHP) that followed spike trains. Cells that showed full spike frequency adaptation had large amplitude slow AHPs, whereas cells that discharged tonically had slow AHPs of much smaller amplitude. During spike trains, all cells showed a similar broadening of their action potentials. Biocytin-filled neurons showed a range of pyramidal-like morphologies, differed in dendritic complexity, had spiny dendrites, and differed in the degree to which they clearly exhibited apical versus basal dendrites. Quantitative analysis revealed no association between cell morphology and firing properties. We conclude that the discharge properties of neurons in the lateral nucleus, in response to somatic current injections, are determined by the differential distribution of ionic conductances rather than through mechanisms that rely on cell morphology.

Increasingly inequitable distribution of general practitioners in Australia, 1986-96

Objective: To document trends in the distribution of general practitioners (GPs) in Australia between 1986 and 1996, adjusted for community need. Methods: Data on the location of GPs, population size and crude mortality in statistical divisions (SD) were obtained from the Australian Bureau of Statistics Census of Population and Housing in 1986 and 1996. From these data, we calculated measures of distribution equality (number of people sharing each GP in each SD) and distribution equity (number of people sharing each GP divided by the crude mortality rate; the Robin Hood Index), and analysed temporal changes in the distribution of GPs. Results: Nationally the number of people sharing each GP fell 11% from 1,038 in 1986 to 921 in 1996. However, in 41 of 57 SDs (72%, p=0.01) the number of people sharing a GP actually increased over this time, and the average Robin Hood Index across SDs fell from 0.943 to 0.783 (p=0.004), indicating increasingly inequitable distribution. Comparing the Robin Hood index values of all SDs ranked in pairs, the value fell in 53 of 57 (93%, p<0.001) paired SDs over the decade. These patterns demonstrate increasing inequity over the decade. The number of people sharing each GP was consistently and substantially lower in the capital city SDs and the Robin Hood Index values were consistently and substantially higher (overserved) compared with country SDs. Conclusions: Despite there being more GPs per capita in Australia, their distribution became increasingly unequal and inequitable between 1986 and 1996, such that rural and remote areas became increasingly poorly served.

A landscape-scale test of the predictive ability of a spatially explicit model for population viability analysis

1. Although population viability analysis (PVA) is widely employed, forecasts from PVA models are rarely tested. This study in a fragmented forest in southern Australia contrasted field data on patch occupancy and abundance for the arboreal marsupial greater glider Petauroides volans with predictions from a generic spatially explicit PVA model. This work represents one of the first landscape-scale tests of its type. 2. Initially we contrasted field data from a set of eucalypt forest patches totalling 437 ha with a naive null model in which forecasts of patch occupancy were made, assuming no fragmentation effects and based simply on remnant area and measured densities derived from nearby unfragmented forest. The naive null model predicted an average total of approximately 170 greater gliders, considerably greater than the true count (n = 81). 3. Congruence was examined between field data and predictions from PVA under several metapopulation modelling scenarios. The metapopulation models performed better than the naive null model. Logistic regression showed highly significant positive relationships between predicted and actual patch occupancy for the four scenarios (P = 0.001-0.006). When the model-derived probability of patch occupancy was high (0.50-0.75, 0.75-1.00), there was greater congruence between actual patch occupancy and the predicted probability of occupancy. 4. For many patches, probability distribution functions indicated that model predictions for animal abundance in a given patch were not outside those expected by chance. However, for some patches the model either substantially over-predicted or under-predicted actual abundance. Some important processes, such as inter-patch dispersal, that influence the distribution and abundance of the greater glider may not have been adequately modelled. 5. Additional landscape-scale tests of PVA models, on a wider range of species, are required to assess further predictions made using these tools. This will help determine those taxa for which predictions are and are not accurate and give insights for improving models for applied conservation management.

Quasars from a complete spectroscopic survey

With the advent of multi-fibre spectrographs such as the 'Two-Degree Field' (2dF) instrument at the Angle-Australian Telescope, quasar surveys that are free of any preselection of candidates and any biases this implies have become possible for the first time. The first of these is that which is being undertaken as part of the Fornax Spectroscopic Survey, a survey of the area around the Fornax Cluster of galaxies, and aims to obtain the spectra of all objects in the magnitude range 16.5 < b(j) < 19.7. To date, 3679 objects in the central pi -deg(2) area have been successfully identified from their spectral characteristics. Of these, 71 are found to be quasars, 61 with redshifts 0.3 < z < 2.2 and 10 with redshifts z > 2.2. Using this complete quasar sample, a new determination of quasar number counts is made, enabling an independent check of existing quasars surveys. Cumulative counts per square degree at a magnitude limit of b(j) < 19.5 are found to be 11.5 +/- 2.2 for 0.3 < z < 2.2, 2.22 +/- 0.93 for z > 2.2 and 13.7 +/- 3.1 for z > 0.3. Given the likely detection of extra quasars in the Fornax survey, we make a more detailed examination of existing quasar selection techniques. First, looking at the use of a stellar criterion, four of the 71 quasars are 'non-stellar' on the basis of the automated plate measuring facility (APM) b(j) classification, however inspection shows all are consistent with stellar, but misclassified due to image confusion. Examining the ultraviolet excess and multicolour selection techniques, for the selection criteria investigated, ultraviolet excess would find 69 +/- 6 per cent of our 0.3 < z < 2.2 quasars and only 50(-18)(+14), per cent of our z > 2.2 quasars, while the completeness level for multicolour selection is found to be 90(-4)(+3) per cent for 0.3 < z < 2.2 quasars and 80(-12)(+14) per cent for z > 2.2 quasars. The extra quasars detected by our all-object survey thus have unusually red star-like colours, and this appears to be a result of the continuum shape rather than any emission features. An intrinsic dust extinction model may, at least partly, account for the red colours.

The evolution and star formation of dwarf galaxies in the Fornax Cluster

We present the results of a spectroscopic survey of 675 bright (16.5 < b(J) < 18) galaxies in a 6 degrees field centred on the Fornax cluster with the FLAIR-II spectrograph on the UK Schmidt Telescope. Three galaxy samples were observed: compact galaxies to search for new blue compact dwarfs, candidate M 32-like compact dwarf ellipticals, and a subset of the brightest known cluster members in order to study the cluster dynamics. We measured redshifts for 516 galaxies, of which 108 were members of the Fornax Cluster. Defining dwarf galaxies to be those with b(J) greater than or equal to 15 (M-B greater than or equal to - 16.5), there are a total of 62 dwarf cluster galaxies in our sample. Nine of these are new cluster members previously misidentified as background galaxies. The cluster dynamics show that the dwarf galaxies are still falling into the cluster whereas the giants are virialized. We classified the observed galaxies as late-type if we detected H alpha emission at an equivalent width greater than 1 Angstrom. The spectra were obtained through fixed apertures, so they reflect activity in the galaxy cores, but this does not significantly bias the classifications of the compact dwarfs in our sample. The new classifications reveal a higher rate of star formation among the dwarf galaxies than suggested by morphological classification: 35 per cent have significant H alpha emission indicative of star formations but only 19 per cent were morphologically classified as late-types. The star-forming dwarf galaxies span the full range of physical sizes and we find no evidence in our data for a distinct class of star-forming blue compact dwarf (BCD) galaxy. The distribution of scale sizes is consistent with evolutionary processes which transform late-type dwarfs to early-type dwarfs. The fraction of dwarfs with active star formation drops rapidly towards the cluster centre: this is the usual density-morphology relation confirmed here for dwarf galaxies. The star-forming dwarfs are concentrated in the outer regions of the cluster, the most extreme in an infalling subcluster. We estimate gas depletion time-scales for five dwarfs with detected Hi emission: these are long (of order 10(10) yr), indicating that an active gas removal process must be involved if they are transformed into gas-poor dwarfs as they fall further into the cluster. Finally, in agreement with our previous results, we find no compact dwarf elliptical (M 32-like) galaxies in the Fornax Cluster.

The spectra of red quasars

We measure the spectral properties of a representative sub-sample of 187 quasars, drawn from the Parkes Half-Jansky, Flat-radio-spectrum Sample (PHFS). Quasars with a wide range of rest-frame optical/UV continuum slopes are included in the analysis: their colours range over 2 < B-K < 7. We present composite spectra of red and blue sub-samples of the PHFS quasars. and tabulate their emission line properties. The median Hbeta and [0 111] emission line equivalent widths of the red quasar sub-sample are a factor of ten weaker than those of the blue quasar sub-sample. No significant differences are seen between the equivalent width distributions of the C IV, C III] and Mg 11 lines. Both the colours and the emission line equivalent widths of the red quasars can be explained by the addition of a featureless red synchrotron continuum component to an otherwise normal blue quasar spectrum. The red synchrotron component must have a spectrum at least as red as a power-law of the form F-nu proportional to nu(-2.8). The relative strengths of the blue and red components span two orders of magnitude at rest-frame 500 nm. The blue component is weaker relative to the red component in low optical luminosity sources. This suggests that the fraction of accretion energy going into optical emission from the jet is greater in low luminosity quasars. This correlation between colour and luminosity may be of use in cosmological distance scale work. This synchrotron model does not, however, fit similar to10% of the quasars, which have both red colours and high equivalent width emission lines. We hypothesise that these red, strong-lined quasars have intrinsically weak Big Blue Bumps. There is no discontinuity in spectral properties between the BL Lac objects in our sample and the other quasars. BL Lac objects appear to be the red, low equivalent width tail of a continuous distribution. The synchrotron emission component only dominates the spectrum at longer wavelengths, so existing BL Lac surveys will be biased against high redshift objects. This will affect measurements of BL Lac evolution. The blue PHFS quasars have significantly higher equivalent width C IV, Hbeta and [0 111] emission than a matched sample of optically selected QSOs.

Spatial and temporal variation in distribution of mangroves in Moreton Bay, subtropical Australia: a comparison of pattern metrics and change detection analyses based on aerial photographs

An assessment of the changes in the distribution and extent of mangroves within Moreton Bay, southeast Queensland, Australia, was carried out. Two assessment methods were evaluated: spatial and temporal pattern metrics analysis, and change detection analysis. Currently, about 15,000 ha of mangroves are present in Moreton Bay. These mangroves are important ecosystems, but are subject to disturbance from a number of sources. Over the past 25 years, there has been a loss of more than 3800 ha, as a result of natural losses and mangrove clearing (e.g. for urban and industrial development, agriculture and aquaculture). However, areas of new mangroves have become established over the same time period, offsetting these losses to create a net loss of about 200 ha. These new mangroves have mainly appeared in the southern bay region and the bay islands, particularly on the landward edge of existing mangroves. In addition, spatial patterns and species composition of mangrove patches have changed. The pattern metrics analysis provided an overview of mangrove distribution and change in the form of single metric values, while the change detection analysis gave a more detailed and spatially explicit description of change. An analysis of the effects of spatial scales on the pattern metrics indicated that they were relatively insensitive to scale at spatial resolutions less than 50 m, but that most metrics became sensitive at coarser resolutions, a finding which has implications for mapping of mangroves based on remotely sensed data. (C) 2003 Elsevier Science B.V. All rights reserved.

A technological and organisational explanation for the size distribution of firms

This paper combines insights from the literature on the economics of organisation with traditional models of market structure to construct a theory of equilibrium firm size heterogeneity under the assumption of a homogenous product industry. It is possible that configurations consisting entirely of small firms (run by entrepreneurs with limited attention) and with larger firms (using managerial techniques to substitute away these limits to allow increasing returns technologies to become profitable) can arise in equilibrium. However, there also exist equilibrium configurations with the co-existence of large and small firms. The efficiency properties of these respective equilibria are discussed. Finally, the implications of an expanding market size are considered.

Micromorphological study of slumping in a hardsetting seedbed under various wetting conditions

Slumping of hardsetting seedbeds upon wetting has not been extensively studied despite the likelihood that it determines the physical properties after drying. Slumping results from processes similar to those involved in crusting except that overburden pressure can dominate rather than rainfall kinetic energy. Only a few studies have dealt with the morphological description of slumping. To simulate different climatic and management conditions, repacked seedbeds of a hardsetting sandy-loam soil were subjected to a range of wetting conditions, e.g. capillary rise, immersion, and rainfall simulation. Slumping processes were characterized using qualitative and quantitative micromorphological observations of polished blocks and thin sections from resin-impregnated samples. A morphogenetical framework was proposed to help description of the complex associations of processes which can lead to structural collapse (crusting and slumping) on wetting. Three main stages were considered, i.e. aggregate disruption or abrasion, relocation of the released material, and compaction. In the hardsetting material studied here, structural collapse under slow wetting occurred at the bottom of cores due to aggregate coalescence under overburden pressure. Coalescence required aggregate cohesion being reduced by microcracking; therefore, it differed from the coalescence previously described in unstable silty loam soils where microcracking was not necessary for aggregates to coalesce. Macroporosity decreased most strongly under fast wetting due to physical dispersion and aggregate breakdown. Under simulated rainfall, compaction by raindrops could not be distinguish from aggregate breakdown. The role of overburden pressure and of rainfall kinetic energy remains to be stated; new data are required including measurement of total porosity in the initial, wet, and dry states. (C) 2003 Elsevier B.V. All rights reserved.

The secondary vascular system of Actinopterygii: interspecific variation in origins and investment

Vascular casts of 3 species of Chondrichthyes, 1 of Dipnoi, 1 of Chondrostei and 14 species of the Teleostei were examined by light and scanning electron microscopy in order to give a qualitative and quantitative analysis of interarterial anastomoses (iaas) that indicate the presence (or absence) of a secondary vascular system (SVS). Anastomoses were found to originate from a variety of different primary blood vessels, many of which have not been previously identified as giving rise to secondary vessels. Segmental arteries derived from the dorsal aorta and supplying body musculature were major sites of origin of the SVS, although there was considerable variation in where, in the hierarchy of arterial branching, the anastomoses occurred. The degree of investment in a SVS was species specific, with more active species having a higher degree of secondary vascularisation. This difference was quantified using an absolute count of iaas between Anguilla reinhardtii and Trachinotus baillonii. A range of general features of the SVS is also described. No evidence of iaas was found on the coeliac, mesenteric or renal circulation in any species. Evidence of iaas was lacking in the dipnoan and chondrichthyan species examined, suggesting that a SVS is restricted to Actinopterygii. The presence and distribution of a SVS does not appear to be exclusively linked to phylogenetic position, but rather to the physiological adaptation of the species.