sig.score: Function to compute signature scores as linear combination of...
In genefu: Computation of Gene Expression-Based Signatures in Breast Cancer

Description Usage Arguments Value Author(s) References Examples

This function computes a signature score from a gene list (aka gene signature), i.e. a signed average as published in Sotiriou et al. 2006 and Haibe-Kains et al. 2009.

1 2	sig.score(x, data, annot, do.mapping = FALSE, mapping, size = 0, cutoff = NA, signed = TRUE, verbose = FALSE)

`x`	Matrix containing the gene(s) in the gene list in rows and at least three columns: "probe", "EntrezGene.ID" and "coefficient" standing for the name of the probe, the NCBI Entrez Gene id and the coefficient giving the direction and the strength of the association of each gene in the gene list.
`data`	Matrix of gene expressions with samples in rows and probes in columns, dimnames being properly defined.
`annot`	Matrix of annotations with at least one column named "EntrezGene.ID", dimnames being properly defined.
`do.mapping`	`TRUE` if the mapping through Entrez Gene ids must be performed (in case of ambiguities, the most variant probe is kept for each gene), `FALSE` otherwise.
`mapping`	Matrix with columns "EntrezGene.ID" and "probe" used to force the mapping such that the probes are not selected based on their variance.
`size`	Integer specifying the number of probes to be considered in signature computation. The probes will be sorted by absolute value of coefficients.
`cutoff`	Only the probes with coefficient greater than `cutoff` will be considered in signature computation.
`signed`	`TRUE` if only the sign of the coefficient must be considered in signature computation, `FALSE` otherwise.
`verbose`	`TRUE` to print informative messages, `FALSE` otherwise.

`score`	Signature score.
`mapping`	Mapping used if necessary.
`probe`	If mapping is performed, this matrix contains the correspondence between the gene list (aka signature) and gene expression data.

Benjamin Haibe-Kains

Sotiriou C, Wirapati P, Loi S, Harris A, Bergh J, Smeds J, Farmer P, Praz V, Haibe-Kains B, Lallemand F, Buyse M, Piccart MJ and Delorenzi M (2006) "Gene expression profiling in breast cancer: Understanding the molecular basis of histologic grade to improve prognosis", Journal of National Cancer Institute, 98:262-272

Haibe-Kains B (2009) "Identification and Assessment of Gene Signatures in Human Breast Cancer", PhD thesis at Universite Libre de Bruxelles, http://theses.ulb.ac.be/ETD-db/collection/available/ULBetd-02182009-083101/

## load NKI data
data(nkis)
## load GGI signature
data(sig.ggi)
## make of ggi signature a gene list
ggi.gl <- cbind(sig.ggi[ ,c("probe", "EntrezGene.ID")],
  "coefficient"=ifelse(sig.ggi[ ,"grade"] == 1, -1, 1))
## computation of signature scores
ggi.score <- sig.score(x=ggi.gl, data=data.nkis, annot=annot.nkis,
  do.mapping=TRUE, signed=TRUE, verbose=TRUE)
str(ggi.score)

Loading required package: survcomp
Loading required package: survival
Loading required package: prodlim
Loading required package: mclust
Package 'mclust' version 5.4.2
Type 'citation("mclust")' for citing this R package in publications.
Loading required package: limma
Loading required package: biomaRt
Loading required package: iC10
Loading required package: pamr
Loading required package: cluster
Loading required package: iC10TrainingData
Loading required package: AIMS
Loading required package: e1071
Loading required package: Biobase
Loading required package: BiocGenerics
Loading required package: parallel

Attaching package: 'BiocGenerics'

The following objects are masked from 'package:parallel':

    clusterApply, clusterApplyLB, clusterCall, clusterEvalQ,
    clusterExport, clusterMap, parApply, parCapply, parLapply,
    parLapplyLB, parRapply, parSapply, parSapplyLB

The following object is masked from 'package:limma':

    plotMA

The following objects are masked from 'package:stats':

    IQR, mad, sd, var, xtabs

The following objects are masked from 'package:base':

    Filter, Find, Map, Position, Reduce, anyDuplicated, append,
    as.data.frame, basename, cbind, colMeans, colSums, colnames,
    dirname, do.call, duplicated, eval, evalq, get, grep, grepl,
    intersect, is.unsorted, lapply, lengths, mapply, match, mget,
    order, paste, pmax, pmax.int, pmin, pmin.int, rank, rbind,
    rowMeans, rowSums, rownames, sapply, setdiff, sort, table, tapply,
    union, unique, unsplit, which, which.max, which.min

Welcome to Bioconductor

    Vignettes contain introductory material; view with
    'browseVignettes()'. To cite Bioconductor, see
    'citation("Biobase")', and for packages 'citation("pkgname")'.

probe candidates: 54/128
List of 3
 $ score  : Named num [1:150] -0.188 0.169 -0.043 -0.239 -0.234 ...
  ..- attr(*, "names")= chr [1:150] "NKI_123" "NKI_327" "NKI_291" "NKI_370" ...
 $ mapping: Named int [1:2] 54 128
  ..- attr(*, "names")= chr [1:2] "mapped" "total"
 $ probe  : chr [1:54, 1:3] "201584_s_at" "201710_at" "202094_at" "202107_s_at" ...
  ..- attr(*, "dimnames")=List of 2
  .. ..$ : chr [1:54] "201584_s_at" "201710_at" "202094_at" "202107_s_at" ...
  .. ..$ : chr [1:3] "probe" "EntrezGene.ID" "new.probe"