Computationally Efficient Methods for MDL-Optimal Density Estimation and Data Clustering

The Minimum Description Length (MDL) principle is a general, well-founded theoretical formalization of statistical modeling. The most important notion of MDL is the stochastic complexity, which can be interpreted as the shortest description length of a given sample of data relative to a model class. The exact definition of the stochastic complexity has gone through several evolutionary steps. The latest instantation is based on the so-called Normalized Maximum Likelihood (NML) distribution which has been shown to possess several important theoretical properties. However, the applications of this modern version of the MDL have been quite rare because of computational complexity problems, i.e., for discrete data, the definition of NML involves an exponential sum, and in the case of continuous data, a multi-dimensional integral usually infeasible to evaluate or even approximate accurately. In this doctoral dissertation, we present mathematical techniques for computing NML efficiently for some model families involving discrete data. We also show how these techniques can be used to apply MDL in two practical applications: histogram density estimation and clustering of multi-dimensional data.

Tools for compositional data with a total

Compositional data analysis usually deals with relative information between parts where the total (abundances, mass, amount, etc.) is unknown or uninformative. This article addresses the question of what to do when the total is known and is of interest. Tools used in this case are reviewed and analysed, in particular the relationship between the positive orthant of D-dimensional real space, the product space of the real line times the D-part simplex, and their Euclidean space structures. The first alternative corresponds to data analysis taking logarithms on each component, and the second one to treat a log-transformed total jointly with a composition describing the distribution of component amounts. Real data about total abundances of phytoplankton in an Australian river motivated the present study and are used for illustration.

Co-distribution and co-infection of chikungunya and dengue viruses

BACKGROUND Chikungunya and dengue infections are spatio-temporally related. The current review aims to determine the geographic limits of chikungunya, dengue and the principal mosquito vectors for both viruses and to synthesise current epidemiological understanding of their co-distribution. METHODS Three biomedical databases (PubMed, Scopus and Web of Science) were searched from their inception until May 2015 for studies that reported concurrent detection of chikungunya and dengue viruses in the same patient. Additionally, data from WHO, CDC and Healthmap alerts were extracted to create up-to-date global distribution maps for both dengue and chikungunya. RESULTS Evidence for chikungunya-dengue co-infection has been found in Angola, Gabon, India, Madagascar, Malaysia, Myanmar, Nigeria, Saint Martin, Singapore, Sri Lanka, Tanzania, Thailand and Yemen; these constitute only 13 out of the 98 countries/territories where both chikungunya and dengue epidemic/endemic transmission have been reported. CONCLUSIONS Understanding the true extent of chikungunya-dengue co-infection is hampered by current diagnosis largely based on their similar symptoms. Heightened awareness of chikungunya among the public and public health practitioners in the advent of the ongoing outbreak in the Americas can be expected to improve diagnostic rigour. Maps generated from the newly compiled lists of the geographic distribution of both pathogens and vectors represent the current geographical limits of chikungunya and dengue, as well as the countries/territories at risk of future incursion by both viruses. These describe regions of co-endemicity in which lab-based diagnosis of suspected cases is of higher priority.

What are the driving factors influencing the size distribution of airborne synthetic clay particles emitted from a jet milling process?

In the field of workplace air quality, measuring and analyzing the size distribution of airborne particles to identify their sources and apportion their contribution has become widely accepted, however, the driving factors that influence this parameter, particularly for nanoparticles (< 100 nm), have not been thoroughly determined. Identification of driving factors, and in turn, general trends in size distribution of emitted particles would facilitate the prediction of nanoparticles’ emission behavior and significantly contribute to their exposure assessment. In this study, a comprehensive analysis of the particle number size distribution data, with a particular focus on the ultrafine size range of synthetic clay particles emitted from a jet milling machine was conducted using the multi-lognormal fitting method. The results showed relatively high contribution of nanoparticles to the emissions in many of the tested cases, and also, that both surface treatment and feed rate of the machine are significant factors influencing the size distribution of the emitted particles of this size. In particular, applying surface treatments and increasing the machine feed rate have the similar effect of reducing the size of the particles, however, no general trend was found in variations of size distribution across different surface treatments and feed rates. The findings of our study demonstrate that for this process and other activities, where no general trend is found in the size distribution of the emitted airborne particles due to dissimilar effects of the driving factors, each case must be treated separately in terms of workplace exposure assessment and regulations.

Reputation model based on rating data and application in recommender systems

This thesis introduced two novel reputation models to generate accurate item reputation scores using ratings data and the statistics of the dataset. It also presented an innovative method that incorporates reputation awareness in recommender systems by employing voting system methods to produce more accurate top-N item recommendations. Additionally, this thesis introduced a personalisation method for generating reputation scores based on users' interests, where a single item can have different reputation scores for different users. The personalised reputation scores are then used in the proposed reputation-aware recommender systems to enhance the recommendation quality.

LRP5 regulates human body fat distribution by modulating adipose progenitor biology in a dose- and depot-specific fashion

Summary Common variants in WNT pathway genes have been associated with bone mass and fat distribution, the latter predicting diabetes and cardiovascular disease risk. Rare mutations in the WNT co-receptors LRP5 and LRP6 are similarly associated with bone and cardiometabolic disorders. We investigated the role of LRP5 in human adipose tissue. Subjects with gain-of-function LRP5 mutations and high bone mass had enhanced lower-body fat accumulation. Reciprocally, a low bone mineral density-associated common LRP5 allele correlated with increased abdominal adiposity. Ex vivo LRP5 expression was higher in abdominal versus gluteal adipocyte progenitors. Equivalent knockdown of LRP5 in both progenitor types dose-dependently impaired β-catenin signaling and led to distinct biological outcomes: diminished gluteal and enhanced abdominal adipogenesis. These data highlight how depot differences in WNT/β-catenin pathway activity modulate human fat distribution via effects on adipocyte progenitor biology. They also identify LRP5 as a potential pharmacologic target for the treatment of cardiometabolic disorders. © 2015 The Authors.

Phospholipids of lipid-containing bacteriophages and their transbilayer distribution

In this study we used electro-spray ionization mass-spectrometry to determine phospholipid class and molecular species compositions in bacteriophages PM2, PRD1, Bam35 and phi6 as well as their hosts. To obtain compositional data of the individual leaflets, phospholipid transbilayer distribution in the viral membranes was studied. We found that 1) the membranes of all studied bacteriophage are enriched in PG as compared to the host membranes, 2) molecular species compositions in the phage and host membranes are similar, and 3) phospholipids in the viral membranes are distributed asymmetrically with phosphatidylglycerol enriched in the outer leaflet and phosphatidylethanolamine in the inner one (except Bam35). Alternative models for selective incorporation of phospholipids to phages and for the origins of the asymmetric phospholipid transbilayer distribution are discussed. Notably, the present data are also useful when constructing high resolution structural models of bacteriophages, since diffraction methods cannot provide a detailed structure of the membrane due to high motility of the lipids and lack of symmetric organization of membrane proteins.

Distribution of non-biting midges (Diptera, Chironomidae) in subarctic lakes of Finnish Lapland applications in lake classification and palaeolimnology

Climate change contributes directly or indirectly to changes in species distributions, and there is very high confidence that recent climate warming is already affecting ecosystems. The Arctic has already experienced the greatest regional warming in recent decades, and the trend is continuing. However, studies on the northern ecosystems are scarce compared to more southerly regions. Better understanding of the past and present environmental change is needed to be able to forecast the future. Multivariate methods were used to explore the distributional patterns of chironomids in 50 shallow (≤ 10m) lakes in relation to 24 variables determined in northern Fennoscandia at the ecotonal area from the boreal forest in the south to the orohemiarctic zone in the north. Highest taxon richness was noted at middle elevations around 400 m a.s.l. Significantly lower values were observed from cold lakes situated in the tundra zone. Lake water alkalinity had the strongest positive correlation with the taxon richness. Many taxa had preference for lakes either on tundra area or forested area. The variation in the chironomid abundance data was best correlated with sediment organic content (LOI), lake water total organic carbon content, pH and air temperature, with LOI being the strongest variable. Three major lake groups were separated on the basis of their chironomid assemblages: (i) small and shallow organic-rich lakes, (ii) large and base-rich lakes, and (iii) cold and clear oligotrophic tundra lakes. Environmental variables best discriminating the lake groups were LOI, taxon richness, and Mg. When repeated, this kind of an approach could be useful and efficient in monitoring the effects of global change on species ranges. Many species of fast spreading insects, including chironomids, show a remarkable ability to track environmental changes. Based on this ability, past environmental conditions have been reconstructed using their chitinous remains in the lake sediment profiles. In order to study the Holocene environmental history of subarctic aquatic systems, and quantitatively reconstruct the past temperatures at or near the treeline, long sediment cores covering the last 10000 years (the Holocene) were collected from three lakes. Lower temperature values than expected based on the presence of pine in the catchment during the mid-Holocene were reconstructed from a lake with great water volume and depth. The lake provided thermal refuge for profundal, cold adapted taxa during the warm period. In a shallow lake, the decrease in the reconstructed temperatures during the late Holocene may reflect the indirect response of the midges to climate change through, e.g., pH change. The results from three lakes indicated that the response of chironomids to climate have been more or less indirect. However, concurrent shifts in assemblages of chironomids and vegetation in two lakes during the Holocene time period indicated that the midges together with the terrestrial vegetation had responded to the same ultimate cause, which most likely was the Holocene climate change. This was also supported by the similarity in the long-term trends in faunal succession for the chironomid assemblages in several lakes in the area. In northern Finnish Lapland the distribution of chironomids were significantly correlated with physical and limnological factors that are most likely to change as a result of future climate change. The indirect and individualistic response of aquatic systems, as reconstructed using the chironomid assemblages, to the climate change in the past suggests that in the future, the lake ecosystems in the north do not respond in one predictable way to the global climate change. Lakes in the north may respond to global climate change in various ways that are dependent on the initial characters of the catchment area and the lake.

Vertical distribution of aerosols over the east coast of India inferred from airborne LIDAR measurements

The information on altitude distribution of aerosols in the atmosphere is essential in assessing the impact of aerosol warming on thermal structure and stability of the atmosphere.In addition, aerosol altitude distribution is needed to address complex problems such as the radiative interaction of aerosols in the presence of clouds. With this objective,an extensive, multi-institutional and multi-platform field experiment (ICARB-Integrated Campaign for Aerosols, gases and Radiation Budget) was carried out under the Geosphere Biosphere Programme of the Indian Space Research Organization (ISRO-GBP) over continental India and adjoining oceans during March to May 2006. Here, we present airborne LIDAR measurements carried out over the east Coast of the India during the ICARB field campaign. An increase in aerosol extinction (scattering + absorption) was observed from the surface upwards with a maximum around 2 to 4 km. Aerosol extinction at higher atmospheric layers (>2 km) was two to three times larger compared to that of the surface. A large fraction (75-85%) of aerosol column optical depth was contributed by aerosols located above 1 km. The aerosol layer heights (defined in this paper as the height at which the gradient in extinction coefficient changes sign) showed a gradual decrease with an increase in the offshore distance. A large fraction (60-75%) of aerosol was found located above clouds indicating enhanced aerosol absorption above clouds. Our study implies that a detailed statistical evaluation of the temporal frequency and spatial extent of elevated aerosol layers is necessary to assess their significance to the climate. This is feasible using data from space-borne lidars such as CALIPSO,which fly in formation with other satellites like MODIS AQUA and MISR, as part of the A-Train constellation.

Continuous Distribution Kinetics for Photopolymerization of Alkyl Methacrylates

The photopolymerization of methyl,ethyl,butyl, and hexyl methacrylates in solution was studied. The effect of initial initiator and monomer concentrations on the time evolution of polymer concentration (M) over bar (n) and PDI was examined. The reversible chain addition and beta-scission, and primary radical termination steps were included in the mechanism along with the classical steps. The rate equations were derived using continuous distribution kinetics and solved numerically to fit the experimental data. The regressed rate coefficients compared well with the literature data. The model predicted the instantaneous increase in (M) over bar (n) and PDI to steady state values. The rate coefficients exhibited a linear increase with the size of alkyl chain of the alkyl methacrylates.

The reproducibility of acquiring three dimensional gait and plantar pressure data using established protocols in participants with and without type 2 diabetes and foot ulcers

Background Several prospective studies have suggested that gait and plantar pressure abnormalities secondary to diabetic peripheral neuropathy contributes to foot ulceration. There are many different methods by which gait and plantar pressures are assessed and currently there is no agreed standardised approach. This study aimed to describe the methods and reproducibility of three-dimensional gait and plantar pressure assessments in a small subset of participants using pre-existing protocols. Methods Fourteen participants were conveniently sampled prior to a planned longitudinal study; four patients with diabetes and plantar foot ulcers, five patients with diabetes but no foot ulcers and five healthy controls. The repeatability of measuring key biomechanical data was assessed including the identification of 16 key anatomical landmarks, the measurement of seven leg dimensions, the processing of 22 three-dimensional gait parameters and the analysis of four different plantar pressures measures at 20 foot regions. Results The mean inter-observer differences were within the pre-defined acceptable level (<7 mm) for 100 % (16 of 16) of key anatomical landmarks measured for gait analysis. The intra-observer assessment concordance correlation coefficients were > 0.9 for 100 % (7 of 7) of leg dimensions. The coefficients of variations (CVs) were within the pre-defined acceptable level (<10 %) for 100 % (22 of 22) of gait parameters. The CVs were within the pre-defined acceptable level (<30 %) for 95 % (19 of 20) of the contact area measures, 85 % (17 of 20) of mean plantar pressures, 70 % (14 of 20) of pressure time integrals and 55 % (11 of 20) of maximum sensor plantar pressure measures. Conclusion Overall, the findings of this study suggest that important gait and plantar pressure measurements can be reliably acquired. Nearly all measures contributing to three-dimensional gait parameter assessments were within predefined acceptable limits. Most plantar pressure measurements were also within predefined acceptable limits; however, reproducibility was not as good for assessment of the maximum sensor pressure. To our knowledge, this is the first study to investigate the reproducibility of several biomechanical methods in a heterogeneous cohort.

Adsorption of Nitrogen on Activated Carbon-Refit of Experimental Data and Derivation of Properties Required for Design of Equipment

Volumetric method based adsorption measurements of nitrogen on two specimens of activated carbon (Fluka and Sarabhai) reported by us are refitted to two popular isotherms, namely, Dubunin−Astakhov (D−A) and Toth, in light of improved fitting methods derived recently. Those isotherms have been used to derive other data of relevance in design of engineering equipment such as the concentration dependence of heat of adsorption and Henry’s law coefficients. The present fits provide a better representation of experimental measurements than before because the temperature dependence of adsorbed phase volume and structural heterogeneity of micropore distribution have been accounted for in the D−A equation. A new correlation to the Toth equation is a further contribution. The heat of adsorption in the limiting uptake condition is correlated with the Henry’s law coefficients at the near zero uptake condition.

Data-aware task scheduling for all-to-all comparison problems in heterogeneous distributed systems

Solving large-scale all-to-all comparison problems using distributed computing is increasingly significant for various applications. Previous efforts to implement distributed all-to-all comparison frameworks have treated the two phases of data distribution and comparison task scheduling separately. This leads to high storage demands as well as poor data locality for the comparison tasks, thus creating a need to redistribute the data at runtime. Furthermore, most previous methods have been developed for homogeneous computing environments, so their overall performance is degraded even further when they are used in heterogeneous distributed systems. To tackle these challenges, this paper presents a data-aware task scheduling approach for solving all-to-all comparison problems in heterogeneous distributed systems. The approach formulates the requirements for data distribution and comparison task scheduling simultaneously as a constrained optimization problem. Then, metaheuristic data pre-scheduling and dynamic task scheduling strategies are developed along with an algorithmic implementation to solve the problem. The approach provides perfect data locality for all comparison tasks, avoiding rearrangement of data at runtime. It achieves load balancing among heterogeneous computing nodes, thus enhancing the overall computation time. It also reduces data storage requirements across the network. The effectiveness of the approach is demonstrated through experimental studies.

Aerosol size distribution and its connection to cloud droplet number concentration

In order to predict the current state and future development of Earth s climate, detailed information on atmospheric aerosols and aerosol-cloud-interactions is required. Furthermore, these interactions need to be expressed in such a way that they can be represented in large-scale climate models. The largest uncertainties in the estimate of radiative forcing on the present day climate are related to the direct and indirect effects of aerosol. In this work aerosol properties were studied at Pallas and Utö in Finland, and at Mount Waliguan in Western China. Approximately two years of data from each site were analyzed. In addition to this, data from two intensive measurement campaigns at Pallas were used. The measurements at Mount Waliguan were the first long term aerosol particle number concentration and size distribution measurements conducted in this region. They revealed that the number concentration of aerosol particles at Mount Waliguan were much higher than those measured at similar altitudes in other parts of the world. The particles were concentrated in the Aitken size range indicating that they were produced within a couple of days prior to reaching the site, rather than being transported over thousands of kilometers. Aerosol partitioning between cloud droplets and cloud interstitial particles was studied at Pallas during the two measurement campaigns, First Pallas Cloud Experiment (First PaCE) and Second Pallas Cloud Experiment (Second PaCE). The method of using two differential mobility particle sizers (DMPS) to calculate the number concentration of activated particles was found to agree well with direct measurements of cloud droplet. Several parameters important in cloud droplet activation were found to depend strongly on the air mass history. The effects of these parameters partially cancelled out each other. Aerosol number-to-volume concentration ratio was studied at all three sites using data sets with long time-series. The ratio was found to vary more than in earlier studies, but less than either aerosol particle number concentration or volume concentration alone. Both air mass dependency and seasonal pattern were found at Pallas and Utö, but only seasonal pattern at Mount Waliguan. The number-to-volume concentration ratio was found to follow the seasonal temperature pattern well at all three sites. A new parameterization for partitioning between cloud droplets and cloud interstitial particles was developed. The parameterization uses aerosol particle number-to-volume concentration ratio and aerosol particle volume concentration as the only information on the aerosol number and size distribution. The new parameterization is computationally more efficient than the more detailed parameterizations currently in use, but the accuracy of the new parameterization was slightly lower. The new parameterization was also compared to directly observed cloud droplet number concentration data, and a good agreement was found.

Prevalence, codetection and seasonal distribution of upper airway viruses and bacteria in children with acute respiratory illnesses with cough as a symptom

Most studies exploring the role of upper airway viruses and bacteria in paediatric acute respiratory infections (ARI) focus on specific clinicaldiagnoses and/or do not account for virus–bacteria interactions. We aimed to describe the frequency and predictors of virus and bacteria codetection in children with ARI and cough, irrespective of clinical diagnosis. Bilateral nasal swabs, demographic, clinical and risk factor data were collected at enrollment in children aged <15 years presenting to an emergency department with an ARI and where cough was a symptom. Swabs were tested by polymerase chain reaction for 17 respiratory viruses and seven respiratory bacteria. Logistic regression was used to investigate associations between child characteristics and codetection of the organisms of interest. Between December 2011 and August 2014, swabs were collected from 817 (93.3%) of 876 enrolled children, median age 27.7 months (interquartile range13.9–60.3 months). Overall, 740 (90.6%) of 817 specimens were positive for any organism. Both viruses and bacteria were detected in 423 specimens (51.8%). Factors associated with codetection were age (adjusted odds ratio (aOR) for age <12 months = 4.9, 95% confidence interval (CI) 3.0, 7.9; age 12 to <24 months = 6.0, 95% CI 3.7, 9.8; age 24 to <60 months = 2.4, 95% CI 1.5, 3.9), male gender (aOR 1.46; 95% CI 1.1, 2.0), child care attendance (aOR 2.0; 95% CI 1.4, 2.8) and winter enrollment (aOR 2.0; 95% CI 1.3, 3.0). Haemophilus influenzae dominated the virus–bacteria pairs. Virus–H. influenzae interactions in ARI should be investigated further, especially as the contribution of nontypeable H. influenzae to acute and chronic respiratory diseases is being increasingly recognized.