posteriorProbability: Calculate the posterior probability for integer copy numbers.
In benilton/crlmm: Genotype Calling (CRLMM) and Copy Number Analysis tool for Affymetrix SNP 5.0 and 6.0 and Illumina arrays

Description Usage Arguments Details Value Note Author(s) See Also Examples

Calculate the posterior probability for integer copy numbers using the bivariate normal prediction regions.

1	posteriorProbability(object, predictRegion, copyNumber = 0:4, w)

`object`	A `CNSet` object.
`predictRegion`	A list containing the bivariate normal prediction region for each of the possible genotypes.
`copyNumber`	Integer vector.
`w`	numeric vector of prior probabilities for each of the copy number states. Must be the same length as `copyNumber` and sum to 1.

This is currently under development.

An array (features x samples x copy number)

This is under development. Use at your own risk.

R. Scharpf

predictionRegion, genotypes

data(cnSetExample)
pr <- predictionRegion(cnSetExample, copyNumber=0:4)
pp <- posteriorProbability(cnSetExample, predictRegion=pr)
dim(pp)

## multiple batches
data(cnSetExample2)
pr <- predictionRegion(cnSetExample2, copyNumber=0:4)
pp <- posteriorProbability(cnSetExample2, predictRegion=pr)