<!DOCTYPE HTML PUBLIC "-//W3C//DTD HTML 3.2 Final//EN"> <!--Converted with LaTeX2HTML 98.1p1 release (March 2nd, 1998) originally by Nikos Drakos (nikos@cbl.leeds.ac.uk), CBLU, University of Leeds * revised and updated by: Marcus Hennecke, Ross Moore, Herb Swan * with significant contributions from: Jens Lippmann, Marek Rouchal, Martin Wilck and others --> <HTML> <HEAD> <TITLE>Power of test statistics</TITLE> <META NAME="description" CONTENT="Power of test statistics"> <META NAME="keywords" CONTENT="vol2"> <META NAME="resource-type" CONTENT="document"> <META NAME="distribution" CONTENT="global"> <META HTTP-EQUIV="Content-Type" CONTENT="text/html; charset=iso-8859-1"> <LINK REL="STYLESHEET" HREF="vol2.css"> <LINK REL="next" HREF="node228.html"> <LINK REL="previous" HREF="node226.html"> <LINK REL="up" HREF="node222.html"> <LINK REL="next" HREF="node228.html"> </HEAD> <BODY > <!--Navigation Panel--> <A NAME="tex2html4248" HREF="node228.html"> <IMG WIDTH="37" HEIGHT="24" ALIGN="BOTTOM" BORDER="0" ALT="next" SRC="icons.gif/next_motif.gif"></A> <A NAME="tex2html4245" HREF="node222.html"> <IMG WIDTH="26" HEIGHT="24" ALIGN="BOTTOM" BORDER="0" ALT="up" SRC="icons.gif/up_motif.gif"></A> <A NAME="tex2html4239" HREF="node226.html"> <IMG WIDTH="63" HEIGHT="24" ALIGN="BOTTOM" BORDER="0" ALT="previous" SRC="icons.gif/previous_motif.gif"></A> <A NAME="tex2html4247" HREF="node1.html"> <IMG WIDTH="65" HEIGHT="24" ALIGN="BOTTOM" BORDER="0" ALT="contents" SRC="icons.gif/contents_motif.gif"></A> <BR> <B> Next:</B> <A NAME="tex2html4249" HREF="node228.html">Time domain analysis</A> <B> Up:</B> <A NAME="tex2html4246" HREF="node222.html">Basic principles of time</A> <B> Previous:</B> <A NAME="tex2html4240" HREF="node226.html">Corrections to the probability</A> <BR> <BR> <!--End of Navigation Panel--> <H2><A NAME="SECTION001725000000000000000"> Power of test statistics</A> </H2> <P> The question which method of the detection of features in the test statistic is the most sensitive, is of considerable importance for all practical TSA applications. It directly translates into the comparison of the power of different statistics. Let <I>Y</I> be a signal of some physical meaning, i.e. different from the pure noise. A power <IMG WIDTH="44" HEIGHT="72" ALIGN="MIDDLE" BORDER="0" SRC="img469.gif" ALT="$\beta$">of the test statistic is the probability that <I>H</I><SUB><I>o</I></SUB> is accepted for the given <I>Y</I>: <!-- MATH: $\beta = 1 - p_Y(S)$ --> <IMG WIDTH="98" HEIGHT="44" ALIGN="MIDDLE" BORDER="0" SRC="img470.gif" ALT="$\beta = 1 - p_Y(S)$">. In other words, this is the probability of the false rejection of the no-noise signal <I>Y</I> as pure noise. Generally, the relative power of various statistics depends on the type of the signal <I>Y</I>. An important practical consequence of this fact is that there is no such thing as the universally best method for time series analysis. A method that is good for one type of signals may be poor for another one. An extensive list of signal types and recommended methods is clearly out of scope of the present manual. However, in order to guide the user's intuition we quote below the results of a comparison of the powers of various test statistics for two types of oscillations, both of small amplitude and observed at more or less even intervals (Schwarzenberg-Czerny, 1989): <P> Let us first consider a sinusoidal oscillation with an amplitude not exceeding the noise. Then, all statistics based on models with 3 parameters have similar values of the test power: power spectrum, Scargle statistics, <IMG WIDTH="12" HEIGHT="43" ALIGN="MIDDLE" BORDER="0" SRC="img471.gif" ALT="$\chi^2(3)$"> fit of a sinusoid, ORT(1), AOV(3) and corrected PDM(3) statistics. (Please note that we depart here from the conventional notation by indicating in brackets the number <I>n</I><SUB><I>m</I></SUB> of the degrees of freedom of the model - e.g. the number of series terms or phase bins - instead of the number of the residual degrees of freedom <I>n</I><SUB><I>r</I></SUB>.) Statistics with more than 3 parameters, e.g. ORT(<I>n</I>/2+1), AOV(<I>n</I>), PDM(<I>n</I>) and <IMG WIDTH="14" HEIGHT="48" ALIGN="MIDDLE" BORDER="0" SRC="img472.gif" ALT="$\chi^2(n)$"> with <I>n</I>>3 and an extended Fourier series have less power. Our final choice is guided by the availability of the analytical probability distribution of the test statistics. Summing up, we recommend tu use for the detection of sinusoidal and other smooth oscillations of small amplitude the statistics with a coarse phase resolution, e.g. ORT(1), Scargle and AOV(3), . <P> For a narrow gaussian pulse or eclipse of width <I>w</I> repeating with period <I>P</I> the most powerful statistics are these with the matching resolution: ORT(<I>P</I>/2<I>w</I>), AOV(<I>P</I>/<I>w</I>) and <IMG WIDTH="62" HEIGHT="78" ALIGN="MIDDLE" BORDER="0" SRC="img473.gif" ALT="$\chi^2($">P/w). Power spectrum, Scargle, ORT(1), AOV(3), <IMG WIDTH="10" HEIGHT="48" ALIGN="MIDDLE" BORDER="0" SRC="img474.gif" ALT="$\chi^2(3)$">, ORT(n/2), AOV(n) and <IMG WIDTH="14" HEIGHT="43" ALIGN="MIDDLE" BORDER="0" SRC="img475.gif" ALT="$\chi^2(n)$">, <IMG WIDTH="50" HEIGHT="44" ALIGN="MIDDLE" BORDER="0" SRC="img476.gif" ALT="$n\gg P/w$"> all have less power. Note the equivalence of the <IMG WIDTH="12" HEIGHT="43" ALIGN="MIDDLE" BORDER="0" SRC="img477.gif" ALT="$\chi^2(3)$"> and Scargle's statistic (Lomb, 1976, Scargle, 1982) and the near-equivalence of the power spectrum and Scargle's statistics in the case of nearly uniformly sampled observations. Considering both test power and computational convenience we recommend for signals with sharp features, e.g. narrow pulses or eclipses, to use the ORT and AOV with the resolution matched to the width of these features. <P> <HR> <!--Navigation Panel--> <A NAME="tex2html4248" HREF="node228.html"> <IMG WIDTH="37" HEIGHT="24" ALIGN="BOTTOM" BORDER="0" ALT="next" SRC="icons.gif/next_motif.gif"></A> <A NAME="tex2html4245" HREF="node222.html"> <IMG WIDTH="26" HEIGHT="24" ALIGN="BOTTOM" BORDER="0" ALT="up" SRC="icons.gif/up_motif.gif"></A> <A NAME="tex2html4239" HREF="node226.html"> <IMG WIDTH="63" HEIGHT="24" ALIGN="BOTTOM" BORDER="0" ALT="previous" SRC="icons.gif/previous_motif.gif"></A> <A NAME="tex2html4247" HREF="node1.html"> <IMG WIDTH="65" HEIGHT="24" ALIGN="BOTTOM" BORDER="0" ALT="contents" SRC="icons.gif/contents_motif.gif"></A> <BR> <B> Next:</B> <A NAME="tex2html4249" HREF="node228.html">Time domain analysis</A> <B> Up:</B> <A NAME="tex2html4246" HREF="node222.html">Basic principles of time</A> <B> Previous:</B> <A NAME="tex2html4240" HREF="node226.html">Corrections to the probability</A> <!--End of Navigation Panel--> <ADDRESS> <I>Petra Nass</I> <BR><I>1999-06-15</I> </ADDRESS> </BODY> </HTML>