Calculation of copy numbers

Share:

Description

Starting from a matrix of coverages, the corresponding matrix of copy numbers is calculated. A null model for significance assessment of copy number alterations is fitted. Each amplicon in each sample is assessed for significance. P-values are reported in three ways: without multiple testing correction, correted for the number of samples and corrected for both, the number of samples and the number of amplicons.

Usage

1
2
calculate.CN(coverage, method.pooled = "amplicon", method.mt = "Bonferroni",
thres.cov = 100, thres.p = 0.05)

Arguments

coverage

A numeric matrix containing the coverages of each amlicon (rows) in each of samples (columns).

method.pooled

Method used for the estimation of the null model. Either one common null model for all amplicons ("pooled") or individual null models for each of the amplicons ("amplicon") are fitted. ("amplicon")

method.mt

Method for multiple testing correction: must be equal to "Bonferroni", "Holm", "Hochberg", "SidakSS", "SidakSD", "BH", "BY", "ABH" or "TSBH" as described in the package multtest.

thres.cov

Theshold for the minimal mean coverage of an amplicon to be included in the analysis.

thres.p

Significance level for calling copy number alterations.

Value

Matrix of copy numbers with the estimated null model ("model") and the p-values for each amplicon in each sample ("P") as attributes.

Examples

1
2
3
4
5
## Not run: 
data(coverage)
CN <- calculate.CN(coverage)

## End(Not run)