Description Usage Arguments Author(s) Examples
Raw MA objects always contains probes with the same gene,or lots of control probes which are considered not differentially expressed in every arrays,and all of this need to deal with before MA go to the lmFit pipeline.I believe MApipeline() give a fast way to do this job. MApipeline2 is the plus version of MApipeline().It annotates MA object via sequence enhanced annotation strategy derectly..
1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 | MApipeline(MA,
MA.probeid,
MA.colnames,
control.probes,
array.annotation,
eset,
enhanced.annotation=T,
gpl.path=NULL,
control.weight = F)
MApipeline2(MA,
MA.probeid,
MA.colnames,
control.probes,
array.annotation,
eset,
gpl.path=NULL,
control.weight = F)
|
MA |
a MA List |
MA.probeid |
the probeid of M matrix or A matrix in MA List.Note that it was the same compose with eset feature data. |
MA.colnames |
the new colnames of M matrix or A matrix in MA List.Note that it was the subset of colnames of eset exprs data. |
control.probes |
a list containing control probes.See AddDesignList() |
array.annotation |
a list containing array annotation information.See AddDesignList() |
eset |
a project from getElist() |
control.weight |
whether to gather the expression of control probes and weight them with 2.In most time we recommand control.weight.F.You can also compared two results and decide which one is the best choice. |
Weibin Huang<654751191@qq.com>
1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 | ## example
library(lucky)
load("E:/RCloud/database/DataDownload/GEOqueryDownload/GSE17154/GSE17154_rawEset.rda")
load("E:/RCloud/database/DataDownload/GEOqueryDownload/GSE17154/GSE17154_MAList.rda")
load("E:/RCloud/database/DataDownload/GEOqueryDownload/GSE17154/targets.rda")
## parameters
MA.probeid=1:44544
MA.colnames = Fastextra(targets$FileName,"_",1)
control.probes <- list(
l1=list(control.col="Reporter Group[role]",control.symbol = "Control"))
array.annotation = list(
probeid.col="ID",
symbol.cols=c("Symbol"),
sequence.col="SEQUENCE",
anno.cols=c("SYMBOL"),
db.anno=common.annot)
## Runing
MA1 <- MApipeline(MA,
MA.probeid,
MA.colnames,
control.probes,
array.annotation,
eset,
control.weight = F)
MA2 <- MApipeline(MA,
MA.probeid,
MA.colnames,
control.probes,
array.annotation,
eset,
control.weight = T)
ma1 <- MApipeline2(MA,
MA.probeid,
MA.colnames,
control.probes,
array.annotation,
eset,
gpl.path=NULL,
control.weight = T)
ma1 <- MApipeline2(MA,
MA.probeid,
MA.colnames,
control.probes,
array.annotation,
eset,
gpl.path=NULL,
control.weight = F)
|
Add the following code to your website.
For more information on customizing the embed code, read Embedding Snippets.