8,927 research outputs found
Fast model-fitting of Bayesian variable selection regression using the iterative complex factorization algorithm
Bayesian variable selection regression (BVSR) is able to jointly analyze
genome-wide genetic datasets, but the slow computation via Markov chain Monte
Carlo (MCMC) hampered its wide-spread usage. Here we present a novel iterative
method to solve a special class of linear systems, which can increase the speed
of the BVSR model-fitting tenfold. The iterative method hinges on the complex
factorization of the sum of two matrices and the solution path resides in the
complex domain (instead of the real domain). Compared to the Gauss-Seidel
method, the complex factorization converges almost instantaneously and its
error is several magnitude smaller than that of the Gauss-Seidel method. More
importantly, the error is always within the pre-specified precision while the
Gauss-Seidel method is not. For large problems with thousands of covariates,
the complex factorization is 10 -- 100 times faster than either the
Gauss-Seidel method or the direct method via the Cholesky decomposition. In
BVSR, one needs to repetitively solve large penalized regression systems whose
design matrices only change slightly between adjacent MCMC steps. This slight
change in design matrix enables the adaptation of the iterative complex
factorization method. The computational innovation will facilitate the
wide-spread use of BVSR in reanalyzing genome-wide association datasets.Comment: Accepted versio
A similarity-based community detection method with multiple prototype representation
Communities are of great importance for understanding graph structures in
social networks. Some existing community detection algorithms use a single
prototype to represent each group. In real applications, this may not
adequately model the different types of communities and hence limits the
clustering performance on social networks. To address this problem, a
Similarity-based Multi-Prototype (SMP) community detection approach is proposed
in this paper. In SMP, vertices in each community carry various weights to
describe their degree of representativeness. This mechanism enables each
community to be represented by more than one node. The centrality of nodes is
used to calculate prototype weights, while similarity is utilized to guide us
to partitioning the graph. Experimental results on computer generated and
real-world networks clearly show that SMP performs well for detecting
communities. Moreover, the method could provide richer information for the
inner structure of the detected communities with the help of prototype weights
compared with the existing community detection models
The belief noisy-or model applied to network reliability analysis
One difficulty faced in knowledge engineering for Bayesian Network (BN) is
the quan-tification step where the Conditional Probability Tables (CPTs) are
determined. The number of parameters included in CPTs increases exponentially
with the number of parent variables. The most common solution is the
application of the so-called canonical gates. The Noisy-OR (NOR) gate, which
takes advantage of the independence of causal interactions, provides a
logarithmic reduction of the number of parameters required to specify a CPT. In
this paper, an extension of NOR model based on the theory of belief functions,
named Belief Noisy-OR (BNOR), is proposed. BNOR is capable of dealing with both
aleatory and epistemic uncertainty of the network. Compared with NOR, more rich
information which is of great value for making decisions can be got when the
available knowledge is uncertain. Specially, when there is no epistemic
uncertainty, BNOR degrades into NOR. Additionally, different structures of BNOR
are presented in this paper in order to meet various needs of engineers. The
application of BNOR model on the reliability evaluation problem of networked
systems demonstrates its effectiveness
Optimal stopping under probability distortion
We formulate an optimal stopping problem for a geometric Brownian motion
where the probability scale is distorted by a general nonlinear function. The
problem is inherently time inconsistent due to the Choquet integration
involved. We develop a new approach, based on a reformulation of the problem
where one optimally chooses the probability distribution or quantile function
of the stopped state. An optimal stopping time can then be recovered from the
obtained distribution/quantile function, either in a straightforward way for
several important cases or in general via the Skorokhod embedding. This
approach enables us to solve the problem in a fairly general manner with
different shapes of the payoff and probability distortion functions. We also
discuss economical interpretations of the results. In particular, we justify
several liquidation strategies widely adopted in stock trading, including those
of "buy and hold", "cut loss or take profit", "cut loss and let profit run" and
"sell on a percentage of historical high".Comment: Published in at http://dx.doi.org/10.1214/11-AAP838 the Annals of
Applied Probability (http://www.imstat.org/aap/) by the Institute of
Mathematical Statistics (http://www.imstat.org
- âŠ