Extract private sequences — getPrivate • AnalyzAIRR

This function extract private sequences from the whole dataset, i.e. sequences found exclusively in one sample.

Usage

getPrivate(
  x,
  level = c("ntCDR3", "aaCDR3", "ntClone", "aaClone"),
  singletons = FALSE
)

Arguments

x: an object of class RepSeqExperiment
level: a character specifying the level of the repertoire to be taken into account. Should be one of aaClone","ntClone", "ntCDR3" or "aaCDR3".
singletons: a boolean indicating whether or not private sequences with a count of 1 should be extracted. Default is FALSE.

Value

an object of class RepSeqExperiment composed exclusively of private sequences.

Examples


data(RepSeqData)

privateclones <- getPrivate(RepSeqData,
                            level = "ntClone",
                            singletons = FALSE)