The design of hazard risk assessment matrices for ranking occupational health risks and their application in mining and minerals processing

Two hazard risk assessment matrices for the ranking of occupational health risks are described. The qualitative matrix uses qualitative measures of probability and consequence to determine risk assessment codes for hazard-disease combinations. A walk-through survey of an underground metalliferous mine and concentrator is used to demonstrate how the qualitative matrix can be applied to determine priorities for the control of occupational health hazards. The semi-quantitative matrix uses attributable risk as a quantitative measure of probability and uses qualitative measures of consequence. A practical application of this matrix is the determination of occupational health priorities using existing epidemiological studies. Calculated attributable risks from epidemiological studies of hazard-disease combinations in mining and minerals processing are used as examples. These historic response data do not reflect the risks associated with current exposures. A method using current exposure data, known exposure-response relationships and the semi-quantitative matrix is proposed for more accurate and current risk rankings.

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.

Path integrals for continuous-time Markov chains

This note presents a method of evaluating the distribution of a path integral for Markov chains on a countable state space.

A zero-augmented gamma mixed model for longitudinal data with many zeros

In many occupational safety interventions, the objective is to reduce the injury incidence as well as the mean claims cost once injury has occurred. The claims cost data within a period typically contain a large proportion of zero observations (no claim). The distribution thus comprises a point mass at 0 mixed with a non-degenerate parametric component. Essentially, the likelihood function can be factorized into two orthogonal components. These two components relate respectively to the effect of covariates on the incidence of claims and the magnitude of claims, given that claims are made. Furthermore, the longitudinal nature of the intervention inherently imposes some correlation among the observations. This paper introduces a zero-augmented gamma random effects model for analysing longitudinal data with many zeros. Adopting the generalized linear mixed model (GLMM) approach reduces the original problem to the fitting of two independent GLMMs. The method is applied to evaluate the effectiveness of a workplace risk assessment teams program, trialled within the cleaning services of a Western Australian public hospital.

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.

A simulated annealing algorithm for finding consensus sequences

Motivation: A consensus sequence for a family of related sequences is, as the name suggests, a sequence that captures the features common to most members of the family. Consensus sequences are important in various DNA sequencing applications and are a convenient way to characterize a family of molecules. Results: This paper describes a new algorithm for finding a consensus sequence, using the popular optimization method known as simulated annealing. Unlike the conventional approach of finding a consensus sequence by first forming a multiple sequence alignment, this algorithm searches for a sequence that minimises the sum of pairwise distances to each of the input sequences. The resulting consensus sequence can then be used to induce a multiple sequence alignment. The time required by the algorithm scales linearly with the number of input sequences and quadratically with the length of the consensus sequence. We present results demonstrating the high quality of the consensus sequences and alignments produced by the new algorithm. For comparison, we also present similar results obtained using ClustalW. The new algorithm outperforms ClustalW in many cases.

Efficient simulation of a Tandem Jackson Network

The two-node tandem Jackson network serves as a convenient reference model for the analysis and testing of different methodologies and techniques in rare event simulation. In this paper we consider a new approach to efficiently estimate the probability that the content of the second buffer exceeds some high level L before it becomes empty, starting from a given state. The approach is based on a Markov additive process representation of the buffer processes, leading to an exponential change of measure to be used in an importance sampling procedure. Unlike changes of measures proposed and studied in recent literature, the one derived here is a function of the content of the first buffer. We prove that when the first buffer is finite, this method yields asymptotically efficient simulation for any set of arrival and service rates. In fact, the relative error is bounded independent of the level L; a new result which is not established for any other known method. When the first buffer is infinite, we propose a natural extension of the exponential change of measure for the finite buffer case. In this case, the relative error is shown to be bounded (independent of L) only when the second server is the bottleneck; a result which is known to hold for some other methods derived through large deviations analysis. When the first server is the bottleneck, experimental results using our method seem to suggest that the relative error is bounded linearly in L.

Laying date and laying sequence influence the sex ratio of crimson rosella broods

We examine the patterns of sex allocation in crimson rosellas Platycercus elegans, a socially monogamous Australian parrot. Overall, 41.8% of nestlings were male, a significant female bias. However underlying this population-level bias were non-random patterns of sex allocation within broods. Broods produced early in the season were female-biased, but the proportion of males in a brood increased as the breeding season progressed. Female rosellas may obtain greater fitness benefits from early-fledging daughters than sons because daughters can breed as 1-year-olds whereas sons do not breed until they are at least 2 years old. Laying date and laying sequence also interacted to influence the sex ratio of eggs. The sex of early-laid eggs strongly followed the brood level pattern, whereas the sex of middle- and late-laid eggs did not change significantly as the season progressed. Nevertheless, late-laid eggs were very unlikely to be male at the end of the season. We argue these differing seasonal patterns reflect the relative costs and benefits to producing early-hatched males and females at different times of the season. Female rosellas appear to maximise the probability that daughters are able to breed early but to minimise competitive asymmetries within the brood. In particular, late-hatched male chicks are disadvantaged if their oldest sibling is male, explaining the dearth of broods containing late-hatched males at the end of the breeding season.

Acute symptoms, not rectally administered sucralfate, predict for late radiation proctitis: Longer term follow-up of a phase III trial - trans-Tasman radiation oncology group

Image : To assess the potential for sucralfate administered rectally to reduce the risk of late rectal morbidity in patients undergoing nonconformal radiotherapy (RT) for carcinoma of the prostate and to study the variables potentially contributing to late rectal morbidity and particularly to explore the relationship between acute and late toxicity. Image : Eighty-six patients with localized prostate carcinoma were randomized in a double-blind, placebo-controlled study to a daily enema of 3 g of sucralfate in a 15-mL suspension or the same suspension without sucralfate. The enema began the first day of RT and was continued for 2 weeks after treatment completion. The primary end point of the study was acute Radiation Therapy Oncology Group (RTOG)/European Organization for Research and Treatment of Cancer (EORTC) toxicity; however, the patients were followed for an additional 5 years on a 6-month basis. The evaluation included late RTOG/EORTC toxicity and a patient self-assessment questionnaire. Image : With a median follow-up of 5 years, the Kaplan-Meier probability of late Grade 2 RTOG/EORTC toxicity was 12% (95% confidence interval [CI] 2–22%) for placebo and 5% (95% CI 0–12%) for sucralfate (p = 0.26). The probability of late rectal bleeding was 59% (95% CI 45–73%) for placebo and 54% (95% CI 40–68%) for sucralfate. No statistically significant difference was found between the treatment arms for the peak incidence of any of the other patient self-assessment variables. Cox proportional hazards modeling indicated acute RTOG/EORTC toxicity of Grade 2 or greater was associated with a hazard ratio of 2.74 (95% CI 1.31–5.73) for the development of late toxicity of Grade 1 or greater. Substituting the patient self-assessment variables for acute RTOG/EORTC toxicity revealed that rectal pain of a moderate or severe grade during RT was the best predictor of the subsequent development of late toxicity, with a hazard ratio of 3.44 (95% CI 1.68–7). Image : The results of this study do not support the use of sucralfate administered rectally as a method for reducing the late toxicity of nonconformal RT for prostate cancer. There appears to be an association between the development of acute and subsequent late toxicity, although the nature of this association remains to be determined

Two-dimensional balanced sampling plans excluding contiguous units

A balanced sampling plan excluding contiguous units (or BSEC for short) was first introduced by Hedayat, Rao and Stufken in 1988. These designs can be used for survey sampling when the units are arranged in one-dimensional ordering and the contiguous units in this ordering provide similar information. In this paper, we generalize the concept of a BSEC to the two-dimensional situation and give constructions of two-dimensional BSECs with block size 3. The existence problem is completely solved in the case where lambda = 1.

Do female rainbowfish (Melanotaenia spp.) prefer to shoal with familiar individuals under predation pressure?

Shoaling with familiar individuals may have many benefits including enhanced escape responses or increased foraging efficiency. This study describes the results of two complimentary experiments. The first utilised a simple binary choice experiment to determine if rainbowfish (Melanotaenia spp.) preferred to shoal with familiar individuals or with strangers. The second experiment used a free range situation where familiar and unfamiliar individuals were free to intermingle and were then exposed to a predator threat. Like many other small species of fish, rainbowfish were capable of identifying and distinguishing between individuals and choose to preferentially associate with familiar individuals as opposed to strangers. Contrary to expectations. however. rainbowrish did not significantly increase their preference for familiar individuals in the presence of a stationary predator model. Griffiths [J Fish Biol (1997) 51:489-4951 conducted similar studies under semi-natural conditions examining, the shoaling preferences of European minnows and showed similar results. Both the current study and that of Griffiths were conducted using predator wary populations of fish. It is suggested that, in predator sympatric populations, the benefits of shoaling with familiar individuals are such that it always pays to stay close to familiar individuals even when the probability If predator attack is remote.

Model specification tests in nonparametric stochastic regression models

In this paper, we consider testing for additivity in a class of nonparametric stochastic regression models. Two test statistics are constructed and their asymptotic distributions are established. We also conduct a small sample study for one of the test statistics through a simulated example. (C) 2002 Elsevier Science (USA).

Cross-national comparison of internal migration: issues and measures

Our objectives are to identify the issues that researchers encounter when measuring internal migration in different countries and to propose key indicators that analysts can use to compare internal migration at the 'national' level. We establish the benefits to be gained by a rigorous approach to cross-national comparisons of internal migration and discuss issues that affect such comparisons. We then distinguish four dimensions of internal migration on which countries can be compared and, for each dimension, identify a series of summary measures. We illustrate the issues and measures proposed by comparing migration in Australia and Great Britain.

Estimating and evaluating the statistics of gapped local-alignment scores

We present a novel maximum-likelihood-based algorithm for estimating the distribution of alignment scores from the scores of unrelated sequences in a database search. Using a new method for measuring the accuracy of p-values, we show that our maximum-likelihood-based algorithm is more accurate than existing regression-based and lookup table methods. We explore a more sophisticated way of modeling and estimating the score distributions (using a two-component mixture model and expectation maximization), but conclude that this does not improve significantly over simply ignoring scores with small E-values during estimation. Finally, we measure the classification accuracy of p-values estimated in different ways and observe that inaccurate p-values can, somewhat paradoxically, lead to higher classification accuracy. We explain this paradox and argue that statistical accuracy, not classification accuracy, should be the primary criterion in comparisons of similarity search methods that return p-values that adjust for target sequence length.

Determinants of tissue Doppler measures of regional diastolic function during dobutamine stress echocardiography

Background Diastolic dysfunction induced by ischemia may alter transmitral blood flow, but this reflects global ventricular function, and pseudonormalization may occur with increased preload. Tissue Doppler may assess regional diastolic function and is relatively load-independent, but limited data exist regarding its application to stress testing. We sought to examine the stress response of regional diastolic parameters to dobutomine echocardiography (DbE). Methods Sixty-three patients underwent study with DbE: 20 with low probability of coronary artery disease (CAD) and 43 with CAD who underwent angiography. A standard DbE protocol was used, and segments were categorized as ischemic, scar, or normal. Color tissue Doppler was acquired at baseline and peak stress, and waveforms in the basal and mid segments were used to measure early filling (Em), late filling (Am), and E deceleration time. Significant CAD was defined by stenoses >50% vessel diameter. Results Diastolic parameters had limited feasibility because of merging of Em and Am waves at high heart rates and limited reproducibility. Nonetheless, compared with normal segments, segments subtended with significant stenoses showed a lower Em velocity at rest (6.2 +/- 2.6 cm/s vs 4.8 +/- 2.2 cm/s, P < .0001) and peak (7.5 +/- 4.2 cm/s vs 5.1 +/- 3.6 cm/s, P < .0001), Abnormal segments also showed a shorter E deceleration time (51 +/- 27 ms vs 41 +/- 27 ms, P = .0001) at base and peak. No changes were documented in Am. The same pattern was seen with segments identified as ischemic with wall motion score. However, in the absence of ischemia, segments of patients with left ventricular hypertrophy showed a lower Em velocity, with blunted Em responses to stress. Conclusion Regional diastolic function is sensitive to ischemia. However, a number of practical limitations limit the applicability of diastolic parameters for the quantification of stress echocardiography.