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.

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

`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. |

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

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

Questions? Problems? Suggestions? Tweet to @rdrrHQ or email at ian@mutexlabs.com.

Please suggest features or report bugs with the GitHub issue tracker.

All documentation is copyright its authors; we didn't write any of that.