Fast Shortest Path Distance Estimation in Large Networks

We study the problem of preprocessing a large graph so that point-to-point shortest-path queries can be answered very fast. Computing shortest paths is a well studied problem, but exact algorithms do not scale to huge graphs encountered on the web, social networks, and other applications. In this paper we focus on approximate methods for distance estimation, in particular using landmark-based distance indexing. This approach involves selecting a subset of nodes as landmarks and computing (offline) the distances from each node in the graph to those landmarks. At runtime, when the distance between a pair of nodes is needed, we can estimate it quickly by combining the precomputed distances of the two nodes to the landmarks. We prove that selecting the optimal set of landmarks is an NP-hard problem, and thus heuristic solutions need to be employed. Given a budget of memory for the index, which translates directly into a budget of landmarks, different landmark selection strategies can yield dramatically different results in terms of accuracy. A number of simple methods that scale well to large graphs are therefore developed and experimentally compared. The simplest methods choose central nodes of the graph, while the more elaborate ones select central nodes that are also far away from one another. The efficiency of the suggested techniques is tested experimentally using five different real world graphs with millions of edges; for a given accuracy, they require as much as 250 times less space than the current approach in the literature which considers selecting landmarks at random. Finally, we study applications of our method in two problems arising naturally in large-scale networks, namely, social search and community detection.

On classical and other methods of discriminant analysis and estimation of log-odds

For two multinormal populations with equal covariance matrices the likelihood ratio discriminant function, an alternative allocation rule to the sample linear discriminant function when n1 ≠ n2 ,is studied analytically. With the assumption of a known covariance matrix its distribution is derived and the expectation of its actual and apparent error rates evaluated and compared with those of the sample linear discriminant function. This comparison indicates that the likelihood ratio allocation rule is robust to unequal sample sizes. The quadratic discriminant function is studied, its distribution reviewed and evaluation of its probabilities of misclassification discussed. For known covariance matrices the distribution of the sample quadratic discriminant function is derived. When the known covariance matrices are proportional exact expressions for the expectation of its actual and apparent error rates are obtained and evaluated. The effectiveness of the sample linear discriminant function for this case is also considered. Estimation of true log-odds for two multinormal populations with equal or unequal covariance matrices is studied. The estimative, Bayesian predictive and a kernel method are compared by evaluating their biases and mean square errors. Some algebraic expressions for these quantities are derived. With equal covariance matrices the predictive method is preferable. Where it derives this superiority is investigated by considering its performance for various levels of fixed true log-odds. It is also shown that the predictive method is sensitive to n1 ≠ n2. For unequal but proportional covariance matrices the unbiased estimative method is preferred. Product Normal kernel density estimates are used to give a kernel estimator of true log-odds. The effect of correlation in the variables with product kernels is considered. With equal covariance matrices the kernel and parametric estimators are compared by simulation. For moderately correlated variables and large dimension sizes the product kernel method is a good estimator of true log-odds.

An EORTC phase I study of epirubicin in combination with fixed doses of cyclophosphamide and infusional 5-fu (CEF-infu) as primary treatment of large operable or locally advanced/inflammatory breast cancer.

PURPOSE: The association of continuous infusion 5-fluorouracil, epirubicin (50 mg/m2 q 3 weeks) and a platinum compound (cisplatin or carboplatin) was found to be very active in patients with either locally advanced/inflammatory (LA/I) [1, 2] or large operable (LO) breast cancer (BC) [3]. The same rate of activity in terms of response rate (RR) and response duration was observed in LA/I BC patients when cisplatin was replaced by cyclophosphamide [4]. The dose of epirubicin was either 50 mg/m2 [ 1, 2, 3] or 60 mg/m2/cycle [4]. The main objective of this study was to determine the maximum tolerated dose (MTD) of epirubicin when given in combination with fixed doses of cyclophosphamide and infusional 5-fluorouracil (CEF-infu) as neoadjuvant therapy in patients with LO or LA/I BC for a maximum of 6 cycles. PATIENTS AND METHODS: Eligible patients had LO or LA/I BC, a performance status 0-1, adequate organ function and were <65 years old. Cyclophosphamide was administered at the dose of 400 mg/m2 day 1 and 8, q 4 weeks and infusional 5-fluorouracil 200 mg/m2/day was given day 1-28, q 4 weeks. Epirubicin was escalated from 30 to 45 and to 60 mg/m2 day 1 and 8; dose escalation was permitted if 0/3 or 1/6 patients experienced dose limiting toxicity (DLT) during the first 2 cycles of therapy. DLT for epirubicin was defined as febrile neutropenia, grade 4 neutropenia lasting for >7 days, grade 4 thrombocytopenia, or any non-haematological toxicity of CTC grade > or =3, excluding alopecia and plantar-palmar erythrodysesthesia (this toxicity was attributable to infusional 5-fluorouracil and was not considered a DLT of epirubicin). RESULTS: A total of 21 patients, median age 44 years (range 29-63) have been treated. 107 courses have been delivered, with a median number of 5 cycles per patient (range 4-6). DLTs on cycles I and 2 on level 1, 2, 3: grade 3 (G3) mucositis occurred in 1/10 patients treated at the third dose level. An interim analysis showed that G3 PPE occurred in 5/16 pts treated with the 28-day infusional 5-FU schedule at the 3 dose levels. The protocol was subsequently amended to limit the duration of infusional 5-fluorouracil infusion from 4 to 3 weeks. No G3 PPE was detected in 5 patients treated with this new schedule. CONCLUSIONS: This study establishes that epirubicin 60mg/m2 day 1 and 8, cyclophosphamide 400mg/m2 day 1 and 8 and infusional 5-fluorouracil 200 mg/m2/day day 1-21. q 4 weeks is the recommended dose level. Given the encouraging activity of this regimen (15/21 clinical responses) we have replaced infusional 5-fluorouracil by oral capecitabine in a recently activated study.

Uncovering compounds by synergy of cluster expansion and high-throughput methods.

Predicting from first-principles calculations whether mixed metallic elements phase-separate or form ordered structures is a major challenge of current materials research. It can be partially addressed in cases where experiments suggest the underlying lattice is conserved, using cluster expansion (CE) and a variety of exhaustive evaluation or genetic search algorithms. Evolutionary algorithms have been recently introduced to search for stable off-lattice structures at fixed mixture compositions. The general off-lattice problem is still unsolved. We present an integrated approach of CE and high-throughput ab initio calculations (HT) applicable to the full range of compositions in binary systems where the constituent elements or the intermediate ordered structures have different lattice types. The HT method replaces the search algorithms by direct calculation of a moderate number of naturally occurring prototypes representing all crystal systems and guides CE calculations of derivative structures. This synergy achieves the precision of the CE and the guiding strengths of the HT. Its application to poorly characterized binary Hf systems, believed to be phase-separating, defines three classes of alloys where CE and HT complement each other to uncover new ordered structures.

Elucidating solvent contributions to solution reactions with ab initio QM/MM methods.

Computer simulations of reaction processes in solution in general rely on the definition of a reaction coordinate and the determination of the thermodynamic changes of the system along the reaction coordinate. The reaction coordinate often is constituted of characteristic geometrical properties of the reactive solute species, while the contributions of solvent molecules are implicitly included in the thermodynamics of the solute degrees of freedoms. However, solvent dynamics can provide the driving force for the reaction process, and in such cases explicit description of the solvent contribution in the free energy of the reaction process becomes necessary. We report here a method that can be used to analyze the solvent contributions to the reaction activation free energies from the combined QM/MM minimum free-energy path simulations. The method was applied to the self-exchange S(N)2 reaction of CH(3)Cl + Cl(-), showing that the importance of solvent-solute interactions to the reaction process. The results were further discussed in the context of coupling between solvent and solute molecules in reaction processes.

Adaptive Mixture Modelling Metropolis Methods for Bayesian Analysis of Non-linear State-Space Models.

We describe a strategy for Markov chain Monte Carlo analysis of non-linear, non-Gaussian state-space models involving batch analysis for inference on dynamic, latent state variables and fixed model parameters. The key innovation is a Metropolis-Hastings method for the time series of state variables based on sequential approximation of filtering and smoothing densities using normal mixtures. These mixtures are propagated through the non-linearities using an accurate, local mixture approximation method, and we use a regenerating procedure to deal with potential degeneracy of mixture components. This provides accurate, direct approximations to sequential filtering and retrospective smoothing distributions, and hence a useful construction of global Metropolis proposal distributions for simulation of posteriors for the set of states. This analysis is embedded within a Gibbs sampler to include uncertain fixed parameters. We give an example motivated by an application in systems biology. Supplemental materials provide an example based on a stochastic volatility model as well as MATLAB code.

Validation of Thin Layer Chromatographic Methods

Thin-layer and high-performance thin-layer chromatography (TLC/HPTLC) methods for assaying compound(s) in a sample must be validated to ensure that they are fit for their intended purpose and, where applicable, meet the strict regulatory requirements for controlled products. Two validation approaches are identified in the literature, i.e. the classic and the alternative, which is using accuracy profiles.Detailed procedures of the two approaches are discussed based on the validation of methods for pharmaceutical analysis, which is an area considered having more strict requirements. Estimation of the measurement uncertainty from the validation approach using accuracy profiles is also described.Examples of HPTLC methods, developed and validated to assay sulfamethoxazole and trimethoprim on the one hand and lamivudine, stavudine, and nevirapine on the other, in their fixed-dose combination tablets, are further elaborated.

Nested heuristic methods for the location-routeing problem

The concept of 'nested methods' is adopted to solve the location-routeing problem. Unlike the sequential and iterative approaches, in this method we treat the routeing element as a sub-problem within the larger problem of location. Efficient techniques that take into account the above concept and which use a neighbourhood structure inspired from computational geometry are presented. A simple version of tabu search is also embedded into our methods to improve the solutions further. Computational testing is carried out on five sets of problems of 400 customers with five levels of depot fixed costs, and the results obtained are encouraging.

Assessing marine plankton community structure from long-term monitoring data with multivariate autoregressive (MAR) models: a comparison of fixed station versus spatially distributed sampling data

We examined how marine plankton interaction networks, as inferred by multivariate autoregressive (MAR) analysis of time-series, differ based on data collected at a fixed sampling location (L4 station in the Western English Channel) and four similar time-series prepared by averaging Continuous Plankton Recorder (CPR) datapoints in the region surrounding the fixed station. None of the plankton community structures suggested by the MAR models generated from the CPR datasets were well correlated with the MAR model for L4, but of the four CPR models, the one most closely resembling the L4 model was that for the CPR region nearest to L4. We infer that observation error and spatial variation in plankton community dynamics influenced the model performance for the CPR datasets. A modified MAR framework in which observation error and spatial variation are explicitly incorporated could allow the analysis to better handle the diverse time-series data collected in marine environments.

RNA expression analysis from formalin fixed paraffin embedded tissues

Formalin fixation and paraffin embedding (FFPE) is the most commonly used method worldwide for tissue storage. This method preserves the tissue integrity but causes extensive damage to nucleic acids stored within the tissue. As methods for measuring gene expression such as RT-PCR and microarray are adopted into clinical practice there is an increasing necessity to access the wealth of information locked in the Formalin fixation and paraffin embedding archives. This paper reviews the progress in this field and discusses the unique opportunities that exist for the application of these techniques in the development of personalized medicine.

A Wavelet-Prony Method for Modeling of Fixed-Speed Wind Farm Low-Frequency Power Pulsations

The increasing penetration of wind generation on the Island of Ireland has been accompanied by close investigation of low-frequency pulsations contained within active power flow. A primary concern is excitation of low-frequency oscillation modes already present on the system, particularly the 0.75 Hz mode as a consequence of interconnection between the Northern and Southern power system networks. In order to determine whether the prevalence of wind generation has a negative effect (excites modes) or positive impact (damping of modes) on the power system, oscillations must be measured and characterised. Using time – frequency methods, this paper presents work that has been conducted to extract features from low-frequency active power pulsations to determine the composition of oscillatory modes which may impact on dynamic stability. The paper proposes a combined wavelet-Prony method to extract modal components and determine damping factors. The method is exemplified using real data obtained from wind farm measurements.

Long-term occlusal and soft-tissue profile outcomes after treatment of Class II Division 1 malocclusion with fixed appliances

Introduction: Our objective was to determine which factors were predictive of good long-term outcomes after fixed appliance treatment of Class II Division 1 malocclusion. Methods: Two hundred seven patients with Class II Division 1 malocclusion were examined in early adulthood at a mean of 4.6 years after treatment with fixed appliances. The peer assessment rating index was used to evaluate dental alignment and occlusal relationships. The soft-tissue profile was assessed with the Holdaway angle. Results: Logistic regression identified 3 pretreatment variables that were predictive of a good facial profile (Holdaway angle) at recall: the lower lip to E-plane distance (P

Development and independent validation of a prognostic assay for stage ii colon cancer using formalin-fixed paraffin-embedded tissue

Purpose: Current prognostic factors are poor at identifying patients at risk of disease recurrence after surgery for stage II colon cancer. Here we describe a DNA microarray-based prognostic assay using clinically relevant formalin-fixed paraffin-embedded (FFPE) samples. Patients and Methods: A gene signature was developed from a balanced set of 73 patients with recurrent disease (high risk) and 142 patients with no recurrence (low risk) within 5 years of surgery. Results: The 634-probe set signature identified high-risk patients with a hazard ratio (HR) of 2.62 (P <.001) during cross validation of the training set. In an independent validation set of 144 samples, the signature identified high-risk patients with an HR of 2.53 (P <.001) for recurrence and an HR of 2.21 (P = .0084) for cancer-related death. Additionally, the signature was shown to perform independently from known prognostic factors (P <.001). Conclusion: This gene signature represents a novel prognostic biomarker for patients with stage II colon cancer that can be applied to FFPE tumor samples. © 2011 by American Society of Clinical Oncology.

Optimized bit level architectures for IIR filtering

Optimized circuits for implementing high-performance bit-parallel IIR filters are presented. Circuits constructed mainly from simple carry save adders and based on most-significant-bit (MSB) first arithmetic are described. Two methods resulting in systems which are 100% efficient in that they are capable of sampling data every cycle are presented. In the first approach the basic circuit is modified so that the level of pipelining used is compatible with the small, but fixed, latency associated with the computation in question. This is achieved through insertion of pipeline delays (half latches) on every second row of cells. This produces an area-efficient solution in which the throughput rate is determined by a critical path of 76 gate delays. A second approach combines the MSB first arithmetic methods with the scattered look-ahead methods. Important design issues are addressed, including wordlength truncation, overflow detection, and saturation.

Analysis of Fixed-Speed Wind Farm Low-Frequency Power Pulsations using a Wavelet-Prony Method

The increasing penetration of wind generation on the Island of Ireland has been accompanied by close investigation of low-frequency periodic pulsations contained within the active power flow from different wind farms. A primary concern is excitation of existing low-frequency oscillation modes already present on the system, particularly the 0.75 Hz mode as a consequence of the interconnected Northern and Southern power system networks. Recently grid code requirements on the Northern Ireland power system have been updated stipulating that wind farms connected after 2005 must be able to control the magnitude of oscillations in the range of 0.25 - 1.75 Hz to within 1% of the wind farm's registered output. In order to determine whether wind farm low-frequency oscillations have a negative effect (excite other modes) or possibly a positive impact (damping of existing modes) on the power system, the oscillations at the point of connection must be measured and characterised. Using time - frequency methods, research presented in this paper has been conducted to extract signal features from measured low-frequency active power pulsations produced by wind farms to determine the effective composition of possible oscillatory modes which may have a detrimental effect on system dynamic stability. The paper proposes a combined wavelet-Prony method to extract modal components and determine damping factors. The method is exemplified using real data obtained from wind farm measurements.