• ### A Deep Active Survival Analysis Approach for Precision Treatment Recommendations: Application of Prostate Cancer(1804.03280)

April 10, 2018 cs.LG, cs.CY, stat.ML
Survival analysis has been developed and applied in the number of areas including manufacturing, finance, economics and healthcare. In healthcare domain, usually clinical data are high-dimensional, sparse and complex and sometimes there exists few amount of time-to-event (labeled) instances. Therefore building an accurate survival model from electronic health records is challenging. With this motivation, we address this issue and provide a new survival analysis framework using deep learning and active learning with a novel sampling strategy. First, our approach provides better representation with lower dimensions from clinical features using labeled (time-to-event) and unlabeled (censored) instances and then actively trains the survival model by labeling the censored data using an oracle. As a clinical assistive tool, we introduce a simple effective treatment recommendation approach based on our survival model. In the experimental study, we apply our approach on SEER-Medicare data related to prostate cancer among African-Americans and white patients. The results indicate that our approach outperforms significantly than baseline models.
• ### Deep Stock Representation Learning: From Candlestick Charts to Investment Decisions(1709.03803)

Feb. 18, 2018 q-fin.CP
We propose a novel investment decision strategy (IDS) based on deep learning. The performance of many IDSs is affected by stock similarity. Most existing stock similarity measurements have the problems: (a) The linear nature of many measurements cannot capture nonlinear stock dynamics; (b) The estimation of many similarity metrics (e.g. covariance) needs very long period historic data (e.g. 3K days) which cannot represent current market effectively; (c) They cannot capture translation-invariance. To solve these problems, we apply Convolutional AutoEncoder to learn a stock representation, based on which we propose a novel portfolio construction strategy by: (i) using the deeply learned representation and modularity optimisation to cluster stocks and identify diverse sectors, (ii) picking stocks within each cluster according to their Sharpe ratio (Sharpe 1994). Overall this strategy provides low-risk high-return portfolios. We use the Financial Times Stock Exchange 100 Index (FTSE 100) data for evaluation. Results show our portfolio outperforms FTSE 100 index and many well known funds in terms of total return in 2000 trading days.
• ### Cross sections for inelastic meson-meson scattering(1708.03062)

Feb. 17, 2018 hep-ph, nucl-th
We study two kinds of inelastic meson-meson scattering. The first kind is inelastic 2-to-2 meson-meson scattering that is governed by quark interchange as well as quark-antiquark annihilation and creation. Cross-section formulas are provided to get unpolarized cross sections for $\pi K \to \rho K^\ast$ for $I=1/2$, $\pi K^\ast \to \rho K$ for $I=1/2$, $\pi K^\ast \to \rho K^\ast$ for $I=1/2$, and $\rho K \to \rho K^\ast$ for $I=1/2$. Near threshold, quark interchange dominates the reactions near the critical temperature. The second kind is 2-to-1 meson-meson scattering with the process that a quark in an initial meson and an antiquark in another initial meson annihilate into a gluon and subsequently the gluon is absorbed by the other antiquark or quark. The transition potential for the process is derived. Four Feynman diagrams at tree level contribute to the 2-to-1 meson-meson scattering. Starting from the $S$-matrix element, the isospin-averaged unpolarized cross section with transition amplitudes is derived. The cross sections for $\pi \pi \to \rho$ and $\pi K \to K^*$ decrease with increasing temperature.
• ### Blind Demixing for Low-Latency Communication(1801.02158)

Jan. 7, 2018 cs.IT, math.IT
In the next generation wireless networks, lowlatency communication is critical to support emerging diversified applications, e.g., Tactile Internet and Virtual Reality. In this paper, a novel blind demixing approach is developed to reduce the channel signaling overhead, thereby supporting low-latency communication. Specifically, we develop a low-rank approach to recover the original information only based on a single observed vector without any channel estimation. Unfortunately, this problem turns out to be a highly intractable non-convex optimization problem due to the multiple non-convex rankone constraints. To address the unique challenges, the quotient manifold geometry of product of complex asymmetric rankone matrices is exploited by equivalently reformulating original complex asymmetric matrices to the Hermitian positive semidefinite matrices. We further generalize the geometric concepts of the complex product manifolds via element-wise extension of the geometric concepts of the individual manifolds. A scalable Riemannian trust-region algorithm is then developed to solve the blind demixing problem efficiently with fast convergence rates and low iteration cost. Numerical results will demonstrate the algorithmic advantages and admirable performance of the proposed algorithm compared with the state-of-art methods.
• ### A Predictive Approach Using Deep Feature Learning for Electronic Medical Records: A Comparative Study(1801.02961)

Jan. 6, 2018 cs.LG, stat.ML
Massive amount of electronic medical records accumulating from patients and populations motivates clinicians and data scientists to collaborate for the advanced analytics to extract knowledge that is essential to address the extensive personalized insights needed for patients, clinicians, providers, scientists, and health policy makers. In this paper, we propose a new predictive approach based on feature representation using deep feature learning and word embedding techniques. Our method uses different deep architectures for feature representation in higher-level abstraction to obtain effective and more robust features from EMRs, and then build prediction models on the top of them. Our approach is particularly useful when the unlabeled data is abundant whereas labeled one is scarce. We investigate the performance of representation learning through a supervised approach. First, we apply our method on a small dataset related to a specific precision medicine problem, which focuses on prediction of left ventricular mass indexed to body surface area (LVMI) as an indicator of heart damage risk in a vulnerable demographic subgroup (African-Americans). Then we use two large datasets from eICU collaborative research database to predict the length of stay in Cardiac-ICU and Neuro-ICU based on high dimensional features. Finally we provide a comparative study and show that our predictive approach leads to better results in comparison with others.
• ### SUBIC: A Supervised Bi-Clustering Approach for Precision Medicine(1709.09929)

Sept. 26, 2017 cs.LG, stat.ML
Traditional medicine typically applies one-size-fits-all treatment for the entire patient population whereas precision medicine develops tailored treatment schemes for different patient subgroups. The fact that some factors may be more significant for a specific patient subgroup motivates clinicians and medical researchers to develop new approaches to subgroup detection and analysis, which is an effective strategy to personalize treatment. In this study, we propose a novel patient subgroup detection method, called Supervised Biclustring (SUBIC) using convex optimization and apply our approach to detect patient subgroups and prioritize risk factors for hypertension (HTN) in a vulnerable demographic subgroup (African-American). Our approach not only finds patient subgroups with guidance of a clinically relevant target variable but also identifies and prioritizes risk factors by pursuing sparsity of the input variables and encouraging similarity among the input variables and between the input and target variables
• ### Millimeter Wave Communications for Future Mobile Networks(1705.06072)

May 17, 2017 cs.IT, math.IT
Millimeter wave (mmWave) communications have recently attracted large research interest, since the huge available bandwidth can potentially lead to rates of multiple Gbps (gigabit per second) per user. Though mmWave can be readily used in stationary scenarios such as indoor hotspots or backhaul, it is challenging to use mmWave in mobile networks, where the transmitting/receiving nodes may be moving, channels may have a complicated structure, and the coordination among multiple nodes is difficult. To fully exploit the high potential rates of mmWave in mobile networks, lots of technical problems must be addressed. This paper presents a comprehensive survey of mmWave communications for future mobile networks (5G and beyond). We first summarize the recent channel measurement campaigns and modeling results. Then, we discuss in detail recent progresses in multiple input multiple output (MIMO) transceiver design for mmWave communications. After that, we provide an overview of the solution for multiple access and backhauling, followed by analysis of coverage and connectivity. Finally, the progresses in the standardization and deployment of mmWave for mobile networks are discussed.
• ### Cooling-Rate Effects in Sodium Silicate Glasses: Bridging the Gap between Molecular Dynamics Simulations and Experiments(1704.08209)

Although molecular dynamics (MD) simulations are commonly used to predict the structure and properties of glasses, they are intrinsically limited to short time scales, necessitating the use of fast cooling rates. It is therefore challenging to compare results from MD simulations to experimental results for glasses cooled on typical laboratory time scales. Based on MD simulations of a sodium silicate glass with varying cooling rate (from 0.01 to 100 K/ps), here we show that thermal history primarily affects the medium-range order structure, while the short-range order is largely unaffected over the range of cooling rates simulated. This results in a decoupling between the enthalpy and volume relaxation functions, where the enthalpy quickly plateaus as the cooling rate decreases, whereas density exhibits a slower relaxation. Finally, we demonstrate that the outcomes of MD simulations can be meaningfully compared to experimental values if properly extrapolated to slower cooling rates.
• ### SAFS: A Deep Feature Selection Approach for Precision Medicine(1704.05960)

April 20, 2017 cs.LG, stat.ML
In this paper, we propose a new deep feature selection method based on deep architecture. Our method uses stacked auto-encoders for feature representation in higher-level abstraction. We developed and applied a novel feature learning approach to a specific precision medicine problem, which focuses on assessing and prioritizing risk factors for hypertension (HTN) in a vulnerable demographic subgroup (African-American). Our approach is to use deep learning to identify significant risk factors affecting left ventricular mass indexed to body surface area (LVMI) as an indicator of heart damage risk. The results show that our feature learning and representation approach leads to better results in comparison with others.
• ### Sparse Hierarchical Solvers with Guaranteed Convergence(1611.03189)

March 13, 2017 math.NA
Solving sparse linear systems from discretized PDEs is challenging. Direct solvers have in many cases quadratic complexity (depending on geometry), while iterative solvers require problem dependent preconditioners to be robust and efficient. Approximate factorization preconditioners, such as incomplete LU factorization, provide cheap approximations to the system matrix. However, even a highly accurate preconditioner may have deteriorating performance when the condition number of the system matrix increases. By increasing the accuracy on low-frequency errors, we propose a novel hierarchical solver with improved robustness with respect to the condition number of the linear system. This solver retains the linear computational cost and memory footprint of the original algorithm.
• ### Interaction of two filaments in a long filament channel associated with twin coronal mass ejections(1701.05122)

Jan. 18, 2017 astro-ph.SR
Using the high-quality observations of the Solar Dynamics Observatory, we present the interaction of two filaments (F1 and F2) in a long filament channel associated with twin coronal mass ejections (CMEs) on 2016 January 26. Before the eruption, a sequence of rapid cancellation and emergence of the magnetic flux has been observed, which likely triggered the ascending of the west filament (F1). The east footpoints of rising F1 moved toward the east far end of the filament channel, accompanying with post-eruption loops and flare ribbons. It likely indicated a large-scale eruption involving the long filament channel, resulted from the interaction between F1 and the east filament (F2). Some bright plasma flew over F2, and F2 stayed at rest during the eruption, likely due to the confinement of its overlying lower magnetic field. Interestingly, the impulsive F1 pushed its overlying magnetic arcades to form the first CME, and F1 finally evolved into the second CME after the collision with the nearby coronal hole. We suggest that the interaction of F1 and the overlying magnetic field of F2 led to the merging reconnection that form a longer eruptive filament loop. Our results also provide a possible picture of the origin of twin CMEs, and show the large-scale magnetic topology of the coronal hole is important for the eventual propagation direction of CMEs.
• ### Autoencoder Regularized Network For Driving Style Representation Learning(1701.01272)

Jan. 5, 2017 cs.AI, cs.NE, cs.CV
In this paper, we study learning generalized driving style representations from automobile GPS trip data. We propose a novel Autoencoder Regularized deep neural Network (ARNet) and a trip encoding framework trip2vec to learn drivers' driving styles directly from GPS records, by combining supervised and unsupervised feature learning in a unified architecture. Experiments on a challenging driver number estimation problem and the driver identification problem show that ARNet can learn a good generalized driving style representation: It significantly outperforms existing methods and alternative architectures by reaching the least estimation error on average (0.68, less than one driver) and the highest identification accuracy (by at least 3% improvement) compared with traditional supervised learning methods.
• ### Dynamic Spectrum Leasing with Two Sellers(1612.05702)

Dec. 17, 2016 cs.GT
This paper studies dynamic spectrum leasing in a cognitive radio network. There are two spectrum sellers, who are two primary networks, each with an amount of licensed spectrum bandwidth. When a seller has some unused spectrum, it would like to lease the unused spectrum to secondary users. A coordinator helps to perform the spectrum leasing stage-by-stage. As the two sellers may have different leasing period, there are three epochs, in which seller 1 has spectrum to lease in Epochs II and III, while seller 2 has spectrum to lease in Epochs I and II. Each seller needs to decide how much spectrum it should lease to secondary users in each stage of its leasing period, with a target at revenue maximization. It is shown that, when the two sellers both have spectrum to lease (i.e., in Epoch II), the spectrum leasing can be formulated as a non-cooperative game. Nash equilibria of the game are found in closed form. Solutions of the two users in the three epochs are derived.
• ### Low-Rank Matrix Completion for Mobile Edge Caching in Fog-RAN via Riemannian Optimization(1608.07800)

Sept. 4, 2016 cs.IT, math.IT
The upcoming big data era is likely to demand tremendous computation and storage resources for communications. By pushing computation and storage to network edges, fog radio access networks (Fog-RAN) can effectively increase network throughput and reduce transmission latency. Furthermore, we can exploit the benefits of cache enabled architecture in Fog-RAN to deliver contents with low latency. Radio access units (RAUs) need content delivery from fog servers through wireline links whereas multiple mobile devices acquire contents from RAUs wirelessly. This work proposes a unified low-rank matrix completion (LRMC) approach to solving the content delivery problem in both wireline and wireless parts of Fog-RAN. To attain a low caching latency, we present a high precision approach with Riemannian trust-region method to solve the challenging LRMC problem by exploiting the quotient manifold geometry of fixed-rank matrices. Numerical results show that the new approach has a faster convergence rate, is able to achieve optimal results, and outperforms other state-of-art algorithms.
• ### Slipping Magnetic Reconnection of Flux Rope Structures as a Precursor to an Eruptive X-class Solar Flare(1608.02057)

Aug. 6, 2016 astro-ph.SR
We present the quasi-periodic slipping motion of flux rope structures prior to the onset of an eruptive X-class flare on 2015 March 11, obtained by the \emph{Interface Region Imaging Spectrograph} (\emph{IRIS}) and the \emph{Solar Dynamics Observatory} (\emph{SDO}). The slipping motion occurred at the north part of the flux rope and seemed to successively peel off the flux rope. The speed of the slippage was 30$-$40 km s$^{-1}$, with an average period of 130$\pm$30 s. The Si {\sc iv} 1402.77 {\AA} line showed a redshift of 10$-$30 km s$^{-1}$ and a line width of 50$-$120 km s$^{-1}$ at the west legs of slipping structures, indicative of reconnection downflow. The slipping motion lasted about 40 min and the flux rope started to rise up slowly at the late stage of the slippage. Then an X2.1 flare was initiated and the flux rope was impulsively accelerated. One of the flare ribbons swept across a negative-polarity sunspot and the penumbral segments of the sunspot decayed rapidly after the flare. We studied the magnetic topology at the flaring region and the results showed the existence of a twisted flux rope, together with quasi-separatrix layers (QSLs) structures binding the flux rope. Our observations imply that quasi-periodic slipping magnetic reconnection occurs along the flux-rope-related QSLs in the preflare stage, which drives the later eruption of the flux rope and the associated flare.
• ### Reading and Writing Single-Atom Magnets(1607.03977)

July 14, 2016 cond-mat.mes-hall
The highest-density magnetic storage media will code data in single-atom bits. To date, the smallest individually addressable bistable magnetic bits on surfaces consist of 5-12 atoms. Long magnetic relaxation times were demonstrated in molecular magnets containing one lanthanide atom, and recently in ensembles of single holmium (Ho) atoms supported on magnesium oxide (MgO). Those experiments indicated the possibility for data storage at the fundamental limit, but it remained unclear how to access the individual magnetic centers. Here we demonstrate the reading and writing of individual Ho atoms on MgO, and show that they independently retain their magnetic information over many hours. We read the Ho states by tunnel magnetoresistance and write with current pulses using a scanning tunneling microscope. The magnetic origin of the long-lived states is confirmed by single-atom electron paramagnetic resonance (EPR) on a nearby Fe sensor atom, which shows that Ho has a large out-of-plane moment of $(10.1 \pm 0.1)$ $\mu_{\rm B}$ on this surface. In order to demonstrate independent reading and writing, we built an atomic scale structure with two Ho bits to which we write the four possible states and which we read out remotely by EPR. The high magnetic stability combined with electrical reading and writing shows that single-atom magnetic memory is possible.
• ### Quantifying the Topology and Evolution of a Magnetic Flux Rope Associated with Multi-flare Activities(1604.07502)

April 26, 2016 astro-ph.SR
Magnetic flux rope (MFR) plays an important role in solar activities. A quantitative assessment of the topology of an MFR and its evolution is crucial for a better understanding of the relationship between the MFR and the associated activities. In this paper, we investigate the magnetic field of active region 12017 from 2014 March 28 to 29, where 12 flares were triggered by the intermittent eruptions of a filament (either successful or confined). Using the vector magnetic field data from the Helioseismic and Magnetic Imager on board the \textit{Solar Dynamics Observatory}, we calculate the magnetic energy and helicity injection in the active region, and extrapolate the 3D magnetic field with a nonlinear force-free field model. From the extrapolations, we find an MFR that is cospatial with the filament. We further determine the configuration of this MFR by a closed quasi-separatrix layer (QSL) around it. Then, we calculate the twist number and the magnetic helicity for the field lines composing the MFR. The results show that the closed QSL structure surrounding the MFR gets smaller as a consequence of the flare occurrence. We also find that the flares in our sample are mainly triggered by kink instability. Moreover, the twist number varies more sensitively than other parameters to the occurrence of flares.
• ### Bidirectional outflows as evidence of magnetic reconnection leading to a solar microflare(1603.00941)

March 3, 2016 astro-ph.SR
Magnetic reconnection is a rapid energy release process that is believed to be responsible for flares on the Sun and stars. Nevertheless, such flare-related reconnection is mostly detected to occur in the corona, while there have been few studies concerning the reconnection in the chromosphere or photosphere. Here we present both spectroscopic and imaging observations of magnetic reconnection in the chromosphere leading to a microflare. During the flare peak time, chromospheric line profiles show significant blueshifted/redshifted components on the two sides of the flaring site, corresponding to upflows and downflows with velocities of $\pm$(70--80) km s$^{-1}$, comparable with the local Alfv\'{e}n speed as expected by the reconnection in the chromosphere. The three-dimensional nonlinear force-free field configuration further discloses twisted field lines (a flux rope) at a low altitude, cospatial with the dark threads in He I 10830 \r{A} images. The instability of the flux rope may initiate the flare-related reconnection. These observations provide clear evidence of magnetic reconnection in the chromosphere and show the similar mechanisms of a microflare to those of major flares.
• ### Modeling and Simulation for Fluid-Rotating Structure Interaction(1510.05152)

Nov. 5, 2015 math.NA
In this paper, we study a dynamic fluid-structure interaction (FSI) model for an elastic structure that is immersed and spinning in the fluid. We develop a linear constitutive model to describe the motion of a rotational elastic structure which is suitable for the application of arbitrary Lagrangian-Eulerian (ALE) method in FSI simulation. Additionally, a novel ALE mapping method is designed to generate the moving fluid mesh while the deformable structure spins in a non-axisymmetric fluid channel. The structure velocity is adopted as the principle unknown to form a monolithic saddle-point system together with fluid velocity and pressure. We discretize the nonlinear saddle-point system with mixed finite element method and Newton's linearization, and prove that the derived saddle-point problem is well-posed. The developed methodology is applied to a self-defined elastic structure and a realistic hydro-turbine under a prescribed angular velocity. Both illustrate the satisfactory numerical results of an elastic structure that is deforming and rotating while interacting with the fluid. The numerical validation is also conducted to demonstrate the modeling consistency.
• ### On the 2012 October 23 circular ribbon flare: emission features and magnetic topology(1505.02914)

May 12, 2015 astro-ph.SR
Circular ribbon flares are usually related to spine-fan type magnetic topology containing null-points. In this paper, we investigate an X-class circular ribbon flare on 2012 October 23, using the multi-wavelength data from the \textit{Solar Dynamics Observatory}, \textit{Hinode}, and the \textit{Ramaty High Energy Solar Spectroscopic Imager}. In \ion{Ca}{2} H emission, the flare showed three ribbons with two highly elongated ones inside and outside a quasi-circular one, respectively. A hot channel was displayed in the extreme ultraviolet (EUV) emissions that infers the existence of a magnetic flux rope. Two hard X-ray (HXR) sources in the 12--25 keV energy band were located at the footpoints of this hot channel. Using a nonlinear force-free magnetic field extrapolation, we identify three topological structures: (1) a 3D null-point, (2) a flux rope below the fan of the null-point, and (3) a large-scale quasi-separatrix layers (QSL) induced by the quadrupolar-like magnetic field of the active region. We find that the null-point is embedded within the large-scale QSL. In our case, all three identified topological structures must be considered to explain all the emission features associated with the observed flare. Besides, the HXR sources are regarded as the consequence of the reconnection within or near the border of the flux rope.
• ### A fast clustering algorithm for mining social network data(1403.1214)

Aug. 29, 2014 physics.soc-ph, cs.SI
Many groups with diverse convictions are interacting online. Interactions in online communities help people to engage each other and enhance understanding across groups. Online communities include multiple sub-communities whose members are similar due to social ties, characteristics, or ideas on a topic. In this research, we are interested in understanding the changes in the relative size and activity of these sub-communities, their merging or splitting patterns, and the changes in the perspectives of the members of these sub-communities due to endogenous dynamics inside the community.
• ### On modeling nonhomogeneous Poisson process for stochastic simulation input analysis(1402.7112)

Aug. 29, 2014 stat.AP
A validated simulation model primarily requires performing an appropriate input analysis mainly by determining the behavior of real-world processes using probability distributions. In many practical cases, probability distributions of the random inputs vary over time in such a way that the functional forms of the distributions and/or their parameters depend on time. This paper answers the question whether a sequence of observations from a process follows the same statistical distribution, and if not, where the exact change points are. We propose a Likelihood Ratio Test (LRT) based method to detect multiple change points when observations follow non-stationary Poisson process with diverse occurrence rates over time. Results from a comprehensive Monte Carlo study indicate satisfactory performance for the proposed method.
• ### Well-posedness and Robust Preconditioners for the Discretized Fluid-Structure Interaction Systems(1403.0046)

April 24, 2014 math.NA
In this paper we develop a family of preconditioners for the linear algebraic systems arising from the arbitrary Lagrangian-Eulerian discretization of some fluid-structure interaction models. After the time discretization, we formulate the fluid-structure interaction equations as saddle point problems and prove the uniform well-posedness. Then we discretize the space dimension by finite element methods and prove their uniform well-posedness by two different approaches under appropriate assumptions. The uniform well-posedness makes it possible to design robust preconditioners for the discretized fluid-structure interaction systems. Numerical examples are presented to show the robustness and efficiency of these preconditioners.
• ### A predictive analytics approach to reducing avoidable hospital readmission(1402.5991)

March 12, 2014 cs.AI, stat.AP