Fixed a bug in rsem-find-DE, which will crash the script if only one DE gene/isoform...

diff --git a/rsem-find-DE b/rsem-find-DE
index a69ebd427e67d1381691e6fde9b1ff552326bb75..98db6d8b91e141f148f3274fd305da107e182bb5 100755 (executable)
--- a/rsem-find-DE
+++ b/rsem-find-DE
@@ -58,7 +58,7 @@ if (is.null(ngvector)) {
 }
 stopifnot(!is.null(EBOut))
 
-PP <- GetPPMat(EBOut)
+PP <- as.data.frame(GetPPMat(EBOut))
 fc_res <- PostFC(EBOut)
 
 # soft threshold, default output