data("STATegRa_S3")
require(ggplot2)
B1 <- createOmicsExpressionSet(Data=Block1.PCA,pData=ed.PCA,
pDataDescr=c("classname"))
B2 <- createOmicsExpressionSet(Data=Block2.PCA,
pData=ed.PCA,pDataDescr=c("classname"))
# Omics components analysis
discoRes <- omicsCompAnalysis(Input=list(B1,B2),Names=c("expr","mirna"),
method="DISCOSCA",Rcommon=2,Rspecific=c(2,2),
center=TRUE,scale=TRUE,weight=TRUE)
jiveRes <- omicsCompAnalysis(Input=list(B1,B2),Names=c("expr","mirna"),
method="JIVE",Rcommon=2,Rspecific=c(2,2),
center=TRUE,scale=TRUE,weight=TRUE)
# DISCO-SCA plotVAF
plotVAF(discoRes)
# JIVE plotVAF
plotVAF(jiveRes)
Run the code above in your browser using DataLab