Gene Network Reconstruction using Global-Local Shrinkage Priors

Reconstructing a gene network from high-throughput molecular data is an important but challenging task, as the number of parameters to estimate easily is much larger than the sample size. A conventional remedy is to regularize or penalize the model likelihood. In network models, this is often done $\textit{locally}$ in the neighborhood of each node or gene. However, estimation of the many regularization parameters is often difficult and can result in large statistical uncertainties. In this paper we propose to combine local regularization with $\textit{global}$ shrinkage of the regularization parameters to borrow strength between genes and improve inference. We employ a simple Bayesian model with nonsparse, conjugate priors to facilitate the use of fast variational approximations to posteriors. We discuss empirical Bayes estimation of hyperparameters of the priors, and propose a novel approach to rank-based posterior thresholding. Using extensive model- and data-based simulations, we demonstrate that the proposed inference strategy outperforms popular (sparse) methods, yields more stable edges, and is more reproducible. The proposed method, termed $\texttt{ShrinkNet}$, is then applied to Glioblastoma to investigate the interactions between genes associated with patient survival.This work was supported by the Center for Medical Systems Biology (CMSB), and the European Union Grant EpiRadBio, established by the Netherlands Genomics Initiative/Netherlands Organization for Scientific Research (NGI/NWO), nr. FP7-269553

First release of Apertif imaging survey data

Context. Apertif is a phased-array feed system for the Westerbork Synthesis Radio Telescope, providing forty instantaneous beams over 300 MHz of bandwidth. A dedicated survey program utilizing this upgrade started on 1 July 2019, with the last observations taken on 28 February 2022. The imaging survey component provides radio continuum, polarization, and spectral line data. Aims. Public release of data is critical for maximizing the legacy of a survey. Toward that end, we describe the release of data products from the first year of survey operations, through 30 June 2020. In particular, we focus on defining quality control metrics for the processed data products. Methods. The Apertif imaging pipeline, Apercal, automatically produces non-primary beam corrected continuum images, polarization images and cubes, and uncleaned spectral line and dirty beam cubes for each beam of an Apertif imaging observation. For this release, processed data products are considered on a beam-by-beam basis within an observation. We validate the continuum images by using metrics that identify deviations from Gaussian noise in the residual images. If the continuum image passes validation, we release all processed data products for a given beam. We apply further validation to the polarization and line data products and provide flags indicating the quality of those data products. Results. We release all raw observational data from the first year of survey observations, for a total of 221 observations of 160 independent target fields, covering approximately one thousand square degrees of sky. Images and cubes are released on a per beam basis, and 3374 beams (of 7640 considered) are released. The median noise in the continuum images is 41.4 uJy beam−1, with a slightly lower median noise of 36.9 uJy beam−1 in the Stokes V polarization image. The median angular resolution is 11.6″/sin δ. The median noise for all line cubes, with a spectral resolution of 36.6 kHz, is 1.6 mJy beam−1, corresponding to a 3-σ H I column density sensitivity of 1.8 × 1020 atoms cm−2 over 20 km s−1 (for a median angular resolution of 24″ × 15″). Line cubes at lower frequency have slightly higher noise values, consistent with the global RFI environment and overall Apertif system performance. We also provide primary beam images for each individual Apertif compound beam. The data are made accessible using a Virtual Observatory interface and can be queried using a variety of standard tools

A LOFAR observation of ionospheric scintillation from two simultaneous travelling ionospheric disturbances

This paper presents the results from one of the first observations of ionospheric scintillation taken using the Low-Frequency Array (LOFAR). The observation was of the strong natural radio source Cassiopeia A, taken overnight on 18–19 August 2013, and exhibited moderately strong scattering effects in dynamic spectra of intensity received across an observing bandwidth of 10–80 MHz. Delay-Doppler spectra (the 2-D FFT of the dynamic spectrum) from the first hour of observation showed two discrete parabolic arcs, one with a steep curvature and the other shallow, which can be used to provide estimates of the distance to, and velocity of, the scattering plasma. A cross-correlation analysis of data received by the dense array of stations in the LOFAR “core” reveals two different velocities in the scintillation pattern: a primary velocity of ~20–40 ms−1 with a north-west to south-east direction, associated with the steep parabolic arc and a scattering altitude in the F-region or higher, and a secondary velocity of ~110 ms−1 with a north-east to south-west direction, associated with the shallow arc and a scattering altitude in the D-region. Geomagnetic activity was low in the mid-latitudes at the time, but a weak sub-storm at high latitudes reached its peak at the start of the observation. An analysis of Global Navigation Satellite Systems (GNSS) and ionosonde data from the time reveals a larger-scale travelling ionospheric disturbance (TID), possibly the result of the high-latitude activity, travelling in the north-west to south-east direction, and, simultaneously, a smaller-scale TID travelling in a north-east to south-west direction, which could be associated with atmospheric gravity wave activity. The LOFAR observation shows scattering from both TIDs, at different altitudes and propagating in different directions. To the best of our knowledge this is the first time that such a phenomenon has been reported

Tracking of an electron beam through the solar corona with LOFAR

© ESO 2018. The Sun's activity leads to bursts of radio emission, among other phenomena. An example is type-III radio bursts. They occur frequently and appear as short-lived structures rapidly drifting from high to low frequencies in dynamic radio spectra. They are usually interpreted as signatures of beams of energetic electrons propagating along coronal magnetic field lines. Here we present novel interferometric LOFAR (LOw Frequency ARray) observations of three solar type-III radio bursts and their reverse bursts with high spectral, spatial, and temporal resolution. They are consistent with a propagation of the radio sources along the coronal magnetic field lines with nonuniform speed. Hence, the type-III radio bursts cannot be generated by a monoenergetic electron beam, but by an ensemble of energetic electrons with a spread distribution in velocity and energy. Additionally, the density profile along the propagation path is derived in the corona. It agrees well with three-fold coronal density model by (1961, ApJ, 133, 983)

LOFAR 144-MHz follow-up observations of GW170817

ABSTRACT We present low-radio-frequency follow-up observations of AT 2017gfo, the electromagnetic counterpart of GW170817, which was the first binary neutron star merger to be detected by Advanced LIGO–Virgo. These data, with a central frequency of 144 MHz, were obtained with LOFAR, the Low-Frequency Array. The maximum elevation of the target is just 13${_{.}^{\circ}}$7 when observed with LOFAR, making our observations particularly challenging to calibrate and significantly limiting the achievable sensitivity. On time-scales of 130–138 and 371–374 d after the merger event, we obtain 3σ upper limits for the afterglow component of 6.6 and 19.5 mJy beam−1, respectively. Using our best upper limit and previously published, contemporaneous higher frequency radio data, we place a limit on any potential steepening of the radio spectrum between 610 and 144 MHz: the two-point spectral index $\alpha ^{610}_{144} \gtrsim$ −2.5. We also show that LOFAR can detect the afterglows of future binary neutron star merger events occurring at more favourable elevations.</jats:p