4d748181 |
\name{posteriorProbability}
\alias{posteriorProbability}
\alias{posteriorProbability,CNSet-method}
%- Also NEED an '\alias' for EACH other topic documented here.
\title{Calculate the posterior probability for integer copy numbers.}
\description{
Calculate the posterior probability for integer copy numbers using the
bivariate normal prediction regions.
}
\usage{
|
1821febe |
posteriorProbability(object, predictRegion, copyNumber = 0:4, w)
|
4d748181 |
}
%- maybe also 'usage' for other objects documented here.
\arguments{
\item{object}{
A \code{CNSet} object.
}
\item{predictRegion}{
A list containing the bivariate normal prediction region for each of
the possible genotypes.
}
\item{copyNumber}{
Integer vector.
|
1821febe |
}
\item{w}{ numeric vector of prior probabilities for each of the copy
number states. Must be the same length as \code{copyNumber} and sum
|
6dc7b945 |
to 1.}
|
1821febe |
|
4d748181 |
}
|
3b0b2109 |
\details{
This is currently under development.
|
4d748181 |
}
|
3b0b2109 |
|
4d748181 |
\value{
An array (features x samples x copy number)
}
\author{
R. Scharpf
}
|
3b0b2109 |
|
4d748181 |
\note{
|
3b0b2109 |
This is under development. Use at your own risk.
|
4d748181 |
}
%% ~Make other sections like Warning with \section{Warning }{....} ~
\seealso{
|
3b0b2109 |
\code{\link{predictionRegion}}, \code{\link{genotypes}} %\code{\link{calculatePosteriorMean}}
|
4d748181 |
}
\examples{
|
3b0b2109 |
data(cnSetExample)
pr <- predictionRegion(cnSetExample, copyNumber=0:4)
pp <- posteriorProbability(cnSetExample, predictRegion=pr)
|
aaef0d80 |
dim(pp)
|
8e21c728 |
## multiple batches
|
3b0b2109 |
data(cnSetExample2)
pr <- predictionRegion(cnSetExample2, copyNumber=0:4)
pp <- posteriorProbability(cnSetExample2, predictRegion=pr)
|
4d748181 |
}
% Add one or more standard keywords, see file 'KEYWORDS' in the
% R documentation directory.
\keyword{array}
|