classify: Apply the binary predictor to FSA peaks

In FSAtools: Fragment Analysis and Capillary Sequencing Tool Kit

Apply the binary predictor to FSA peaks

Description

Predict to which class the sample is most likely to belong, using a modified LPS model.

Usage

  classify(x, plot = TRUE)

Arguments

x	A fsa object with peaks and model attributes.
plot	Single logical value, whether to plot a visual representation of the prediction or not.

Value

Returns a list :

score	The raw score used to make the prediction.
p	The probability to belong to each of the two groups.
class	The final prediction, as a group name. May be NA if no probability passes the model threshold.

Author(s)

Sylvain Mareschal

References

Radmacher MD, McShane LM, Simon R. A paradigm for class prediction using gene expression profiles. J Comput Biol. 2002;9(3):505-11.

Wright G, Tan B, Rosenwald A, Hurt EH, Wiestner A, Staudt LM. A gene expression-based method to diagnose clinically distinct subgroups of diffuse large B cell lymphoma. Proc Natl Acad Sci U S A. 2003 Aug 19;100(17):9991-6.

Bohers E, Mareschal S, Bouzelfen A, Marchand V, Ruminy P, Maingonnat C, Menard AL, Etancelin P, Bertrand P, Dubois S, Alcantara M, Bastard C, Tilly H, Jardin F. Targetable activating mutations are very frequent in GCB and ABC diffuse large B-cell lymphoma. Genes Chromosomes Cancer. 2014 Feb;53(2):144-53.

See Also

read.fsa, peaks.fsa, add.model, generic.process

Examples

  # Example FSA file provided
  fsa <- read.fsa(system.file("extdata/fsa_GEP/A5918.fsa", package="FSAtools"))
  fsa <- align.fsa(fsa)
  
  # Add model from design file
  design <- designFile(system.file("extdata/design_GEP.conf", package="FSAtools"))
  fsa <- add.model(fsa, model=design$GLOBALS$MODEL)
  
  # Add peak heights
  fsa <- peaks.fsa(fsa, peaks=design$GLOBALS$PEAKS)
  
  # Classify sample
  fsa <- classify(fsa, plot=TRUE)
  print(attr(fsa, "classification"))

FSAtools documentation built on Aug. 19, 2023, 1:06 a.m.