doseresNMA antidep/drnma.plot.hist.R

# Histogram of DE parameters: beta1, beta2 and regression coef: g
histDRparam <- function(x=x,metareg=F){ # ,  drug.lab=drug.lab
  param.lab <-c('beta.1.','beta.2.')
  if(metareg){
    param.lab0 <- 'g.'
    param.lab <- append(param.lab,param.lab0)
  }
  sims <- x$BUGSoutput$sims.array
  sims.mat <- apply(sims,3L,c) # merge the 3 chains (columns into one)
  sims.param <- sims.mat%>%data.frame()%>%select(starts_with(param.lab))
  plotdata <- data.frame(sim=unlist(c(sims.param)),param=rep(colnames(sims.param),each=nrow(sims.param)))
  
  g_hist <- ggplot(plotdata, aes(x=sim)) +
    geom_histogram(alpha=0.6, binwidth = 0.001,color="darkblue", fill="lightblue")
  g_split <- g_hist + ggplot2::facet_wrap(~plotdata$param,scales = 'free')
  g_labs <- g_split + ggplot2::labs(y="Frequency", x="iterations")
  g_labs+theme(axis.line = element_line(colour = "black"),
           panel.grid.major = element_blank(),
           panel.grid.minor = element_blank())+theme_bw()
}
htx-r/doseresNMA documentation built on Jan. 28, 2021, 5:32 a.m.