alakazam (version 0.2.2)

countClones: Tabulates clones sizes

Description

countClones determines the number of sequences and total copy number of clonal groups.

Usage

countClones(data, groups = NULL, copy = NULL, clone = "CLONE")

Arguments

data

data.frame with Change-O style columns containing clonal assignments.

groups

character vector defining data columns containing grouping variables. If group=NULL, then do not group data.

copy

name of the data column containing copy numbers for each sequence. If this value is specified, then total copy abundance is determined by the sum of copy numbers within each clonal group.

clone

name of the data column containing clone identifiers.

Value

• A data.frame summarizing clone counts and frequencies with columns:
  • CLONE: clone identifier.
  • SEQ_COUNT: total number of sequences for the clone.
  • SEQ_FREQ: frequency of the clone as a fraction of the total number of sequences within each group.
  • COPY_COUNT: sum of the copy counts in thecopycolumn. Only present if thecopyargument is specified.
  • COPY_FREQ: frequency of the clone as a fraction of the total copy number within each group. Only present if thecopyargument is specified.
  Also includes additional columns specified in the groups argument.

Examples

# Load example data
file <- system.file("extdata", "ExampleDb.gz", package="alakazam")
df <- readChangeoDb(file)

# Without copy numbers
clones <- countClones(df, groups="SAMPLE")

# With copy numbers and multiple groups
clones <- countClones(df, groups=c("SAMPLE", "ISOTYPE"), copy="DUPCOUNT")