14753e6b |
\name{annotFeaturesComp}
\alias{annotFeaturesComp}
\title{Features annotation}
\usage{
annotFeaturesComp(icaSet, params,
type = toupper(typeID(icaSet)["geneID_annotation"]),
featureId = typeID(icaSet)["featureID_biomart"],
geneId = typeID(icaSet)["geneID_biomart"])
}
\arguments{
\item{icaSet}{An object of class \code{\link{IcaSet}}
whose features have to be annotated. The attribute
\code{annotation} of this object contains the annotation
package to be used.}
\item{params}{An object of class
\code{\link{MineICAParams}} containing the parameters of
the analysis.}
\item{type}{The ID of the object of the annotation
package to be used for the annotation, must be available
in \code{ls("package:package_name")}}
\item{featureId}{The type of the feature IDs, in the
\code{biomaRt} way (type \code{listFilters(mart)} to
choose one). Used when \code{annotation(icaSet)} is of
length 0.}
\item{geneId}{The type of the gene IDs, in the
\code{biomaRt} way (type \code{listAttributes(mart)} to
choose one). Used when \code{annotation(icaSet)} is of
length 0.}
}
\value{
This function returns the argument \code{icaSet} with
attributes \code{SByGene} and \code{datByGene} filled.
}
\description{
##' This function annotates the features of an object of
class \code{\link{IcaSet}}, and fills its attributes
\code{SByGene} and \code{datByGene}.
}
\details{
This function is called by function
\code{\link{annotInGene}} which will check the validity
of the attributes \code{annotation, typeID, chipManu} and
eventually \code{chipVersion} of \code{icaSet}. If
available, the attribute \code{annotation} of argument
\code{icaSet} must be an annotation package and will be
used to annotate the \code{featureNames} of
\code{icaSet}. If attribute \code{annotation} of argument
\code{icaSet} is not available (of length 0),
\code{biomaRt} is used to annotate the features.
This function fills the attributes \code{SByGene} and
\code{datByGene} of the argument \code{icaSet}. When
several feature IDs are available for a same gene ID, the
median value of the corresponding features IDs is
attributed to the gene (the median of projection values
is used for attribute \code{SByGene}, and the median of
expression values is used for attribute
\code{datByGene}).
When attribute \code{chipManu} of the argument
\code{icaSet} is "illumina", the features are first
converted into nuID using the package 'lumi*Mapping' and
then annotated into genes. In that case, features can
only be annotated in ENTREZID or SYMBOL. It means that
\code{typeID(icaSet)['geneID_annotation']} must be either
'ENTREZID' or 'SYMBOL'. You will need to annotate
yourself the \code{\link{IcaSet}} object if you want to
use different IDs.
}
\examples{
## load an example of IcaSet
data(icaSetCarbayo)
params <- buildMineICAParams()
require(hgu133a.db)
####===================================================
## Use of annotation package contained in annotation(icaSet)
####====================================================
## annotation in SYMBOL
icaSetCarbayo_annot <- annotFeaturesComp(icaSet=icaSetCarbayo, params=params, type="SYMBOL")
# arg 'type' is optional since the function uses contents of typeID(icaSet) as the defaults,
# it is specified in these examples for pedagogy views
## annotation in Entrez Gene
icaSetCarbayo_annot <- annotFeaturesComp(icaSet=icaSetCarbayo, params=params, type="ENTREZID")
\dontrun{
####===================================================
## Use of biomaRt, when annotation(icaSet) is of length 0
####====================================================
## empty attribute 'annotation' of the IcaSet object
# when this attribute is not specified, biomaRt is used for annotation
annotation(icaSetCarbayo) <- character()
# make sure the mart attribute is correctly defined
mart(icaSetCarbayo) <- useMart(biomart="ensembl", dataset="hsapiens_gene_ensembl")
## make sure elements "featureID_biomaRt" and "geneID_biomaRt" of typeID(icaSet) are correctly filled
# they will be used by function 'annotFeaturesComp' through biomaRt to query the database
typeID(icaSetCarbayo)
## run annotation of HG-U133A probe set IDs into Gene Symbols using biomaRt
icaSetCarbayo_annot <- annotFeaturesComp(icaSet=icaSetCarbayo, params=params)
}
}
\author{
Anne Biton
}
\seealso{
\code{\link{annotFeatures}},
\code{\link{annotFeaturesWithBiomaRt}},
\code{\link{annotInGene}}
}
|