data(simudata2)
n<-nrow(simudata2)
m<-ncol(simudata2)
#the expression data is from the second column to m
exprs.simu2<-as.matrix(simudata2[,2:m])
#prepare the group ID number for each sample array
data(simu2_groupid)
#prepare the design matrix for each group of samples
data(simu2_compgroup)
#fit 2 correlation motif to the data
motif.fitted<-cormotiffit(exprs.simu2, simu2_groupid,simu2_compgroup,K=2)
plotMotif(motif.fitted)
Run the code above in your browser using DataLab