champ.CNA {ChAMP} | R Documentation

## Inference of Copy Number Abberrations from intensity values.

### Description

This function enables CNA profiles to be built using methylation data from Illumina HumanMethylation450K and HumanMethylationEPIC BeadChips. This function provide options to find Copy Number Abberrations between two phenotype (.e.g. Cancer & Normal), or the function would take the average value of your dataset as control and detect if some value are out of average status. For user want to detect abberrations between phenotypes, they can specify controlGroup in parameter, or they can simply used packaged dataset as control. Two kinds of plot would be returned, the abberrations of each sample, and the abberrations of each phenotype. The older version of ChAMP provide batchcorrect for intensity dataset, but it's nolonger provided here, user may use champ.runCombat() function to correct batch effect just like they correct beta matrix.

### Usage

        champ.CNA(intensity=myLoad$intensity,
                  pheno=myLoad$pd$Sample_Group,
                  control=TRUE,
                  controlGroup="champCtls",
                  sampleCNA=TRUE,
                  groupFreqPlots=TRUE,
                  Rplot=FALSE,
                  PDFplot=TRUE,
                  freqThreshold=0.3,
                  resultsDir="./CHAMP_CNA",
                  arraytype="450K")

### Arguments

`intensity` | A matrix of intensity values for each sample. (default = myLoad$intensity)
`pheno` | This is a categorical vector representing phenotype of factor wish to be analysed, for example "Cancer", "Normal"...Tow or even more phenotypes are allowed. (default = myLoad$pd$Sample_Group)
`control` | If champ.CNA() should calculate the difference between groups(controls and case) of not(with average). (default = TRUE)
`controlGroup` | which phenotype in your pheno parameter shall be treated as control type is you want to comparision between two groups. If this value was missing or invalid, the function would automatically use packaged Blood sample(champCtls) as control. (default = "champCtls")
`sampleCNA` | If sampleCNA=TRUE, then each sample's Copy Number Abberrations would be calculated and plotted. (default = TRUE)
`groupFreqPlots` | If groupFreqPlots=TRUE, then each group's Copy Number Abberrations Frequence would be calculated and plotted. (default = TRUE)
`freqThreshold` | If groupFreqPlots=T, then freqThreshold will be used as the cutoff for calling a gain or loss. (default = 0.3)
`PDFplot` | If PDFplot would be generated and save in resultsDir. (default = TRUE)
`Rplot` | If Rplot would be generated and save in resultsDir. Note if you are doing analysis on a server remotely, please make sure the server could connect your local graph applications. (For example X11 for linux.) (default = TRUE)
`arraytype` | Choose microarray type is 450K or EPIC.
`resultsDir` | The directory where PDF files would be saved. (default = "./CHAMP_CNA/")

### Value

`sampleResult` | The Copy Number Abberrations result calculated and ploted for each Sample.
`groupResult` | The Copy Number Abberrations result calculated and ploted for each Group.

### Author(s)

Feber, A
adapted by Yuan Tian

### References

Feber, A et. al. (2014). CNA profiling using high density DNA methylation arrays. Genome Biology.

### Examples

        ## Not run:
            myLoad <- champ.load(directory=system.file("extdata",package="ChAMPdata"))
            myCNA <- champ.CNA()

    ## End(Not run)

[Package ChAMP version 2.12.4 Index]