Faster k-Medoids Clustering: Improving the PAM, CLARA, and CLARANS Algorithms

Clustering non-Euclidean data is difficult, and one of the most used algorithms besides hierarchical clustering is the popular algorithm Partitioning Around Medoids (PAM), also simply referred to as k-medoids. In Euclidean geometry the mean-as used in k-means-is a good estimator for the cluster center, but this does not hold for arbitrary dissimilarities. PAM uses the medoid instead, the object with the smallest dissimilarity to all others in the cluster. This notion of centrality can be used with any (dis-)similarity, and thus is of high relevance to many domains such as biology that require the use of Jaccard, Gower, or more complex distances. A key issue with PAM is its high run time cost. We propose modifications to the PAM algorithm to achieve an O(k)-fold speedup in the second SWAP phase of the algorithm, but will still find the same results as the original PAM algorithm. If we slightly relax the choice of swaps performed (at comparable quality), we can further accelerate the algorithm by performing up to k swaps in each iteration. With the substantially faster SWAP, we can now also explore alternative strategies for choosing the initial medoids. We also show how the CLARA and CLARANS algorithms benefit from these modifications. It can easily be combined with earlier approaches to use PAM and CLARA on big data (some of which use PAM as a subroutine, hence can immediately benefit from these improvements), where the performance with high k becomes increasingly important. In experiments on real data with k=100, we observed a 200-fold speedup compared to the original PAM SWAP algorithm, making PAM applicable to larger data sets as long as we can afford to compute a distance matrix, and in particular to higher k (at k=2, the new SWAP was only 1.5 times faster, as the speedup is expected to increase with k)

Reconstructions of deltaic environments from Holocene palynological records in the Volga delta, northern Caspian Sea

This article was made available through open access by the Brunel Open Access Publishing Fund.New palynological and ostracod data are presented from the Holocene Volga delta, obtained from short cores and surface samples collected in the Damchik region, near Astrakhan, Russian Federation in the northern Caspian Sea. Four phases of delta deposition are recognized and constrained by accelerated mass spectrometry (AMS) radiocarbon ages. Palynological records show that erosive channels, dunes (Baer hills) and inter-dune lakes were present during the period 11,500–8900 cal. BP at the time of the Mangyshlak Caspian lowstand. The period 8900–3770 cal. BP was characterized regionally by extensive steppe vegetation, with forest present at times with warmer, more humid climates, and with halophytic and xerophytic vegetation present at times of drought. The period 3770–2080 cal. BP was a time of active delta deposition, with forest or woodland close to the delta, indicating relatively warm and humid climates and variable Caspian Sea levels. From 2080 cal. BP to the present-day, aquatic pollen is frequent in highstand intervals and herbaceous pollen and fungal hyphae frequent in lowstand intervals. Soils and incised valley sediments are associated with the regional Derbent regression and may be time-equivalent with the ‘Medieval Warm Period’. Fungal spores are an indicator of erosional or aeolian processes, whereas fungal hyphae are associated with soil formation. Freshwater algae, ostracods and dinocysts indicate mainly freshwater conditions during the Holocene with minor brackish influences. Dinocysts present include Spiniferites cruciformis, Caspidinium rugosum, Impagidinium caspienense and Pterocysta cruciformis, the latter a new record for the Caspian Sea. The Holocene Volga delta is a partial analogue for the much larger oil and gas bearing Mio-Pliocene palaeo-Volga delta.Funding for the data collection and field work was provided from the following sources: 1 – IGCP-UNESCO 2003–2008 (Project 481 CASPAGE, Dating Caspian Sea Level Change); 2 – NWO, Netherlands Science Foundation and RFFI, Russian Science Foundation 2005–2008 (Programme: ‘VHR Seismic Stratigraphy and Paleoecology of the Holocene Volga Delta’); and 3 – BP Exploration (Caspian Sea) Sea Ltd. (Azeri-Chirag-Gunashli) 2005–2008 (‘Unravelling the Small-Scale Stratigraphy and Sediment Dynamics of the Modern Volga Delta Using VHR Marine Geophysics’). The palynological work was funded jointly by BP Exploration (Caspian Sea) Ltd., Delft University of Technology and KrA Stratigraphic Ltd. Ostracod analyses were funded by StrataData Ltd. and funding for two additional radiocarbon dates provided by Deltares

Prognostic potential of body composition indices in detecting risk of musculoskeletal injury in army officer cadet profiles.

OBJECTIVES: High values in most of the body composition indices have been related to musculoskeletal injuries, but limited data exists on the accuracy of these diagnoses when detecting musculoskeletal injuries in military populations. METHODS: The suitability of body fat percentage, body mass index, fat mass index and fat free mass index to identify injury risk was examined in a group of army officer recruits. All body composition diagnoses were measured in 268 male army officer recruits prior to the commencement of basic combat training. Musculoskeletal injury was identified using codes from the International Classification of Diseases. The area under the curve, in the receiver operating characteristic curve, was used to quantify the overall ability to discriminate between those who were injured and those who were not. RESULTS: The statistics indicated that all indices, apart from body mass index, had a significant possibility to detect musculoskeletal injury potential (p 22, for fat mass index >6.5 and for fat free mass index <16.5. CONCLUSION: Body mass index values can not similarly detect the possibility of occurrence of musculoskeletal injuries in army officer recruits, just as other body composition diagnoses related to fat mass or/and free fat mass. However, the cut off-points related to the overall diagnostic performance of each body composition index should be used with caution and in accordance with the aims of each experimental setting

The influence of the atmospheric boundary layer on nocturnal layers of noctuids and other moths migrating over southern Britain

Insects migrating at high altitude over southern Britain have been continuously monitored by automatically-operating, vertical-looking radars over a period of several years. During some occasions in the summer months, the migrants were observed to form well-defined layer concentrations, typically at heights of 200-400 m, in the stable night-time atmosphere. Under these conditions, insects are likely to have control over their vertical movements and are selecting flight heights which are favourable for long-range migration. We therefore investigated the factors influencing the formation of these insect layers by comparing radar measurements of the vertical distribution of insect density with meteorological profiles generated by the UK Met. Office’s Unified Model (UM). Radar-derived measurements of mass and displacement speed, along with data from Rothamsted Insect Survey light traps provided information on the identity of the migrants. We present here three case studies where noctuid and pyralid moths contributed substantially to the observed layers. The major meteorological factors influencing the layer concentrations appeared to be: (a) the altitude of the warmest air, (b) heights corresponding to temperature preferences or thresholds for sustained migration and (c), on nights when air temperatures are relatively high, wind-speed maxima associated with the nocturnal jet. Back-trajectories indicated that layer duration may have been determined by the distance to the coast. Overall, the unique combination of meteorological data from the UM and insect data from entomological radar described here show considerable promise for systematic studies of high-altitude insect layering

The rise of fully turbulent flow

Over a century of research into the origin of turbulence in wallbounded shear flows has resulted in a puzzling picture in which turbulence appears in a variety of different states competing with laminar background flow. At slightly higher speeds the situation changes distinctly and the entire flow is turbulent. Neither the origin of the different states encountered during transition, nor their front dynamics, let alone the transformation to full turbulence could be explained to date. Combining experiments, theory and computer simulations here we uncover the bifurcation scenario organising the route to fully turbulent pipe flow and explain the front dynamics of the different states encountered in the process. Key to resolving this problem is the interpretation of the flow as a bistable system with nonlinear propagation (advection) of turbulent fronts. These findings bridge the gap between our understanding of the onset of turbulence and fully turbulent flows.Comment: 31 pages, 9 figure

Arthroscopic Treatment of Acetabular Retroversion With Acetabuloplasty and Subspine Decompression: A Matched Comparison With Patients Undergoing Arthroscopic Treatment for Focal Pincer-Type Femoroacetabular Impingement.

BackgroundGlobal acetabular retroversion is classically treated with open reverse periacetabular osteotomy. Given the low morbidity and recent success associated with the arthroscopic treatment of femoroacetabular impingement (FAI), there may also be a role for arthroscopic treatment of acetabular retroversion. However, the safety and outcomes after hip arthroscopic surgery for retroversion need further study, and the effect of impingement from the anterior inferior iliac spine (subspine) in patients with retroversion is currently unknown.HypothesisArthroscopic treatment for global acetabular retroversion will be safe, and patients will have similar outcomes compared with a matched group undergoing arthroscopic treatment for focal pincer-type FAI.Study designCohort study; Level of evidence, 2.MethodsPatients undergoing hip arthroscopic surgery for symptomatic global acetabular retroversion were prospectively enrolled and compared with a matched group of patients undergoing arthroscopic surgery for focal pincer-type FAI. Both groups underwent the same arthroscopic treatment protocol. All patients were administered patient-reported outcome (PRO) measures, including the 12-item Short-Form Health Survey (SF-12) Physical Component Summary (PCS) and a Mental Component Summary (MCS), modified Harris Hip Score (mHHS), Hip disability and Osteoarthritis Outcome Score (HOOS), and visual analog scale (VAS) for pain preoperatively and at 1 year postoperatively.ResultsThere were no differences in age, sex, or body mass index between 39 hips treated for global acetabular retroversion and 39 hips treated for focal pincer-type FAI. There were no major or minor complications in either group. Patients who underwent arthroscopic treatment for global acetabular retroversion demonstrated similar significant improvements in postoperative PRO scores (scores increased by 17 to 43 points) as patients who underwent arthroscopic treatment for focal pincer-type FAI. Patients treated for retroversion who also underwent subspine decompression had greater improvement than patients who did not undergo subspine decompression for the HOOS-Pain (33.7 ± 15.3 vs 22.5 ± 17.6, respectively; P = .046) and HOOS-Quality of Life (49.7 ± 18.8 vs 34.6 ± 22.0, respectively; P = .030) scores.ConclusionArthroscopic treatment for acetabular retroversion is safe and provides significant clinical improvement similar to arthroscopic treatment for pincer-type FAI. Patients with acetabular retroversion who also underwent arthroscopic subspine decompression demonstrated greater improvements in pain and quality of life outcomes than those who underwent arthroscopic treatment without subspine decompression

Tour guiding, organisational culture and learning: lessons from an entrepreneurial company

This paper examines the impacts of organisational culture on the learning and development of tour guides. Drawing on a case study of a small entrepreneurial tour company, the paper considers the nature of the organisation's culture, the tours it provides, including their narrative contents and the processes of organisational learning and socialisation. The paper suggests that the development of a learning culture within such an organisation may benefit from the provision of appropriate learning opportunities among the guides and facilitators who coordinate guide development

HIV Prevention in Care and Treatment Settings: Baseline Risk Behaviors among HIV Patients in Kenya, Namibia, and Tanzania.

HIV care and treatment settings provide an opportunity to reach people living with HIV/AIDS (PLHIV) with prevention messages and services. Population-based surveys in sub-Saharan Africa have identified HIV risk behaviors among PLHIV, yet data are limited regarding HIV risk behaviors of PLHIV in clinical care. This paper describes the baseline sociodemographic, HIV transmission risk behaviors, and clinical data of a study evaluating an HIV prevention intervention package for HIV care and treatment clinics in Africa. The study was a longitudinal group-randomized trial in 9 intervention clinics and 9 comparison clinics in Kenya, Namibia, and Tanzania (N = 3538). Baseline participants were mostly female, married, had less than a primary education, and were relatively recently diagnosed with HIV. Fifty-two percent of participants had a partner of negative or unknown status, 24% were not using condoms consistently, and 11% reported STI symptoms in the last 6 months. There were differences in demographic and HIV transmission risk variables by country, indicating the need to consider local context in designing studies and using caution when generalizing findings across African countries. Baseline data from this study indicate that participants were often engaging in HIV transmission risk behaviors, which supports the need for prevention with PLHIV (PwP). TRIAL REGISTRATION: ClinicalTrials.gov NCT01256463

Segregation of Mn, Si, Al, and oxygen during the friction stir welding of DH36 steel

This work investigates the role of welding speed in elemental segregation of Mn, Si, Al, and oxygen during friction stir welding (FSW) in DH36 steel. The experimental work undertaken showed that when the speed of the FSW process exceeds 500 RPM with a traverse speed of 400 mm/min, then elemental segregation of Mn, Si, Al, and O occurred. The mechanism of this segregation is not fully understood; additionally, the presence of oxygen within these segregated elements needs investigation. This work examines the elemental segregation within DH36 steel by conducting heat treatment experiments on unwelded samples incrementally in the range of 1200–1500 °C and at cooling rates similar to that in FSW process. The results of heat treatments were compared with samples welded under two extremes of weld tool speeds, namely W1 low tool speeds (200 RPM with traverse speed of 100 mm/min) and W2 high tool speeds (550 RPM with traverse speed of 400 mm/min). The results from the heat treatment trials showed that segregation commences when the temperature exceeds 1400 °C and Mn, Si, Al, and oxygen segregation progress occurs at 1450 °C and at a cooling rate associated with acicular ferrite formation. It was also found that high rotational speeds exceeding 500 RPM caused localized melting at the advancing-trailing side of the friction stir-welded samples. The study aims to estimate peak temperature limits at which elemental segregation does not occur and hence prevent their occurrence in practice by applying the findings to the tool’s rotational and traverse speed that correspond to the defined temperature

Complexity and Expressivity of Branching- and Alternating-Time Temporal Logics with Finitely Many Variables

We show that Branching-time temporal logics CTL and CTL*, as well as Alternating-time temporal logics ATL and ATL*, are as semantically expressive in the language with a single propositional variable as they are in the full language, i.e., with an unlimited supply of propositional variables. It follows that satisfiability for CTL, as well as for ATL, with a single variable is EXPTIME-complete, while satisfiability for CTL*, as well as for ATL*, with a single variable is 2EXPTIME-complete,--i.e., for these logics, the satisfiability for formulas with only one variable is as hard as satisfiability for arbitrary formulas.Comment: Prefinal version of the published pape