mergeExternalData: Merge external data
In FRASER: Find RAre Splicing Events in RNA-Seq Data

Description Usage Arguments Details Value Examples

To boost its own sequencing data, one can download existing and precounted data. This function merges the existing FraserDataSet with external count data.

1	mergeExternalData(fds, countFiles, sampleIDs, annotation = NULL)

`fds`	A `FraserDataSet`
`countFiles`	A character vector of file names pointing to the external count data. The vector has to be names or the files have to start with `k_j`, `k_theta`, `n_psi3`, `n_psi5`, `n_theta`.
`sampleIDs`	The samples to be merged from the external data.
`annotation`	A sample annotation of the external data (optional).

For more details on existing datasets have a look at: <https://github.com/gagneurlab/drop#datasets>

Since FRASER can not hand NA values, the merge will return only the intersecting regions and will drop any non overlapping features. This has to be kept in mind when analysing rare disease samples.

Merged FraserDataSet object.

anno <- fread(system.file("extdata", "externalCounts", 
        "annotation.tsv.gz", package="FRASER"))
ctsFiles <- list.files(full.names = TRUE, pattern="counts",
        system.file("extdata", "externalCounts", package="FRASER"))

fds <- createTestFraserDataSet()
fds_merged <- mergeExternalData(fds, ctsFiles, anno[,sampleID], anno)

K(fds, "psi5")
K(fds_merged, "psi5")