Skip to content

Commit

Permalink
Update Software.R
Browse files Browse the repository at this point in the history
  • Loading branch information
jfeldman396 authored Jan 4, 2023
1 parent 3008e61 commit 728bcf0
Showing 1 changed file with 2 additions and 2 deletions.
4 changes: 2 additions & 2 deletions R/Software.R
Original file line number Diff line number Diff line change
Expand Up @@ -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]])
Expand Down

0 comments on commit 728bcf0

Please sign in to comment.