Fixing the relationship among a set of experimental quantities is a fundamental issue in many scientific disciplines. In the two-dimensional case, the classical approach is to compute the linear correlation coefficient from a scatterplot. This method, however, implicitly assumes a linear relationship between the variables. Such assumption is not always correct. With the use of the partial correlation coefficients, an extension to the multi-dimensional case is possible. However, the problem of the assumed mutual linear relationship among the variables still remains. A relatively recent approach which permits to avoid this problem is modeling the joint probability density function (PDF) of the data with copulas. These are functions which contain all the information on the relationship between two random variables. Although in principle this approach can work also with multi-dimensional data, theoretical as well computational difficulties often limit its use to the two-dimensional case. In this paper, we consider an approach, based on so-called vine copulas, which overcomes this limitation and at the same time is amenable to a theoretical treatment and feasible from the computational point of view. We apply this method to published data on the near-IR and far-IR luminosities and atomic and molecular masses of the Herschel Reference Sample. We determine the relationship among the luminosities and gas masses and show that the far-IR luminosity can be considered as the key parameter which relates all the other three galaxy properties. Once removed from the 4D relation, the residual relation among the other three is negligible. This may be interpreted as that the correlation between the gas masses and near-IR luminosity is driven by the far-IR luminosity, likely by the star-formation activity of the galaxy.

The detection of signals hidden in noise is one of the oldest and common problems in astronomy. Various solutions have been proposed in the past such as the parametric approaches based on the least-squares fit of theoretical templates or the non-parametric techniques as the phase-folding method. Most of them, however, are suited only for signals with specific time evolution. For generic signals the spectral approach based on the periodogram is potentially the most effective. In astronomy the main problem in working with the periodogram is that often the sampling of the signals is irregular. This complicates its efficient computation (the fast Fourier transform cannot be directly used) but overall the determination of its statistical characteristics. The Lomb-Scargle periodogram (LSP) provides a solution to this last important issue, but its main drawback is the assumption of a very specific model of the datawhich is not correct for most of the practical applications. These issues are not always considered in literature with theoretical and practical consequences of no easy solution. Moreover, apart from pathological samplings, it is common believe that the LSP and the classical periodogram (CP) usually provide almost identical results. In general, this is true but here it is shown that there are situations where the LSP is less effective than the CP in the detection of signals in noise. There are no compelling reasons, therefore, to use the LSP instead of the CP which is directly connected to the correlation function of the observed signal with the sinusoidal functions at the various frequencies of interest.

The matched filter (MF) represents one of the main tools to detect signals from known sources embedded in the noise. In the Gaussian case the noise is assumed to be the realization of a Gaussian random field (GRF). The most important property of the MF, the maximization of the probability of detection subject to a constant probability of false detection or false alarm (PFA), makes it one of the most popular techniques. However, the MF technique relies upon the a priori knowledge of the number and the position of the searched signals in the GRF which usually are not available. A typical way out is to assume that the position of a signal coincides with one of the peaks in the matched filtered data. A detection is claimed when the probability that a given peak is due only to the noise (i.e. the PFA) is smaller than a prefixed threshold. In this case the probability density function (PDF) of the amplitudes has to be used for the computation of the PFA, which is different from the Gaussian. Moreover, the probability that a detection is false depends on the number of peaks present in the filtered GRF, the greater the number of peaks in a GRF, the higher the probability of peaks due to the noise that exceed the detection threshold. If not taken into account, the PFA can be severely underestimated. Many solutions proposed to this problem are non-parametric hence not able to exploit all the available information. This limitation has been overcome by means of two efficient parametric approaches, one based on the PDF of the peak amplitudes of a smooth and isotropic GRF whereas the other uses the Gumbel distribution (the asymptotic PDF of the corresponding extreme). Simulations and ALMA maps show that, although the two methods produce almost identical results, the first is more flexible and allows us to check the reliability of the detection procedure.

Why do we write this note? It is erroneous to pretend to extract physical information from the experimental light curves (time series) of astrophysical systems by means of linear stochastic differential equations (LSDE). In general, the time evolution of these systems is governed by a set of nonlinear differential equations. Hence, the LSDEs are not suitable to model their dynamics. In spite of this, recently the LSDEs have been proposed as tools for the analysis of AGN light curves. Their use in this context seems to be dictated by their simplicity rather than by a real physical argument. We stress in this note that the correct approach to the analysis of signals coming from systems with nonlinear dynamics is to tackle the problem using methodologies in well defined physical contexts.

We discuss the results of the relationships between the K-band and stellar mass, far-infrared luminosities, star formation rate, dust and gas masses of nearby galaxies computing the bivariate K-band Luminosity Function (BLF) and bivariate K-band Mass Function (BMF) of the Herschel Reference Survey (HRS), a volume-limited sample with full wavelength coverage. We derive the BLFs and BMFs from the K-band and stellar mass, far-infrared luminosities, star formation rate, dust and gas masses cumulative distributions using a copula method which is outlined in detail. The use of the bivariate computed taking into account the upper limits allows us to derive on a more solid statistical ground the relationship between the observed physical quantities. The analysis shows that the behaviour of the morphological (optically selected) subsamples is quite different. A statistically meaningful result can be obtained over the whole HRS sample only from the relationship between the K-band and the stellar mass, while for the remaining physical quantities (dust and gas masses, far-IR luminosity and star formation rate), the analysis is distinct for late-type (LT) and early-type galaxies (ETG). However, the number of ETGs is small to perform a robust statistical analysis, and in most of the case results are discussed only for the LTG subsample. The Luminosity and Mass Functions (LFs, MFs) of LTGs are generally dependent on the K-band and the various dependencies are discussed in detail. We are able to derive the corresponding LFs and MFs and compare them with those computed with other samples. Our statistical analysis allows us to characterise the HRS, that, although non homogeneously selected and partially biased towards low IR luminosities, may be considered as representative of the local LT galaxy population.

The matched filter (MF) is widely used to detect signals hidden within the noise. If the noise is Gaussian, its performances are well-known and describable in an elegant analytical form. The treatment of non-Gaussian noises is often cumbersome as in most cases there is no analytical framework. This is true also for Poisson noise which, especially in the low-number count regime, presents the additional difficulty to be discrete. For this reason in the past methods have been proposed based on heuristic or semi-heuristic arguments. Recently, an analytical form of the MF has been introduced but the computation of the probability of false detection or false alarm (PFA) is based on numerical simulations. To overcome this inefficient and time consuming approach we propose here an effective method to compute the PFA based on the saddle point approximation (SA). We provide the theoretical framework and support our findings by means of numerical simulations. We discuss also the limitations of the MF in practical applications.

The matched filter (MF) is one of the most popular and reliable techniques to the detect signals of known structure and amplitude smaller than the level of the contaminating noise. Under the assumption of stationary Gaussian noise, MF maximizes the probability of detection subject to a constant probability of false detection or false alarm (PFA). This property relies upon a priori knowledge of the position of the searched signals, which is usually not available. Recently, it has been shown that when applied in its standard form, MF may severely underestimate the PFA. As a consequence the statistical significance of features that belong to noise is overestimated and the resulting detections are actually spurious. For this reason, an alternative method of computing the PFA has been proposed that is based on the probability density function (PDF) of the peaks of an isotropic Gaussian random field. In this paper we further develop this method. In particular, we discuss the statistical meaning of the PFA and show that, although useful as a preliminary step in a detection procedure, it is not able to quantify the actual reliability of a specific detection. For this reason, a new quantity is introduced called the specific probability of false alarm (SPFA), which is able to carry out this computation. We show how this method works in targeted simulations and apply it to a few interferometric maps taken with the Atacama Large Millimeter/submillimeter Array (ALMA) and the Australia Telescope Compact Array (ATCA). We select a few potential new point sources and assign an accurate detection reliability to these sources.

The detection reliability of weak signals is a critical issue in many astronomical contexts and may have severe consequences for determining number counts and luminosity functions, but also for optimising the use of telescope time in follow-up observations. Because of its optimal properties, one of the most popular and widely-used detection technique is the matched filter (MF). This is a linear filter designed to maximise the detectability of a signal of known structure that is buried in additive Gaussian random noise. In this work we show that in the very common situation where the number and position of the searched signals within a data sequence (e.g. an emission line in a spectrum) or an image (e.g. a point-source in an interferometric map) are unknown, this technique, when applied in its standard form, may severely underestimate the probability of false detection. This is because the correct use of the MF relies upon a-priori knowledge of the position of the signal of interest. In the absence of this information, the statistical significance of features that are actually noise is overestimated and detections claimed that are actually spurious. For this reason, we present an alternative method of computing the probability of false detection that is based on the probability density function (PDF) of the peaks of a random field. It is able to provide a correct estimate of the probability of false detection for the one-, two- and three-dimensional case. We apply this technique to a real two-dimensional interferometric map obtained with ALMA.

We study the relationship between the K-band and the sub-millimetre (submm) emissions of nearby galaxies by computing the bivariate K-band-submm luminosity function (BLF) of the Herschel Reference Survey (HRS), a volume-limited sample observed in submm with Hersche/SPIRE. We derive the BLF from the K-band and submm cumulative distributions using a copula method. Using the BLF allows us to derive the relationship between the luminosities on more solid statistical ground. The analysis shows that over the whole HRS sample, no statistically meaningful conclusion can be derived for any relationship between the K-band and the submm luminosity. However, a very tight relationship between these luminosities is highlighted, by restricting our analysis to late-type galaxies. The luminosity function of late-type galaxies computed in the K-band and in the submm are dependent and the dependence is caused by the link, between the stellar mass and the cold dust mass, which has been already observed

This paper deals with the detection problem of extragalactic point-sources in multi-frequency, microwave sky maps that will be obtainable in future cosmic microwave background radiation (CMB) experiments with instruments capable of very high spatial resolution. With spatial resolutions that can be of order of 0.1-1.0 arcsec or better, the extragalactic point-sources will appear isolated. The same holds also for the compact structures due to the Sunyaev-Zeldovich (SZ) effect (both thermal and kinetic). This situation is different from the maps obtainable with instruments as WMAP or PLANCK where, because of the smaller spatial resolution (approximately 5-30 arcmin), the point-sources and the compact structures due to the SZ effect form a uniform noisy background (the "confusion noise"). Hence, the point-source detection techniques developed in the past are based on the assumption that all the emissions that contribute to the microwave background can be modeled with homogeneous and isotropic (often Gaussian) random fields and make use of the corresponding spatial power-spectra. In the case of very high resolution observations such an assumption cannot be adopted since it still holds only for the CMB. Here, we propose an approach based on the assumption that the diffuse emissions that contribute to the microwave background can be locally approximated by two-dimensional low order polynomials. In particular, two sets of numerical techniques are presented containing two different algorithms each. The performance of the algorithms is tested with numerical experiments that mimic the physical scenario expected for high Galactic latitude observations with the Atacama Large Millimeter/Submillimeter Array (ALMA).

Detection of a signal hidden by noise within a time series is an important problem in many astronomical searches, i.e. for light curves containing the contributions of periodic/semi-periodic components due to rotating objects and all other astrophysical time-dependent phenomena. One of the most popular tools for use in such studies is the "periodogram", whose use in an astronomical context is often not trivial. The "optimal" statistical properties of the periodogram are lost in the case of irregular sampling of signals, which is a common situation in astronomical experiments. Parts of these properties are recovered by the "Lomb-Scargle" (LS) technique, but at the price of theoretical difficulties, that can make its use unclear, and of algorithms that require the development of dedicated software if a fast implementation is necessary. Such problems would be irrelevant if the LS periodogram could be used to significantly improve the results obtained by approximated but simpler techniques. In this work we show that in many astronomical applications simpler techniques provide results similar to those obtainable with the LS periodogram. The meaning of the "Nyquist frequency" is also discussed in the case of irregular sampling.

We propose two detection techniques that take advantage of a small sky area approximation and are based on modifications of the "internal linear combination" (ILC) method, an approach widely used in Cosmology for the separation of the various components that contribute to the microwave background. The main advantage of the proposed approach, especially in handling multi-frequency maps of the same region, is that it does not require the "a priori" knowledge of the spatial power-spectrum of either the CMB and/or the Galactic foreground. Hence, it is more robust, easier and more intuitive to use. The performance of the proposed algorithms is tested with numerical experiments that mimic the physical scenario expected for high Galactic latitude observations with the Atacama Large Millimeter/submillimeter Array (ALMA).

In experimental microwave maps, point-sources can strongly affect the estimation of the power-spectrum and/or the test of Gaussianity of the Cosmic Microwave Background (CMB) component. As a consequence, their removal from the sky maps represents a critical step in the analysis of the CMB data. Before removing a source, however, it is necessary to detect it and source extraction consists of a delicate preliminary operation. In the literature, various techniques have been presented to detect point-sources in the sky maps. The most sophisticated ones exploit the multi-frequency nature of the observations that is typical of the CMB experiments. These techniques have "optimal" theoretical properties and, at least in principle, are capable of remarkable performances. Actually, they are rather difficult to use and this deteriorates the quality of the obtainable results. In this paper, we present a new technique, the "weighted matched filter" (WMF), that is quite simple to use and hence more robust in practical applications. Such technique shows particular efficiency in the detection of sources whose spectra have a slope different from zero. We apply this method to three Southern Hemisphere sky regions - each with an area of 400 square degrees - of the seven years Wilkinson Microwave Anisotropy Probe (WMAP) maps and compare the resulting sources with those of the two seven-year WMAP point-sources catalogues. In these selected regions we find seven additional sources not previously listed in WMAP catalogues and discuss their most likely identification and spectral properties.

The periodogram is a popular tool that tests whether a signal consists only of noise or if it also includes other components. The main issue of this method is to define a critical detection threshold that allows identification of a component other than noise, when a peak in the periodogram exceeds it. In the case of signals sampled on a regular time grid, determination of such a threshold is relatively simple. When the sampling is uneven, however, things are more complicated. The most popular solution in this case is to use the "Lomb-Scargle" periodogram, but this method can be used only when the noise is the realization of a zero-mean, white (i.e. flat-spectrum) random process. In this paper, we present a general formalism based on matrix algebra, which permits analysis of the statistical properties of a periodogram independently of the characteristics of noise (e.g. colored and/or non-stationary), as well as the characteristics of sampling.

Although the ``internal linear'' combination method (ILC) is a technique widely used for the separation of the Cosmic Microwave Background signal from the Galactic foregrounds, its characteristics are not yet well defined. This can lead to misleading conclusions about the actual potentialities and limits of such approach in real applications. Here we discuss briefly some facts about ILC that to our knowledge are not fully worked out in literature and yet have deep effects in the interpretation of the results.

The detection of point-sources in microwave maps is a critical step in the analysis of the Cosmic Microwave Background (CMB) data. If not properly removed, these sources have adverse effects on the estimation of the power-spectrum and/or the test of Gaussianity of the CMB component. The techniques to extract point sources from an observed image used up to now do not show a general consensus about their real performance and properties. Their characteristics have been studied essentially through numerical simulations based on semi-empirical models of the CMB and the Galactic foreground. Such models often have different levels of sophistication and/or are based on different physical assumptions (e.g. the number of Galactic components and level of the noise). Moreover, the application of a given technique to a set of data (either simulated or experimental) requires the tuning of one or more parameters that unavoidably is a subjective operation. Hence, a reliable comparison is difficult. What is missing is a statistical analysis of the properties of the proposed methodologies. This is the aim of the present paper. The statistical properties of the detection techniques in the context of two different criteria, i.e. the Neyman-Pearson criterion and the maximization of the signal-to-noise ratio, are analyzed through an analytical approach. Various examples, techniques and the conditions under which they work satisfactorily are shown together with their limits and implementation details. Exploiting some a priori information, it is possible to develop simple algorithms with performances similar to those of more sophisticated but complex techniques. In this respect, a detection algorithm, useful in future ground-based experiments as, for example the Atacama Large Millimeter/submillimeter Array, is presented.

Foreground contamination is the fundamental hindrance to the cosmic microwave background (CMB) signals and its separation from it represents a fundamental question in Cosmology. One of the most popular algorithm used to disentangle foregrounds from the CMB signals is the "internal linear combination" method (ILC). In its original version, this technique is applied directly to the observed maps. In recent literature, however, it is suggested that in the harmonic (Fourier) domain it is possible to obtain better results since a separation can be attempted where the various Fourier frequencies are given different weights. This is seen as a useful characteristic in the case of noisy data. Here, we argue that the benefits of using such an approach are overestimated. Better results can be obtained if a classic procedure is adopted where data are filtered before the separation is carried out.

AIMS: The separation of foreground contamination from cosmic microwave background (CMB) observations is one of the most challenging and important problem of digital signal processing in Cosmology. In literature, various techniques have been presented, but no general consensus about their real performances and properties has been reached. This is due to the characteristics of these techniques that have been studied essentially through numerical simulations based on semi-empirical models of the CMB and the Galactic foregrounds. Such models often have different level of sophistication and/or are based on different physical assumptions (e.g., the number of the Galactic components and the level of the noise). Hence, a reliable comparison is difficult. What actually is missing is a statistical analysis of the properties of the proposed methodologies. Here, we consider the "Internal Linear Combination" method (ILC) which, among the separation techniques, requires the smallest number of "a priori" assumptions. This feature is of particular interest in the context of the CMB polarization measurements at small angular scales where the lack of knowledge of the polarized backgrounds represents a serious limit. METHODS: The statistical characteristics of ILC are examined through an analytical approach and the basic conditions are fixed in a way to work satisfactorily. RESULTS: ILC provides satisfactory results only under rather restrictive conditions. This is a critical fact to take into consideration in planning the future ground-based observations (e.g., with ALMA) where, contrary to the satellite experiments, there is the possibility to have a certain control of the experimental conditions.

AIMS: One of the most challenging and important problem of digital signal processing in Cosmology is the separation of foreground contamination from cosmic microwave background (CMB). This problem becomes even more difficult in situations, as the CMB polarization observations, where the amount of available "a priori" information is limited. In this case, it is necessary to resort to the "blind separation" methods. One important member of this class is represented by the "Independent Components Analysis" (ICA). In its original formulation, this method has various interesting characteristics, but also some limits. One of the most serious is the difficulty to take into account any information available in advance. In particular, ICA is not able to exploit the fact that emission of CMB is the same at all the frequencies of observations. Here, we show how to deal with this question. The connection of the proposed methodology with the "Internal Linear Composition" (ILC) technique is also illustrated. METHODS: A modification of the classic ICA approach is presented and its characteristics are analyzed both analytically and by means of numerical experiments. RESULTS: The modified version of ICA appears to provide more stable results and of better quality.

Context. One important step in the statistical analysis of the Ly-alpha forest data is the study of their second order properties. Usually, this is accomplished by means of the two-point correlation function or, alternatively, the K-function. In the computation of these functions it is necessary to take into account the presence of strong metal line complexes and strong Ly-alpha lines that can hidden part of the Ly-alpha forest and represent a non negligible source of bias. Aims. In this work, we show quantitatively what are the effects of the gaps introduced in the spectrum by the strong lines if they are not properly accounted for in the computation of the correlation properties. We propose a geometric method which is able to solve this problem and is computationally more efficient than the Monte Carlo (MC) technique that is typically adopted in Cosmology studies. The method is implemented in two different algorithms. The first one permits to obtain exact results, whereas the second one provides approximated results but is computationally very efficient. The proposed approach can be easily extended to deal with the case of two or more lists of lines that have to be analyzed at the same time. Methods. Numerical experiments are presented that illustrate the consequences to neglect the effects due to the strong lines and the excellent performances of the proposed approach. Results. The proposed method is able to remarkably improve the estimates of both the two-point correlation function and the K-function.

The arguments presented by Gonzalez-Nuevo et. al (2006) in favour of the Mexican Hat Wavelet Family (MHWF) are critically discussed here. These authors allege the optimal properties of this new class of filters in the detection of point sources embedded in a noise background but their claim is not based upon a solid mathematical foundation and proof.

The Kluzniak & Abramowicz model explains high frequency, double peak, "3:2" QPOs observed in neutron star and black hole sources in terms of a non-linear parametric resonance between radial and vertical epicyclic oscillations of an almost Keplerian accretion disk. The 3:2 ratio of epicyclic frequencies occurs only in strong gravity. Rebusco (2004) and Horak (2004) studied the model analytically: they proved that a small forcing may indeed excite the parametric 3:2 resonance, but they have not explained the physical nature of the forcing. Here we integrate their equations numerically, dropping the ad hoc forcing, and adding instead a stochastic term to mimic the action of the very complex processes that occur in disks as, for example, MRI turbulence. We demonstrate that the presence of the stochastic term triggers the resonance in epicyclic oscillations of nearly Keplerian disks, and influences their pattern.

In this note we re-propose the arguments presented in Vio (2005) examining the superiority of the 'bi-parametric scale adaptive filter' (BSAF) when dealing with source detection as claimed by (Lopez-Caniego 2005a) and confirmed in Lopez-Caniego (2005b). We suggest here that the dispute can be easily solved if these authors provide the community with a detailed derivation of a basic equation.

In this note we stress the necessity of a careful check of both the theoretical arguments and the numerical experiments used by Lopez-Caniego et al. (2005) to support the superior performances of the biparametric scale adaptive filter (BSAF) with respect to the classic matched filter (MF) in the detection of sources on a random Gaussian background.

In this paper we consider the problem of the limits concerning the physical information that can be extracted from the analysis of one or more time series (light curves) typical of astrophysical objects. On the basis of theoretical considerations and numerical simulations, we show that with no a priori physical model there are not so many possibilities to obtain interpretable results. For this reason, the practice to develop more and more sophisticated statistical methods of time series analysis is not very productive. Only techniques of data analysis developed in a specific physical context can be expected to provide useful results. The field of {\it stochastic dynamics} appears to be an useful framework for such an approach. In particular, it is shown that modelling the experimental time series by means of the stochastic differential equations (SDE) represents a valuable tool of analysis. For example, the use of SDE permits to make the analysis of a continuous signal independent from the frequency sampling with which the experimental time series have been obtained. In this respect, an efficient approach based on the extended Kalman-filter technique is presented. Freely downloadable software is made available.

It is well-known that the noise associated with the collection of an astronomical image by a CCD camera is, in large part, Poissonian. One would expect, therefore, that computational approaches that incorporate this a priori information will be more effective than those that do not. The Richardson-Lucy (RL) algorithm, for example, can be viewed as a maximum-likelihood (ML) method for image deblurring when the data noise is assumed to be Poissonian. Least-squares (LS) approaches, on the other hand, arises from the assumption that the noise is Gaussian with fixed variance across pixels, which is rarely accurate. Given this, it is surprising that in many cases results obtained using LS techniques are relatively insensitive to whether the noise is Poissonian or Gaussian. Furthermore, in the presence of Poisson noise, results obtained using LS techniques are often comparable with those obtained by the RL algorithm. We seek an explanation of these phenomena via an examination of the regularization properties of particular LS algorithms. In addition, a careful analysis of the RL algorithm yields an explanation as to why it is more effective than LS approaches for star-like objects, and why it provides similar reconstructions for extended objects. We finish with a convergence analysis of the RL algorithm. Numerical results are presented throughout the paper. It is important to stress that the subject treated in this paper is not academic. In fact, in comparison with many ML algorithms, the LS algorithms are much easier to use and to implement, often provide faster convergence rates, and are much more flexible regarding the incorporation of constraints on the solution. Consequently, if little to no improvement is gained in the use of an ML approach over an LS algorithm, the latter will often be the preferred approach.

We consider the estimation of the regularization parameter for the simultaneous deblurring of multiple noisy images via Tikhonov regularization. We approach the problem in three ways. We first reduce the problem to a single-image deblurring for which the regularization parameter can be estimated through a classic generalized cross-validation (GCV) method. A modification of this function is used for correcting the undersmoothing typical of the original technique. With a second method, we minimize an average least-squares fit to the images and define a new GCV function. In the last approach, we use the classical $GCV$ on a single higher-dimensional image obtained by concatanating all the images into a single vector. With a reliable estimator of the regularization parameter, one can fully exploit the excellent computational characteristics typical of direct deblurring methods, which, especially for large images, makes them competitive with the more flexible but much slower iterative algorithms. The performance of the techniques is analyzed through numerical experiments. We find that under the independent homoscedastic and Gaussian assumptions made on the noise, the three approaches provide almost identical results with the first single image providing the practical advantage that no new software is required and the same image can be used with other deblurring algorithms.

We consider the simultaneous deblurring of a set of noisy images whose point spread functions are different but known and spatially invariant, and the noise is Gaussian. Currently available iterative algorithms that are typically used for this type of problem are computationally expensive, which makes their application for very large images impractical. We present a simple extension of a classical least-squares (LS) method where the multi-image deblurring is efficiently reduced to a computationally efficient single-image deblurring. In particular, we show that it is possible to remarkably improve the ill-conditioning of the LS problem by means of stable operations on the corresponding normal equations, which in turn speed up the convergence rate of the iterative algorithms. The performance and limitations of the method are analyzed through numerical simulations. Its connection with a column weighted least-squares approach is also considered in an appendix.

In this draft we comment on the results concerning the performances of matched filters, scale adaptive filters and Mexican hat wavelet that recently appeared in literature in the context of point source detection in Cosmic Microwave Background maps. In particular, we show that, contrary to what has been claimed, the use of the matched filters still appear to be the most reliable and efficient method to disantangle point sources from the backgrounds, even when using detection criterion that, differently from the classic $n\sigma$ thresholding rule, takes into account not only the height of the peaks in the signal corresponding to the candidate sources but also their curvature.

In this second paper in a series dedicated to developing efficient numerical techniques for the deblurring Cosmic Microwave Background (CMB) maps, we consider the case of asymmetric point spread functions (PSF). Although conceptually this problem is not different from the symmetric case, there are important differences from the computational point of view because it is no longer possible to use some of the efficient numerical techniques that work with symmetric PSFs. We present procedures that permit the use of efficient techniques even when this condition is not met. In particular, two methods are considered: a procedure based on a Kronecker approximation technique that can be implemented with the numerical methods used with symmetric PSFs but that has the limitation of requiring only mildly asymmetric PSFs. The second is a variant of the classic Tikhonov technique that works even with very asymmetric PSFs but that requires discarding the edges of the maps. We provide details for efficient implementations of the algorithms. Their performance is tested on simulated CMB maps.

Digital deblurring of images is an important problem that arises in multifrequency observations of the Cosmic Microwave Background (CMB) where, because of the width of the point spread functions (PSF), maps at different frequencies suffer a different loss of spatial resolution. Deblurring is useful for various reasons: first, it helps to restore high frequency components lost through the smoothing effect of the instrument's PSF; second, emissions at various frequencies observed with different resolutions can be better studied on a comparable resolution; third, some map-based component separation algorithms require maps with similar level of degradation. Because of computational efficiency, deblurring is usually done in the frequency domain. But this approach has some limitations as it requires spatial invariance of the PSF, stationarity of the noise, and is not flexible in the selection of more appropriate boundary conditions. Deblurring in real space is more flexible but usually not used because of its high computational cost. In this paper (the first in a series on the subject) we present new algorithms that allow the use of real space deblurring techniques even for very large images. In particular, we consider the use of Tikhonov deblurring of noisy maps with applications to {\it PLANCK}. We provide details for efficient implementations of the algorithms. Their performance is tested on Gaussian and non-Gaussian simulated CMB maps, and PSFs with both circular and elliptical symmetry. Matlab code is made available.

We provide theoretical procedures and practical recipes to simulate non-Gaussian correlated, homogeneous random fields with prescribed marginal distributions and cross-correlation structure, either in a N-dimensional Cartesian space or on the celestial sphere. We illustrate our methods using far-infrared maps obtained with the Infrared Space Observatory. However, the methodology presented here can be used in other astrophysical applications that require modeling correlated features in sky maps, for example, the simulation of multifrequency sky maps where backgrounds, sources and noise are correlated and can be modeled by random fields.

It is the purpose of the paper to describe the virtues of time-frequency methods for signal processing applications, having astronomical time series in mind. Different methods are considered and their potential usefulness respectively drawbacks are discussed and illustrated by examples. As areas where one can hope for a successful application of joint time-frequency analysis (JTFA), we describe specifically the problem of signal denoising as well as the question of signal separation which allows to separate signals (possibly overlapping in time or frequency, but) which are living on disjoint parts of the time-frequency plane. Some recipes for practical use of the algorithms are also provided.

Point-source contamination in high-precision Cosmic Microwave Background (CMB) maps severely affects the precision of cosmological parameter estimates. Among the methods that have been proposed for source detection, wavelet techniques based on ``optimal'' filters have been proposed.In this paper we show that these filters are in fact only restrictive cases of a more general class of matched filters that optimize signal-to-noise ratio and that have, in general, better source detection capabilities, especially for lower amplitude sources. These conclusions are confirmed by some numerical experiments. \keywords{Methods: data analysis -- Methods: statistical

In this paper we will consider the problem of the numerical simulation of non-Gaussian, scalar random fields with a prescribed correlation structure provided either by a theoretical model or computed on a set of observational data. Although, the numerical generation of a generic, non-Gaussian random field is a trivial operation, the task becomes tough when constraining the field with a prefixed correlation structure. At this regards, three numerical methods, useful for astronomical applications, are presented. The limits and capabilities of each method are discussed and the pseudo-codes describing the numerical implementation are provided for two of them.

We present an analysis of the ultraviolet and X-Ray variability of the Broad- Line-Radio Galaxy 3C390.3 over 15 years. The UV continuum showed large variations with amplitudes of up to a factor of 10. We find: (1) The variations of CIV and Ly_alpha are highly correlated with the UV continuum, and are delayed with respect to the continuum variations by 50-110 days with the red wing of both CIV and Ly_alpha , leading the blue wing; (2) The CIV/Ly_alpha ratio is positively correlated with both the continuum flux and UV line strength, a behavior different from other AGNs studied so far; (3) The blue sides of the Ly_alpha and CIV profiles are similar to the blue side of the Balmer lines, while the red sides are different, suggesting a different origin for the red peak in the Balmer lines. Our results suggest: (1) The broad CIV and Ly_alpha emitting gas is infalling towards the central object; (2) The overall behaviour of the CIV/Ly_alpha ratio and the absence of a big blue bump, strongly indicate the coexistence of optically thick as well as optically thin BLR clouds; (3) Assuming circular symmetry and predominantly circular motion, the BLR gas is situated at 83$\pm$25 lightdays from the central source; (4) Under these assumptions and with the derived circular velocity of $v_{rot} \simeq$ 2850 km s$^{-1}$, the central mass inside this radius is confined to $1.3~10^8 M\odot < M_{CM}< 4.0~10^8 M\odot$; (6) Comparing our results with those obtained from VLBI and observations of the Fe $K\alpha$ line, suggests the association of the BLR with a disk, inclined at 98$\pm$12 degrees with respect to the direction of superluminal motion of the radio blobs.

The long-term variability of a sample of 486 optically selected QSOs in the fields of the SA94, SA57 and SGP has been studied. The relations of variability with luminosity and redshift have been investigated by means of statistical estimators that are ``robust'' and allow at the same time to eliminate the influence of the measurement errors. The analysis in the QSOs rest frame of both the ensemble structure function ($SF$) and individual variability indices show: 1) a negative correlation between variability and luminosity; 2) a significant positive correlation between variability and $z$; 3) such correlations may be equally well parameterized either with a model in which the timescale of the variability is fixed for all the QSOs ($\tau \sim 2.4$ yr), while the amplitude linearly increases with the absolute magnitude and $z$, or with a model in which the timescale of the variability linearly depends on the absolute magnitude and the amplitude is only a function of $z$. The same analysis carried out in the observer's frame shows: 1) a negative correlation between variability and luminosity; 2) the timescale of variability does not depend significantly either on the absolute magnitude or on $z$; 3) the ensemble structure function is well represented by a parameterization in which, with a fixed timescale of about 5.5 yr, the amplitude linearly increases with the absolute magnitude; 4) although the general behaviour of the $SF$ does not show a systematic variation of the timescale and/or amplitude with $z$, if we examine the average variability index for objects with $-25>M_B>-27$, we find that below $z=1$ QSOs are significantly less variable than at higher redshift. Implications in terms of the black-hole, starburst and microlensing models are briefly discussed.

This is the first paper in a series aimed at defining a statistically significant sample of QSOs in the range $ 15 < B < 18.75$ and $ 0.3 < z < 2.2$. The selection is carried out using direct plates obtained at the ESO and UK Schmidt Telescopes, scanned with the COSMOS facility and searched for objects with an ultraviolet excess. Follow-up spectroscopy, carried out at ESO La Silla, is used to classify each candidate. In this initial paper, we describe the scientific objectives of the survey; the selection and observing techniques used. We present the first sample of 285 QSOs ($M_B < -23$) in a 153 deg$^2$ area, covered by the six ``deep'' fields, intended to obtain significant statistics down $B \simeq 18.75$ with unprecedented photometric accuracy. From this database, QSO counts are determined in the magnitude range $ 17 < B < 18.75$.

The issue of the X-ray to optical luminosity relationship ($L_o-L_x$) is addressed for both optically and X-ray selected quasar samples. We have applied a generalized regression algorithm for the case of samples involving censored data, with errors on both the dependent and the independent variable. Contrary to some previous results, we find that such relationship is consistent with being a linear one ($L_x\propto L_o$). We argue that previous reports of non-linear relationships (i.e. $L_x\propto L_o^e$ with $e < 1$) are due to the neglect of the influence of the photometric errors, the precise knowledge of which strongly influences the reliability of the results. Further progresses in the determination of the $L_o-L_x$ relationship can be achieved with ROSAT observations of the new generation of large bright quasar surveys.