\name{RankingBstat}
\alias{RankingBstat}
\title{Ranking based on the B-statistic}
\description{
 The B-statistic was motivated in a bayesian framework described by
 Lonnstedt and Speed (2002). It is implemented in the package \code{sma},
 the function is just a wrapper.\cr
 For \code{S4} method information, see \link{RankingBstat-methods}.    
}
\usage{
RankingBstat(x, y, type = c("paired", "onesample"), gene.names = NULL, ...)
}

\arguments{
  \item{x}{A \code{matrix} of gene expression values with rows
           corresponding to genes and columns corresponding to observations or alternatively an object of class \code{ExpressionSet}.\cr
           If \code{type = paired}, the first half of the columns corresponds to 
           the first measurements and the second half to the second ones. 
           For instance, if there are 10 observations, each measured twice,
           stored in an expression matrix \code{expr}, 
           then \code{expr[,1]} is paired with \code{expr[,11]}, \code{expr[,2]}
           with \code{expr[,12]}, and so on.}
  \item{y}{If \code{x} is a matrix, then \code{y} may be
           a \code{numeric} vector or a factor with at most two levels.\cr
           If \code{x} is an \code{ExpressionSet}, then \code{y}
           is a character specifying the phenotype variable in
           the output from \code{pData}.\cr
           If \code{type = paired}, take care that the coding is
           analogously to the requirement concerning \code{x}}
  \item{type}{\describe{
               \item{"paired":}{paired test. Take care that the coding of \code{y}
                                 is correct (s. above)}
                \item{"onesample":}{\code{y} has only one level. 
                                    Test whether the true mean is different
                                    from zero.}
                }
                \code{"unpaired"} is \emph{not} possible.}
  \item{gene.names}{An optional vector of gene names.}
  \item{\dots}{Furher arguments passed to \code{stat.bayesian}
               from the package \code{sma}}
}

\value{An object of class \link{GeneRanking}.}
\references{Lonnstedt, I., Speed, T. (2002). \cr
            Replicated microarray data.
            \emph{Statistica sinica, 12, 31-46}}
\author{Martin Slawski \email{martin.slawski@campus.lmu.de} \cr
        Anne-Laure Boulesteix \url{http://www.slcmsr.net/boulesteix}}
\seealso{
         \link{GetRepeatRanking}, \link{RankingTstat}, \link{RankingFC}, \link{RankingWelchT}, \link{RankingWilcoxon},
          \link{RankingBaldiLong}, \link{RankingFoxDimmic}, \link{RankingLimma}, 
          \link{RankingEbam}, \link{RankingWilcEbam}, \link{RankingSam}, 
      \link{RankingShrinkageT}, \link{RankingSoftthresholdT}, 
          \link{RankingPermutation}, \link{RankingGap}}
\keyword{univar}
\examples{
## Load toy gene expression data
data(toydata)
### class labels
yy <- toydata[1,]
### gene expression
xx <- toydata[-1,]
### run RankingBstat
Bstat <- RankingBstat(xx, yy, type="paired")
}