R/DTK.plot.R

In DTK: Dunnett-Tukey-Kramer Pairwise Multiple Comparison Test Adjusted for Unequal Variances and Unequal Sample Sizes

DTK.plot <-
function (x = "DTK.test output") {
    out = x[[2]]
    a = x[[1]]
    n = nrow(out)
    plot(c(max(out[, 3]), min(out[, 2])), c(1, n), type = "n",
        xlab = "Mean Difference", ylab = "Mean Comparison", yaxt = "n")
    axis(2, at = seq(1, n), labels = rownames(out))
    title(main = paste(paste(((1 - a) * 100), "%", sep = ""),
        "Confidence Intervals"))
    for (i in 1:n) {
        lines(x=c(-1e100,1e100),y=c(i,i),col="gray",lty=2)
        if(out[i,2]>0 || out[i,3]<0){
          lines(out[i, 2:3], y = c(i, i),col="red",lwd=2)
          points(x = out[i, 1], y = i,col="red")
        } else{
          lines(out[i, 2:3], y = c(i, i))
          points(x = out[i, 1], y = i)
        }  
     }
  }