R/plotC.R

Defines functions plotC

plotC <- function (object, columns = NULL) 
{
     if (is.null(object[["C"]])) stop("No random parameters to plot.")
  
     C <- object[["C"]]
     gNIV <- ncol(C) - 2
     numIDs <- nrow(C)
     graphics.off()
     par(ask = TRUE)
     
     if (is.null(columns)) columns <- 1:gNIV + 2
     
     for (i in columns) {
          plot(density(C[, i]), type = "l", main = paste(colnames(C)[i], "\n",
               sum(C[,i] >= 0), " (", round(sum(C[,i] >= 0) / length(C[,i]) * 100, 1),
               "%) >= 0", sep = ""), 
               xlab = "Utility", ylab = "Density")
               dev.flush()
     }
     par(ask = FALSE)
}

Try the RSGHB package in your browser

Any scripts or data that you put into this service are public.

RSGHB documentation built on May 29, 2017, 6:23 p.m.