diff --git a/R/Software.R b/R/Software.R
index 0d6e934..0b98e15 100644
--- a/R/Software.R
+++ b/R/Software.R
@@ -745,9 +745,9 @@ get_predictive_Y<-function(mcmc,
     cs<- sample(unique_z, n, replace = T, prob =  pi_h[unique_z])
     etas = NULL
     #need to fix this
-    etas = as.matrix(do.call(rbind, data.frame(t(sapply(unique(cs),function(x) return(rbind(etas,
+    etas = do.call(rbind, data.frame(t(sapply(unique(cs),function(x) return(rbind(etas,
                                                                        mvrnorm(sum(cs == x),mu[[x]],
-                                                                               Delta[[x]]))))))))
+                                                                               Delta[[x]])))))))
 
   }else{
     etas = mvrnorm(n,mu[[unique_z]], Delta[[unique_z]])