Automatic ICD-10 classification of cancers from free-text death certificates

Objective Death certificates provide an invaluable source for cancer mortality statistics; however, this value can only be realised if accurate, quantitative data can be extracted from certificates – an aim hampered by both the volume and variable nature of certificates written in natural language. This paper proposes an automatic classification system for identifying cancer related causes of death from death certificates. Methods Detailed features, including terms, n-grams and SNOMED CT concepts were extracted from a collection of 447,336 death certificates. These features were used to train Support Vector Machine classifiers (one classifier for each cancer type). The classifiers were deployed in a cascaded architecture: the first level identified the presence of cancer (i.e., binary cancer/nocancer) and the second level identified the type of cancer (according to the ICD-10 classification system). A held-out test set was used to evaluate the effectiveness of the classifiers according to precision, recall and F-measure. In addition, detailed feature analysis was performed to reveal the characteristics of a successful cancer classification model. Results The system was highly effective at identifying cancer as the underlying cause of death (F-measure 0.94). The system was also effective at determining the type of cancer for common cancers (F-measure 0.7). Rare cancers, for which there was little training data, were difficult to classify accurately (F-measure 0.12). Factors influencing performance were the amount of training data and certain ambiguous cancers (e.g., those in the stomach region). The feature analysis revealed a combination of features were important for cancer type classification, with SNOMED CT concept and oncology specific morphology features proving the most valuable. Conclusion The system proposed in this study provides automatic identification and characterisation of cancers from large collections of free-text death certificates. This allows organisations such as Cancer Registries to monitor and report on cancer mortality in a timely and accurate manner. In addition, the methods and findings are generally applicable beyond cancer classification and to other sources of medical text besides death certificates.

Using harvest statistics to monitor temporal variation in kangaroo density and harvest rate

Management of the commercial harvest of kangaroos relies on quotas set annually as a proportion of regular estimates of population size. Surveys to generate these estimates are expensive and, in the larger states, logistically difficult; a cheaper alternative is desirable. Rainfall is a disappointingly poor predictor of kangaroo rate of increase in many areas, but harvest statistics (sex ratio, carcass weight, skin size and animals shot per unit time) potentially offer cost-effective indirect monitoring of population abundance (and therefore trend) and status (i.e. under-or overharvest). Furthermore, because harvest data are collected continuously and throughout the harvested areas, they offer the promise of more intensive and more representative coverage of harvest areas than aerial surveys do. To be useful, harvest statistics would need to have a close and known relationship with either population size or harvest rate. We assessed this using longterm (11-22 years) data for three kangaroo species (Macropus rufus, M. giganteus and M. fuliginosus) and common wallaroos (M. robustus) across South Australia, New South Wales and Queensland. Regional variation in kangaroo body size, population composition, shooter efficiency and selectivity required separate analyses in different regions. Two approaches were taken. First, monthly harvest statistics were modelled as a function of a number of explanatory variables, including kangaroo density, harvest rate and rainfall. Second, density and harvest rate were modelled as a function of harvest statistics. Both approaches incorporated a correlated error structure. Many but not all regions had relationships with sufficient precision to be useful for indirect monitoring. However, there was no single relationship that could be applied across an entire state or across species. Combined with rainfall-driven population models and applied at a regional level, these relationships could be used to reduce the frequency of aerial surveys without compromising decisions about harvest management.

Travelling with medications and medical equipment across international borders

The international traveller needs to plan ahead to ensure medicines are available and used as directed for optimal therapeutic outcome. The planning needs to take account of legal and customs requirements for travelling with medicines for personal use. The standard advice by travel health providers is that travellers should check with the country of destination for requirements when travelling into the country with medicines for personal use. This is akin to introducing a barrier to care for this category of travellers. Innovative method of care for this group of traveller is needed.

Non-medical prescribing in Australia

The announcement in the 2009 federal budget to allow nurse practitioners and midwives access to the Pharmaceutical Benefits Scheme (PBS) and the Medicare Benefits Scheme,1and the subsequent announcement of a November 2010 start date,2has brought non-medical prescribing into the public arena. Non-medical prescribing is not a new concept in Australia as nurse practitioners, podiatrists and optometrists have been authorised to prescribe under various state legislations for some time. However, state legislation is not uniform in relation to authorisation or formulary. Midwives are currently seeking prescribing rights,3and other groups such as physiotherapists and pharmacists are likely to seek them in the future.

Association between endometriosis and the interleukin 1A (IL1A) locus

STUDY QUESTION Are single-nucleotide polymorphisms (SNPs) at the interleukin 1A (IL1A) gene locus associated with endometriosis risk? SUMMARY ANSWER We found evidence for strong association between IL1A SNPs and endometriosis risk. WHAT IS KNOWN ALREADY Genetic factors contribute substantially to the complex aetiology of endometriosis and the disease has an estimated heritability of ∼51%. We, and others, have conducted genome-wide association (GWA) studies for endometriosis, which identified a total of nine independent risk loci. Recently, two small Japanese studies reported eight SNPs (rs6542095, rs11677416, rs3783550, rs3783525, rs3783553, rs2856836, rs1304037 and rs17561) at the IL1A gene locus as suggestively associated with endometriosis risk. There is also evidence of a link between inflammation and endometriosis. STUDY DESIGN, SIZE, DURATION We sought to further investigate the eight IL1A SNPs for association with endometriosis using an independent sample of 3908 endometriosis cases and 8568 controls of European and Japanese ancestry. The study was conducted between October 2013 and July 2014. PARTICIPANTS/MATERIALS, SETTING, METHODS By leveraging GWA data from our previous multi-ethnic GWA meta-analysis for endometriosis, we imputed variants in the IL1A region, using a recent 1000 Genomes reference panel. After combining summary statistics for the eight SNPs from our European and Japanese imputed data with the published results, a fixed-effect meta-analysis was performed. An additional meta-analysis restricted to endometriosis cases with moderate-to-severe (revised American Fertility Society stage 3 or 4) disease versus controls was also performed. MAIN RESULTS AND THE ROLE OF CHANCE All eight IL1A SNPs successfully replicated at P < 0.014 in the European imputed data with concordant direction and similar size to the effects reported in the original Japanese studies. Of these, three SNPs (rs6542095, rs3783550 and rs3783525) also showed association with endometriosis at a nominal P < 0.05 in our independent Japanese sample. Fixed-effect meta-analysis of the eight SNPs for moderate-to-severe endometriosis produced a genome-wide significant association for rs6542095 (odds ratio = 1.21; 95% confidence interval = 1.13–1.29; P = 3.43 × 10−8). LIMITATIONS, REASONS FOR CAUTION The meta-analysis for moderate-to-severe endometriosis included results of moderate-to-severe endometriosis cases from our European data sets and all endometriosis cases from the Japanese data sets, as disease stage information was not available for endometriosis cases in the Japanese data sets. WIDER IMPLICATIONS OF THE FINDINGS SNP rs6542095 is located ∼2.3 kb downstream of the IL1A gene and ∼6.9 kb upstream of cytoskeleton-associated protein 2-like (CKAP2L) gene. The IL1A gene encodes the IL1a protein, a member of the interleukin 1 cytokine family which is involved in various immune responses and inflammatory processes. These results provide important replication in an independent Japanese sample and, for the first time, association of the IL1A locus in endometriosis patients of European ancestry. SNPs within the IL1A locus may regulate other genes, but if IL1A is the target, our results provide supporting evidence for a link between inflammatory responses and the pathogenesis of endometriosis. STUDY FUNDING/COMPETING INTEREST(S) The research was funded by grants from the Australian National Health and Medical Research Council and Wellcome Trust. None of the authors has competing interests for the study.

Noise sensitivity medical, psychological and genetic aspects

Noise can be defined as unwanted sound. It may adversely affect the health and well-being of individuals. Noise sensitivity is a personality trait covering attitudes towards noise in general and a predictor of noise annoyance. Noise sensitive individuals are more affected by noise than less sensitive individuals. The determinants and characteristics related to noise sensitivity are rather poorly known. The risk of health effects caused by noise can be hypothesized to be higher for noise sensitive individuals compared to those who are not noise sensitive. A cardiovascular disease may be an example of outcomes. The general aim of the present study was to investigate the association of noise sensitivity with specific somatic and psychological factors, including the genetic component of noise sensitivity, and the association of noise sensitivity with mortality. The study was based on the Finnish Twin Cohort of same-sex twin pairs born before 1958. In 1988 a questionnaire was sent to twin pairs discordant for hypertension. 1495 individuals (688 men, 807 women) aged 31 88 years replied, including 573 twin pairs. 218 of the subjects lived in the Helsinki Metropolitan Area. Self-reported noise sensitivity, lifetime noise exposure and hypertension were obtained from the questionnaire study in 1988 and other somatic and psychological factors from the questionnaire study in 1981 for the same individuals. In addition, noise map information (1988 1992) from the Helsinki Metropolitan Area and mortality follow-up 1989 2003 were used. To evaluate the stability and validity of noise sensitivity, a new questionnaire was sent in 2002 to a sample of the subjects who had replied to the 1988 questionnaire. Of all subjects who had answered the question on noise sensitivity, 38 % were noise sensitive. Noise sensitivity was independent of noise exposure levels indicated in noise maps. Subjects with high noise sensitivity reported more transportation noise exposure than subjects with low noise sensitivity. Noise sensitive subjects reported transportation noise exposure outside the environmental noise map areas almost twice as often as non-sensitive subjects. Noise sensitivity was associated with hypertension, emphysema, use of psychotropic drugs, smoking, stress and hostility, even when lifetime noise exposure was adjusted for. Monozygotic twin pairs were more similar with regards noise sensitivity than dizygotic twin pairs, and quantitative genetic modelling indicated significant familiality. The best fitting genetic model provided an estimate of heritability of 36 %. Follow-up of subjects in the case-control study showed that cardiovascular mortality was significantly increased among noise sensitive women, but not among men. For coronary heart mortality the interaction of noise sensitivity and lifetime noise exposure was statistically significant in women. In conclusion, noise sensitivity has both somatic and psychological components. It does aggregate in families and probably has a genetic component. Noise sensitivity may be a risk factor for cardiovascular mortality in women.

Amber teething necklaces – medical marvel or maternal myth?

BACKGROUND: Baltic amber-bead necklaces or bracelets are commonly used for managing teething symptoms in infants. The effectiveness of these beads is claimed to be from succinic acid release (a compound with analgesic and anti-inflammatory properties), which is then absorbed through the skin. AIM: To investigate whether Baltic amber teething necklaces purchased in Australia contained succinic acid, and to quantify succinic acid release from the beads. METHODS: Infrared spectroscopy was used to confirm that the teething necklaces were made of Baltic amber. The amount of succinic acid contained within the beads was quantified, and succinic acid release from intact beads was measured in phosphate buffered saline (PBS) pH 5.5 or octanol to simulate aqueous or oily skin environments. RESULTS: Each necklace (33 beads in length) contained 19.17±4.89 mg of succinic acid (mean±se). Over a 6-month period, no succinic acid was detected in PBS, while 0.13±0.09 mg of succinic acid per necklace was released in octanol. Only one replicate of amber beads in octanol released succinic acid, and they had fragmented, with shards free-floating in the solvent. DISCUSSION: It is likely succinic acid was only detected because the beads were breaking down in octanol, which does not occur when worn around the neck of a child. Furthermore, the hydrophilic properties of succinic acid would not favour its absorption across hydrophobic layers of the skin and into the bloodstream. CONCLUSION: While the teething necklaces do contain small quantities of succinic acid, it is highly unlikely to be released from intact beads.

Susceptibility locus on chromosome 1q23-25 for a schizophrenia subtype resembling deficit schizophrenia identified by latent class analysis

Context: Identifying susceptibility genes for schizophrenia may be complicated by phenotypic heterogeneity, with some evidence suggesting that phenotypic heterogeneity reflects genetic heterogeneity. Objective: To evaluate the heritability and conduct genetic linkage analyses of empirically derived, clinically homogeneous schizophrenia subtypes. Design: Latent class and linkage analysis. Setting: Taiwanese field research centers. Participants: The latent class analysis included 1236 Han Chinese individuals with DSM-IV schizophrenia. These individuals were members of a large affected-sibling-pair sample of schizophrenia (606 ascertained families), original linkage analyses of which detected a maximum logarithm of odds (LOD) of 1.8 (z = 2.88) on chromosome 10q22.3. Main Outcome Measures: Multipoint exponential LOD scores by latent class assignment and parametric heterogeneity LOD scores. Results: Latent class analyses identified 4 classes, with 2 demonstrating familial aggregation. The first (LC2) described a group with severe negative symptoms, disorganization, and pronounced functional impairment, resembling “deficit schizophrenia.” The second (LC3) described a group with minimal functional impairment, mild or absent negative symptoms, and low disorganization. Using the negative/deficit subtype, we detected genome-wide significant linkage to 1q23-25 (LOD = 3.78, empiric genome-wide P = .01). This region was not detected using the DSM-IV schizophrenia diagnosis, but has been strongly implicated in schizophrenia pathogenesis by previous linkage and association studies.Variants in the 1q region may specifically increase risk for a negative/deficit schizophrenia subtype. Alternatively, these results may reflect increased familiality/heritability of the negative class, the presence of multiple 1q schizophrenia risk genes, or a pleiotropic 1q risk locus or loci, with stronger genotype-phenotype correlation with negative/deficit symptoms. Using the second familial latent class, we identified nominally significant linkage to the original 10q peak region. Conclusion: Genetic analyses of heritable, homogeneous phenotypes may improve the power of linkage and association studies of schizophrenia and thus have relevance to the design and analysis of genome-wide association studies.

Optimal simultaneous maximum a posterioriestimation of states, noise statistics and parameters I. Algorithm

The simultaneous state and parameter estimation problem for a linear discrete-time system with unknown noise statistics is treated as a large-scale optimization problem. The a posterioriprobability density function is maximized directly with respect to the states and parameters subject to the constraint of the system dynamics. The resulting optimization problem is too large for any of the standard non-linear programming techniques and hence an hierarchical optimization approach is proposed. It turns out that the states can be computed at the first levelfor given noise and system parameters. These, in turn, are to be modified at the second level.The states are to be computed from a large system of linear equations and two solution methods are considered for solving these equations, limiting the horizon to a suitable length. The resulting algorithm is a filter-smoother, suitable for off-line as well as on-line state estimation for given noise and system parameters. The second level problem is split up into two, one for modifying the noise statistics and the other for modifying the system parameters. An adaptive relaxation technique is proposed for modifying the noise statistics and a modified Gauss-Newton technique is used to adjust the system parameters.

Optimal simultaneous maximum a posterioriestimation of states, noise statistics and parameters II. Numerical performance

A very general and numerically quite robust algorithm has been proposed by Sastry and Gauvrit (1980) for system identification. The present paper takes it up and examines its performance on a real test example. The example considered is the lateral dynamics of an aircraft. This is used as a vehicle for demonstrating the performance of various aspects of the algorithm in several possible modes.

The Spin-Statistics Relation and Noncommutative Quantum Field Theory

The efforts of combining quantum theory with general relativity have been great and marked by several successes. One field where progress has lately been made is the study of noncommutative quantum field theories that arise as a low energy limit in certain string theories. The idea of noncommutativity comes naturally when combining these two extremes and has profound implications on results widely accepted in traditional, commutative, theories. In this work I review the status of one of the most important connections in physics, the spin-statistics relation. The relation is deeply ingrained in our reality in that it gives us the structure for the periodic table and is of crucial importance for the stability of all matter. The dramatic effects of noncommutativity of space-time coordinates, mainly the loss of Lorentz invariance, call the spin-statistics relation into question. The spin-statistics theorem is first presented in its traditional setting, giving a clarifying proof starting from minimal requirements. Next the notion of noncommutativity is introduced and its implications studied. The discussion is essentially based on twisted Poincaré symmetry, the space-time symmetry of noncommutative quantum field theory. The controversial issue of microcausality in noncommutative quantum field theory is settled by showing for the first time that the light wedge microcausality condition is compatible with the twisted Poincaré symmetry. The spin-statistics relation is considered both from the point of view of braided statistics, and in the traditional Lagrangian formulation of Pauli, with the conclusion that Pauli's age-old theorem stands even this test so dramatic for the whole structure of space-time.

Joseph A. Rosen Papers 1921-1938

With few exceptions, the bulk of the collection pertains to the work of the Agro-Joint. Records of the Agro-Joint Director General. Agreements of the American Relief Administration (ARA) and the Joint Distribution Committee with the Soviet government, 1922-1923. Agreements between the Agro-Joint and the Soviet government, 1924, 1927, 1928. Agreements of the Agro-Joint and the American Society for Jewish Farm Settlements (ASJFS) with the Soviet government, 1929, 1930, 1933, 1938. Materials relating to relief work of the JDC within the framework of the American Relief Administration, 1922, including the appointment of J. Rosen as the JDC representative at the ARA. Statistics, reports, miscellaneous correspondence relating to JDC activities in Russia. Minutes, memos, reports, legal documents, certificate of incorporation, and general correspondence relating to the ASJFS, its formation, fund-raising activities, 1927-1939. Records of the Agro-Joint Main Office, Moscow. Annual and periodi c reports of the Agro-Joint including statistics, financial estimates, financial reports, analyses of expenditures, relating to Agro-Joint work, 1924-1937. General correspondence files: incoming and outgoing letters, reports, and memoranda. Materials relating to land surveys and allocations in the Crimea: statistics, surveys, memos, correspondence, relating to the Salsk district, Chernomor district, Changar peninsula, Azov, Kuban, Odessa district, Samara district, Povolzhe, Krivoy Rog, Kherson, The Far East, Siberia. Materials relating to contacts with KOMZET. Correspondence, minutes of KOMZET meetings, statistical information, reports. By-laws of the OZET (Obshchestvo po Zemleustroystvu Trudyachtchikhsya Evreev - Association For the Settlement of Toiling Jews On Land) and AGRO-KUSTBANK (Evreysky Agrarno-Kustarny Bank - Jewish Agricultural and House Workers Bank). Register of Agro-Joint assets transferred to KOMZET. Records of the Agro-Joint Agricultural Department. Materials

Geospatial environmental data modelling applications using remote sensing, GIS and spatial statistics

This thesis presents novel modelling applications for environmental geospatial data using remote sensing, GIS and statistical modelling techniques. The studied themes can be classified into four main themes: (i) to develop advanced geospatial databases. Paper (I) demonstrates the creation of a geospatial database for the Glanville fritillary butterfly (Melitaea cinxia) in the Åland Islands, south-western Finland; (ii) to analyse species diversity and distribution using GIS techniques. Paper (II) presents a diversity and geographical distribution analysis for Scopulini moths at a world-wide scale; (iii) to study spatiotemporal forest cover change. Paper (III) presents a study of exotic and indigenous tree cover change detection in Taita Hills Kenya using airborne imagery and GIS analysis techniques; (iv) to explore predictive modelling techniques using geospatial data. In Paper (IV) human population occurrence and abundance in the Taita Hills highlands was predicted using the generalized additive modelling (GAM) technique. Paper (V) presents techniques to enhance fire prediction and burned area estimation at a regional scale in East Caprivi Namibia. Paper (VI) compares eight state-of-the-art predictive modelling methods to improve fire prediction, burned area estimation and fire risk mapping in East Caprivi Namibia. The results in Paper (I) showed that geospatial data can be managed effectively using advanced relational database management systems. Metapopulation data for Melitaea cinxia butterfly was successfully combined with GPS-delimited habitat patch information and climatic data. Using the geospatial database, spatial analyses were successfully conducted at habitat patch level or at more coarse analysis scales. Moreover, this study showed it appears evident that at a large-scale spatially correlated weather conditions are one of the primary causes of spatially correlated changes in Melitaea cinxia population sizes. In Paper (II) spatiotemporal characteristics of Socupulini moths description, diversity and distribution were analysed at a world-wide scale and for the first time GIS techniques were used for Scopulini moth geographical distribution analysis. This study revealed that Scopulini moths have a cosmopolitan distribution. The majority of the species have been described from the low latitudes, sub-Saharan Africa being the hot spot of species diversity. However, the taxonomical effort has been uneven among biogeographical regions. Paper III showed that forest cover change can be analysed in great detail using modern airborne imagery techniques and historical aerial photographs. However, when spatiotemporal forest cover change is studied care has to be taken in co-registration and image interpretation when historical black and white aerial photography is used. In Paper (IV) human population distribution and abundance could be modelled with fairly good results using geospatial predictors and non-Gaussian predictive modelling techniques. Moreover, land cover layer is not necessary needed as a predictor because first and second-order image texture measurements derived from satellite imagery had more power to explain the variation in dwelling unit occurrence and abundance. Paper V showed that generalized linear model (GLM) is a suitable technique for fire occurrence prediction and for burned area estimation. GLM based burned area estimations were found to be more superior than the existing MODIS burned area product (MCD45A1). However, spatial autocorrelation of fires has to be taken into account when using the GLM technique for fire occurrence prediction. Paper VI showed that novel statistical predictive modelling techniques can be used to improve fire prediction, burned area estimation and fire risk mapping at a regional scale. However, some noticeable variation between different predictive modelling techniques for fire occurrence prediction and burned area estimation existed.

Hazard Function Models to Estimate Mortality Rates Affecting Fish Populations with Application to the Sea Mullet (Mugil cephalus) Fishery on the Queensland Coast (Australia)

Fisheries management agencies around the world collect age data for the purpose of assessing the status of natural resources in their jurisdiction. Estimates of mortality rates represent a key information to assess the sustainability of fish stocks exploitation. Contrary to medical research or manufacturing where survival analysis is routinely applied to estimate failure rates, survival analysis has seldom been applied in fisheries stock assessment despite similar purposes between these fields of applied statistics. In this paper, we developed hazard functions to model the dynamic of an exploited fish population. These functions were used to estimate all parameters necessary for stock assessment (including natural and fishing mortality rates as well as gear selectivity) by maximum likelihood using age data from a sample of catch. This novel application of survival analysis to fisheries stock assessment was tested by Monte Carlo simulations to assert that it provided unbiased estimations of relevant quantities. The method was applied to the data from the Queensland (Australia) sea mullet (Mugil cephalus) commercial fishery collected between 2007 and 2014. It provided, for the first time, an estimate of natural mortality affecting this stock: 0.22±0.08 year −1 .