共查询到20条相似文献，搜索用时 41 毫秒
1.
针对包含饱和样本数据的频数幂律分布拟合,提出一个新的幂律分布指数的极大似然估计方法的修正公式.对比研究显示,修正公式适用于剔除异常饱和值的幂律频数拟合.如果不剔除饱和值,幂律指数的估计只能使用修正前的公式,其误差随幂律指数变化,指数较小时误差较大.由此建议,对于包含饱和样本的频数分布拟合,首先剔除异常的饱和值,然后对剩余不含饱和值的子集使用修正公式进行参数估计. 相似文献
2.
3.
T. D. Kitching A. Amara F. B. Abdalla B. Joachimi A. Refregier 《Monthly notices of the Royal Astronomical Society》2009,399(4):2107-2128
In the absence of any compelling physical model, cosmological systematics are often misrepresented as statistical effects and the approach of marginalizing over extra nuisance systematic parameters is used to gauge the effect of the systematic. In this article, we argue that such an approach is risky at best since the key choice of function can have a large effect on the resultant cosmological errors.
As an alternative we present a functional form-filling technique in which an unknown, residual, systematic is treated as such. Since the underlying function is unknown, we evaluate the effect of every functional form allowed by the information available (either a hard boundary or some data). Using a simple toy model, we introduce the formalism of functional form filling. We show that parameter errors can be dramatically affected by the choice of function in the case of marginalizing over a systematic, but that in contrast the functional form-filling approach is independent of the choice of basis set.
We then apply the technique to cosmic shear shape measurement systematics and show that a shear calibration bias of | m ( z )| ≲ 10−3 (1 + z )0.7 is required for a future all-sky photometric survey to yield unbiased cosmological parameter constraints to per cent accuracy.
A module associated with the work in this paper is available through the open source icosmo code available at http://www.icosmo.org . 相似文献
As an alternative we present a functional form-filling technique in which an unknown, residual, systematic is treated as such. Since the underlying function is unknown, we evaluate the effect of every functional form allowed by the information available (either a hard boundary or some data). Using a simple toy model, we introduce the formalism of functional form filling. We show that parameter errors can be dramatically affected by the choice of function in the case of marginalizing over a systematic, but that in contrast the functional form-filling approach is independent of the choice of basis set.
We then apply the technique to cosmic shear shape measurement systematics and show that a shear calibration bias of | m ( z )| ≲ 10
A module associated with the work in this paper is available through the open source icosmo code available at http://www.icosmo.org . 相似文献
4.
A. Pagliaro V. Antonuccio-Delogu U. Becciani M. Gambera 《Monthly notices of the Royal Astronomical Society》1999,310(3):835-841
We present and discuss a method to identify substructures in combined angular-redshift samples of galaxies within clusters. The method relies on the use of the discrete wavelet transform (hereafter DWT) and has already been applied to the analysis of the Coma cluster. The main new ingredient of our method with respect to previous studies lies in the fact that we make use of a 3D data set rather than a 2D one. We test the method on mock cluster catalogues with spatially localized substructures and on a N -body simulation. Our main conclusion is that our method is able to identify the existing substructures provided that: (a) the subclumps are detached in part or all of the phase space, (b) one has a statistically significant number of redshifts, increasing as the distance decreases due to redshift distortions; (c) one knows a priori the scale on which substructures are to be expected. We have found that to allow an accurate recovery we must have both a significant number of galaxies (≈200 for clusters at z ≥0.4 or about 800 at z ≤0.4) and a limiting magnitude for completeness m B =16.
The only true limitation to our method seems to be the necessity of knowing a priori the scale on which the substructure is to be found. This is an intrinsic drawback of the method and no improvement in numerical codes based on this technique could make up for it. 相似文献
The only true limitation to our method seems to be the necessity of knowing a priori the scale on which the substructure is to be found. This is an intrinsic drawback of the method and no improvement in numerical codes based on this technique could make up for it. 相似文献
5.
Stéphane Colombi rew Jaffe Dmitri Novikov Christophe Pichon 《Monthly notices of the Royal Astronomical Society》2009,393(2):511-526
A method to rapidly estimate the Fourier power spectrum of a point distribution is presented. This method relies on a Taylor expansion of the trigonometric functions. It yields the Fourier modes from a number of fast Fourier transforms (FFTs), which is controlled by the order N of the expansion and by the dimension D of the system. In three dimensions, for the practical value N = 3 , the number of FFTs required is 20.
We apply the method to the measurement of the power spectrum of a periodic point distribution that is a local Poisson realization of an underlying stationary field. We derive an explicit analytic expression for the spectrum, which allows us to quantify – and correct for – the biases induced by discreteness and by the truncation of the Taylor expansion, and to bound the unknown effects of aliasing of the power spectrum. We show that these aliasing effects decrease rapidly with the order N . For N = 3 , they are expected to be, respectively, smaller than ∼10−4 and 0.02 at half the Nyquist frequency and at the Nyquist frequency of the grid used to perform the FFTs. The only remaining significant source of errors is reduced to the unavoidable cosmic/sample variance due to the finite size of the sample.
The analytical calculations are successfully checked against a cosmological N -body experiment. We also consider the initial conditions of this simulation, which correspond to a perturbed grid. This allows us to test a case where the local Poisson assumption is incorrect. Even in that extreme situation, the third-order Fourier–Taylor estimator behaves well, with aliasing effects restrained to at most the per cent level at half the Nyquist frequency.
We also show how to reach arbitrarily large dynamic range in Fourier space (i.e. high wavenumber), while keeping statistical errors in control, by appropriately 'folding' the particle distribution. 相似文献
We apply the method to the measurement of the power spectrum of a periodic point distribution that is a local Poisson realization of an underlying stationary field. We derive an explicit analytic expression for the spectrum, which allows us to quantify – and correct for – the biases induced by discreteness and by the truncation of the Taylor expansion, and to bound the unknown effects of aliasing of the power spectrum. We show that these aliasing effects decrease rapidly with the order N . For N = 3 , they are expected to be, respectively, smaller than ∼10
The analytical calculations are successfully checked against a cosmological N -body experiment. We also consider the initial conditions of this simulation, which correspond to a perturbed grid. This allows us to test a case where the local Poisson assumption is incorrect. Even in that extreme situation, the third-order Fourier–Taylor estimator behaves well, with aliasing effects restrained to at most the per cent level at half the Nyquist frequency.
We also show how to reach arbitrarily large dynamic range in Fourier space (i.e. high wavenumber), while keeping statistical errors in control, by appropriately 'folding' the particle distribution. 相似文献
6.
Dae-Won Kim Pavlos Protopapas Charles Alcock Yong-Ik Byun Federica B. Bianco 《Monthly notices of the Royal Astronomical Society》2009,397(1):558-568
We present a detrending algorithm for the removal of trends in time series. Trends in time series could be caused by various systematic and random noise sources such as cloud passages, changes of airmass, telescope vibration, CCD noise or defects of photometry. Those trends undermine the intrinsic signals of stars and should be removed. We determine the trends from subsets of stars that are highly correlated among themselves. These subsets are selected based on a hierarchical tree clustering algorithm. A bottom-up merging algorithm based on the departure from normal distribution in the correlation is developed to identify subsets, which we call clusters. After identification of clusters, we determine a trend per cluster by weighted sum of normalized light curves. We then use quadratic programming to detrend all individual light curves based on these determined trends. Experimental results with synthetic light curves containing artificial trends and events are presented. Results from other detrending methods are also compared. The developed algorithm can be applied to time series for trend removal in both narrow and wide field astronomy. 相似文献
7.
8.
9.
Chris Koen 《Monthly notices of the Royal Astronomical Society》2005,361(3):887-896
The theory of low-order linear stochastic differential equations is reviewed. Solutions to these equations give the continuous time analogues of discrete time autoregressive time-series. Explicit forms for the power spectra and covariance functions of first- and second-order forms are given. A conceptually simple method is described for fitting continuous time autoregressive models to data. Formulae giving the standard errors of the parameter estimates are derived. Simulated data are used to verify the performance of the methods. Irregularly spaced observations of the two hydrogen-deficient stars FQ Aqr and NO Ser are analysed. In the case of FQ Aqr the best-fitting model is of second order, and describes a quasi-periodicity of about 20 d with an e-folding time of 3.7 d. The NO Ser data are best fitted by a first-order model with an e-folding time of 7.2 d. 相似文献
10.
11.
12.
Coryn A. L. Bailer-Jones Mike Irwin & Ted von Hippel 《Monthly notices of the Royal Astronomical Society》1998,298(2):361-377
We investigate the application of neural networks to the automation of MK spectral classification. The data set for this project consists of a set of over 5000 optical (3800–5200 Å) spectra obtained from objective prism plates from the Michigan Spectral Survey. These spectra, along with their two-dimensional MK classifications listed in the Michigan Henry Draper Catalogue, were used to develop supervised neural network classifiers. We show that neural networks can give accurate spectral type classifications (σ68 = 0.82 subtypes, σrms = 1.09 subtypes) across the full range of spectral types present in the data set (B2–M7). We show also that the networks yield correct luminosity classes for over 95 per cent of both dwarfs and giants with a high degree of confidence. Stellar spectra generally contain a large amount of redundant information. We investigate the application of principal components analysis (PCA) to the optimal compression of spectra. We show that PCA can compress the spectra by a factor of over 30 while retaining essentially all of the useful information in the data set. Furthermore, it is shown that this compression optimally removes noise and can be used to identify unusual spectra. This paper is a continuation of the work carried out by von Hippel et al. (Paper I). 相似文献
13.
14.
15.
In this article we describe a case study of how NOAO is considering improving its management of Target‐of‐Opportunity (ToO) observations by integrating VOEvent into the flow of activities. We believe that using VOEvent to help document and track the use of ToO time will improve the user experience of ToOs at NOAO. It will also greatly aid in the management of the process and of the resulting data, allowing us to better track the ownership and provenance of the data and any resulting data products. Finally, it will provide an important method of archival access to the data and data “collections,” which might include not only processed data from a single VOEvent triggered observation but could also include multiple observations traceable to a single (or set of related) VOEvents. (© 2008 WILEY‐VCH Verlag GmbH & Co. KGaA, Weinheim) 相似文献
16.
17.
多任务学习方法在机器学习、计算机视觉、人工智能领域已得到广泛关注,利用任务间的相关性,将多个任务同时学习的效果优于每个任务单独学习的情况.采用多任务Lasso回归法(Multi-task Lasso Regression)用于恒星光谱物理参量的估计,不仅可以获取不同物理参量间的共同的特征信息,而且也可以很好地保留不同物理参量的特有的补充信息.使用恒星大气模拟模型合成光谱库ELODIE中的光谱数据和美国大型巡天项目Sloan发布的SDSS实测光谱数据进行实验,模型估算精度优于相关文献中的方法,特别是对重力加速度(lg g)和化学丰度([Fe/H])的估计.实验中通过改变光谱的分辨率,施加不同信噪比(SNR)的噪声,来说明模型的稳定性强.结果表明,模型精度受光谱分辨率和噪声的影响,但噪声对其影响更大,可见,多任务Lasso回归法不仅操作简便,稳定性强,而且也提高了模型的整体预测精度. 相似文献
18.
The auroras on Jupiter and Saturn can be studied with a high sensitivity and resolution by the Hubble Space Telescope ( HST ) ultraviolet (UV) and far-ultraviolet Space Telescope Imaging Spectrograph (STIS) and Advanced Camera for Surveys (ACS) instruments. We present results of automatic detection and segmentation of Jupiter's auroral emissions as observed by the HST ACS instrument with the VOronoi Image SEgmentation (VOISE). VOISE is a dynamic algorithm for partitioning the underlying pixel grid of an image into regions according to a prescribed homogeneity criterion. The algorithm consists of an iterative procedure that dynamically constructs a tessellation of the image plane based on a Voronoi diagram, until the intensity of the underlying image within each region is classified as homogeneous. The computed tessellations allow the extraction of quantitative information about the auroral features, such as mean intensity, latitudinal and longitudinal extents and length-scales. These outputs thus represent a more automated and objective method of characterizing auroral emissions than manual inspection. 相似文献
19.
F. Feroz M. P. Hobson M. Bridges 《Monthly notices of the Royal Astronomical Society》2009,398(4):1601-1614
We present further development and the first public release of our multimodal nested sampling algorithm, called M ulti N est . This Bayesian inference tool calculates the evidence, with an associated error estimate, and produces posterior samples from distributions that may contain multiple modes and pronounced (curving) degeneracies in high dimensions. The developments presented here lead to further substantial improvements in sampling efficiency and robustness, as compared to the original algorithm presented in Feroz & Hobson, which itself significantly outperformed existing Markov chain Monte Carlo techniques in a wide range of astrophysical inference problems. The accuracy and economy of the M ulti N est algorithm are demonstrated by application to two toy problems and to a cosmological inference problem focusing on the extension of the vanilla Λ cold dark matter model to include spatial curvature and a varying equation of state for dark energy. The M ulti N est software, which is fully parallelized using MPI and includes an interface to C osmo MC, is available at http://www.mrao.cam.ac.uk/software/multinest/ . It will also be released as part of the SuperBayeS package, for the analysis of supersymmetric theories of particle physics, at http://www.superbayes.org . 相似文献
20.
S.L. Allen 《Astronomische Nachrichten》2008,329(3):298-300
A trivial modification to the XML schema of VOEvent v1.1 allows the inclusion of W3C digital signatures. Signatures enable identification, identification enables trust, and trust enables authorization. Such changes would inhibit abuse of the VOEvent networks. (© 2008 WILEY‐VCH Verlag GmbH & Co. KGaA, Weinheim) 相似文献