\HeaderA{twilight.pval}{A function to compute p-values from expression sets}{twilight.pval}
\keyword{nonparametric}{twilight.pval}
\begin{Description}\relax
A function to compute two-sample t, Z and fold change equivalent test statistics (paired or unpaired) and correlation coefficients. Based on permutations, expected test statistics as given in Tusher et al. (2001) and empirical \eqn{p}{}-values are computed. Additional output are \eqn{q}{}-values computed as given in Storey and Tibshirani (2003). The resulting object is of class \code{twilight} and can be passed to functions \code{twilight} or \code{plot.twilight}.
\end{Description}
\begin{Usage}
\begin{verbatim}
twilight.pval(xin, yin, method = "fc", paired = FALSE, B = 1000, yperm = NULL, balance = FALSE, quant.ci = 0.95, s0=NULL, verbose = TRUE, filtering = FALSE)
\end{verbatim}
\end{Usage}
\begin{Arguments}
\begin{ldescription}
\item[\code{xin}] Either an expression set (\code{exprSet}) or a data matrix with rows corresponding to features and columns corresponding to samples. 
\item[\code{yin}] A numerical vector containing class labels. The higher label denotes the case, the lower label the control samples to test case vs. control. For correlation scores, \code{yin} can be any numerical vector of length equal to the number of samples. 
\item[\code{method}] Character string: \code{"fc"} for fold change equivalent test (that is log ratio test), \code{"t"} for t-test, and \code{"z"} for Z-test. With \code{"pearson"} or \code{"spearman"}, the test statistic is either Pearson's correlation coefficient or Spearman's rank correlation coefficient.
\item[\code{paired}] Logical value. Depends on whether the samples are paired. Ignored if \code{method="pearson"} or \code{method="spearman"}.
\item[\code{B}] Numerical value specifying the number of permutations. 
\item[\code{yperm}] Optional matrix containing in each row a permutation of the class labels in binary(!) format for two-sample testing. For computation of correlation scores, the rows of \code{yperm} have to contain the appropriate values or ranks. Use this argument carefully! If \code{yperm} is specified, no other permutation will be done. Please note that the first row of \code{yperm} MUST be the input vector \code{yin}. Otherwise, the \eqn{p}{}-value calculation will be incorrect.
\item[\code{balance}] Logical value. Depends on whether balanced or unbalanced permutations should be done. Ignored if \code{method="pearson"} or \code{method="spearman"}.
\item[\code{quant.ci}] Probability value for confidence lines. Lines are symmetric and denote the \code{quant.ci}-quantile of maximal absolute differences between each permutatin and the expected scores. 
\item[\code{s0}] Fudge factor for variance correction in the Z-test. Takes effect only if \code{method="z"}. If \code{s0=NULL}: The fudge factor is set to the median of the pooled standard deviations.
\item[\code{verbose}] Logical value for message printing. 
\item[\code{filtering}] Logical value for filtering for permutations of class labels that produce a set of complete null scores. Invokes function \code{twilight.filtering}. If \code{yperm} is specified, no filtering will be done. Note that the filtering is done on unbalanced permutations even if \code{balance=TRUE}. 
\end{ldescription}
\end{Arguments}
\begin{Details}\relax
See vignette.
\end{Details}
\begin{Value}
Returns a \code{twilight} object consisting of a \code{data.frame} named \code{result} with variables
\begin{ldescription}
\item[\code{observed}] Observed test statistics.
\item[\code{expected}] Mean of order statistics of the permutation statistics.
\item[\code{candidate}] Binary vector. \code{"1"} for observations exceeding the confidence lines.
\item[\code{pvalue}] Empirical \eqn{p}{}-values from two-sided hypothesis tests.
\item[\code{qvalue}] \eqn{q}{}-values computed as described in Storey and Tibshirani (2003).
\item[\code{index}] Index of the original ordering.
\end{ldescription}

Values are sorted by absolute \code{observed} scores.
\end{Value}
\begin{Note}\relax
Additional output consists of

\Tabular{ll}{
\code{ci.line} & Quantile corresponding to \code{quant.ci}, passed for plotting. \\
\code{pi0} & Estimated prior probability. \\
\code{call} & Character string of function arguments. \\
\code{quant.ci} & Passes \code{quant.ci} for plotting. \\
}

The remaining slots are left free for function \code{twilight}.
\end{Note}
\begin{Author}\relax
Stefanie Scheid \url{http://www.molgen.mpg.de/~scheid}
\end{Author}
\begin{References}\relax
Scheid S and Spang R (2004): A stochastic downhill search algorithm for estimating the local false discovery rate, \emph{IEEE TCBB} \bold{1(3)}, 98--108.

Storey JD and Tibshirani R (2003): Statistical significance for genomewide studies, \emph{PNAS} \bold{100(16)}, 9440--9445.

Tusher VG, Tibshirani R and Chu G (2001): Significance analysis of mircroarrays applied to the ionizing response, \emph{PNAS} \bold{98(9)}, 5116--5121.
\end{References}
\begin{SeeAlso}\relax
\code{\LinkA{twilight}{twilight}}, \code{\LinkA{plot.twilight}{plot.twilight}}, \code{\LinkA{twilight.combi}{twilight.combi}}, \code{\LinkA{twilight.filtering}{twilight.filtering}}, \code{\LinkA{expval}{expval}}
\end{SeeAlso}
\begin{Examples}
\begin{ExampleCode}
data(expval)  # twilight object created by twilight.pval on Golub data set
expval
plot(expval,which="scores")
plot(expval,which="qvalues")
\end{ExampleCode}
\end{Examples}

