Typing mineral deposits using their grades and tonnages in an artificial neural network

A test of the ability of a probabilistic neural network to classify deposits into types on the basis of deposit tonnage and average Cu, Mo, Ag, Au, Zn, and Pb grades is conducted. The purpose is to examine whether this type of system might serve as a basis for integrating geoscience information available in large mineral databases to classify sites by deposit type. Benefits of proper classification of many sites in large regions are relatively rapid identification of terranes permissive for deposit types and recognition of specific sites perhaps worthy of exploring further. Total tonnages and average grades of 1,137 well-explored deposits identified in published grade and tonnage models representing 13 deposit types were used to train and test the network. Tonnages were transformed by logarithms and grades by square roots to reduce effects of skewness. All values were scaled by subtracting the variable's mean and dividing by its standard deviation. Half of the deposits were selected randomly to be used in training the probabilistic neural network and the other half were used for independent testing. Tests were performed with a probabilistic neural network employing a Gaussian kernel and separate sigma weights for each class (type) and each variable (grade or tonnage). Deposit types were selected to challenge the neural network. For many types, tonnages or average grades are significantly different from other types, but individual deposits may plot in the grade and tonnage space of more than one type. Porphyry Cu, porphyry Cu-Au, and porphyry Cu-Mo types have similar tonnages and relatively small differences in grades. Redbed Cu deposits typically have tonnages that could be confused with porphyry Cu deposits, also contain Cu and, in some situations, Ag. Cyprus and kuroko massive sulfide types have about the same tonnages. Cu, Zn, Ag, and Au grades. Polymetallic vein, sedimentary exhalative Zn-Pb, and Zn-Pb skarn types contain many of the same metals. Sediment-hosted Au, Comstock Au-Ag, and low-sulfide Au-quartz vein types are principally Au deposits with differing amounts of Ag. Given the intent to test the neural network under the most difficult conditions, an overall 75% agreement between the experts and the neural network is considered excellent. Among the largestclassification errors are skarn Zn-Pb and Cyprus massive sulfide deposits classed by the neuralnetwork as kuroko massive sulfides—24 and 63% error respectively. Other large errors are the classification of 92% of porphyry Cu-Mo as porphyry Cu deposits. Most of the larger classification errors involve 25 or fewer training deposits, suggesting that some errors might be the result of small sample size. About 91% of the gold deposit types were classed properly and 98% of porphyry Cu deposits were classes as some type of porphyry Cu deposit. An experienced economic geologist would not make many of the classification errors that were made by the neural network because the geologic settings of deposits would be used to reduce errors. In a separate test, the probabilistic neural network correctly classed 93% of 336 deposits in eight deposit types when trained with presence or absence of 58 minerals and six generalized rock types. The overall success rate of the probabilistic neural network when trained on tonnage and average grades would probably be more than 90% with additional information on the presence of a few rock types.

Improving precision and reducing bias in biological surveys: estimating false-negative error rates

The use of presence/absence data in wildlife management and biological surveys is widespread. There is a growing interest in quantifying the sources of error associated with these data. We show that false-negative errors (failure to record a species when in fact it is present) can have a significant impact on statistical estimation of habitat models using simulated data. Then we introduce an extension of logistic modeling, the zero-inflated binomial (ZIB) model that permits the estimation of the rate of false-negative errors and the correction of estimates of the probability of occurrence for false-negative errors by using repeated. visits to the same site. Our simulations show that even relatively low rates of false negatives bias statistical estimates of habitat effects. The method with three repeated visits eliminates the bias, but estimates are relatively imprecise. Six repeated visits improve precision of estimates to levels comparable to that achieved with conventional statistics in the absence of false-negative errors In general, when error rates are less than or equal to50% greater efficiency is gained by adding more sites, whereas when error rates are >50% it is better to increase the number of repeated visits. We highlight the flexibility of the method with three case studies, clearly demonstrating the effect of false-negative errors for a range of commonly used survey methods.

Statistical modeling and projections of lung cancer mortality in 4 industrialized countries

The purpose of this work was to model lung cancer mortality as a function of past exposure to tobacco and to forecast age-sex-specific lung cancer mortality rates. A 3-factor age-period-cohort (APC) model, in which the period variable is replaced by the product of average tar content and adult tobacco consumption per capita, was estimated for the US, UK, Canada and Australia by the maximum likelihood method. Age- and sex-specific tobacco consumption was estimated from historical data on smoking prevalence and total tobacco consumption. Lung cancer mortality was derived from vital registration records. Future tobacco consumption, tar content and the cohort parameter were projected by autoregressive moving average (ARIMA) estimation. The optimal exposure variable was found to be the product of average tar content and adult cigarette consumption per capita, lagged for 2530 years for both males and females in all 4 countries. The coefficient of the product of average tar content and tobacco consumption per capita differs by age and sex. In all models, there was a statistically significant difference in the coefficient of the period variable by sex. In all countries, male age-standardized lung cancer mortality rates peaked in the 1980s and declined thereafter. Female mortality rates are projected to peak in the first decade of this century. The multiplicative models of age, tobacco exposure and cohort fit the observed data between 1950 and 1999 reasonably well, and time-series models yield plausible past trends of relevant variables. Despite a significant reduction in tobacco consumption and average tar content of cigarettes sold over the past few decades, the effect on lung cancer mortality is affected by the time lag between exposure and established disease. As a result, the burden of lung cancer among females is only just reaching, or soon will reach, its peak but has been declining for I to 2 decades in men. Future sex differences in lung cancer mortality are likely to be greater in North America than Australia and the UK due to differences in exposure patterns between the sexes. (c) 2005 Wiley-Liss, Inc.

The importance of wet-powder dynamic mechanical properties in understanding granulation

Granule impact deformation has long been recognised as important in determining whether or not two colliding granules will coalesce. Work in the last 10 years has highlighted the fact that viscous effects are significant in granulation. The relative strengths of different formulations can vary with strain rate. Therefore, traditional strength measurements made at pseudo-static conditions give no indication, even qualitatively, of how materials will behave at high strain rates, and hence are actually misleading when used to model granule coalescence. This means that new standard methods need to be developed for determining the strain rates encountered by granules inside industrial equipment and also for measuring the mechanical properties of granules at these strain rates. The constitutive equations used in theoretical models of granule coalescence also need to be extended to include strain-rate dependent components.

Analysis and application of an equilibrium model for in vitro bioassay systems with three components: Receptor, hormone and hormone-binding-protein

A simple theoretical framework is presented for bioassay studies using three component in vitro systems. An equilibrium model is used to derive equations useful for predicting changes in biological response after addition of hormone-binding-protein or as a consequence of increased hormone affinity. Sets of possible solutions for receptor occupancy and binding protein occupancy are found for typical values of receptor and binding protein affinity constants. Unique equilibrium solutions are dictated by the initial condition of total hormone concentration. According to the occupancy theory of drug action, increasing the affinity of a hormone for its receptor will result in a proportional increase in biological potency. However, the three component model predicts that the magnitude of increase in biological potency will be a small fraction of the proportional increase in affinity. With typical initial conditions a two-fold increase in hormone affinity for its receptor is predicted to result in only a 33% increase in biological response. Under the same conditions an Ii-fold increase in hormone affinity for receptor would be needed to produce a two-fold increase in biological potency. Some currently used bioassay systems may be unrecognized three component systems and gross errors in biopotency estimates will result if the effect of binding protein is not calculated. An algorithm derived from the three component model is used to predict changes in biological response after addition of binding protein to in vitro systems. The algorithm is tested by application to a published data set from an experimental study in an in vitro system (Lim et al., 1990, Endocrinology 127, 1287-1291). Predicted changes show good agreement (within 8%) with experimental observations. (C) 1998 Academic Press Limited.

Evolution of three-dimensional folds for a non-Newtonian plate in a viscous medium

We derive analytical solutions for the three-dimensional time-dependent buckling of a non-Newtonian viscous plate in a less viscous medium. For the plate we assume a power-law rheology. The principal, axes of the stretching D-ij in the homogeneously deformed ground state are parallel and orthogonal to the bounding surfaces of the plate in the flat state. In the model formulation the action of the less viscous medium is replaced by equivalent reaction forces. The reaction forces are assumed to be parallel to the normal vector of the deformed plate surfaces. As a consequence, the buckling process is driven by the differences between the in-plane stresses and out of plane stress, and not by the in-plane stresses alone as assumed in previous models. The governing differential equation is essentially an orthotropic plate equation for rate dependent material, under biaxial pre-stress, supported by a viscous medium. The differential problem is solved by means of Fourier transformation and largest growth coefficients and corresponding wavenumbers are evaluated. We discuss in detail fold evolutions for isotropic in-plane stretching (D-11 = D-22), uniaxial plane straining (D-22 = 0) and in-plane flattening (D-11 = -2D(22)). Three-dimensional plots illustrate the stages of fold evolution for random initial perturbations or initial embryonic folds with axes non-parallel to the maximum compression axis. For all situations, one dominant set of folds develops normal to D-11, although the dominant wavelength differs from the Biot dominant wavelength except when the plate has a purely Newtonian viscosity. However, in the direction parallel to D-22, there exist infinitely many modes in the vicinity of the dominant wavelength which grow only marginally slower than the one corresponding to the dominant wavelength. This means that, except for very special initial conditions, the appearance of a three-dimensional fold will always be governed by at least two wavelengths. The wavelength in the direction parallel to D-11 is the dominant wavelength, and the wavelength(s) in the direction parallel to D-22 is determined essentially by the statistics of the initial state. A comparable sensitivity to the initial geometry does not exist in the classic two-dimensional folding models. In conformity with tradition we have applied Kirchhoff's hypothesis to constrain the cross-sectional rotations of the plate. We investigate the validity of this hypothesis within the framework of Reissner's plate theory. We also include a discussion of the effects of adding elasticity into the constitutive relations and show that there exist critical ratios of the relaxation times of the plate and the embedding medium for which two dominant wavelengths develop, one at ca. 2.5 of the classical Biot dominant wavelength and the other at ca. 0.45 of this wavelength. We propose that herein lies the origin of parasitic folds well known in natural examples.

Sustained gene expression in transplanted skin fibroblasts in rats

Retrovirus-mediated gene transfer into adult skin fibroblasts has provided measurable amounts of therapeutic proteins in animal models. However, the major problem emerging from these experiments was a limited time of vector encoded gene expression once transduced cells were engrafted We hypothesized that sustained transduced gene expression in quiescent fibroblasts in vivo might be obtained by using a fibronectin (Fn) promoter. Fibronectin plays a key role in cell adhesion, migration and wound healing and is up-regulated in quiescent fibroblasts. Retroviral vectors containing human adenosine deaminase (ADA) cDNA linked to rat fibronectin promoter (LNFnA) or viral LTR promoter (LASN) were compared for their ability to express ADA from transduced primary rat skin fibroblasts in vivo. Skin grafts formed from fibroblasts transduced with LNFnA showed strong human ADA enzyme activity from 1 week to 3 months. In contrast, skin grafts containing LASN-transduced fibroblasts tested positive for human ADA for weeks 1 and 2, were faintly positive at week 3 and showed no human ADA expression at 1, 2 and 3 months. Thus, a fibronectin promoter provided sustained transduced gene expression at high levels for at least 3 months in transplanted rat skin fibroblasts, perhaps permitting the targeting of this tissue for human gene therapy.

Trends in the incidence, severity, and short-term outcome of stroke in Perth, Western Australia

Background and Purpose-This report describes trends in the key indices of cerebrovascular disease over 6 years from the end of the 1980s in a geographically defined segment of the city of Perth, Western Australia. Methods-Identical methods were used to find and assess all cases of suspected stroke in a population of approximately 134 000 residents in a triangular area of the northern suburbs of Perth. Case fatality was measured as vital status at 28 days after the onset of symptoms. Data for first-ever strokes and for all strokes for equivalent periods of 12 months in 1989-1990 and 1995-1996 were compared by age-standardized rates and proportions and Poisson regression. Results-There were 355 strokes in 328 patients and 251 first-ever strokes (71%) for 1989-1990 and 290 events in 281 patients and 213 first-ever strokes (73%) for 1995-1996. In Poisson models including age and period, overall trends in the incidence of both first-ever strokes (rate ratio = 0.75; 95% confidence limits, 0.63, 0.90) and all strokes (rate ratio = 0.73; 95% confidence limits, 0.62, 0.85) were obviously significant, but only the changes in men were independently significant. Case fatality did not change, and the balance between hemorrhagic and occlusive strokes in 1995-1996 was almost indistinguishable from that observed in 1989-1990. Conclusions-Our results, which are the only longitudinal population-based data available for Australia for key indices of stroke, suggest that it is a change in the frequency of stroke, rather than its outcome, that is chiefly responsible nationally for the fall in mortality from cerebrovascular disease.

Predictors of death and vascular events in the elderly - The Perth Community Stroke Study

Background and Purpose-The goal of the present study was to identify risk factors for vascular disease in the elderly. Methods-We conducted a prospective study of control subjects from a population-based study of stroke in Perth, Western Australia, that was completed in 1989 to 1990 and used record linkage and a survey of survivors to identify deaths and nonfatal vascular events. Data validated through reference to medical records were analyzed with the use of Cox proportional hazards models. Results-Follow-up for the 931 subjects was 88% complete. By June 24, 1994, 198 (24%) of the subjects had died (96 from vascular disease), and there had been 45 nonfatal strokes or myocardial infarctions. The hazard ratio for diabetes exceeded 2.0 for all end points, whereas the consumption of meat >4 times weekly was associated with a reduction in risk of less than or equal to 30%. In most models, female sex and consumption of alcohol were associated with reduced risks, whereas previous myocardial infarction was linked to an increase in risk. Conclusions-There are only limited associations between lifestyle and major vascular illness in old age. Effective health promotion activities in early and middle life may be the key to a longer and healthier old age.

Accelerating seismic energy release and evolution of event time and size statistics: Results from two heterogeneous cellular automaton models

The evolution of event time and size statistics in two heterogeneous cellular automaton models of earthquake behavior are studied and compared to the evolution of these quantities during observed periods of accelerating seismic energy release Drier to large earthquakes. The two automata have different nearest neighbor laws, one of which produces self-organized critical (SOC) behavior (PSD model) and the other which produces quasi-periodic large events (crack model). In the PSD model periods of accelerating energy release before large events are rare. In the crack model, many large events are preceded by periods of accelerating energy release. When compared to randomized event catalogs, accelerating energy release before large events occurs more often than random in the crack model but less often than random in the PSD model; it is easier to tell the crack and PSD model results apart from each other than to tell either model apart from a random catalog. The evolution of event sizes during the accelerating energy release sequences in all models is compared to that of observed sequences. The accelerating energy release sequences in the crack model consist of an increase in the rate of events of all sizes, consistent with observations from a small number of natural cases, however inconsistent with a larger number of cases in which there is an increase in the rate of only moderate-sized events. On average, no increase in the rate of events of any size is seen before large events in the PSD model.

Estimating HIV incidence rates from age prevalence data in epidemic situations

We present a method of estimating HIV incidence rates in epidemic situations from data on age-specific prevalence and changes in the overall prevalence over time. The method is applied to women attending antenatal clinics in Hlabisa, a rural district of KwaZulu/Natal, South Africa, where transmission of HIV is overwhelmingly through heterosexual contact. A model which gives age-specific prevalence rates in the presence of a progressing epidemic is fitted to prevalence data for 1998 using maximum likelihood methods and used to derive the age-specific incidence. Error estimates are obtained using a Monte Carlo procedure. Although the method is quite general some simplifying assumptions are made concerning the form of the risk function and sensitivity analyses are performed to explore the importance of these assumptions. The analysis shows that in 1998 the annual incidence of infection per susceptible woman increased from 5.4 per cent (3.3-8.5 per cent; here and elsewhere ranges give 95 per cent confidence limits) at age 15 years to 24.5 per cent (20.6-29.1 per cent) at age 22 years and declined to 1.3 per cent (0.5-2.9 per cent) at age 50 years; standardized to a uniform age distribution, the overall incidence per susceptible woman aged 15 to 59 was 11.4 per cent (10.0-13.1 per cent); per women in the population it was 8.4 per cent (7.3-9.5 per cent). Standardized to the age distribution of the female population the average incidence per woman was 9.6 per cent (8.4-11.0 per cent); standardized to the age distribution of women attending antenatal clinics, it was 11.3 per cent (9.8-13.3 per cent). The estimated incidence depends on the values used for the epidemic growth rate and the AIDS related mortality. To ensure that, for this population, errors in these two parameters change the age specific estimates of the annual incidence by less than the standard deviation of the estimates of the age specific incidence, the AIDS related mortality should be known to within +/-50 per cent and the epidemic growth rate to within +/-25 per cent, both of which conditions are met. In the absence of cohort studies to measure the incidence of HIV infection directly, useful estimates of the age-specific incidence can be obtained from cross-sectional, age-specific prevalence data and repeat cross-sectional data on the overall prevalence of HIV infection. Several assumptions were made because of the lack of data but sensitivity analyses show that they are unlikely to affect the overall estimates significantly. These estimates are important in assessing the magnitude of the public health problem, for designing vaccine trials and for evaluating the impact of interventions. Copyright (C) 2001 John Wiley & Sons, Ltd.

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.

Analytical solution for tidal propagation in a coupled semi-confined/phreatic coastal aquifer

An analytical solution is derived for tidal fluctuations in a coupled coastal aquifer system consisting of a semi-confined aquifer, a thin semi-permeable layer and a phreatic aquifer. Based on the solution, we study the interactions (via leakage) between the confined and unconfined aquifers in response to tides. The results show that, under certain conditions, leakage from the confined aquifer can affect considerably the tidal water table fluctuation in the phreatic aquifer and vice versa. Ignoring these effects could lead to errors in estimating aquifer properties based on tidal signals. (C) 2002 Elsevier Science Ltd. All rights reserved.

Entanglement of indistinguishable particles in condensed-matter physics

The concept of entanglement in systems where the particles are indistinguishable has been the subject of much recent interest and controversy. In this paper we study the notion of entanglement of particles introduced by Wiseman and Vaccaro [Phys. Rev. Lett. 91, 097902 (2003)] in several specific physical systems, including some that occur in condensed-matter physics. The entanglement of particles is relevant when the identical particles are itinerant and so not distinguished by their position as in spin models. We show that entanglement of particles can behave differently than other approaches that have been used previously, such as entanglement of modes (occupation-number entanglement) and the entanglement in the two-spin reduced density matrix. We argue that the entanglement of particles is what could actually be measured in most experimental scenarios and thus its physical significance is clear. This suggests that entanglement of particles may be useful in connecting theoretical and experimental studies of entanglement in condensed-matter systems.

The role of vaccines in the control of STDs: HPV vaccines

Prophylactic vaccines for genital human papillomavirus (HPV) infection have been shown to be feasible in animal models, and suitable vaccine material based on virus-like particles can be produced in bulk at reasonable cost. Initiation of phase III clinical trials will follow definition of trial outcome measures through further epidemiological studies, and development-of assays of host protective immunity. Vaccines could in principle eliminate HPV-related disease, as the human race is the only natural host for the relevant papillomaviruses (PVs). Therapeutic vaccines for genital HPV infection are also possible, but have not yet been demonstrated as feasible in practice because the choice of vaccine antigens is difficult, the method of their optimal delivery is uncertain, and the nature of the relevant antiviral immunity is unknown. PV species specificity will require trials to be conducted in man, which will slow definition of an ideal vaccine.