Maximum Likelihood Estimation of Mixture Densities for Binned and Truncated Multivariate Data

Binning and truncation of data are common in data analysis and machine learning. This paper addresses the problem of fitting mixture densities to multivariate binned and truncated data. The EM approach proposed by McLachlan and Jones (Biometrics, 44: 2, 571-578, 1988) for the univariate case is generalized to multivariate measurements. The multivariate solution requires the evaluation of multidimensional integrals over each bin at each iteration of the EM procedure. Naive implementation of the procedure can lead to computationally inefficient results. To reduce the computational cost a number of straightforward numerical techniques are proposed. Results on simulated data indicate that the proposed methods can achieve significant computational gains with no loss in the accuracy of the final parameter estimates. Furthermore, experimental results suggest that with a sufficient number of bins and data points it is possible to estimate the true underlying density almost as well as if the data were not binned. The paper concludes with a brief description of an application of this approach to diagnosis of iron deficiency anemia, in the context of binned and truncated bivariate measurements of volume and hemoglobin concentration from an individual's red blood cells.

Transmission dynamics of the Echinococcus granulosus sheep-dog strain (G1 genotype) in camels in Tunisia

Cystic echinococcosis, caused by Echinococcus grantilosus, is highly endemic in North Africa and the Middle East. This paper examines the abundance and prevalence of infection of E. granulosus in camels in Tunisia. No cysts were found in 103 camels from Kebili, whilst 19 of 188 camels from Benguerden (10.1%) were infected. Of the cysts found 95% were considered fertile with the presence of protoscolices and 80% of protoscolices were considered viable by their ability to exclude aqueous eosin. Molecular techniques were used on cyst material from camels and this demonstrated that the study animals were infected with the G1 sheep strain of E. granulosus. Observed data were fitted to a mathematical model by maximum likelihood techniques to define the parameters and their confidence limits and the negative binomial distribution was used to define the error variance in the observed data. The infection pressure to camels was somewhat lower in comparison to sheep reported in an earlier study. However, because camels are much longer-lived animals, the results of the model fit suggested that older camels have a relatively high prevalence rate, reaching a most likely value of 32% at age 15 years. This could represent an important source of transmission to dogs and hence indirectly to man of this zonotic strain. In common with similar studies on other species, there was no evidence of parasite-induced immunity in camels. (C) 2004 Elsevier B.V. All rights reserved.

Anthropometric measures in relation to Basal Cell Carcinoma: a longitudinal study

Background: The relationship between anthropometric indices and risk of basal cell carcinoma ( BCC) is largely unknown. We aimed to examine the association between anthropometric measures and development of BCC and to demonstrate whether adherence to World Health Organisation guidelines for body mass index, waist circumference, and waist/ hip ratio was associated with risk of BCC, independent of sun exposure. Methods: Study participants were participants in a community- based skin cancer prevention trial in Nambour, a town in southeast Queensland ( latitude 26 degrees S). In 1992, height, weight, and waist and hip circumferences were measured for all 1621 participants and weight was remeasured at the end of the trial in 1996. Prevalence proportion ratios were calculated using a log- binomial model to estimate the risk of BCC prior to or prevalent in 1992, while Poisson regression with robust error variances was used to estimate the relative risk of BCC during the follow- up period. Results: At baseline, 94 participants had a current BCC, and 202 had a history of BCC. During the 5- year follow- up period, 179 participants developed one or more new BCCs. We found no significant association between any of the anthropometric measures or indices and risk of BCC after controlling for potential confounding factors including sun exposure. There was a suggestion that short- term weight gain may increase the risk of developing BCC for women only. Conclusion: Adherence to World Health Organisation guidelines for body mass index, waist circumference and waist/ hip ratio is not significantly associated with occurrence of basal cell carcinomas of the skin.

Indirect evidence of density-dependent population regulation in Aponomma hydrosauri (Acari : Ixodidae), an ectoparasite of reptiles

The extent to which density-dependent processes regulate natural populations is the subject of an ongoing debate. We contribute evidence to this debate showing that density-dependent processes influence the population dynamics of the ectoparasite Aponomma hydrosauri (Acari: Ixodidae), a tick species that infests reptiles in Australia. The first piece of evidence comes from an unusually long-term dataset on the distribution of ticks among individual hosts. If density-dependent processes are influencing either host mortality or vital rates of the parasite population, and those distributions can be approximated with negative binomial distributions, then general host-parasite models predict that the aggregation coefficient of the parasite distribution will increase with the average intensity of infections. We fit negative binomial distributions to the frequency distributions of ticks on hosts, and find that the estimated aggregation coefficient k increases with increasing average tick density. This pattern indirectly implies that one or more vital rates of the tick population must be changing with increasing tick density, because mortality rates of the tick's main host, the sleepy lizard, Tiliqua rugosa, are unaffected by changes in tick burdens. Our second piece of evidence is a re-analysis of experimental data on the attachment success of individual ticks to lizard hosts using generalized linear modelling. The probability of successful engorgement decreases with increasing numbers of ticks attached to a host. This is direct evidence of a density-dependent process that could lead to an increase in the aggregation coefficient of tick distributions described earlier. The population-scale increase in the aggregation coefficient is indirect evidence of a density-dependent process or processes sufficiently strong to produce a population-wide pattern, and thus also likely to influence population regulation. The direct observation of a density-dependent process is evidence of at least part of the responsible mechanism.

Polymorphisms in glutathione S-transferases and non-melanoma skin cancer risk in Australian renal transplant recipients

Caucasian renal transplant recipients from Queensland, Australia have the highest non-melanoma skin cancer (NMSC) risk worldwide. Although ultraviolet light (UVR) exposure is critical, genetic factors also appear important. We and others have shown that polymorphism in the glutathione S-transferases (GST) is associated with NMSC in UK recipients. However, the effect of high UVR exposure and differences in immunosuppressive regimen on these associations is unknown. In this study, we examined allelism in GSTM1, GSTM3, GSTT1 and GSTP1 in 361 Queensland renal transplant recipients. Data on squamous (SCC) and basal cell carcinoma (BCC), UVR/tobacco exposure and genotype were obtained. Associations with both NMSC risk and numbers were examined using logistic and negative binomial regression, respectively. In the total group, GSTM1 AB [P = 0.049, rate ratio (RR) = 0.23] and GSTM3 AA (P = 0.015, RR = 0.50) were associated with fewer SCC. Recipients were then stratified by prednisolone dose (less than or equal to7 versus >7 mg/day). In the low-dose group, GSTT1 null (P = 0.006, RR = 0.20) and GSTP1 Val/Val (P = 0.021, RR = 0.20) were associated with SCC numbers. In contrast, in the high-dose group, GSTM1 AB (P = 0.009, RR = 0.05), GSTM3 AB (P = 0.042, RR = 2.29) and BB (P = 0.014, RR = 5.31) and GSTP1 Val/Val (P = 0.036, RR = 2.98) were associated with SCC numbers. GSTM1 AB (P = 0.016) and GSTP1 Val/Val (P = 0.046) were also associated with fewer BCC in this group. GSTP1 associations were strongest in recipients with lower UVR/tobacco exposure. The data confirm our UK findings, suggesting that protection against UVR-induced oxidative stress is important in NMSC development in recipients, but that this effect depends on the immunosuppressant regimen.

Food intake and risk of squamous cell carcinoma of the skin in a community: The Nambour skin cancer cohort study

There is some evidence that dietary factors may modify the risk of squamous cell carcinoma (SCC) of the skin, but the association between food intake and SCC has not been evaluated prospectively. We examined the association between food intake and SCC incidence among 1,056 randomly selected adults living in an Australian sub-tropical community. Measurement-error corrected estimates of intake in 15 food groups were defined from a validated food frequency questionnaire in 1992. Associations with SCC risk were assessed using Poisson and negative binomial regression to the persons affected and tumour counts, respectively, based on incident, histologically confirmed tumours occurring between 1992 and 2002. After multivariable adjustment, none of the food groups was significantly associated with SCC risk. Stratified analysis in participants with a past history of skin cancer showed a decreased risk of SCC tumours for high intakes of green leafy vegetables (RR = 0.45, 95% CI = 0.22-0.91; p for trend = 0.02) and an increased risk for high intake of unmodified dairy products (RR = 2.53, 95% CI: 1.15-5.54; p for trend = 0.03). Food intake was not associated with SCC risk in persons who had no past history of skin cancer. These findings suggest that consumption of green leafy vegetables may help prevent development of subsequent SCCs of the skin among people with previous skin cancer and that consumption of unmodified dairy products, such as whole milk, cheese and yoghurt, may increase SCC risk in susceptible persons.

A Comparison of Bayesian and Sampling Theory Inferences in a Probit Model

HE PROBIT MODEL IS A POPULAR DEVICE for explaining binary choice decisions in econometrics. It has been used to describe choices such as labor force participation, travel mode, home ownership, and type of education. These and many more examples can be found in papers by Amemiya (1981) and Maddala (1983). Given the contribution of economics towards explaining such choices, and given the nature of data that are collected, prior information on the relationship between a choice probability and several explanatory variables frequently exists. Bayesian inference is a convenient vehicle for including such prior information. Given the increasing popularity of Bayesian inference it is useful to ask whether inferences from a probit model are sensitive to a choice between Bayesian and sampling theory techniques. Of interest is the sensitivity of inference on coefficients, probabilities, and elasticities. We consider these issues in a model designed to explain choice between fixed and variable interest rate mortgages. Two Bayesian priors are employed: a uniform prior on the coefficients, designed to be noninformative for the coefficients, and an inequality restricted prior on the signs of the coefficients. We often know, a priori, whether increasing the value of a particular explanatory variable will have a positive or negative effect on a choice probability. This knowledge can be captured by using a prior probability density function (pdf) that is truncated to be positive or negative. Thus, three sets of results are compared:those from maximum likelihood (ML) estimation, those from Bayesian estimation with an unrestricted uniform prior on the coefficients, and those from Bayesian estimation with a uniform prior truncated to accommodate inequality restrictions on the coefficients.

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.

An integrated cognitive model of panic disorder: The role of positive and negative cognitions

One reason for the neglect of the role of positive factors in cognitive-behavioural therapy (CBT) may relate to a failure to develop cognitive models that integrate positive and negative cognitions. Bandura [Psychol. Rev. 84 (1977) 191; Anxiety Res. 1 (1988) 77] proposed that self-efficacy beliefs mediate a range of emotional and behavioural outcomes. However, in panic disorder, cognitively based research to date has largely focused on catastrophic misinterpretation of bodily sensations. Although a number of studies support each of the predictions associated with the account of panic disorder that is based on the role of negative cognitions, a review of the literature indicated that a cognitively based explanation of the disorder may be considerably strengthened by inclusion of positive cognitions that emphasize control or coping. Evidence to support an Integrated Cognitive Model (ICM) of panic disorder was examined and the theoretical implications of this model were discussed in terms of both schema change and compensatory skills accounts of change processes in CBT. (C) 2004 Elsevier Ltd. All rights reserved.

The moderating role of social support in Karasek's job strain model

This paper examines whether social support is a boundary-determining criterion in the job strain model of Karasek (1979). The particular focus is the extent to which different sources of social support, work overload and task control influence job satisfaction, depersonalization and supervisor assessments of work performance. Hypotheses are tested using prospective survey data from 80 clerical staff in a university setting. Results revealed 3-way interactions among levels of support (supervisor, co-worker, non-work), perceived task control and work overload on levels of work performance and employee adjustment (self-report). After controlling for levels of negative affect in all analyses, there was evidence that high levels of supervisor support mitigated against the negative effects of high strain jobs on levels of job satisfaction and reduced reported levels of depersonalization. Moreover, high levels of non-work support and co-worker support also mitigated against the negative effects of high strain jobs on levels of work performance. The results are discussed in terms of the importance of social support networks both at, and beyond, the work context.

Interpretation of negative values of the interaction parameter in the adsorption equation through the effects of surface layer heterogeneity

The problem of the negative values of the interaction parameter in the equation of Frumkin has been analyzed with respect to the adsorption of nonionic molecules on energetically homogeneous surface. For this purpose, the adsorption states of a homologue series of ethoxylated nonionic surfactants on air/water interface have been determined using four different models and literature data (surface tension isotherms). The results obtained with the Frumkin adsorption isotherm imply repulsion between the adsorbed species (corresponding to negative values of the interaction parameter), while the classical lattice theory for energetically homogeneous surface (e.g., water/air) admits attraction alone. It appears that this serious contradiction can be overcome by assuming heterogeneity in the adsorption layer, that is, effects of partial condensation (formation of aggregates) on the surface. Such a phenomenon is suggested in the Fainerman-Lucassen-Reynders-Miller (FLM) 'Aggregation model'. Despite the limitations of the latter model (e.g., monodispersity of the aggregates), we have been able to estimate the sign and the order of magnitude of Frumkin's interaction parameter and the range of the aggregation numbers of the surface species. (C) 2004 Elsevier B.V All rights reserved.

Ferromagnetism, paramagnetism, and a Curie-Weiss metal in an electron-doped Hubbard model on a triangular lattice

Motivated by the unconventional properties and rich phase diagram of NaxCoO2 we consider the electronic and magnetic properties of a two-dimensional Hubbard model on an isotropic triangular lattice doped with electrons away from half-filling. Dynamical mean-field theory (DMFT) calculations predict that for negative intersite hopping amplitudes (t < 0) and an on-site Coulomb repulsion, U, comparable to the bandwidth, the system displays properties typical of a weakly correlated metal. In contrast, for t > 0 a large enhancement of the effective mass, itinerant ferromagnetism, and a metallic phase with a Curie-Weiss magnetic susceptibility are found in a broad electron doping range. The different behavior encountered is a consequence of the larger noninteracting density of states (DOS) at the Fermi level for t > 0 than for t < 0, which effectively enhances the mass and the scattering amplitude of the quasiparticles. The shape of the DOS is crucial for the occurrence of ferromagnetism as for t > 0 the energy cost of polarizing the system is much smaller than for t < 0. Our observation of Nagaoka ferromagnetism is consistent with the A-type antiferromagnetism (i.e., ferromagnetic layers stacked antiferromagnetically) observed in neutron scattering experiments on NaxCoO2. The transport and magnetic properties measured in NaxCoO2 are consistent with DMFT predictions of a metal close to the Mott insulator and we discuss the role of Na ordering in driving the system towards the Mott transition. We propose that the Curie-Weiss metal phase observed in NaxCoO2 is a consequence of the crossover from a bad metal with incoherent quasiparticles at temperatures T > T-* and Fermi liquid behavior with enhanced parameters below T-*, where T-* is a low energy coherence scale induced by strong local Coulomb electron correlations. Our analysis also shows that the one band Hubbard model on a triangular lattice is not enough to describe the unusual properties of NaxCoO2 and is used to identify the simplest relevant model that captures the essential physics in NaxCoO2. We propose a model which allows for the Na ordering phenomena observed in the system which, we propose, drives the system close to the Mott insulating phase even at large dopings.

Expression of a truncated Brca1 protein delays lactational mammary development in transgenic mice

To address the hypothesis that certain disease-associated mutants of the breast-ovarian cancer susceptibility gene BRCA1 have biological activity in vivo, we have expressed a truncated Brca1 protein (trBrca1) in cell-lines and in the mammary gland of transgenic mice. Immunofluorescent analysis of transfected cell-lines indicates that trBRCA1 is a stable protein and that it is localized in the cell cytoplasm. Functional analysis of these cell-lines indicates that expression of trBRCA1 confers an increased radiosensitivity phenotype on mammary epithelial cells, consistent with abrogation of the BRCA1 pathway. MMTV-trBrca1 transgenic mice from two independent lines displayed a delay in lactational mammary gland development, as demonstrated by altered histological profiles of lobuloalveolar structures. Cellular and molecular analyses indicate that this phenotype results from a defect in differentiation, rather than altered rates of proliferation or apoptosis. The results presented in this paper are consistent with trBrca1 possessing dominant-negative activity and playing an important role in regulating normal mammary development. They may also have implications for germline carriers of BRCA1 mutations.

Dynamical properties of a strongly correlated model for quarter-filled layered organic molecular crystals

The dynamical properties of an extended Hubbard model, which is relevant to quarter-filled layered organic molecular crystals, are analyzed. We have computed the dynamical charge correlation function, spectral density, and optical conductivity using Lanczos diagonalization and large-N techniques. As the ratio of the nearest-neighbor Coulomb repulsion, V, to the hopping integral, t, increases there is a transition from a metallic phase to a charge-ordered phase. Dynamical properties close to the ordering transition are found to differ from the ones expected in a conventional metal. Large-N calculations display an enhancement of spectral weight at low frequencies as the system is driven closer to the charge-ordering transition in agreement with Lanczos calculations. As V is increased the charge correlation function displays a collective mode which, for wave vectors close to (pi,pi), increases in amplitude and softens as the charge-ordering transition is approached. We propose that inelastic x-ray scattering be used to detect this mode. Large-N calculations predict superconductivity with d(xy) symmetry close to the ordering transition. We find that this is consistent with Lanczos diagonalization calculations, on lattices of 20 sites, which find that the binding energy of two holes becomes negative close to the charge-ordering transition.