In neuroimaging, hundreds to hundreds of thousands of tests are performed
across a set of brain regions or all locations in an image. Recent studies have
shown that the most common family-wise error (FWE) controlling procedures in
imaging, which rely on classical mathematical inequalities or Gaussian random
field theory, yield FWE rates that are far from the nominal level. Depending on
the approach used, the FWER can be exceedingly small or grossly inflated. Given
the widespread use of neuroimaging as a tool for understanding neurological and
psychiatric disorders, it is imperative that reliable multiple testing
procedures are available. To our knowledge, only permutation joint testing
procedures have been shown to reliably control the FWER at the nominal level.
However, these procedures are computationally intensive due to the increasingly
available large sample sizes and dimensionality of the images, and analyses can
take days to complete. Here, we develop a parametric bootstrap joint testing
procedure. The parametric bootstrap procedure works directly with the test
statistics, which leads to much faster estimation of adjusted \emph{p}-values
than resampling-based procedures while reliably controlling the FWER in sample
sizes available in many neuroimaging studies. We demonstrate that the procedure
controls the FWER in finite samples using simulations, and present region- and
voxel-wise analyses to test for sex differences in developmental trajectories
of cerebral blood flow.
We present the final results from a targeted search for brown dwarfs with
unusual near-infrared colors. From a positional cross-match of SDSS, 2MASS and
WISE, we have identified 144 candidate peculiar L and T dwarfs. Spectroscopy
confirms that 20 of the objects are peculiar or are candidate binaries. Nine of
the 420 objects in our sample are young ($\lesssim$200 Myr; 2.1%) and another 8
(1.9%) are unusually red with no signatures of youth. With a spectroscopic
$J-K_s$ color of 2.58 $\pm$ 0.11 mag, one of the new objects, the L6 dwarf
2MASS J03530419+0418193, is among the reddest field dwarfs currently known and
is one of the reddest objects with no signatures of youth known to date. We
have also discovered another potentially very low gravity object, the L1 dwarf
2MASS J00133470+1109403, and independently identified the young L7 dwarf 2MASS
J00440332+0228112, first reported by Schneider and collaborators. Our results
confirm that signatures of low gravity are no longer discernible in low to
moderate resolution spectra of objects older than $\sim$200 Myr. The 1.9% of
unusually red L dwarfs that do not show other signatures of youth could be
slightly older, up to $\sim$400 Myr. In this case a red $J-K_s$ color may be
more diagnostic of moderate youth than individual spectral features. However,
its is also possible that these objects are relatively metal-rich, and so have
an enhanced atmospheric dust content.
Authors' comments: 27 pages, 7 figures, 5 tables, accepted to AJ
Video deblurring is a challenging problem as the blur is complex and usually
caused by the combination of camera shakes, object motions, and depth
variations. Optical flow can be used for kernel estimation since it predicts
motion trajectories. However, the estimates are often inaccurate in complex
scenes at object boundaries, which are crucial in kernel estimation. In this
paper, we exploit semantic segmentation in each blurry frame to understand the
scene contents and use different motion models for image regions to guide
optical flow estimation. While existing pixel-wise blur models assume that the
blur kernel is the same as optical flow during the exposure time, this
assumption does not hold when the motion blur trajectory at a pixel is
different from the estimated linear optical flow. We analyze the relationship
between motion blur trajectory and optical flow, and present a novel pixel-wise
non-linear kernel model to account for motion blur. The proposed blur model is
based on the non-linear optical flow, which describes complex motion blur more
effectively. Extensive experiments on challenging blurry videos demonstrate the
proposed algorithm performs favorably against the state-of-the-art methods.
Authors' comments: ICCV 2017
In this article, we provide a systematic way of creating generalized Moran
sets using an analogous iterated function system (IFS) procedure. We use a
step-wise adjustable IFS to introduce some variance (such as
non-self-similarity) in the fractal limit sets. The process retains the
computational simplicity of a standard IFS procedure. In our construction of
the generalized Moran sets, we relax the second Moran Structure Condition so
that the limit set is not necessarily self-similar. We also weaken the fourth
Moran Structure Condition so that there are no limitations on the ratios of the
diameters of the sets across a generation. Moreover, we provide upper and lower
bounds for the Hausdorff dimension of the fractals created from this
generalized process. Specific examples (Cantor-like sets, Sierpinski-like
Triangles, etc) with the calculations of their corresponding dimensions are
studied.
Authors' comments: 31 pages, 20 figures
The paper deals with point-wise estimates for the heat kernel of a nonlocal
convolution type operator with a kernel that decays at least exponentially at
infinity. It is shown that the large time behaviour of the heat kernel depends
essentially on whether $|x|\ll t$, or $|x|\sim t$, or $|x|\gg t$. We obtain
sharp point-wise upper bounds for the heat kernel in all these regions. In the
first region the nonlocal heat kernel behaves like the classical one, while in
the other regions we observe an essential difference.
We present an analysis on the potential prevalence of mid infrared excesses
in stars with and without planetary companions. Based on an extended database
of stars detected with the ${\it WISE}$ satellite, we studied two stellar
samples: one with 236 planet hosts and another with 986 objects for which
planets have been searched but not found. We determined the presence of an
excess over the photosphere by comparing the observed flux ratio at 22 $\mu$m
and 12 $\mu$m ($f_{22}/f_{12}$) with the corresponding synthetic value, derived
from results of classical model photospheres. We found a detection rate of
0.85$\%$ at 22 $\mu$m (2 excesses) in the sample of stars with planets and
0.1$\%$ (1 detection) for the stars without planets. The difference of the
detection rate between the two samples is not statistically significant, a
result that is independent of the different approaches found in the literature
to define an excess in the wavelength range covered by ${\it WISE}$
observations. As an additional result, we found that the ${\it WISE}$ fluxes
required a normalisation procedure to make them compatible with synthetic data,
probably pointing out a revision of the ${\it WISE}$ data calibration.
Authors' comments: 10 pages, 6 figures, 3 tables. Accepted for publication in MNRAS
Wide-angle photometric surveys of previously uncharted sky areas or
wavelength regimes will always bring in unexpected sources whose existence and
properties cannot be easily predicted from earlier observations: novelties or
even anomalies. Such objects can be efficiently sought for with novelty
detection algorithms. Here we present an application of such a method, called
one-class support vector machines (OCSVM), to search for anomalous patterns
among sources preselected from the mid-infrared AllWISE catalogue covering the
whole sky. To create a model of expected data we train the algorithm on a set
of objects with spectroscopic identifications from the SDSS DR13 database,
present also in AllWISE. OCSVM detects as anomalous those sources whose
patterns - WISE photometric measurements in this case - are inconsistent with
the model. Among the detected anomalies we find artefacts, such as objects with
spurious photometry due to blending, but most importantly also real sources of
genuine astrophysical interest. Among the latter, OCSVM has identified a sample
of heavily reddened AGN/quasar candidates distributed uniformly over the sky
and in a large part absent from other WISE-based AGN catalogues. It also
allowed us to find a specific group of sources of mixed types, mostly stars and
compact galaxies. By combining the semi-supervised OCSVM algorithm with
standard classification methods it will be possible to improve the latter by
accounting for sources which are not present in the training sample but are
otherwise well-represented in the target set. Anomaly detection adds
flexibility to automated source separation procedures and helps verify the
reliability and representativeness of the training samples. It should be thus
considered as an essential step in supervised classification schemes to ensure
completeness and purity of produced catalogues.
Authors' comments: 14 pages, 15 figures
We present the first results of an ongoing variability monitoring program of
active galactic nuclei (AGNs) using the 46\,cm telescope of the Wise
observatory in Israel. The telescope has a field of view of $1.25^{\circ}
\times 0.84^{\circ}$ and is specially equipped with five narrow band filters at
4300, 5200, 5700, 6200 and 7000\,\AA\, in order to perform photometric
reverberation mapping studies of the central engine of AGNs. The program aims
to observe a sample of 27 AGNs (V $<17$ mag) selected according to tentative
continuum and line time delay measurements obtained in previous works. We
describe the autonomous operation of the telescope together with the fully
automatic pipeline used to achieve high-performance unassisted observations,
data reduction, and light curves extraction using different photometric
methods. The science verification data presented here, demonstrates the
performance of the monitoring program in particular for efficiently photometric
reverberation mapping of AGNs with additional capabilities to carry out
complementary studies of other transient and variable phenomena such as
variable stars studies.
Authors' comments: 14 pages, 13 figures, accepted for publication in PASP
The interplay of rotation and manifested magnetic activity on ultracool
dwarfs (UCDs) is of key importance for gathering clues as to the operation of
the dynamos within these objects. A number of magnetized UCDs host kG-strength
magnetic fields. It was recently reported that the L8 dwarf WISEP
J060738.65+242953.4 is a radio-emitting UCD that is likely observed pole-on,
due to its lack of photometric variability and narrow spectral lines. Follow-up
radio observations at Arecibo Observatory, together with careful analysis of
previously published details, however, suggest that the scientific and
statistical significance of the radio and spectroscopic data have been
overstated. If the UCD is observed along its aligned spin/magnetic axis, the
absence of observed H$\alpha$ activity may present challenges to the auroral
model of UCD magnetism, although short or long-term cyclic magnetic activity
may explain this behavior. Monte Carlo simulations presented here suggest that
the source probably rotates with v sin i=6-12 km s$^{-1}$, indicating that its
inclination angle and rotational velocity are unexceptional and that its
angular momentum has evolved as expected for brown dwarfs observed in ~1
Myr-old clusters. The discovery and verification of the most rapidly and
slowest rotating brown dwarfs place valuable constraints on the angular
momentum evolution and magnetic activity history of these objects.
Authors' comments: Published in ApJ; 7 pages, 1 figure
The two-component broad-band spectral energy distributions of blazars were
suggested to form a sequence in which (1) the peak frequency of the low-energy
(synchrotron) component nu_syn is anticorrelated with the synchrotron
luminosity L_syn, and (2) the luminosity ratio of the high-energy (inverse
Compton) to synchrotron components q = L_IC / L_syn (Compton dominance)
increases with L_syn from the BL Lac objects (BL Lacs) to the flat spectrum
radio quasars (FSRQs). The Compton dominance parameter is an important probe of
plasma magnetisation in the blazar zones within relativistic jets. We
investigate a sample of blazars detected by WISE in the mid-infrared (MIR) band
and by Fermi/LAT in the GeV gamma-ray band, with the focus on the distribution
of luminosities and photon indices. Our findings are the following: (1) the MIR
photon index Gamma_W12 is a useful probe of the blazar sequence, with the
exception of low-luminosity BL Lacs that are most likely contaminated by their
host galaxies (L_W1 ~ 10^44 erg/s and Gamma_W12 < 1); (2) Gamma_W12 is
correlated with the gamma-ray photon index Gamma_1-100GeV, with the MIR
luminosity L_W1, and with the Fermi/WISE Compton dominance q_FW = L_1GeV/L_W1;
(3) a clean separation between FSRQs and BL Lacs can be seen in the parameter
space of Gamma_W12 and q_FW; (4) the observed distribution of MIR luminosity
L_W1 vs. Compton dominance q_FW for the entire sample of blazars can be modeled
as a sequence of lepto-magnetic jet powers in the range (log_10 P_eB) in
[42:45] with the preference for sub-equipartition magnetic fields P_B/P_e in
[0.05:1], assuming fixed bulk Lorentz factor Gamma_j = 15, fixed jet opening
angle Gamma_j Theta_j = 0.3, fixed radiative efficiency of jet electrons
epsilon_em = 50%, and that external radiation luminosity scales like L_ext
propto (P_eB)^(1.6) (parameter degeneracies are discussed).
Authors' comments: 9 pages, 5 figures, submitted to A&A
Context. Mars-crossing asteroids (MCs) are a dynamically unstable group
between the main belt and the near-Earth populations. Characterising the
physical properties of a large sample of MCs can help to understand the
original sources of many near-Earth asteroids, some of which may produce
meteorites on Earth. Aims. Our aim is to provide diameters and albedos of MCs
with available WISE/NEOWISE data. Methods. We used the near-Earth asteroid
thermal model to find the best-fitting values of equivalent diameter and,
whenever possible, the infrared beaming parameter. With the diameter and
tabulated asteroid absolute magnitudes we also computed the visible geometric
albedos. Results. We determined the diameters and beaming parameters of 404
objects observed during the fully cryogenic phase of the WISE mission, most of
which have not been published elsewhere. We also obtained 1572 diameters from
data from the 3-Band and posterior non-cryogenic phases using a default value
of beaming parameter. The average beaming parameter is 1.2 +/- 0.2 for objects
smaller than 10 km, which constitute most of our sample. This is higher than
the typical value of 1.0 found for the whole main belt and is possibly related
to the fact that WISE is able to observe many more small objects at shorter
heliocentric distances, i.e. at higher phase angles. We argue that this is a
better default value for modelling Mars-crossing asteroids from the
WISE/NEOWISE catalogue and discuss the effects of this choice on the diameter
and albedo distributions. We find a double-peaked distribution for the visible
geometric albedos, which is expected since this population is compositionally
diverse and includes objects in the major spectral complexes. However, the
distribution of beaming parameters is homogeneous for both low- and high-albedo
objects.
Authors' comments: 8 pages, 6 figures, accepted for publication in Astronomy &
Astrophysics
Many real datasets contain values missing not at random (MNAR). In this
scenario, investigators often perform list-wise deletion, or delete samples
with any missing values, before applying causal discovery algorithms. List-wise
deletion is a sound and general strategy when paired with algorithms such as
FCI and RFCI, but the deletion procedure also eliminates otherwise good samples
that contain only a few missing values. In this report, we show that we can
more efficiently utilize the observed values with test-wise deletion while
still maintaining algorithmic soundness. Here, test-wise deletion refers to the
process of list-wise deleting samples only among the variables required for
each conditional independence (CI) test used in constraint-based searches.
Test-wise deletion therefore often saves more samples than list-wise deletion
for each CI test, especially when we have a sparse underlying graph. Our
theoretical results show that test-wise deletion is sound under the justifiable
assumption that none of the missingness mechanisms causally affect each other
in the underlying causal graph. We also find that FCI and RFCI with test-wise
deletion outperform their list-wise deletion and imputation counterparts on
average when MNAR holds in both synthetic and real data.
How to develop slim and accurate deep neural networks has become crucial for
real- world applications, especially for those employed in embedded systems.
Though previous work along this research line has shown some promising results,
most existing methods either fail to significantly compress a well-trained deep
network or require a heavy retraining process for the pruned deep network to
re-boost its prediction performance. In this paper, we propose a new layer-wise
pruning method for deep neural networks. In our proposed method, parameters of
each individual layer are pruned independently based on second order
derivatives of a layer-wise error function with respect to the corresponding
parameters. We prove that the final prediction performance drop after pruning
is bounded by a linear combination of the reconstructed errors caused at each
layer. Therefore, there is a guarantee that one only needs to perform a light
retraining process on the pruned network to resume its original prediction
performance. We conduct extensive experiments on benchmark datasets to
demonstrate the effectiveness of our pruning method compared with several
state-of-the-art baseline methods.
We have reprocessed over 100 terabytes of single-exposure WISE/NEOWISE images
to create the deepest ever full-sky maps at 3-5 microns. We incorporate all
publicly available W1 and W2 imaging - a total of ~8 million exposures in each
band - from ~37 months of observations spanning 2010 January to 2015 December.
Our coadds preserve the native WISE resolution and feature depth of coverage ~3
times greater than that of the AllWISE Atlas stacks. Our coadds are designed to
enable deep forced photometry, in particular for the Dark Energy Camera Legacy
Survey (DECaLS) and Mayall z-Band Legacy Survey (MzLS), both of which are being
used to select targets for the Dark Energy Spectroscopic Instrument (DESI). We
describe newly introduced processing steps aimed at leveraging added redundancy
to remove artifacts, with the intent of facilitating uniform target selection
and searches for rare/exotic objects (e.g. high-redshift quasars and distant
galaxy clusters). Forced photometry depths achieved with these coadds extend
0.56 (0.46) magnitudes deeper in W1 (W2) than is possible with only
pre-hibernation WISE imaging.
Authors' comments: data release available at http://unwise.me
Recent work has found evidence for a difference in the bias and dark matter
halo masses of WISE-selected obscured and unobscured quasars, implying a
distinction between these populations beyond random line-of-sight effects.
However, the significance of this difference in the most up-to-date
measurements is relatively weak, at $\sim$2$\sigma$ for individual measurements
but bolstered by agreement from different techniques, including angular
clustering and cross-correlations with cosmic microwave background (CMB)
lensing maps. Here, we expand the footprint of previous work, aiming to improve
the precision of both methods. In this larger area we correct for position
dependent selection effects, in particular fluctuations of the WISE-selected
quasar density as a function of Galactic latitude. We also measure the
cross-correlation of the obscured and unobscured samples and confirm that they
are well-matched in redshift, both centred at $z=1$. Combined with very similar
detection fractions and magnitude distributions in the long-wavelength WISE
bands, this redshift match strongly supports the fact that IR selection
identifies obscured and unobscured quasars of similar bolometric luminosity.
Finally, we perform cross-correlations with confirmed spectroscopic quasars,
again confirming the results from other methods --- obscured quasars reside in
haloes a factor of 3 times more massive than unobscured quasars. This
difference is significant at the $\sim$5$\sigma$ level when the measurements
are combined, strong support for the idea that obscuration in at least some
quasars is tied to the larger environment, and may have an evolutionary
component.
Authors' comments: 15 pages, 12 figures, accepted for publication in MNRAS
We investigate extremely luminous dusty galaxies in the environments around
WISE-selected hot dust obscured galaxies (Hot DOGs) and WISE/radio-selected
active galactic nuclei (AGNs) at average redshifts of z = 2.7 and z = 1.7,
respectively. Previous observations have detected overdensities of companion
submillimetre-selected sources around 10 Hot DOGs and 30 WISE/radio AGNs, with
overdensities of ~ 2 - 3 and ~ 5 - 6 , respectively. We find that the space
densities in both samples to be overdense compared to normal star-forming
galaxies and submillimetre galaxies (SMGs) in the SCUBA-2 Cosmology Legacy
Survey (S2CLS). Both samples of companion sources have consistent mid-IR
colours and mid-IR to submm ratios as SMGs. The brighter population around
WISE/radio AGNs could be responsible for the higher overdensity reported. We
also find the star formation rate density (SFRDs) are higher than the field,
but consistent with clusters of dusty galaxies. WISE-selected AGNs appear to be
good signposts for protoclusters at high redshift on arcmin scales. The results
reported here provide an upper limit to the strength of angular clustering
using the two-point correlation function. Monte Carlo simulations show no
angular correlation, which could indicate protoclusters on scales larger than
the SCUBA-2 1.5arcmin scale maps.
Authors' comments: 10 pages, 7 figures, 6 tables, MNRAS accepted
CTR prediction in real-world business is a difficult machine learning problem
with large scale nonlinear sparse data. In this paper, we introduce an
industrial strength solution with model named Large Scale Piece-wise Linear
Model (LS-PLM). We formulate the learning problem with $L_1$ and $L_{2,1}$
regularizers, leading to a non-convex and non-smooth optimization problem.
Then, we propose a novel algorithm to solve it efficiently, based on
directional derivatives and quasi-Newton method. In addition, we design a
distributed system which can run on hundreds of machines parallel and provides
us with the industrial scalability. LS-PLM model can capture nonlinear patterns
from massive sparse data, saving us from heavy feature engineering jobs. Since
2012, LS-PLM has become the main CTR prediction model in Alibaba's online
display advertising system, serving hundreds of millions users every day.
The eShel, an off-the-shelf, fiber-fed echelle spectrograph ($R \approx
10,000$), was installed on the 1m telescope at the Wise observatory in Israel.
We report the installation of the multi-order spectrograph, and describe our
pipeline to extract stellar radial velocity from the obtained spectra. We also
introduce a new algorithm---UNICOR, to remove radial-velocity systematics that
can appear in some of the observed orders. We show that the system performance
is close to the photon-noise limit for exposures with more than $10^7$ counts,
with a precision that can get better than 200 m/s for F--K stars, for which the
eShel spectral response is optimal. This makes the eShel at Wise a useful tool
for studying spectroscopic binaries brighter than $m_V=11$. We demonstrate this
capability with orbital solutions of two binaries from projects being performed
at Wise.
Authors' comments: 30 pages, 12 figures. Accepted for publication in PASP
Hot, Dust-Obscured Galaxies, or "Hot DOGs", are a rare, dusty, hyperluminous
galaxy population discovered by the WISE mission. Predominantly at redshifts
2-3, they include the most luminous known galaxies in the universe. Their high
luminosities likely come from accretion onto highly obscured super massive
black holes (SMBHs). We have conducted a pilot survey to measure the SMBH
masses of five z~2 Hot DOGs via broad H_alpha emission lines, using
Keck/MOSFIRE and Gemini/FLAMINGOS-2. We detect broad H_alpha emission in all
five Hot DOGs. We find substantial corresponding SMBH masses for these Hot DOGs
(~ 10^{9} M_sun), and their derived Eddington ratios are close to unity. These
z~2 Hot DOGs are the most luminous AGNs at given BH masses, suggesting they are
accreting at the maximum rates for their BHs. A similar property is found for
known z~6 quasars. Our results are consistent with scenarios in which Hot DOGs
represent a transitional, high-accretion phase between obscured and unobscured
quasars. Hot DOGs may mark a special evolutionary stage before the red quasar
and optical quasar phases, and they may be present at other cosmic epochs.
Authors' comments: 15 pages, 9 figures. Accepted by ApJ
The stochastic block model is widely used for detecting community structures
in network data. How to test the goodness-of-fit of the model is one of the
fundamental problems and has gained growing interests in recent years. In this
article, we propose a novel goodness-of-fit test based on the maximum entry of
the centered and re-scaled adjacency matrix for the stochastic block model. One
noticeable advantage of the proposed test is that the number of communities can
be allowed to grow linearly with the number of nodes ignoring a logarithmic
factor. We prove that the null distribution of the test statistic converges in
distribution to a Gumbel distribution, and we show that both the number of
communities and the membership vector can be tested via the proposed method.
Further, we show that the proposed test has asymptotic power guarantee against
a class of alternatives. We also demonstrate that the proposed method can be
extended to the degree-corrected stochastic block model. Both simulation
studies and real-world data examples indicate that the proposed method works
well.
Authors' comments: 31 pages
