Modelling life history strategies with capture-recapture data: Evolutionary demography of the water skink Eulamprus tympanum

Matrix population models, elasticity analysis and loop analysis can potentially provide powerful techniques for the analysis of life histories. Data from a capture-recapture study on a population of southern highland water skinks (Eulamprus tympanum) were used to construct a matrix population model. Errors in elasticities were calculated by using the parametric bootstrap technique. Elasticity and loop analyses were then conducted to identify the life history stages most important to fitness. The same techniques were used to investigate the relative importance of fast versus slow growth, and rapid versus delayed reproduction. Mature water skinks were long-lived, but there was high immature mortality. The most sensitive life history stage was the subadult stage. It is suggested that life history evolution in E. tympanum may be strongly affected by predation, particularly by birds. Because our population declined over the study, slow growth and delayed reproduction were the optimal life history strategies over this period. Although the techniques of evolutionary demography provide a powerful approach for the analysis of life histories, there are formidable logistical obstacles in gathering enough high-quality data for robust estimates of the critical parameters.

Design options for sample-based mortality surveillance

Background Reliable information on causes of death is a fundamental component of health development strategies, yet globally only about one-third of countries have access to such information. For countries currently without adequate mortality reporting systems there are useful models other than resource-intensive population-wide medical certification. Sample-based mortality surveillance is one such approach. This paper provides methods for addressing appropriate sample size considerations in relation to mortality surveillance, with particular reference to situations in which prior information on mortality is lacking. Methods The feasibility of model-based approaches for predicting the expected mortality structure and cause composition is demonstrated for populations in which only limited empirical data is available. An algorithm approach is then provided to derive the minimum person-years of observation needed to generate robust estimates for the rarest cause of interest in three hypothetical populations, each representing different levels of health development. Results Modelled life expectancies at birth and cause of death structures were within expected ranges based on published estimates for countries at comparable levels of health development. Total person-years of observation required in each population could be more than halved by limiting the set of age, sex, and cause groups regarded as 'of interest'. Discussion The methods proposed are consistent with the philosophy of establishing priorities across broad clusters of causes for which the public health response implications are similar. The examples provided illustrate the options available when considering the design of mortality surveillance for population health monitoring purposes.

Absolute and relative abundance estimates of Australian east cost humpback whales (Megaptera novaeangliae)

The humpback whales that migrate along the east coast of Australia were hunted to near-extinction in the 1950s and early 1960s. Two independent series of land-based surveys conducted over the last 25 years during the whales’ northward migration along the Australian coastline have demonstrated a rapid increase in the size of the population. In 2004 we conducted a survey of the migratory population as a continuation of these series of surveys. Two methods of data analysis were used in line with the previous surveys, both for calculation of absolute and relative abundance. We consider the best estimates for 2004 to be 7,090 ± 660 (95% CI) whales with an annual rate of increase of 10.6 ± 0.5% (95% CI) for 1987 – 2004. The rate of increase agrees with those previously obtained for this population and demonstrates the continuation of a strong post-exploitation recovery. While there are still some uncertainties concerning the absolute abundance estimate and structure of this population, the rate of annual increase should be independent of these and highly robust.

Measuring the accumulated hazards of smoking: global and regional estimates for 2000

Objective: Current prevalence of smoking, even where data are available, is a poor proxy for cumulative hazards of smoking, which depend on several factors including the age at which smoking began, duration of smoking, number of cigarettes smoked per day, degree of inhalation, and cigarette characteristics such as tar and nicotine content or filter type. Methods: We extended the Peto-Lopez smoking impact ratio method to estimate accumulated hazards of smoking for different regions of the world. Lung cancer mortality data were obtained from the Global Burden of Disease mortality database. The American Cancer Society Cancer Prevention Study, phase 11 (CPS-II) with follow up for the years 1982 to 1988 was the reference population. For the global application of the method, never-smoker lung cancer mortality rates were chosen based on the estimated use of coal for household energy in each region. Results: Men in industrialised countries of Europe, North America, and the Western Pacific had the largest accumulated hazards of smoking. Young and middle age males in many regions of the developing world also had large smoking risks. The accumulated hazards of smoking for women were highest in North America followed by Europe. Conclusions: In the absence of detailed data on smoking prevalence and history, lung cancer mortality provides a robust indicator of the accumulated hazards of smoking. These hazards in developing countries are currently more concentrated among young and middle aged males.

The return to schooling: Estimates from a sample of young Australian twins

This study uses a sample of young Australian twins to examine whether the findings reported in [Ashenfelter, Orley and Krueger, Alan, (1994). 'Estimates of the Economic Return to Schooling from a New Sample of Twins', American Economic Review, Vol. 84, No. 5, pp.1157-73] and [Miller, P.W., Mulvey, C and Martin, N., (1994). 'What Do Twins Studies Tell Us About the Economic Returns to Education?: A Comparison of Australian and US Findings', Western Australian Labour Market Research Centre Discussion Paper 94/4] are robust to choice of sample and dependent variable. The economic return to schooling in Australia is between 5 and 7 percent when account is taken of genetic and family effects using either fixed-effects models or the selection effects model of Ashenfelter and Krueger. Given the similarity of the findings in this and in related studies, it would appear that the models applied by [Ashenfelter, Orley and Krueger, Alan, (1994). 'Estimates of the Economic Return to Schooling from a New Sample of Twins', American Economic Review, Vol. 84, No. 5, pp. 1157-73] are robust. Moreover, viewing the OLS and IV estimators as lower and upper bounds in the manner of [Black, Dan A., Berger, Mark C., and Scott, Frank C., (2000). 'Bounding Parameter Estimates with Nonclassical Measurement Error', Journal of the American Statistical Association, Vol. 95, No.451, pp.739-748], it is shown that the bounds on the return to schooling in Australia are much tighter than in [Ashenfelter, Orley and Krueger, Alan, (1994). 'Estimates of the Economic Return to Schooling from a New Sample of Twins', American Economic Review, Vol. 84, No. 5, pp. 1157-73], and the return is bounded at a much lower level than in the US. (c) 2005 Elsevier B.V. All rights reserved.

How many glomerular profiles must be measured to obtain reliable estimates of mean glomerular areas in human renal biopsies?

The objective of this study was to investigate the number of glomerular profiles that are required for accurate estimates of mean profile area in a renal biopsy series. Slides from 384 renal biopsies from one center were reviewed. They contained a median of seven glomerular profiles or of four profiles without sclerosis. Profile areas were measured using stereologic point counting. The true individual mean for each biopsy was calculated and the true population mean for groups of biopsies derived. Individual and population random sample means then were calculated from a random sampling of profiles in each biopsy and were compared with true means for the same biopsies. The effect on the true population means of the entire group of biopsies was also assessed, as the minimum number of glomerular profiles that were required for inclusion was changed. In a single biopsy, random sampling of >= 10 profiles without exclusions and of eight profiles or more without sclerosis reliably estimated the true mean areas. In a group of 30 biopsies, random sampling of five or more glomeruli per biopsy reliably estimated the true population mean. In the aggregate series, inclusion of all 384 biopsies produced the most robust true population mean; the reliability of the estimates decreased as the numbers of eligible biopsies diminished with increasing requisite minimum numbers of profiles per biopsy. We conclude that, while >= 10 profiles might be needed for reliable area estimates in a single biopsy, far fewer profiles per biopsy can suffice when groups of biopsies are studied. In analyses of groups of biopsies, all available biopsies should be used without consideration of the number of glomerular profiles in each. Stipulation of a specific minimum number of glomeruli in each biopsy for inclusion reduces the power of analyses because fewer biopsies are available for evaluation.

Remotely sensed estimates of vegetation structural characteristics in restored wetlands, Southern California

Traditional field sampling approaches for ecological studies of restored habitat can only cover small areas in detail, con be time consuming, and are often invasive and destructive. Spatially extensive and non-invasive remotely sensed data can make field sampling more focused and efficient. The objective of this work was to investigate the feasibility and accuracy of hand-held and airborne remotely sensed data to estimate vegetation structural parameters for an indicator plant species in a restored wetland. High spatial resolution, digital, multispectral camera images were captured from an aircraft over Sweetwater Marsh (San Diego County, California) during each growing season between 1992-1996. Field data were collected concurrently, which included plant heights, proportional ground cover and canopy architecture type, and spectral radiometer measurements. Spartina foliosa (Pacific cordgrass) is the indicator species for the restoration monitoring. A conceptual model summarizing the controls on the spectral reflectance properties of Pacific cordgrass was established. Empirical models were developed relating the stem length, density, and canopy architecture of cordgrass to normalized-difference-vegetation-index values. The most promising results were obtained from empirical estimates of total ground cover using image data that had been stratified into high, middle, and low marsh zones. As part of on-going restoration monitoring activities, this model is being used to provide maps of estimated vegetation cover.

Hemichordates and deuterostome evolution: robust molecular phylogenetic support for a hemichordate plus echinoderm clade

Hemichordates were traditionally allied to the chordates, but recent molecular analyses have suggested that hemichordates are a sister group to the echinoderms, a relationship that has important consequences for the interpretation of the evolution of deuterostome body plans. However, the molecular phylogenetic analyses to date have not provided robust support for the hemichordate + echinoderm clade. We use a maximum likelihood framework, including the parametric bootstrap, to reanalyze DNA data from complete mitochondrial genomes and nuclear 18S rRNA. This approach provides the first statistically significant support for the hemichordate + echinoderm clade from molecular data. This grouping implies that the ancestral deuterostome had features that included an adult with a pharynx and a dorsal nerve cord and an indirectly developing dipleurula-like larva.

The effect of arousals during sleep onset on estimates of sleep onset latency

It is well established that insomniacs overestimate sleep-onset latency. Furthermore, there is evidence that brief arousals from sleep may occur more frequently in insomnia. This study examined the hypothesis that brief arousals from sleep influence the perception of sleep-onset latency. An average of four sleep onsets was obtained from each of 20 normal subjects on each of two nonconsecutive, counterbalanced, experimental nights. The experimental nights consisted of a control night (control condition) and a condition in which a moderate respiratory load was applied to increase the frequency of microarousals during sleep onset (mask condition). Subjective estimation of sleep-onset latency and indices of sleep quality were assessed by self-report inventory. Objective measures of sleep-onset latency and microarousals were assessed using polysomnography. Results showed that sleep-onset latency estimates were longer in the mask condition than in the control condition, an effect not reflected in objective sleep-stage scoring of sleep-onset latency. Furthermore, an increase in the frequency of brief arousals from sleep was detected in the mask condition, and this is a possible source for the sleep-onset latency increase perceived by the subjects. Findings are consistent with the concept of a physiological basis for sleep misperception in insomnia.

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.

The HI Parkes All Sky Survey: southern observations, calibration and robust imaging

The acquisition of HI Parkes All Shy Survey (HIPASS) southern sky data commenced at the Australia Telescope National Facility's Parkes 64-m telescope in 1997 February, and was completed in 2000 March. HIPASS is the deepest HI survey yet of the sky south of declination +2 degrees, and is sensitive to emission out to 170 h(75)(-1) Mpc. The characteristic root mean square noise in the survey images is 13.3 mJy. This paper describes the survey observations, which comprise 23 020 eight-degree scans of 9-min duration, and details the techniques used to calibrate and image the data. The processing algorithms are successfully designed to be statistically robust to the presence of interference signals, and are particular to imaging point (or nearly point) sources. Specifically, a major improvement in image quality is obtained by designing a median-gridding algorithm which uses the median estimator in place of the mean estimator.

Dietary trends: Estimates from food supply and survey data

Objective: To illustrate methodological issues involved in estimating dietary trends in populations using data obtained from various sources in Australia in the 1980s and 1990s. Methods: Estimates of absolute and relative change in consumption of selected food items were calculated using national data published annually on the national food supply for 1982-83 to 1992-93 and responses to food frequency questions in two population based risk factor surveys in 1983 and 1994 in the Hunter Region of New South Wales, Australia. The validity of estimated food quantities obtained from these inexpensive sources at the beginning of the period was assessed by comparison with data from a national dietary survey conducted in 1983 using 24 h recall. Results: Trend estimates from the food supply data and risk factor survey data were in good agreement for increases in consumption of fresh fruit, vegetables and breakfast food and decreases in butter, margarine, sugar and alcohol. Estimates for trends in milk, eggs and bread consumption, however, were inconsistent. Conclusions: Both data sources can be used for monitoring progress towards national nutrition goals based on selected food items provided that some limitations are recognized. While data collection methods should be consistent over time they also need to allow for changes in the food supply (for example the introduction of new varieties such as low-fat dairy products). From time to time the trends derived from these inexpensive data sources should be compared with data derived from more detailed and quantitative estimates of dietary intake.

Two estimates concerning asymptotics of the minimizations of a Ginzburg-Landau functional

We prove two asymptotical estimates for minimizers of a Ginzburg-Landau functional of the form integral(Omega) [1/2 \del u\(2) + 1/4 epsilon(2) (1 - \u\(2))(2) W (x)] dx.