Included EBSeq for downstream differential expression analysis

[rsem.git] / EBSeq / R / GetMultiPP.R

diff --git a/EBSeq/R/GetMultiPP.R b/EBSeq/R/GetMultiPP.R
new file mode 100644 (file)
index 0000000..4ae2882
--- /dev/null
+++ b/EBSeq/R/GetMultiPP.R
@@ -0,0 +1,6 @@
+GetMultiPP <- function(EBout){
+       PP=EBout$PPDE   
+       MAP=colnames(EBout$f)[apply(EBout$f,1,which.max)]
+       AllParti=EBout$AllParti
+       out=list(PP=PP, MAP=MAP,Patterns=AllParti)
+}