• Home
• Variations
• Box Plot
• Quantile Plot
• Scatter Plot
• Y~Y Plot
• Slug Plot
• Data Sets
• Documents & Tips

• Grid
• Admin. Time
• Box'n'Whisker
• Quantile
• a-b-line
• Lines
• Mean Line
• Loess
• 95% CI
• Points

Administration Range

Example

Code Used

#pdf(file="F:/Website/blocks/admin-large.pdf", width = 15, height = 7.5)
### Import Data  
  data <- warfarin.mdvpk.1
  split.column <- "AGE"
    
  data <- data[data$DOSE == 120,]

  data.select <- data[ data$MDV == 0 & data$AMT == 0, ] 
    ### Rows where MDV == 0 and AMT == 0  (Good Rows for plotting x,y)
    
  data.to.plot <- list()
  
  data.to.plot$x <- data.select[,"TIME"]  #x
  data.to.plot$y <- data.select[,"DV"]  #y
  data.to.plot$ids <- data.select[,"ID"]  #ids
  amt <- data[data$AMT==data$DOSE,c("TIME","ID")]
  amt[amt$ID==28,"TIME"] <- 48
  data.to.plot$amt <- amt
    # Done purposefully to show how the multi administration times works.
    # Yes, this picture is producing "false" data.  It is for example purposes only!
  
  data.to.plot$criteria <- data.select[,split.column]

  x.min.lim <- 0
  x.max.lim <- 150
  y.min.lim <- 0
  y.max.lim <- max(data.to.plot[["y"]])
  data.to.plot$x.lim <- c(x.min.lim,x.max.lim)
  data.to.plot$y.lim <- c(y.min.lim,y.max.lim)

  bw <- FALSE

  library(lattice)
### End Import


### Panel Function

  function.of.panel <- function(
    x,
    y,
    groups = groups,
    data.to.plot.in.panel = data.to.plot,
    black.and.white.in.panel = bw,
    subscripts
  ){
  cat("\n\n\nPacket", packet.number(),"\n")
  
  ### Grid
    plot.grid(
      x.divisions = 10,
      y.divisions = 10,
      grid.color = "grey95",
      x.limits = data.to.plot.in.panel[["x.lim"]],
      y.limits = data.to.plot.in.panel[["y.lim"]]
    )

  #### Administration Line
    plot.administration.time(
      data.for.amt = data.to.plot.in.panel,
      ids = groups,
      rows.being.used = subscripts,
      black.and.white = black.and.white.in.panel
    )
		
  }
      
### End Panel



### Make Scatter Plot
  lat.admin <- xyplot(
    y ~ x | as.factor(criteria), 
      # Remove '| as.factor(criteria)' if you do not want to split the data
    data = data.to.plot,
    type = NULL,
      # Taken Care of in Panel
    aspect=.5,
      # X is twice as long as Y
    groups = ids,
    xlab = "Time [units]",
    ylab = "Response [units]",
    xlim = c(x.min.lim - 0.03*(x.max.lim-x.min.lim)*0.5,x.max.lim + 0.03*(x.max.lim-x.min.lim)*0.5), 
      #0.5 is to compensate for the aspect ratio (aspect)
    ylim = c(y.min.lim - 0.03*(y.max.lim-y.min.lim), y.max.lim + 0.03*(y.max.lim-y.min.lim)),
    main = "Administration Time (Range) Example",
    panel = function.of.panel,
	
    key = list(
      title = "Legend",
      text = list(
        c("Administration Range"),
        col="black"
      ),
      points = list(
        col = c(if(bw){c("grey25")}else{c("green")}),
        pch = c(15),
        cex = c(2)*0.75
      ),
      lines = list(
        col = c(if(bw){c("grey25")}else{c("green")}),
        lwd = c(1)*0.75,
        lty = c("solid")
      ),
      type = "l",
      space = "right",
      cex.title = 1,
      cex = 0.75
    )

  )

### Change Color of Trellis Device and Plot Data
  trellis.device(color = !bw,new = FALSE,retain = FALSE)

### Plot Data
  print(lat.admin)

#dev.off()