Description Usage Arguments Details Value Author(s) See Also Examples
View source: R/utility03282012.r
Merge microarray data sets in possibly irregular order.
1 | MetaDE.merge(x,MVperc=0)
|
x |
a list of studies. Each study is a list with components:
|
MVperc |
a threshold to specificy which genes are remained for further analysis. The default is zero. |
The gene expression data sets may be in possibly irregular order with different numbers of genes. This function is used to extact the common genes across studies. The merged data sets have the same genes in the same order.
When we combine a large of number of studies, the number of common genes may be very samll, so we allow to include somes gene appearing in most studies and missing in few studies. The default is zero which means that we only include genes appearing in all the studies.
a list of studies. Each study is a list with components:
x: the gene expression matrix.
y: the outcome.
censoring.status: the censoring status. This only for survival data.
Jia Li and Xingbin Wang
MetaDE.Read
, MetaDE.filter
,ind.analysis
and MetaDE.rawdata
1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 | #================example test MetaDE.merge========================================================#
label1<-rep(0:1,each=5)
label2<-rep(0:1,each=5)
time1=c(4,3,1,1,2,2,3,10,5,4)
event1=c(1,1,1,0,1,1,0,0,0,1)
exp1<-cbind(matrix(rnorm(5*20),20,5),matrix(rnorm(5*20,2),20,5))
exp2<-cbind(matrix(rnorm(5*20),20,5),matrix(rnorm(5*20,1.5),20,5))
rownames(exp1)<-paste("g1",1:20,sep="_")
rownames(exp2)<-paste("g2",1:20,sep="_")
symbol1<-sample(c("SST","VGF","CNP","LPA"),20,replace=TRUE)
symbol2<-sample(c("SST","VGF","CNP","APOE"),20,replace=TRUE)
study1<-cbind(c(NA,NA,symbol1),rbind(rbind(time1,event1),exp1))
study2<-cbind(c(NA,symbol2),rbind(label2,exp2))
setwd(tempdir())
write.table(study1,"study1.txt",sep="\t")
write.table(study2,"study2.txt",sep="\t")
mydata.raw<-MetaDE.Read(c("study1","study2"),via="txt",skip=c(2,1),log=FALSE)
mydata.matched<-MetaDE.match(mydata.raw,"IQR")
mydata.merged<-MetaDE.merge(mydata.matched)
|
Loading required package: survival
Loading required package: impute
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 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, cbind, colMeans, colSums, colnames, 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")'.
Loading required package: combinat
Attaching package: 'combinat'
The following object is masked from 'package:utils':
combn
Loading required package: tools
Add the following code to your website.
For more information on customizing the embed code, read Embedding Snippets.