
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 highdimensional, sparse and complex and
sometimes there exists few amount of timetoevent (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 (timetoevent) 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
SEERMedicare data related to prostate cancer among AfricanAmericans and white
patients. The results indicate that our approach outperforms significantly than
baseline models.

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 translationinvariance. 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 lowrisk
highreturn 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.

We study two kinds of inelastic mesonmeson scattering. The first kind is
inelastic 2to2 mesonmeson scattering that is governed by quark interchange
as well as quarkantiquark annihilation and creation. Crosssection 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 2to1 mesonmeson 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 2to1 mesonmeson scattering. Starting from the
$S$matrix element, the isospinaveraged 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.

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
lowlatency communication. Specifically, we develop a lowrank 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 nonconvex optimization problem due to the multiple nonconvex
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 elementwise extension of the
geometric concepts of the individual manifolds. A scalable Riemannian
trustregion 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 stateofart methods.

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 higherlevel 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
(AfricanAmericans). Then we use two large datasets from eICU collaborative
research database to predict the length of stay in CardiacICU and NeuroICU
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.

Traditional medicine typically applies onesizefitsall 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 (AfricanAmerican). 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 (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.

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 mediumrange order
structure, while the shortrange 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.

In this paper, we propose a new deep feature selection method based on deep
architecture. Our method uses stacked autoencoders for feature representation
in higherlevel 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 (AfricanAmerican). 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.

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
lowfrequency 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.

Using the highquality 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 posteruption loops and flare
ribbons. It likely indicated a largescale 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 largescale magnetic topology of the
coronal hole is important for the eventual propagation direction of CMEs.

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.

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 stagebystage. 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
noncooperative game. Nash equilibria of the game are found in closed form.
Solutions of the two users in the three epochs are derived.

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 (FogRAN) can effectively increase
network throughput and reduce transmission latency. Furthermore, we can exploit
the benefits of cache enabled architecture in FogRAN 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 lowrank matrix completion (LRMC)
approach to solving the content delivery problem in both wireline and wireless
parts of FogRAN. To attain a low caching latency, we present a high precision
approach with Riemannian trustregion method to solve the challenging LRMC
problem by exploiting the quotient manifold geometry of fixedrank matrices.
Numerical results show that the new approach has a faster convergence rate, is
able to achieve optimal results, and outperforms other stateofart algorithms.

We present the quasiperiodic slipping motion of flux rope structures prior
to the onset of an eruptive Xclass 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
negativepolarity 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
quasiseparatrix layers (QSLs) structures binding the flux rope. Our
observations imply that quasiperiodic slipping magnetic reconnection occurs
along the fluxroperelated QSLs in the preflare stage, which drives the later
eruption of the flux rope and the associated flare.

The highestdensity magnetic storage media will code data in singleatom
bits. To date, the smallest individually addressable bistable magnetic bits on
surfaces consist of 512 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 longlived
states is confirmed by singleatom electron paramagnetic resonance (EPR) on a
nearby Fe sensor atom, which shows that Ho has a large outofplane 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 singleatom magnetic memory is possible.

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 forcefree 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 quasiseparatrix 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.

Magnetic reconnection is a rapid energy release process that is believed to
be responsible for flares on the Sun and stars. Nevertheless, such
flarerelated 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$(7080) km
s$^{1}$, comparable with the local Alfv\'{e}n speed as expected by the
reconnection in the chromosphere. The threedimensional nonlinear forcefree
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 flarerelated 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.

In this paper, we study a dynamic fluidstructure 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
LagrangianEulerian (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 nonaxisymmetric fluid channel. The structure
velocity is adopted as the principle unknown to form a monolithic saddlepoint
system together with fluid velocity and pressure. We discretize the nonlinear
saddlepoint system with mixed finite element method and Newton's
linearization, and prove that the derived saddlepoint problem is wellposed.
The developed methodology is applied to a selfdefined elastic structure and a
realistic hydroturbine 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.

Circular ribbon flares are usually related to spinefan type magnetic
topology containing nullpoints. In this paper, we investigate an Xclass
circular ribbon flare on 2012 October 23, using the multiwavelength 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 quasicircular one, respectively. A hot channel was displayed in
the extreme ultraviolet (EUV) emissions that infers the existence of a magnetic
flux rope. Two hard Xray (HXR) sources in the 1225 keV energy band were
located at the footpoints of this hot channel. Using a nonlinear forcefree
magnetic field extrapolation, we identify three topological structures: (1) a
3D nullpoint, (2) a flux rope below the fan of the nullpoint, and (3) a
largescale quasiseparatrix layers (QSL) induced by the quadrupolarlike
magnetic field of the active region. We find that the nullpoint is embedded
within the largescale 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.

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 subcommunities 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 subcommunities, their merging or splitting
patterns, and the changes in the perspectives of the members of these
subcommunities due to endogenous dynamics inside the community.

A validated simulation model primarily requires performing an appropriate
input analysis mainly by determining the behavior of realworld 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 nonstationary Poisson process with diverse
occurrence rates over time. Results from a comprehensive Monte Carlo study
indicate satisfactory performance for the proposed method.

In this paper we develop a family of preconditioners for the linear algebraic
systems arising from the arbitrary LagrangianEulerian discretization of some
fluidstructure interaction models. After the time discretization, we formulate
the fluidstructure interaction equations as saddle point problems and prove
the uniform wellposedness. Then we discretize the space dimension by finite
element methods and prove their uniform wellposedness by two different
approaches under appropriate assumptions. The uniform wellposedness makes it
possible to design robust preconditioners for the discretized fluidstructure
interaction systems. Numerical examples are presented to show the robustness
and efficiency of these preconditioners.

Hospital readmission has become a critical metric of quality and cost of
healthcare. Medicare anticipates that nearly $17 billion is paid out on the 20%
of patients who are readmitted within 30 days of discharge. Although several
interventions such as transition care management and discharge reengineering
have been practiced in recent years, the effectiveness and sustainability
depends on how well they can identify and target patients at high risk of
rehospitalization. Based on the literature, most current risk prediction models
fail to reach an acceptable accuracy level; none of them considers patient's
history of readmission and impacts of patient attribute changes over time; and
they often do not discriminate between planned and unnecessary readmissions.
Tackling such drawbacks, we develop a new readmission metric based on
administrative data that can identify potentially avoidable readmissions from
all other types of readmission. We further propose a tree based classification
method to estimate the predicted probability of readmission that can directly
incorporate patient's history of readmission and risk factors changes over
time. The proposed methods are validated with 201112 Veterans Health
Administration data from inpatients hospitalized for heart failure, acute
myocardial infarction, pneumonia, or chronic obstructive pulmonary disease in
the State of Michigan. Results shows improved discrimination power compared to
the literature (cstatistics>80%) and good calibration.

Reducing potentially preventable readmissions has been identified as an
important issue for decreasing Medicare costs and improving quality of care
provided by hospitals. Based on previous research by medical professionals,
preventable readmissions are caused by such factors as flawed patient
discharging process, inadequate followups after discharging, and noncompliance
of patients on discharging and follow up instructions. It is also found that
the risk of preventable readmission also may relate to some patient's
characteristics, such as age, health condition, diagnosis, and even treatment
specialty. In this study, using both general demographic information and
individual past history of readmission records, we develop a risk prediction
model based on hierarchical nonlinear mixed effect framework to extract
significant prognostic factors associated with patient risk of 30day
readmission. The effectiveness of our proposed approach is validated based on a
real dataset from four VA facilities in the State of Michigan. Simultaneously
explaining both patient and population based variations of readmission process,
such an accurate model can be used to recognize patients with high likelihood
of discharging noncompliances, and then targeted postcare actions can be
designed to reduce further rehospitalization.