git-svn-id: file:///home/git/hedgehog.fhcrc.org/bioconductor/trunk/madman/Rpacks/MineICA@100850 bc3139a8-67e5-0310-9ffc-ced21a209358
... | ... |
@@ -23,10 +23,6 @@ |
23 | 23 |
\alias{organism<-} |
24 | 24 |
\alias{organism,IcaSet-method} |
25 | 25 |
\alias{organism<-,IcaSet-method} |
26 |
-\alias{organism<-,IcaSet,character-method} |
|
27 |
-\alias{setOrganism,IcaSet-method} |
|
28 |
-\alias{getOrganism,IcaSet-method} |
|
29 |
- |
|
30 | 26 |
|
31 | 27 |
\alias{mart} |
32 | 28 |
\alias{mart<-} |
git-svn-id: file:///home/git/hedgehog.fhcrc.org/bioconductor/trunk/madman/Rpacks/MineICA@73179 bc3139a8-67e5-0310-9ffc-ced21a209358
1 | 1 |
new file mode 100644 |
... | ... |
@@ -0,0 +1,341 @@ |
1 |
+\name{IcaSet} |
|
2 |
+\docType{class} |
|
3 |
+\alias{class:IcaSet} |
|
4 |
+\alias{IcaSet} |
|
5 |
+\alias{IcaSet-class} |
|
6 |
+ |
|
7 |
+%\alias{dat} % pour la generique |
|
8 |
+%\alias{dat,IcaSet-method}% # avec la signature de ma classe |
|
9 |
+%\alias{dat<-} %# la generique de remplacement |
|
10 |
+%\alias{dat<-,IcaSet,matrix-method} %# avec la signature de la classe |
|
11 |
+ |
|
12 |
+\alias{[} |
|
13 |
+\alias{[,ANY,ANY,IcaSet-method} |
|
14 |
+\alias{[,IcaSet,ANY-method} |
|
15 |
+\alias{[,IcaSet,ANY,ANY-method} |
|
16 |
+\alias{[,IcaSet,ANY,ANY,ANY-method} |
|
17 |
+\alias{[<-} |
|
18 |
+\alias{[<-,IcaSet,ANY,ANY,ANY,ANY-method} |
|
19 |
+\alias{[<-,IcaSet,ANY,ANY,ANY-method} |
|
20 |
+\alias{[<-,IcaSet,ANY,ANY-method} |
|
21 |
+ |
|
22 |
+\alias{organism} |
|
23 |
+\alias{organism<-} |
|
24 |
+\alias{organism,IcaSet-method} |
|
25 |
+\alias{organism<-,IcaSet-method} |
|
26 |
+\alias{organism<-,IcaSet,character-method} |
|
27 |
+\alias{setOrganism,IcaSet-method} |
|
28 |
+\alias{getOrganism,IcaSet-method} |
|
29 |
+ |
|
30 |
+ |
|
31 |
+\alias{mart} |
|
32 |
+\alias{mart<-} |
|
33 |
+\alias{mart,IcaSet-method} |
|
34 |
+\alias{mart<-,IcaSet-method} |
|
35 |
+\alias{mart<-,IcaSet,character-method} |
|
36 |
+\alias{setMart,IcaSet-method} |
|
37 |
+\alias{getMart,IcaSet-method} |
|
38 |
+ |
|
39 |
+\alias{chipVersion} |
|
40 |
+\alias{chipVersion,IcaSet-method} |
|
41 |
+\alias{chipVersion<-} |
|
42 |
+\alias{chipVersion<-,IcaSet-method} |
|
43 |
+\alias{chipVersion<-,IcaSet,character-method} |
|
44 |
+ |
|
45 |
+%\alias{indComp} |
|
46 |
+%\alias{indComp<-} |
|
47 |
+\alias{indComp,IcaSet-method} |
|
48 |
+\alias{setIndComp,IcaSet-method} |
|
49 |
+\alias{getIndComp,IcaSet-method} |
|
50 |
+\alias{indComp<-,IcaSet,character-method} |
|
51 |
+%\alias{compNames} |
|
52 |
+%\alias{compNames<-} |
|
53 |
+\alias{compNames,IcaSet-method} |
|
54 |
+\alias{setLabelsComp,IcaSet-method} |
|
55 |
+\alias{getLabelsComp,IcaSet-method} |
|
56 |
+\alias{compNames<-,IcaSet,character-method} |
|
57 |
+%\alias{witGenes} |
|
58 |
+%\alias{witGenes<-} |
|
59 |
+\alias{witGenes,IcaSet-method} |
|
60 |
+%\alias{setWitGenes,IcaSet-method} |
|
61 |
+%\alias{getWitGenes,IcaSet-method} |
|
62 |
+\alias{witGenes<-,IcaSet,character-method} |
|
63 |
+\alias{refSamples} |
|
64 |
+\alias{refSamples<-} |
|
65 |
+\alias{refSamples,IcaSet-method} |
|
66 |
+\alias{setRefSamples,IcaSet-method} |
|
67 |
+\alias{getRefSamples,IcaSet-method} |
|
68 |
+\alias{refSamples<-,IcaSet-method} |
|
69 |
+\alias{refSamples<-,IcaSet,character-method} |
|
70 |
+\alias{typeID} |
|
71 |
+\alias{typeID<-} |
|
72 |
+\alias{typeID,IcaSet-method} |
|
73 |
+\alias{typeID<-,IcaSet-method} |
|
74 |
+\alias{setTypeID,IcaSet-method} |
|
75 |
+\alias{getTypeID,IcaSet-method} |
|
76 |
+\alias{typeID<-,IcaSet,list-method} |
|
77 |
+\alias{chipManu} |
|
78 |
+\alias{chipManu<-} |
|
79 |
+\alias{chipManu<-,IcaSet-method} |
|
80 |
+\alias{chipManu<-,IcaSet,character-method} |
|
81 |
+\alias{chipManu,IcaSet-method} |
|
82 |
+\alias{setChipManu,IcaSet-method} |
|
83 |
+\alias{getChipManu,IcaSet-method} |
|
84 |
+ |
|
85 |
+ |
|
86 |
+\alias{Slist,IcaSet-method} |
|
87 |
+\alias{SlistByGene,IcaSet-method} |
|
88 |
+\alias{Alist,IcaSet-method} |
|
89 |
+ |
|
90 |
+ |
|
91 |
+ |
|
92 |
+\title{ |
|
93 |
+ Class to Contain and Describe an ICA decomposition of High-Throughput Data. |
|
94 |
+} |
|
95 |
+ |
|
96 |
+\description{ |
|
97 |
+ Container for high-throughput data and results of ICA decomposition |
|
98 |
+ obtained on these data. \code{IcaSet} class is derived from |
|
99 |
+ \code{\link{eSet}}, and requires a matrix named \code{dat} as |
|
100 |
+ \code{assayData} member. |
|
101 |
+} |
|
102 |
+\section{Extends}{ |
|
103 |
+ Directly extends class \code{\link{eSet}}. |
|
104 |
+} |
|
105 |
+\section{Creating Objects}{ |
|
106 |
+ |
|
107 |
+ \code{new("IcaSet")} |
|
108 |
+ |
|
109 |
+ \code{new("IcaSet", |
|
110 |
+ annotation = character(0), |
|
111 |
+ experimentData = new("MIAME"), |
|
112 |
+ featureData = new("AnnotatedDataFrame"), |
|
113 |
+ phenoData = new("AnnotatedDataFrame"), |
|
114 |
+ protocolData = phenoData[,integer(0)], |
|
115 |
+ dat = new("matrix"), |
|
116 |
+ A=new("data.frame"), |
|
117 |
+ S=new("data.frame"), ...) |
|
118 |
+ } |
|
119 |
+ |
|
120 |
+ This creates an \code{IcaSet} with \code{assayData} implicitly |
|
121 |
+ created to contain \code{dat}. %Additional named matrix arguments |
|
122 |
+ %with the same dimensions as \code{dat} are added to |
|
123 |
+ %\code{assayData}; the row and column names of these additional |
|
124 |
+ %matrices should match those of \code{dat}. |
|
125 |
+ |
|
126 |
+ \code{new("IcaSet", |
|
127 |
+ annotation = character(0), |
|
128 |
+ assayData = assayDataNew(dat=new("matrix")), |
|
129 |
+ experimentData = new("MIAME"), |
|
130 |
+ featureData = new("AnnotatedDataFrame"), |
|
131 |
+ phenoData = new("AnnotatedDataFrame"), |
|
132 |
+ protocolData = phenoData[,integer(0)], |
|
133 |
+ A=new("data.frame"), |
|
134 |
+ S=new("data.frame"), ...) |
|
135 |
+ } |
|
136 |
+ |
|
137 |
+ This creates an \code{IcaSet} with \code{assayData} provided |
|
138 |
+ explicitly. %In this form, the only required named arguments are \code{assayData}. |
|
139 |
+ |
|
140 |
+ |
|
141 |
+ \code{IcaSet} instances are usually created through |
|
142 |
+ \code{new("IcaSet", ...)}. Usually the arguments to \code{new} |
|
143 |
+ include \code{dat} ('features x samples', e.g a matrix of expression |
|
144 |
+ data), \code{phenoData} ('samples x annotations', a |
|
145 |
+ matrix of sample annotations), \code{S} the Source |
|
146 |
+ matrix of the ICA decomposition ('features x comp'), \code{A} the Mixing matrix of the ICA |
|
147 |
+ decomposition ('samples x comp'), \code{annotation} the annotation |
|
148 |
+ package, \code{typeID} the description of the feature and gene IDs. |
|
149 |
+ |
|
150 |
+ The other attributes can be missing, in which case |
|
151 |
+ they are assigned default values. |
|
152 |
+ |
|
153 |
+ The function \code{\link{buildIcaSet}} is a more convenient way to |
|
154 |
+ create \code{IcaSet} instances, and allows to automatically annotate |
|
155 |
+ the features. |
|
156 |
+} |
|
157 |
+\section{Slots}{ |
|
158 |
+ Inherited from \code{eSet}: |
|
159 |
+ \describe{ |
|
160 |
+ \item{\code{annotation}:}{See \code{\link{eSet}}} |
|
161 |
+ \item{\code{assayData}:}{Contains matrices with equal |
|
162 |
+ dimensions, and with column number equal to |
|
163 |
+ \code{nrow(phenoData)}. \code{assayData} must contain a matrix |
|
164 |
+ \code{dat} with rows representing features (e.g., reporters) |
|
165 |
+ and columns representing samples. Class:\code{\link{AssayData-class}}} |
|
166 |
+ \item{\code{experimentData}:}{See \code{\link{eSet}}} |
|
167 |
+ \item{\code{featureData}:}{See \code{\link{eSet}}} |
|
168 |
+ \item{\code{phenoData}:}{See \code{\link{eSet}}} |
|
169 |
+ \item{\code{protocolData}:}{See \code{\link{eSet}}} |
|
170 |
+ |
|
171 |
+ Specific slot: |
|
172 |
+ |
|
173 |
+ \item{\code{organism}:}{Contains the name of the species. Currently |
|
174 |
+ only Human ("Human" or "Homo sapiens") and Mouse ("Mouse" or "Mus |
|
175 |
+ musculus") are supported. Only used when \code{chipManu}="illumina" |
|
176 |
+ } |
|
177 |
+ \item{\code{mart}:}{An output of \code{\link[biomaRt]{useMart}} of package \code{biomaRt}. Only useful if no annotation package is available for argument \code{icaSet}. |
|
178 |
+ } |
|
179 |
+ |
|
180 |
+ \item{\code{datByGene}:}{Data.frame containing the data \code{dat} where |
|
181 |
+ features have been replaced by their annotations (e.g, gene IDs). Rows |
|
182 |
+ represent annotations of the features (e.g., gene IDs) and |
|
183 |
+ columns represent samples.} |
|
184 |
+ \item{\code{A}:}{The mixing matrix of the ICA decomposition, contained |
|
185 |
+ in a data.frame whose |
|
186 |
+ column number equals the number of components and row number equals |
|
187 |
+ \code{nrow(phenoData)} (dimension: 'samples x comp').} |
|
188 |
+ \item{\code{S}:}{The source matrix of the ICA decomposition, contained |
|
189 |
+ in a data.frame whose |
|
190 |
+ column number equals the number of components and row number equals |
|
191 |
+ \code{nrow(assayData)} (dimension: 'features x comp').} |
|
192 |
+ \item{\code{SByGene}:}{The matrix Source of the ICA decomposition, contained |
|
193 |
+ in a data.frame whose |
|
194 |
+ column number equals the number of components and row number equals |
|
195 |
+ \code{nrow(datByGene)} (dimension: 'annotatedFeatures x comp').} |
|
196 |
+ \item{\code{compNames}:}{A vector of component labels with length |
|
197 |
+ equal to the number of component.} |
|
198 |
+ \item{\code{indComp}:}{A vector of component indices with length |
|
199 |
+ equal to the number of component.} |
|
200 |
+ \item{\code{witGenes}:}{A vector of gene IDs with length |
|
201 |
+ equal to the number of component.} |
|
202 |
+ \item{\code{chipManu}:}{The manufacturer of the technology the |
|
203 |
+ data originates from. Useful for the annotation of the features when |
|
204 |
+ data originates from an _illumina_ microarray.} |
|
205 |
+ \item{\code{chipVersion}:}{The version of the chip, only useful for |
|
206 |
+ when \code{chipManu}="illumina"} |
|
207 |
+ \item{\code{refSamples}:}{A vector of sample IDs including the reference samples, e.g the "normal" samples. |
|
208 |
+ Must be included in \code{sampleNames(object)}, i.e in \code{colnames(dat)}.} |
|
209 |
+ \item{\code{typeID}:}{A vector of characters providing the annotation IDs. It includes |
|
210 |
+ three elements: |
|
211 |
+ \describe{ |
|
212 |
+ \item{geneID_annotation}{the IDs from the |
|
213 |
+ package to be used to annotate the features into genes. It will be used to |
|
214 |
+ fill the attributes \code{datByGene} and \code{SByGene} of the \code{icaSet}. |
|
215 |
+ It must match one of the objects the corresponding package supports |
|
216 |
+ (you can access the list of objects by typing ls("package:packagename")). If |
|
217 |
+ no annotation package is provided, this element is not useful.} |
|
218 |
+ \item{geneID_biomart}{the type of gene IDs, as available in |
|
219 |
+ \code{listFilters(mart)}; where mart is specified as described in \code{\link[biomaRt]{useMart}}. |
|
220 |
+ If you have directly built the IcaSet at the |
|
221 |
+ gene level (i.e if no annotation package is used), \code{featureID_biomart} and |
|
222 |
+ \code{geneID_biomart} will be identical.} |
|
223 |
+ \item{featureID_biomart}{the |
|
224 |
+ type of feature IDs, as available in \code{listFilters(mart)}; where |
|
225 |
+ \code{mart} is specified as described in function \code{\link[biomaRt]{useMart}}. |
|
226 |
+ Not useful if you work at the gene level.} }} |
|
227 |
+ |
|
228 |
+} |
|
229 |
+} |
|
230 |
+\section{Methods}{ |
|
231 |
+ |
|
232 |
+ Class-specific methods. |
|
233 |
+ \describe{ |
|
234 |
+ |
|
235 |
+ \item{\code{getComp(IcaSet, ind, |
|
236 |
+ level=c("features","genes"))}}{Given a component index, extract |
|
237 |
+ the corresponding sample contribution values from A, and the |
|
238 |
+ feature (\code{level}="features") or gene (\code{level}="genes") |
|
239 |
+ projections from S. Returns a list with two elements: |
|
240 |
+ \code{contrib} the sample contributions and \code{proj} the |
|
241 |
+ feature or gene projections.} |
|
242 |
+ |
|
243 |
+ Access and set any slot specific to IcaSet: |
|
244 |
+ \item{\code{slotName(IcaSet)}, and |
|
245 |
+ \code{slotName(IcaSet)<-}:}{Accessing and setting any slot |
|
246 |
+ of name \code{slotName} contained in an IcaSet object.} |
|
247 |
+ |
|
248 |
+ \item{\code{IcaSet["slotName"]}, and |
|
249 |
+ \code{IcaSet["slotName"]<-}:}{Accessing and setting any slot |
|
250 |
+ of name \code{slotName} contained in an IcaSet object.} |
|
251 |
+ |
|
252 |
+ Most used accessors and settors: |
|
253 |
+ \item{\code{A(IcaSet)}, and |
|
254 |
+ \code{A(IcaSet)<-}:}{Accessing and setting Mixing matrix \code{A}.} |
|
255 |
+ \item{\code{S(IcaSet)}, and |
|
256 |
+ \code{S(IcaSet)<-}:}{Accessing and setting |
|
257 |
+ the data.frame Source \code{S}.} |
|
258 |
+ \item{\code{Slist(IcaSet)}:}{Accessing |
|
259 |
+ the data.frame Source as a list where names are preserved.} |
|
260 |
+ \item{\code{SByGene(IcaSet)}, and |
|
261 |
+ \code{SByGene(IcaSet)<-}:}{Accessing |
|
262 |
+ and setting the _annotated_ data.frame Source \code{SByGene}.} |
|
263 |
+ \item{\code{SlistByGene(IcaSet)}:}{Accessing |
|
264 |
+ the _annotated_ Source matrix as a list where names are preserved.} |
|
265 |
+ \item{\code{organism(IcaSet)}, \code{organism(IcaSet,characte)<-}}{Access and |
|
266 |
+ set value in the \code{organism} slot.} |
|
267 |
+ \item{\code{dat(IcaSet)}, \code{dat(IcaSet,matrix)<-}}{Access and |
|
268 |
+ set elements named \code{dat} in the \code{AssayData-class} |
|
269 |
+ slot.} |
|
270 |
+ |
|
271 |
+ } |
|
272 |
+ |
|
273 |
+ Derived from \code{\link{eSet}}: |
|
274 |
+ \describe{ |
|
275 |
+ \item{\code{pData(IcaSet)}, \code{pData(IcaSet,value)<-}:}{See \code{\link{eSet}}} |
|
276 |
+ \item{\code{assayData(IcaSet)}:}{See \code{\link{eSet}}} |
|
277 |
+ \item{\code{sampleNames(IcaSet)} and \code{sampleNames(IcaSet)<-}:}{See \code{\link{eSet}}} |
|
278 |
+ \item{\code{featureNames(IcaSet)}, \code{featureNames(IcaSet, value)<-}:}{See \code{\link{eSet}}} |
|
279 |
+ \item{\code{dims(IcaSet)}:}{See \code{\link{eSet}}} |
|
280 |
+ \item{\code{phenoData(IcaSet)}, \code{phenoData(IcaSet,value)<-}:}{See \code{\link{eSet}}} |
|
281 |
+ \item{\code{varLabels(IcaSet)}, \code{varLabels(IcaSet, value)<-}:}{See \code{\link{eSet}}} |
|
282 |
+ \item{\code{varMetadata(IcaSet)}, \code{varMetadata(IcaSet,value)<-}:}{See \code{\link{eSet}}} |
|
283 |
+ \item{\code{varMetadata(IcaSet)}, \code{varMetadata(IcaSet,value)}}{See \code{\link{eSet}}} |
|
284 |
+ \item{\code{experimentData(IcaSet)},\code{experimentData(IcaSet,value)<-}:}{See \code{\link{eSet}}} |
|
285 |
+ \item{\code{pubMedIds(IcaSet)}, \code{pubMedIds(IcaSet,value)}}{See \code{\link{eSet}}} |
|
286 |
+ \item{\code{abstract(IcaSet)}:}{See \code{\link{eSet}}} |
|
287 |
+ \item{\code{annotation(IcaSet)}, \code{annotation(IcaSet,value)<-}}{See \code{\link{eSet}}} |
|
288 |
+ \item{\code{protocolData(IcaSet)}, \code{protocolData(IcaSet,value)<-}}{See \code{\link{eSet}}} |
|
289 |
+ \item{\code{combine(IcaSet,IcaSet)}:}{See \code{\link{eSet}}} |
|
290 |
+ \item{\code{storageMode(IcaSet)}, \code{storageMode(IcaSet,character)<-}:}{See \code{\link{eSet}}} |
|
291 |
+ } |
|
292 |
+ |
|
293 |
+ Standard generic methods: |
|
294 |
+ \describe{ |
|
295 |
+ \item{\code{initialize(IcaSet)}:}{Object instantiation, used |
|
296 |
+ by \code{new}; not to be called directly by the user.} |
|
297 |
+ |
|
298 |
+ \item{\code{validObject(IcaSet)}:}{Validity-checking method, ensuring |
|
299 |
+ that \code{dat} is a member of |
|
300 |
+ \code{assayData}, and that the number of features, genes, samples, |
|
301 |
+ and components are consistent across all the attributes of the |
|
302 |
+ IcaSet object. \code{checkValidity(IcaSet)} imposes this |
|
303 |
+ validity check, and the validity checks of \code{eSet}.} |
|
304 |
+ \item{\code{IcaSet[slotName]}, \code{IcaSet[slotName]<-}:}{Accessing |
|
305 |
+ and setting any slot of name \code{slotName} contained in an |
|
306 |
+ IcaSet object.} |
|
307 |
+ \item{\code{IcaSet[i, j, k]}:}{Extract object of class "IcaSet" |
|
308 |
+ for features or genes with names i, samples with names or |
|
309 |
+ indices j, and components with names or indices k.} |
|
310 |
+ \item{\code{makeDataPackage(object, author, email, packageName, packageVersion, license, biocViews, filePath, description=paste(abstract(object), collapse="\n\n"), ...)}}{ |
|
311 |
+ Create a data package based on an IcaSet object. See |
|
312 |
+ \code{\link{makeDataPackage}}.} |
|
313 |
+ \item{\code{show(IcaSet)}:}{See \code{\link{eSet}}} |
|
314 |
+ \item{\code{dim(IcaSet)}, \code{ncol}:}{See \code{\link{eSet}}} |
|
315 |
+ \item{\code{IcaSet[(index)]}:}{See \code{\link{eSet}}} |
|
316 |
+ \item{\code{IcaSet$}, \code{IcaSet$<-}:}{See \code{\link{eSet}}} |
|
317 |
+ \item{\code{IcaSet[[i]]}, \code{IcaSet[[i]]<-}:}{See \code{\link{eSet}}} |
|
318 |
+ } |
|
319 |
+} |
|
320 |
+ |
|
321 |
+\author{Anne Biton} |
|
322 |
+ |
|
323 |
+\seealso{ |
|
324 |
+ \code{\link{eSet-class}}, \code{\link{buildIcaSet}}, |
|
325 |
+ \code{\link{IcaSet-class}}, \code{\link{MineICAParams-class}}. |
|
326 |
+} |
|
327 |
+ |
|
328 |
+\examples{ |
|
329 |
+# create an instance of IcaSet |
|
330 |
+new("IcaSet") |
|
331 |
+dat <- matrix(runif(100000), nrow=1000, ncol=100) |
|
332 |
+rownames(dat) <- 1:nrow(dat) |
|
333 |
+new("IcaSet", |
|
334 |
+ dat=dat, |
|
335 |
+ A=as.data.frame(matrix(runif(1000), nrow=100, ncol=10)), |
|
336 |
+ S=as.data.frame(matrix(runif(10000), nrow=1000, ncol=10), row.names = 1:nrow(dat))) |
|
337 |
+ |
|
338 |
+ |
|
339 |
+} |
|
340 |
+ |
|
341 |
+\keyword{classes} |